08-图9 关键活动 (30分)

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

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

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

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

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

输入格式:

输入第1行给出两个正整数$N$($\le 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

解析：先做一次拓扑排序（计算入度）算出每个活动的最早完成时间，参照How long does it take， 然后再反过来做一次拓扑排序（计算出度）算出每个活动的最晚完成时间，顺便也算出所有活动的机动时间，最后机动时间等于0的就是关键路径。另一种方法可以用堆栈来存第一次的拓扑序列，这样就不需要算出度了。

#include <cstdio>
#include <cstdlib>
#include <queue>
using namespace std;
#define MAX 105
#define INFINITY 65535
int N, M, A[MAX][MAX], ECT, EarliestTime[MAX] = {0}, LatestTime[MAX], D[MAX][MAX], idx;  //ETC--earliest completion time
int getMax( int arr[] ) {
	int max = 0;
	for(int i = 0; i < N; i++)
		if( max < arr[i] ){
			max = arr[i];
			idx = i;
		}
	return max;
}
int TopSort_Earliest(){
	int V, cnt = 0, Indegree[MAX] = {0};
	queue<int> q;
	//计算各结点的入度
	for( int i = 0; i < N; i++ )
		for( int j = 0; j < N; j++ )
			if( A[i][j] != INFINITY )
				Indegree[j]++;	//对于有向边<i,j>累计终点j的入度
	//入度为0的入队
	for( int i = 0; i < N; i++ )
		if( Indegree[i] == 0 )
			q.push(i);
	while( !q.empty() ) {
		V = q.front();
		q.pop();
		cnt++;
		for( int j = 0; j < N; j++ )
			if( A[V][j] != INFINITY ) {  //<V, j>有有向边
				if( EarliestTime[V] + A[V][j] > EarliestTime[j] ){ //如果 V的最早完成时间 + j所需时间 > j的最早完成时间
					EarliestTime[j] = EarliestTime[V] + A[V][j];
				}
				if( --Indegree[j] == 0 )  //去掉V后，如果j的入度为0
					q.push(j);
			}
	}
	ECT = getMax(EarliestTime); //最早完成时间应是所有元素中最大的
	if( cnt != N ) return 0;	//如果没有取出所有元素，说明图中有回路
	else return 1;
}
void TopSort_Latest() {
	int V, Outdegree[MAX] = {0};
	queue<int> q;
		//计算各结点的出度
	for( int i = 0; i < N; i++ )
		for( int j = 0; j < N; j++ )
			if( A[i][j] != INFINITY )
				Outdegree[i]++;	//对于有向边<i,j>累计起点i的出度
	//出度为0的入队
	for( int i = 0; i < N; i++ )
		if( Outdegree[i] == 0 )
			q.push(i);
	//初始化LatestTime
	for( int i = 0; i < N; i++ )
		LatestTime[i] = INFINITY;
	LatestTime[idx] = ECT;	//将最后一个活动的最晚完成时间设为它的最早完成时间
	while( !q.empty() ) {
		V = q.front();
		q.pop();
		//cnt++;	//不需要再算cnt了
		for( int j = 0; j < N; j++ )
			if( A[j][V] != INFINITY ) {  //<j, V>有有向边
				if( LatestTime[V] - A[j][V] <= LatestTime[j] ) {  //必须用<=,只<的话只能算一条关键路径，<=才能算出所有的关键路径(错误原因)
					LatestTime[j] = LatestTime[V] - A[j][V];
					D[j][V] = LatestTime[V] - EarliestTime[j] - A[j][V];
				}
				if( --Outdegree[j] == 0 )  //去掉V后，如果j的出度为0
					q.push(j);
			}
	}
}
void PrintKeyRoute() {
	for( int i = 0; i < N; i++ )
		for( int j = N - 1; j >= 0; j-- )	//根据题目要求，i相同时要j要逆序输出
			if( D[i][j] == 0 )
				printf("%d->%d\n", i + 1, j + 1);
}
int main(){
	int a, b;
	scanf("%d %d", &N, &M);
	//初始化图的边A，各组的机动时间D
	for( int i = 0; i < N; i++ )
		for( int j = 0; j < N; j++ )
			D[i][j] = A[i][j] = INFINITY;
	//read
	for( int i = 0; i < M; i++ ) {
		scanf("%d %d", &a, &b);
		scanf("%d", &A[--a][--b]);	//题目中编号从1开始（错误原因）
	}
	if( !TopSort_Earliest() )
		printf("0\n");
	else {
		printf("%d\n", ECT);
		TopSort_Latest();
		PrintKeyRoute();
	}
	system("pause");
	return 0;
}