08-图9 关键活动

假定一个工程项目由一组子任务构成，子任务之间有的可以并行执行，有的必须在完成了其它一些子任务后才能执行。“任务调度”包括一组子任务、以及每个子任务可以执行所依赖的子任务集。

比如完成一个专业的所有课程学习和毕业设计可以看成一个本科生要完成的一项工程，各门课程可以看成是子任务。有些课程可以同时开设，比如英语和C程序设计，它们没有必须先修哪门的约束；有些课程则不可以同时开设，因为它们有先后的依赖关系，比如C程序设计和数据结构两门课，必须先学习前者。

但是需要注意的是，对一组子任务，并不是任意的任务调度都是一个可行的方案。比如方案中存在“子任务A依赖于子任务B，子任务B依赖于子任务C，子任务C又依赖于子任务A”，那么这三个任务哪个都不能先执行，这就是一个不可行的方案。

任务调度问题中，如果还给出了完成每个子任务需要的时间，则我们可以算出完成整个工程需要的最短时间。在这些子任务中，有些任务即使推迟几天完成，也不会影响全局的工期；但是有些任务必须准时完成，否则整个项目的工期就要因此延误，这种任务就叫“关键活动”。

请编写程序判定一个给定的工程项目的任务调度是否可行；如果该调度方案可行，则计算完成整个工程项目需要的最短时间，并输出所有的关键活动。

输入格式:

输入第1行给出两个正整数N(≤100)和M，其中N是任务交接点（即衔接相互依赖的两个子任务的节点，例如：若任务2要在任务1完成后才开始，则两任务之间必有一个交接点）的数量。交接点按1_{N编号，M是子任务的数量，依次编号为1}M。随后M行，每行给出了3个正整数，分别是该任务开始和完成涉及的交接点编号以及该任务所需的时间，整数间用空格分隔。

输出格式:

如果任务调度不可行，则输出0；否则第1行输出完成整个工程项目需要的时间，第2行开始输出所有关键活动，每个关键活动占一行，按格式“V->W”输出，其中V和W为该任务开始和完成涉及的交接点编号。关键活动输出的顺序规则是：任务开始的交接点编号小者优先，起点编号相同时，与输入时任务的顺序相反。

输入样例:

7 8
1 2 4
1 3 3
2 4 5
3 4 3
4 5 1
4 6 6
5 7 5
6 7 2

输出样例:

17
1->2
2->4
4->6
6->7

解题思路

先通过拓扑排序判断是任务调度是否可行，即否为有向无环图，同时求的事件i的最早发生事件ve[i]，ve[i]为结点i的所有前继结点与边权之和的最大值。
由于事先不知道汇点编号，所以取ve数组的最大值，由于ve数组的含义是事件的最早开始时间，因此所有事件中ve最大的一定是最后一个事件，也就是汇点
存储ve数组的最大值maxLength，作为vl数组的初始值，vl[i]为结点i的所有后继结点与边权之差的最小值。由于对活动i->j而言，最早开始时间e为ve[i]，最迟开始时间l为vl[j]-w(i->j的边权)，如果e和l相等，则该活动为关键活动

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
const int MAXV = 105;
const int INF = 65535;
struct Node{
	int v,w; //v为后继结点的编号，w为边权 
};
int N,M;
vector<Node> G[MAXV];
stack<int> topOrder; //逆拓扑序列
int inDegree[MAXV] = {0}, ve[MAXV] = {0}, vl[MAXV]; 
bool topSort(){
	queue<int> q;
	for(int i = 1; i <= N; i++)
		if(inDegree[i] == 0)
			q.push(i); //入度为0的入队
	while(!q.empty()){
		int u = q.front();
		q.pop();
		topOrder.push(u);
		for(int i = 0; i < G[u].size(); i++){
			int v = G[u][i].v;
			if(--inDegree[v] == 0)
				q.push(v);
			if(ve[v] < ve[u] + G[u][i].w)
				ve[v] = ve[u] + G[u][i].w;
		}
	}
	if(topOrder.size() == N) return true;
	return false;
}
void criticalPath(int maxLength){ //输出关键路径 
	fill(vl + 1, vl + N + 1, maxLength); //用关键路径长度填充所有的最迟完成时间 
	while(!topOrder.empty()){ //逆拓扑序列求最迟完成时间 
		int u = topOrder.top();
		topOrder.pop();
		for(int i = 0; i < G[u].size(); i++){
			int v = G[u][i].v; //取出后继结点 
			if(vl[u] > vl[v] - G[u][i].w) //最迟完成时间为后继结点的最迟完成时间减去边权中值最小的 
				vl[u] = vl[v] - G[u][i].w; //更新最迟完成时间 
		}
	}
	for(int u = 1; u <= N; u++){ //按起点编号顺序 
		for(int i = G[u].size() - 1; i >= 0; i--){//按输入顺序的倒序 
			int v = G[u][i].v, w = G[u][i].w;
			int e = ve[u], l = vl[v] - w; //活动的最早开始时间e和最迟开始时间l 
			if(e == l)
				printf("%d->%d\n",u,v);
		}
	}
}
int main(){
	int u,v,w;
	cin>>N>>M;
	for(int i = 0; i < M; i++){
		cin>>u>>v>>w;
		Node node = {v,w};
		G[u].push_back(node);
		inDegree[v]++;
	}
	if(!topSort())
		printf("0\n");
	else{
		int ETC = 0;
		for(int i = 1; i <= N; i++)
			if(ETC < ve[i]) ETC = ve[i];
		printf("%d\n", ETC); //关键路径长度
		criticalPath(ETC); 
	}
	return 0;
}