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

自创输入样例

16 20
1 2 6
1 3 4
1 4 5
2 5 1
3 5 1
4 6 2
6 5 0
5 7 9
5 8 7
6 8 4
7 9 2
8 9 4
10 11 4
10 12 3
11 13 5
12 13 3
13 14 1
13 15 6
14 16 5
15 16 2

输出样例

18
1->4
1->2
2->5
4->6
5->8
5->8
5->7
5->7
6->5
7->9
7->9
8->9
8->9
10->11
11->13
13->15
15->16

Note

注意有边的条件是>=0
多边的测试点始终过不了，自己测试的多变可以（将本题与递归排序那题样例结合）

Code

#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
#define MAX 100struct result{int x;int y;
}res[MAX];
int cmp(result a, result b){return a.x == b.x ? a.y > b.y : a.x < b.x;
}
int a[MAX][MAX], num, arcnum, cnt = 0, e[MAX],l[MAX],Indegree[MAX],outdegree[MAX], temp, size = 0;
void Dfs(int i){for(int j = 0; j < num; j++){if(a[i][j] >= 0 && l[j] - a[i][j] - e[i] == 0){res[size].x = i;res[size].y = j;size ++;Dfs(j);}}
}
int main() {cin >> num >>arcnum;int start;for(int i = 0; i < num; i++){       //初始化for(int j = 0; j < num; j++){a[i][j] = -1;}e[i] = 0;l[i] = MAX;Indegree[i] = outdegree[i] = 0;}for(int i = 0; i < arcnum; i++){    //输入int tempa, tempb, tempc;cin >> tempa >> tempb >> tempc;a[tempa - 1][tempb - 1] = tempc;}queue<int> q1, q2, qstart;for(int i = 0; i < num; i++){       //计算初始时的入度出度for(int j = 0; j < num; j++){if(a[j][i] >= 0) Indegree[i]++;if(a[i][j] >= 0) outdegree[i]++;}if(Indegree[i] == 0) { q1.push(i); cnt++; qstart.push(i);}if(outdegree[i] == 0){ q2.push(i); }}while(q1.size() > 0){                //拓扑排序计算e[]temp = q1.front();q1.pop();int visted = 0;for(int i = 0; i < num; i++){if(a[temp][i] >= 0){visted = 1;Indegree[i]--;if(a[temp][i] + e[temp] > e[i]) e[i] = a[temp][i] + e[temp];if(Indegree[i] == 0){q1.push(i); cnt++;}}}if(visted == 0) l[temp] = e[temp];}if(cnt != num){ cout << 0; return 0;}while(q2.size() > 0){                //反向拓扑排序计算l[]temp = q2.front();q2.pop();for(int i = 0; i < num; i++){if(a[i][temp] >= 0){outdegree[i]--;if(l[temp] - a[i][temp] < l[i]) l[i] = -a[i][temp] + l[temp];if(outdegree[i] == 0){q2.push(i);}}}}while(qstart.size() > 0){int temp = qstart.front();qstart.pop();Dfs( temp);}sort(res, res + size, cmp);sort(e, e + num);cout  <<  e[num -1] << endl;for(int i = 0; i < size; i++){cout << res[i].x + 1 << "->" << res[i].y + 1 << endl;}return 0;
}

Note2

利用栈存储拓扑序便可以直接得到逆拓扑序
用邻接矩阵邻接表存储，空间换时间
输出顺序问题只需要改一下循环的方向即可

Code2

#include<iostream>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
int num,arcnum;
vector<int>adj[110]; //邻接表
int a[110][110];    //邻接矩阵
int indegree[110]={0};
int early[110]={0};
int late[110];
int main(){int i,j,temp1,temp2,temp3;scanf("%d %d",&num,&arcnum);for(i=0;i<arcnum;i++){                      //从1开始存scanf("%d %d %d",&temp1,&temp2,&temp3);indegree[temp2]++;adj[temp1].push_back(temp2);a[temp1][temp2]=a[temp2][temp1]=temp3;}queue<int>q;for(i=1;i<=num;i++){if(indegree[i]==0){q.push(i);}}int count=0;stack<int> s;while(!q.empty()){count++;int head=q.front();s.push(head); q.pop();for(i=0;i<adj[head].size();i++){int next=adj[head][i];if(early[next]<early[head]+a[head][next]){early[next]=early[head]+a[head][next];}if(--indegree[next]==0){q.push(next);}}}if(count!=num){printf("0");return 0;}printf("%d\n",early[s.top()]); for(i=1;i<=num;i++){late[i]=early[s.top()];}while(!s.empty()){int head=s.top();s.pop();for(i=0;i<adj[head].size();i++){int next=adj[head][i];if(late[head]>late[next]-a[head][next]){late[head]=late[next]-a[head][next];}}}for(i=1;i<=num;i++){for(j=adj[i].size()-1;j>=0;j--){ int next=adj[i][j];if(late[next]-early[i]-a[i][next]==0){printf("%d->%d\n",i,next);}}}
}