【gzoj4112】新校园【关键路径】

题意

对于一个总工程（DAG），找出不能“拖延”的子工程，也就是说若“拖延”，就会影响总工程的完成时间。

分析

在一个AOE网上面，不能进行“拖延”的工程所连成的链，就叫关键路径，可能不止一条。

求关键路径只需要正向拓扑一遍，求出整个大工程最快完成的时间，建个反图，再拓扑一遍，求出每个子工程最快完成时间，然后用大工程的最快时间减去子工程的最快完成时间（反图上），若刚好等于正图上子工程的最快完成时间，那么这个工程就在关键路径上。

因为反图上做的最快完成时间其实就是这个点后面所有工程完成的最短时间，也就是说他可以“拖延”到 总时间减去后面最短时间 这个时间点。

上代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
using namespace std;

struct node
{
	int to,next;
}e[100010];
struct node2
{
	int to,next;
}e2[100010];

int n,a[210],in[210],in2[210],f[210][2],cnt;
int hd[100010],tot,hd2[100010],tot2;

void add(int x,int y)
{
	e[++tot]=(node){y,hd[x]};
	hd[x]=tot;
}
void add2(int x,int y)
{
	e2[++tot2]=(node2){y,hd2[x]};
	hd2[x]=tot2;
}

void topsort()
{
	queue<int> q;
	for(int i=1;i<=n;i++)
	{
		if(!in[i])
		{
			q.push(i);
			cnt++;
			f[i][0]=a[i];
		}
	}
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		for(int i=hd[x];i;i=e[i].next)
		{
			f[e[i].to][0]=max(f[e[i].to][0],f[x][0]+a[e[i].to]);
			in[e[i].to]--;
			if(!in[e[i].to]) q.push(e[i].to),cnt++; 
		}
	}
}

void retopsort()
{
	queue<int> q;
	for(int i=1;i<=n;i++)
	{
		if(!in2[i])
		{
			q.push(i);
		}
	}
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		for(int i=hd2[x];i;i=e2[i].next)
		{
			f[e2[i].to][1]=max(f[e2[i].to][1],f[x][1]+a[x]);
			in2[e2[i].to]--;
			if(!in2[e2[i].to]) q.push(e2[i].to); 
		}
	}
}

int main()
{
	cin>>n;
	for(int i=1;i<=n;i++) cin>>a[i];
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(i==j) continue;
			int x;
			scanf("%d",&x);
			if(x==1)
			{
				add(j,i),in[i]++;
				add2(i,j),in2[j]++;
			 } 
		}
	}
	topsort();
//	for(int i=1;i<=n;i++) cout<<f[i][0]<<' ';
//	cout<<endl;
	retopsort();
//	for(int i=1;i<=n;i++) cout<<f[i][1]<<' ';
//	cout<<endl;
	if(cnt<n)
	{
		cout<<-1;
		return 0;
	} 
	int mx=0;
	for(int i=1;i<=n;i++) mx=max(f[i][0],mx);
	cout<<mx<<endl; 
	for(int i=1;i<=n;i++)
	{
		if(mx-f[i][1]==f[i][0])
		{
			cout<<i<<' ';
		}
	}
	return 0;
}