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

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

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

Sample

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

Output 

18
1 2
2 5
5 7
7 9

关键路径问题就是绝对不允许延误活动组成的路径。

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#define MAXVEX 50001
 
typedef struct vertex
{
    int pre;
    int rear;
    int weight;
    int out_degree;
    int in_degree;
}Vertex;
 
Vertex vex[MAXVEX];
 
void initvex(int n)
{
    for(int i=1;i<=n;i++){
        vex[i].pre = 0;
        vex[i].rear = 0;
        vex[i].weight = 0;
        vex[i].in_degree = 0;
        vex[i].out_degree = 0;
    }
}
 
void AOE(int n,int m,int pos)
{
    int temp;
    int path[MAXVEX],cost[MAXVEX];
    memset(path,0,sizeof(path));
    memset(cost,0,sizeof(cost));
    for(int i=2;i<=n;i++){
        temp = 0;
        for(int j=1;j<=m;j++){
            if(cost[vex[j].pre]<cost[vex[j].rear]+vex[j].weight ||
(cost[vex[j].pre]==cost[vex[j].rear]+vex[j].weight && vex[j].rear<path[vex[j].pre])){
                cost[vex[j].pre] = cost[vex[j].rear] + vex[j].weight;
                path[vex[j].pre] = vex[j].rear;
                temp = 1;
             }
        }
        if(!temp){
            break;
        }
    }
    printf("%d\n",cost[pos]);
    while(path[pos]){
        printf("%d %d\n",pos,path[pos]);
        pos = path[pos];
    }
}
 
int main()
{
    int n,m,a,b,c,pos;
    while (scanf("%d%d",&n,&m)!=EOF){
        initvex(n);
        for(int i=1;i<=m;i++){
            scanf("%d%d%d",&a,&b,&c);
            vex[i].pre = a;
            vex[i].rear = b;
            vex[i].weight = c;
            vex[b].out_degree++;
            vex[a].in_degree++;
        }
        for(int i=1;i<=n;i++){
            if(vex[i].out_degree==0){
                pos = i;
            }
        }
        AOE(n,m,pos);
    }
    return 0;
}