08-图9 关键活动（测试点2和5）

中国大学MOOC-陈越、何钦铭-数据结构

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

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

比如完成一个专业的所有课程学习和毕业设计可以看成一个本科生要完成的一项工程，各门课程可以看成是子任务。有些课程可以同时开设，比如英语和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

算法思想：

问题为AOE网求关键路径，抽象为有向图，用邻接表存储图，在存储时记录每个结点的入度和出度（出度在这题没用上），然后根据入度进行拓扑排序，记录拓扑序列（很关键！关系到AOE结点的最早和最晚开始时间的计算），如果拓扑序列不存在，即图有回路，直接判0，否则进行关键路径的两个参数（最早和最晚时间）的求解（注：以关键路径终点的最早/晚时间初始化所以结点的最晚时间，可筛除非关键路径），最后对两个相连的A和B结点， 根据B的最晚时间减去A->B这条边的权值是否等于A的最晚时间来输出关键路径。（算法这玩意，想通半小时，想不通半天= =！）

测试点2和5就是拓扑序的问题。

做题过程中参考了  我只有一件白T恤  的算法，十分感谢！链接如下：

https://www.cnblogs.com/zengguoqiang/p/8447386.html

/***************2019.5.6-09:22-12:00（V1.0）  12:47-13:06   ***********/
/***************2019.5.7-09:20-11:50（V2.0）  12:30-12:55(V3.0) ***********/
//08-图9 关键活动V3.0   6h
//关键活动输出的顺序规则是：任务开始的交接点编号小者优先，-------输出时从编号小的开始遍历 
//起点编号相同时，与输入时任务的顺序相反。 -----------插入边表时直接头插法，保证输入顺序相反。 
//一定要用拓扑序来取最早和最晚时间，因为这个时间算法是迭代的，前后间有联系，一定要按序取才能保证正确！                               
#include<stdio.h>
#include<stdlib.h> 
#define MaxVertexNum 108
#include<queue>
using namespace std;

typedef struct AdjVNode *AdjV;
struct AdjVNode{
	int Data;
	int Weight;
	AdjV Next;
};

struct GNode{
	int Data;
	AdjV FirstEdge;
	int Earliest;
	int Lastest; 
	int InD;
	int OutD;
}G[MaxVertexNum];

int main(){
	int N,M;
	int i,j,k;
	int v1,v2,w;
	scanf("%d%d",&N,&M);
	
	//*******************
	//初始化图
	//*******************	
	for(i=1;i<=N;i++){
		G[i].Data=i;
		G[i].FirstEdge=NULL;
		G[i].Earliest=0;
		G[i].Lastest=999999;
		G[i].InD=0;
		G[i].OutD=0;
	} 
	for(i=1;i<=M;i++){
		scanf("%d%d%d",&v1,&v2,&w);
		AdjV X=(AdjV)malloc(sizeof (struct AdjVNode));
		X->Data=v2;
		X->Weight=w;
		X->Next=G[v1].FirstEdge;
		G[v1].FirstEdge=X;
		G[v1].OutD++;
		G[v2].InD++;
	}
	
	//*******************
	//拓扑排序
	//******************* 
	int count=0,CP=0; //CP:cost of Critical Path
	queue<int>Q;
	int Topo[N+1];
	for(i=1;i<=N;i++){
		if(G[i].InD==0){
			Q.push(i);
			count++;
			Topo[count]=i;
		}	
	}
	while(!Q.empty()){
		i=Q.front();
		Q.pop();
		AdjV tmp=G[i].FirstEdge;
		while(tmp){
			j=tmp->Data;
			if(--G[j].InD==0){
				Q.push(j);
				count++;
        	    Topo[count]=j;
			}
			tmp=tmp->Next;
		}
	}
	if(count!=N) printf("0\n");

	else{
		//*******************
		// Get Earliest
		//*******************
		for(i=1;i<=N;i++){
			k=Topo[i];
			AdjV tmp=G[k].FirstEdge;
			while(tmp){
				j=tmp->Data;
				if(G[k].Earliest+tmp->Weight>G[j].Earliest){ //Earliest
					G[j].Earliest=G[k].Earliest+tmp->Weight;
					if(!G[j].FirstEdge){
						G[j].Lastest=G[j].Earliest;
						if(G[j].Lastest>CP) CP=G[j].Lastest;
					} 
				}
				tmp=tmp->Next;
			}
			
		}
		//*******************
		// Get Lastest
		//*******************
		for(i=N;i>0;i--) G[i].Lastest=CP; //very important!!!!remove the non-Critical Path!!!
		for(i=N;i>0;i--){
        	k=Topo[i];
			AdjV tmp=G[k].FirstEdge;
			while(tmp){
				j=tmp->Data;
				if(G[j].Lastest-tmp->Weight<G[k].Lastest){ //Lastest
					G[k].Lastest=G[j].Lastest-tmp->Weight;
				}
				tmp=tmp->Next;
			}
		}
		/*******************/
		//Print
		/*******************/
		printf("%d\n",CP);
		for(i=1;i<=N;i++){
			AdjV tmp=G[i].FirstEdge;
			while(tmp){
				j=tmp->Data;
	    		if(G[i].Earliest==G[j].Lastest-tmp->Weight){
					printf("%d->%d\n",i,j);
				}
				tmp=tmp->Next;
			}
		}
	}
	
    return 0;	
}