一、 关键子工程（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.