关键子工程

http://www.cnblogs.com/hongyang/p/3407666.html

题目描述

²在大型过程的施工前，我们经常把整个工程分为若干个子工程，并把这些子工程编号为 1、2、…、N；这样划分之后，子工程之间就会有一些依赖关系，即一些子工程必须在 某些子工程完成之后才能施工。由于子工程之间有相互依赖关系，因此有两个任务需要 我们去完成：首先，我们需要根据每个子工程的完成时间计算整个工程最少的完成时间；

另一方面，由于一些不可预测的客观因素会使某些子工程延期，因此我们必须知道哪些 子工程的延期会影响整个工程的延期，我们把有这种特性的子工程称为关键子工程；因 此第二个任务就是找出所有的关键子工程，以便集中精力管理好这些子工程，尽量避免 这些子工程延期，达到用最快的速度完成整个工程。为了便于编程，现在我们假设： ²根据预算，每一个子工程都有一个完成时间。 ²子工程之间的依赖关系是：部分子工程必须在一些子工程完成之后才开工。 ²只要满足子工程间的依赖关系，在任何时刻可以有任何多个子工程同时在施工， 也既同时施工的子工程个数不受限制。

整个工程的完成是指：所有子工程的完成。

例如，有5个子工程的工程规划表：

其中，表格中第I+1行J+2列的值如为0表示“子工程I”可以在“子过程J”没完成前施工，为I表示“子工程I”必须在“子过程J”完成后才能施工。

上述工程最快完成时间为 14天，其中子工程1、3、4、5为关键子工程。

又例如，有5个子工程的工程规划表：

上述的子工程规划不合理，因为无法安排子工程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行为按有小到大顺序输出所有关键子过程的编号。

样例输入

样例1:
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

样例2:
5
5 4 12 7 2
0 1 0 0
0 0 0 0
0 0 1 0
1 1 0 0
1 1 1 1
样例输出

样例1:
14
1 3 4 5

样例2:
-1

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <stack>
#define maxn 203
using namespace std;
stack<int>s;
int n,w[maxn],d[maxn],e[maxn],l[maxn],cd[maxn],ans=0,ft[maxn],fst[maxn],cnt=0,ct=0,nxt[maxn*maxn],nt[maxn*maxn],to[maxn*maxn],t[maxn*maxn];
int tot=0;
void add(int x,int y){to[++cnt]=y,nxt[cnt]=fst[x];fst[x]=cnt;}
void add2(int x,int y){t[++ct]=y;nt[ct]=ft[x];ft[x]=ct;}
int main(){
    scanf("%d",&n);
    memset(d,0,sizeof(d));
    memset(e,0,sizeof(e));
    memset(l,0x7f,sizeof(l));
    memset(cd,0,sizeof(cd));
    for(int i=1;i<=n;++i)scanf("%d",&w[i]);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j){
			int c;
		    if(j==i)continue;
		    scanf("%d",&c);
		    if(c==1)add(j,i),add2(i,j),++d[i],++cd[j];
		}
		if(d[i]==0)s.push(i),e[i]=w[i];
	}
	while(s.size()){
	    int x=s.top();s.pop();tot++;
	    for(int i=fst[x];i;i=nxt[i]){
		    e[to[i]]=max(e[to[i]],e[x]+w[to[i]]);
		    --d[to[i]];
		    if(d[to[i]]==0)s.push(to[i]);
		}
		d[x]=-1;
	}
	if(tot!=n){printf("-1\n");return 0;}
	for(int i=1;i<=n;++i)ans=max(ans,e[i]);
	for(int i=1;i<=n;++i)if(cd[i]==0)s.push(i),l[i]=ans;
	while(s.size()){
	    int x=s.top();s.pop();
	    for(int i=ft[x];i;i=nt[i]){
		    if(l[t[i]]>l[x]-w[x])l[t[i]]=l[x]-w[x];
		    --cd[t[i]];
		    if(cd[t[i]]==0)s.push(t[i]);
		}
	}
	printf("%d\n",ans);
	for(int i=1;i<=n;++i)if(e[i]==l[i])printf("%d ",i);
	return 0;
}