数据结构实验之图论：AOE网上的关键路径

数据结构实验之图论十一：AOE网上的关键路径

Time Limit: 1000MS Memory Limit: 65536KB

SubmitStatistic

Problem Description

    一个无环的有向图称为无环图（Directed Acyclic Graph），简称DAG图。
    AOE(Activity On Edge)网：顾名思义，用边表示活动的网，当然它也是DAG。与AOV不同，活动都表示在了边上，如下图所示：
                                     
    如上所示，共有11项活动（11条边），9个事件（9个顶点）。整个工程只有一个开始点和一个完成点。即只有一个入度为零的点（源点）和只有一个出度为零的点（汇点）。
    关键路径：是从开始点到完成点的最长路径的长度。路径的长度是边上活动耗费的时间。如上图所示，1 到2 到 5到7到9是关键路径（关键路径不止一条，请输出字典序最小的），权值的和为18。

Input

    这里有多组数据，保证不超过10组，保证只有一个源点和汇点。输入一个顶点数n(2<=n<=10000),边数m(1<=m <=50000),接下来m行，输入起点sv，终点ev,权值w（1<=sv,ev<=n,sv != ev,1<=w <=20)。数据保证图连通。

Output

    关键路径的权值和，并且从源点输出关键路径上的路径（如果有多条，请输出字典序最小的）。

Example Input

9 11
1 2 6
1 3 4
1 4 5
2 5 1
3 5 1
4 6 2
5 7 9
5 8 7
6 8 4
8 9 4
7 9 2

Example Output

18
1 2
2 5
5 7
7 9

#include <bits/stdc++.h> 
using namespace std; 
typedef struct arcnode//表结点 
{
   int adj;//存储结点;
   arcnode *next; 
   int info;//存储权值； 
} arcnode; 
typedef struct vnode//头结点； 
{
    int data; 
    arcnode *head; 
} adjlist[10001]; 
typedef struct//图的结构 
{ 
   adjlist a; 
   int vexn,arcn; // 结点数 边数 
} MG; 
   int i,j; 
   int indegree[10001];//记录每个点的入度 
   int ve[10001];//最早开始时间； 
   int vl[10001];//最迟开始时间； 
void create(MG &G)//建立有向图邻接表； 
{ 
    int v1,v2,w; 
	arcnode *p;
	for(i=1; i<=G.vexn; i++)  
	G.a[i].head=NULL;//头结点清空；
	for(i=1; i<=G.arcn; i++) {
	  cin >> v1 >> v2 >> w; 
	  p=new arcnode;  
	  p->adj=v2;  
	  p->info=w; 
	  indegree[v2]++; 
	  p->next=G.a[v1].head; 
	  G.a[v1].head=p; 
	 } 
}
 stack<int>s;//用于逆拓扑有序时求vl； 
 int topo1(MG &G)//求ve（最早开始时间即最大路径长度）； 
 { 
    int k; 
	arcnode *p; 
	stack<int>t; 
	for(i=1; i<=G.vexn; i++) 
	if(!indegree[i])//入度为零的结点入栈
	      t.push(i); 
	memset(ve,0,sizeof(ve));//初始化最早开始时间； 
	int count=0;//记录出栈元素个数，判定拓扑排序是否有序； 
	while(!t.empty())
	{ 
	    j=t.top(); 
		t.pop(); 
		s.push(j); 
		count++;  
		for(p=G.a[j].head; p; p=p->next) 
		{ 
		    k=p->adj; 
			if(!(--indegree[k]))//将新入度为零的节点进栈； 
			   t.push(k); 
			if(ve[j]+p->info>ve[k])//更新最大路径长度； 
			   ve[k]=ve[j]+p->info; 
		} 
	} 
	if(count<G.vexn)//拓扑不是有序的，即有环； 
	     return 0; 
	return 1; 
} 
int topo2(MG &G)//求vL（既不影响施工进度的情况下最晚的开始时间； 
{ 
    int k; 
	arcnode *p; 
	if(!topo1(G)) return 0;//拓扑无序，则结束； 
	cout <<  ve[G.vexn] << endl;//关键路径的权值和； 
	for(i=1; i<=G.vexn; i++) //初始化最迟开始时间； 
	vl[i]=ve[G.vexn]; 
	while(!s.empty())//按拓扑逆序求vl； 
	{ 
	    j=s.top(); 
		s.pop(); 
		for(p=G.a[j].head; p; p=p->next)
		{ 
		    k=p->adj;  
			if(vl[k]-p->info<vl[j])//更新最晚开始时间； 
			   vl[j]=vl[k]-p->info; 
		} 
	} 
	int flag=0,a,b;
	for(j=1; j<=G.vexn; j++) //求关键活动； 
	  for(p=G.a[j].head; p; p=p->next) 
	   { 
	       k=p->adj; 
		   int e=ve[j]; 
		   int l=vl[k]-p->info; 
		   if(e==l)//关键路径要求最早开始时间和最晚开始时间相同没有空余时间； 
		   { 
		      if(flag==0)//对输出的特殊处理，要求字典序最小 
			   { 
			       a=j; 
				   b=k; 
				   flag=1;  
				} 
			  else if(a==j&&b>k) 
			     b=k; 
			  else if(b==j)//结束活动的分支时b作为活动的开始时。 
			   { 
			      printf("%d %d\n",a,b); 
				  a=j; 
				  b=k; 
				} 
			}  
		} 
		cout << a << " " << b << endl; 
		return 1; 
} 
int main() { 
     MG G; 
	 while(cin >> G.vexn >> G.arcn) 
	 { 
	     memset(indegree,0,sizeof(indegree)); 
		 create(G);
		 topo2(G); 
	}
      return 0; 
}