7-11 关键活动（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为该任务开始和完成涉及的交接点编号。关键活动输出的顺序规则是：任务开始的交接点编号小者优先，起点编号相同时，与输入时任务的顺序相反。
1 #include<iostream> 2 #include<vector> 3 #include<queue> 4 #include<stack> 5 using namespace std; 6 #define maxNv 101 7 #define maxtime 100000 8 int flag=1; 9 vector<int> earliest(maxNv,0); 10 queue<int> q; 11 vector<int> rudu(maxNv,0); 12 vector<int> rudu2(maxNv,0); 13 vector<int> start(maxNv,maxtime); 14 vector<int> chudu(maxNv,0); 15 int mode[maxNv][maxNv]; 16 struct activity{ 17 int start,end; 18 }; 19 struct graph{ 20 int Nv,Ne; 21 int G[maxNv][maxNv]; 22 }; 23 using Graph=graph*; 24 struct enode{ 25 int v1,v2; 26 int time; 27 }; 28 using edge=enode*; 29 Graph CreateGraph() 30 { 31 Graph gra=new graph(); 32 cin>>gra->Nv>>gra->Ne; 33 for(int i=1;i<=gra->Nv;i++) 34 for(int j=1;j<=gra->Nv;j++) 35 gra->G[i][j]=maxtime; 36 return gra; 37 } 38 void Insert(edge e,Graph gra){ 39 gra->G[e->v1][e->v2]=e->time; 40 rudu[e->v2]++; 41 rudu2[e->v2]++; 42 chudu[e->v1]++; 43 } 44 Graph BuildGraph(){ 45 Graph gra=CreateGraph(); 46 edge e=new enode(); 47 for(int i=1;i<=gra->Ne;i++){ 48 cin>>e->v1>>e->v2>>e->time; 49 Insert(e,gra); 50 } 51 return gra; 52 } 53 void timeset(int v1,Graph gra){ 54 int f=0; 55 for(int i=1;i<=gra->Nv;i++) 56 if(gra->G[i][v1]!=maxtime&&earliest[i]+gra->G[i][v1]>f) 57 f=earliest[i]+gra->G[i][v1]; 58 earliest[v1]=f; 59 } 60 void solve(Graph gra){ 61 int v1=0,v2=0; 62 for(int i=1;i<=gra->Nv;i++) 63 if(rudu[i]==0) { 64 q.push(i); 65 earliest[0]=0; 66 } 67 while(!q.empty()){ 68 v1=q.front(); 69 q.pop(); 70 rudu[v1]--; 71 timeset(v1,gra); 72 for(v2=1;v2<=gra->Nv;v2++) 73 if(gra->G[v1][v2]!=maxtime&&rudu[v2]!=-1){ 74 rudu[v2]--; 75 if(rudu[v2]==0) q.push(v2);} 76 } 77 // for(int i=0;i<gra->Nv;i++) 78 // cout<<rudu[i]<<" "; 79 // cout<<endl; 80 for(int i=1;i<=gra->Nv;i++) 81 if(rudu[i]!=-1) flag=0; 82 if(flag==0) cout<<0; 83 else{ 84 int time=-1; 85 for(int i=1;i<=gra->Nv;i++) 86 if(earliest[i]>time) time=earliest[i]; 87 cout<<time<<endl; 88 queue<int> Q; 89 for(int i=1;i<=gra->Nv;i++) 90 if(chudu[i]==0){ 91 Q.push(i); 92 start[i]=time; 93 } 94 while(!Q.empty()){ 95 v1=Q.front(); 96 Q.pop(); 97 for(v2=1;v2<=gra->Nv;v2++) 98 if(gra->G[v2][v1]!=maxtime&&start[v1]-gra->G[v2][v1]<start[v2]){ 99 start[v2]=start[v1]-gra->G[v2][v1]; Q.push(v2); 100 } 101 } 102 for(int i=1;i<=gra->Nv;i++) 103 for(int j=1;j<=gra->Nv;j++) 104 mode[i][j]=0; 105 for(int i=1;i<=gra->Nv;i++) 106 if(rudu2[i]==0){ 107 Q.push(i); rudu2[i]=-1; 108 } 109 while(!Q.empty()){ 110 v1=Q.front(); 111 Q.pop(); 112 for(v2=1;v2<=gra->Nv;v2++) 113 if(gra->G[v1][v2]!=maxtime){ 114 if(rudu2[v2]!=-1){ 115 Q.push(v2); rudu2[v2]=-1;} 116 if(start[v2]-earliest[v1]==gra->G[v1][v2]) 117 mode[v1][v2]=1; 118 } 119 } 120 for(int i=1;i<=gra->Nv;i++) 121 for(int j=gra->Nv;j>=1;j--) 122 if(mode[i][j]==1) 123 cout<<i<<"->"<<j<<endl; 124 } 125 } 126 int main() 127 { 128 Graph gra=BuildGraph(); 129 solve(gra); 130 return 0; 131 }