AOE网求解关键路径完整代码

泪奔，图终于看到最后啦~
AOE网手动实现比较复杂，学理论的时候也总是记不清（记过就忘），现在回头再看看代码，发现实现起来还是比较方便的~
把书上给的核心代码补齐了，注意以下几点：
1、用栈存拓扑排序序列，出栈的时候就是逆拓扑排序
2、vl初始化的赋值：fill(vl,vl+n,ve[n-1]); （不太理解）
3、e == l 判断是不是关键路径
这一部分强烈建议结合理论来看，理解AOE网模拟求关键路径过程，这个算法也就理解了

输入样例：
第一行为顶点数、边数
7 10
第2-2+边数 行为有向图边的弧头顶点、弧尾顶点、边权
0 1 3
0 2 2
0 3 6
1 3 2
1 4 8
2 3 1
2 5 3
3 4 1
4 6 3
5 6 4
最后一行为每个顶点的入度
0 1 1 3 2 1 2

代码：

#include <stdio.h>
#include <string.h>
#include <queue>
#include <stack>
#include <vector>
#include <algorithm>
using namespace std;

const int maxn = 100;
struct node{
	int v,w;    //记录（有向边的下一个）顶点，边权 
};

//拓扑序列，放在栈中，出栈就是逆拓扑排序 
stack<int>toporder;
vector<node>G[maxn];   //邻接表存储的图 
int n,m,ve[maxn],vl[maxn],indegree[maxn];  //n为顶点数,m为边数 

bool topologicalsort(){
	queue<int>q;
	for(int i = 0;i <n;i++){
		if(indegree[i] == 0){
			q.push(i);
		}
	}
	while(!q.empty()){
		int u = q.front();
		printf("%d ",u);
		q.pop();
		toporder.push(u);    //将u加入拓扑序列 
		for(int i = 0;i <G[u].size();i++){
			int v = G[u][i].v;
			indegree[v]--;
			if(indegree[v] == 0){
				q.push(v);
			}
			if(ve[u] + G[u][i].w >ve[v]){
				ve[v] = ve[u] + G[u][i].w;
			}
		}
	}
	//printf("%d ",toporder.size());
	if(toporder.size() == n) return true;
	else return false;
}

int criticalpath(){
	memset(ve,0,sizeof(ve));
	if(topologicalsort() == false){
		return -1;
	}
	printf("\n");
	fill(vl,vl+n,ve[n-1]);    //vl数组初始化，初始值为汇点的ve值
	while(!toporder.empty()){      //出栈，逆拓扑排序 
		int u = toporder.top();
		toporder.pop();
		for(int i = 0;i <G[u].size();i++){
			int v = G[u][i].v;
			if(vl[v] - G[u][i].w<vl[u]){
				vl[u] = vl[v] - G[u][i].w;
			}
		}
	}
	//遍历邻接表的所有边，计算活动的最早开始时间e和最迟开始时间l
	for(int u = 0;u <n;u++){
		for(int i = 0;i <G[u].size();i++){
			int v = G[u][i].v,w = G[u][i].w;
			int e = ve[u],l = vl[v] - w;
			if(e == l){
				printf("%d->%d\n",u,v);
			}
		}
	}
	return ve[n-1];    //返回关键路径长度 
}

int main(){
	scanf("%d %d",&n,&m);
	int s,t,weight;
	for(int i = 0;i <m;i++){
		scanf("%d %d %d",&s,&t,&weight);
		node pt;   //用一个临时节点先存数据 
		pt.v = t;
		pt.w = weight;
		G[s].push_back(pt);   //邻接表的存法注意下 
	}
	for(int i = 0;i <n;i++){
		scanf("%d",&indegree[i]);    //记录入度 
	}
	//topologicalsort();
	printf("拓扑排序序列为:\n"); 
	int ans = criticalpath();
	printf("%d\n",ans);
} 

/*
9 12
0 1 2
0 2 1
1 3 3
2 4 4
3 6 2
3 5 2
4 5 1
4 7 5
5 6 3
5 7 1
6 8 3
7 8 6
0 1 1 1 1 2 2 2 2

*/

/*
7 10
0 1 3
0 2 2
0 3 6
1 3 2
1 4 8
2 3 1
2 5 3
3 4 1
4 6 3
5 6 4
0 1 1 3 2 1 2

*/