一、 关键子工程（project.c/cpp/pas） 拓扑排序+aov

 一、 关键子工程（project.c/cpp/pas）
在大型工程的施工前，我们把整个工程划分为若干个子工程，并把这些子工程编号为1、
2、……、N；这样划分之后，子工程之间就会有一些依赖关系，即一些子工程必须在某些
子工程完成之后才能施工。由于子工程之间有相互依赖关系，因此有两个任务需要我们去完
成：首先，我们需要计算整个工程最少的完成时间；同时，由于一些不可预测的客观因素会
使某些子工程延期，因此我们必须知道哪些子工程的延期会影响整个工程的延期，我们把有
这种特征的子工程称为关键子工程，因此第二个任务就是找出所有的关键子工程，以便集中
精力管理好这些子工程，尽量避免这些子工程延期，达到用最快的速度完成整个工程。为了
便于编程，现在我们假设：
（1）根据预算，每一个子工程都有一个完成时间。
（2）子工程之间的依赖关系是：部分子工程必须在一些子工程完成之后才开工。
（3）只要满足子工程间的依赖关系，在任何时刻可以有任何多个子工程同时在施工，
也既同时施工的子工程个数不受限制。
（4）整个工程的完成是指：所有子工程的完成。
例如，有五个子工程的工程规划表：
序号  完成时间  子工程 1  子工程 2  子工程 3  子工程 4  子工程5
子工程1  5    0  0  0  0
子工程2  4  0    0  0  0
子工程3  12  0  0    0  0
子工程4  7  1  1  0    0
子工程5  2  1  1  1  1  
其中，表格中第 I+1 行 J+2 列的值如为 0 表示“子工程 I”可以在“子工程 J”没完成
前施工，为 1 表示“子工程 I”必须在“子工程 J”完成后才能施工。上述工程最快完成时
间为14天，其中子工程1、3、4、5为关键子工程。
 
又例如，有五个子工程的工程规划表：
序号  完成时间  子工程1  子工程2  子工程3  子工程4  子工程5
子工程1  5    0  1  0  0
子工程2  4  0    0  0  0
子工程3  12  0  0    1  0
子工程4  7  1  1  0    0
子工程5  2  1  1  1  1  
上述的子工程划分不合理，因为无法安排子工程1，3，4的施工。
 
输入数据：
第1行为N，N是子工程的总个数，N≤200。
第2行为N个正整数，分别代表子工程1、2、……、N的完成时间。
第 3 行到 N+2 行，每行有 N-1 个 0 或 1。其中的第 I+2 行的这些 0，1，分别表示“子
工程I”与子工程1、2、…、I-1、I+1、…N的依赖关系，（I=1、2、……、N）。每行数据之
间均用一个空格分开。
 
输出数据：

如子工程划分不合理，则输出-1；
如子工程划分合理，则用两行输出：第1行为整个工程最少的完成时间。第2行为按由
小到大顺序输出所有关键子工程的编号。
 
样例：
输入文件名：project.in
5                             
5 4 12 7 2                  
0 0 0 0
0 0 0 0
0 0 0 0
1 1 0 0
1 1 1 1
 
输出文件名：project.out
14
1 3 4 5 
 

 

program project;var r:array [0..200,0..200] of longint;    s,d,t,f:array [0..200] of longint;    k:array [0..200] of boolean;    i,j,n,ans,tot:longint;procedure print(ok:boolean);var i:longint;begin assign(output,'project.out'); rewrite(output); if ok then  begin   writeln(ans);   for i:=1 to n do    if k[i] then write(i,' ');   writeln;  end else writeln(-1); close(output); halt;end;begin assign(input,'project.in'); reset(input); readln(n); for i:=1 to n do  read(t[i]); for i:=1 to n do  for j:=1 to n do   if i<>j then    read(r[i,j]); close(input); for i:=1 to n do  for j:=1 to n do   inc(d[i],r[i,j]); tot:=0; for i:=1 to n do  if d[i]=0 then   begin    inc(tot);    s[tot]:=i;   end; i:=0; ans:=0; repeat  inc(i);  for j:=1 to n do   if r[s[i],j]=1 then    if f[j]>f[s[i]] then f[s[i]]:=f[j];  inc(f[s[i]],t[s[i]]);  if f[s[i]]>ans then ans:=f[s[i]];  for j:=1 to n do   if d[j]<>0 then    begin     dec(d[j],r[j,s[i]]);     if d[j]=0 then      begin       inc(tot);       s[tot]:=j;      end;    end; until i=tot; if tot<>n then print(false); for i:=tot downto 1 do  begin   if f[s[i]]=ans then    k[s[i]]:=true;   for j:=1 to n do    if (r[j,s[i]]=1) and (k[j]) and (f[s[i]]+t[j]=f[j]) then     k[s[i]]:=true;  end; print(true);end.