图论-关键路径-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为该任务开始和完成涉及的交接点编号。关键活动输出的顺序规则是：任务开始的交接点编号小者优先，起点编号相同时，与输入时任务的顺序相反。

输入样例:
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

解题思路:
1、建图:增加两个超级节点S、T
2、在g上跑拓扑排序(若不是DAG,输出0)
3、如果是DAG按拓扑排序的顺序依次松弛每个节点(期间记录criPath数组-到这个节点最长路的节点)
4、由每个节点的criPath构成一个图G,在G上跑一遍DFS；
5、输出答案

Code:

#include<iostream>
#include<vector>
#include<queue>
#include<algorithm>
#include<numeric>
using namespace std;
struct T{
	int next;
	int dur;
};
struct D{
	int pre;
	int next;
};
vector<vector<T>>g;
vector<int>est;
vector<int>f,s;
vector<int>Path;
vector<D>resPath;
vector<vector<int>>criPath;
vector<bool>visit;
vector<bool>visited;
int N,M;
int flagCycle=0,ts=0;
void TopSort(int x){
	s[x]=++ts;
	visit[x]=true;
	for(int i=0;i<g[x].size();i++){
		if(s[g[x][i].next]==0)
			TopSort(g[x][i].next);
		else if(f[g[x][i].next]==0&&s[x]>s[g[x][i].next]){
			flagCycle=1;//存在回路 
			return;
		}
	}
	Path.push_back(x);
	f[x]=ts++;
	
}
void criticalAc(int x){
	for(int i=0;i<g[x].size();i++){
		if(est[x]+g[x][i].dur>est[g[x][i].next]){
			est[g[x][i].next]=est[x]+g[x][i].dur;
			criPath[g[x][i].next].clear();//找到比以前更长 
			criPath[g[x][i].next].push_back(x);
		}else if(est[x]+g[x][i].dur==est[g[x][i].next]){//路径长度一样 也需要加进去 也算一条路径 
			criPath[g[x][i].next].push_back(x);
		}
	}
}
void ResP(int x){
	visited[x]=true;
	for(int i=0;i<criPath[x].size();i++){
		if(!visited[criPath[x][i]])
			ResP(criPath[x][i]);
		resPath.push_back({criPath[x][i],x});
	}
}
bool cmp(D e1,D e2){
	return e1.pre==e2.pre?e1.next>e2.next:e1.pre<e2.pre;
}
int main(){
	cin>>N>>M;
	g.resize(N+2);
	numEdg.resize(N+2);
//	RAP.resize(N+2);
	est.resize(N+2);
	s.resize(N+2);
	f.resize(N+2);
	visit.resize(N+2);
	visited.resize(N+2);
	criPath.resize(N+2);	
	for(int i=0;i<M;i++){
		int a,b,c;
		cin>>a>>b>>c;
		g[a].push_back({b,c});
		g[0].push_back({a,0});
		g[b].push_back({N+1,0});
		numEdg[b]++;
		numEdg[a]++;
		numEdg[N+1]++;
	}
	TopSort(0);
	for(int i=0;i<=N+1;i++){
		if(!visit[i])
			TopSort(i);
	}
	if(flagCycle)
		cout<<0;
	else{
		for(int i=N+1;i>=0;i--)
			criticalAc(Path[i]);
		cout<<est[N+1]<<endl;
		ResP(N+1); 
		sort(resPath.begin(),resPath.end(),cmp);
		for(int i=0;i<resPath.size();i++){
			if(resPath[i].pre==0||resPath[i].next==N+1)
				continue;
			cout<<resPath[i].pre<<"->"<<resPath[i].next<<endl;
		}
	}
	
}