关键活动 (30分)【C语言】拓扑排序中找到 关键路径

习题讲解视频

题目：

假定一个工程项目由一组子任务构成，子任务之间有的可以并行执行，有的必须在完成了其它一些子任务后才能执行。“任务调度”包括一组子任务、以及每个子任务可以执行所依赖的子任务集。
比如完成一个专业的所有课程学习和毕业设计可以看成一个本科生要完成的一项工程，各门课程可以看成是子任务。有些课程可以同时开设，比如英语和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

知识点回顾

• AOV网络中，边（弧）代表活动，边上的权值代表活动持续的时间。结点表示事件，即活动开始或者结束的标志。
• 拓扑排序能够得到拓扑序列用来检验活动安排是否合理（即不会出现某事件成环的现象）。
• 两次拓扑排序可以求出关键路径。

算法

拓扑排序

• 正向拓扑排序得到 　earliest （每个事件最早完成时间）
• 反向拓扑排序得到　latest（每个事件最晚完成时间）
• 根据 earliest [ i ]+G->GraphMatrix [ i ] [ j ] == latest [ j ] 判断 <i , j> 是否为关键路径

代码实现

int main()
{	
	int N,E;
	scanf("%d %d",&N,&E);
	Graph G=CreateGraph(N);
	BuildGraph(G,E);
	int EarliestTime=TopSort(G);
	if(EarliestTime!=-1){
	    printf("%d\n",EarliestTime);
	    ReverseTopSort(G,EarliestTime);
	    int i,j;
	    for(i=1;i<=G->VertexNum;i++){
	        for(j=G->VertexNum;j>0;j--){
	             //关键路径的判别条件earliest[i]+G->GraphMatrix[i][j]==latest[j]
	            if((G->GraphMatrix[i][j]!=-1)&&(earliest[i]+G->GraphMatrix[i][j]==latest[j])){
	                printf("%d->%d\n",i,j);
	            }
	        }
	    }
	}else{
	    printf("0");
	}
	return 0;
}
		

TopSort函数

说明：正向拓扑排序得到 earliest（每个事件最早完成时间）

int TopSort(Graph G)
{ 
	int i;
	int *Queue=(int*)malloc(sizeof(int)*(G->VertexNum+1));
	int rear=0,head=0;
	int count=0;//当前收录到集合中元素个数 
	for(i=1;i<=G->VertexNum;i++){
	    if(indegree[i]==0){ 
	        Queue[rear++]=i;
	    }
	}
	int t,MaxVertex;
	while(rear>head){
	    t=Queue[head++];
	    count++;
	    for(i=1;i<=G->VertexNum;i++){
	        if(G->GraphMatrix[t][i]!=-1){
	            if(--indegree[i]==0){
	                Queue[rear++]=i;
	            }
	            if(earliest[i]<earliest[t]+G->GraphMatrix[t][i]){
	                earliest[i]=earliest[t]+G->GraphMatrix[t][i];
	            } 
	        }
	    } 
	}
	free(Queue);
	if(count==G->VertexNum){
	    int MaxSchedule=0;
	    for(i=1;i<=G->VertexNum;i++){
	        if(earliest[i]>MaxSchedule){
	            MaxSchedule=earliest[i];
	        }
	    }
	    return MaxSchedule;
	}else{
	    return -1;
	}
}

ReverseTopSort函数

说明：反向拓扑排序得到latest （每个事件最晚完成时间）

void ReverseTopSort(Graph G,int EarliestTime)//形式上与TopSort相同 
{ 
	int i,j;
	int *Queue=(int*)malloc(sizeof(int)*(G->VertexNum+1));
	int rear=0,head=0;
	for(i=1;i<=G->VertexNum;i++){
	    latest[i]=EarliestTime;//这里的latest要全部初始化为最大值 
	    if(outdegree[i]==0){ 
	        Queue[rear++]=i;
	    }
	}
	int t;
	while(rear>head){
	    t=Queue[head++];
	    for(i=1;i<=G->VertexNum;i++){
	        if((G->GraphMatrix[i][t]!=-1)){//注意这里的结点顺序G->GraphMatrix[i][t] 
	            if(--outdegree[i]==0){
	                Queue[rear++]=i;
	            }
	            if(latest[i]>latest[t]-G->GraphMatrix[i][t]){//注意这里的结点顺序G->GraphMatrix[i][t] 
	                latest[i]=latest[t]-G->GraphMatrix[i][t];
	            } 
	        }
	    } 
	}
	free(Queue);
}

邻接矩阵构造的图

说明：

• void BuildGraph(Graph G,int E)：在图中插入边的同时初始化 indegree（入度）和 outdegree（出度）

#define MAXVERTEXNUM 101
int earliest[MAXVERTEXNUM];//记录每个事件到达的最早时间 
int latest[MAXVERTEXNUM];//记录每个事件到达的最晚时间 
int indegree[MAXVERTEXNUM];//记录每个事件的入度 
int outdegree[MAXVERTEXNUM];//记录每个事件的出度 
typedef struct GNode* Graph;
struct GNode{
	int VertexNum;
	int EdgeNum;
	int GraphMatrix[MAXVERTEXNUM][MAXVERTEXNUM];
};
typedef struct ENode* edge;
struct ENode{
	int V;
	int W;
	int Weight;
};
Graph CreateGraph(int N)
{
	Graph G=(Graph)malloc(sizeof(struct GNode));
	G->VertexNum=N;
	int i,j;
	for(i=1;i<=N;i++){
	    for(j=1;j<=N;j++){
	        G->GraphMatrix[i][j]=-1;
	    }
	}
	return G;
}
void InsertEdge(Graph G,edge L)
{
	G->GraphMatrix[L->V][L->W]=L->Weight;
}
void BuildGraph(Graph G,int E)//在图中插入边的同时初始化入度和出度 
{
	G->EdgeNum=E;
	edge L=(edge)malloc(sizeof(struct ENode));
	int i;
	for(i=0;i<E;i++){
	    scanf("%d %d %d",&(L->V),&(L->W),&(L->Weight));
	    indegree[L->W]++;
	    outdegree[L->V]++; 
	    InsertEdge(G,L);
	}
	free(L);
}