AOE网上的关键路径

AOE网上的关键路径

Time Limit: 1000MS Memory Limit: 65536KB
Submit Statistic Discuss
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
Hint

Author

AC Code：

#include <stdio.h>
#include <stdlib.h>
#include <stack>
#include <string.h>
using namespace std;
typedef struct arcnode//表结点
{
    int adj; //存储节点
    struct arcnode *next;
    int info;//存储权值
} arcnode;
typedef struct vnode//头结点
{
    int data;
    arcnode *first;
} adjlist[10010];
typedef struct
{
    adjlist a;
    int vn, an;
} ALG;
int n, m;
int indegree[10010];//记录每个点的入度
int ve[10010]; //最早开始时间
int vl[10010]; //最迟开始时间
void create(ALG &g)
{
    int i, j, v1, v2, w;
    arcnode *p;
    for(i=1; i<=n; i++)
    {
        g.a[i].first = NULL; //头结点清空
    }
    for(i=1; i<=m; i++)
    {
        scanf("%d %d %d", &v1, &v2, &w);
        p = new arcnode;
        p->adj = v2;
        p->info = w;
        indegree[v2]++;
        p->next = g.a[v1].first;
        g.a[v1].first = p;
    }
}
stack <int> t; //用于逆拓扑有序时求vl
int topo1(ALG &g)//求 ve（最早开始时间即最大路径长度）
{


    int i, j, k;
    arcnode *p;
    stack <int> s;
    for(i=1; i<=n; i++)//入度为0的结点入栈
        if(!indegree[i])
            s.push(i);
    memset(ve, 0, sizeof(ve));
    int count = 0;//记录出栈顶点个数
    while(!s.empty())
    {
        j = s.top();
        s.pop();
        t.push(j);
        count++;
        for(p=g.a[j].first; p; p=p->next) //删去所有已j为起点的出边
        {
            k = p->adj;
            if(!(--indegree[k]))//将新入度为0的结点入栈
                s.push(k);
            if(ve[j]+p->info>ve[k]) //更新最大路径长度
                ve[k] = ve[j]+p->info;
        }
    }
    if(count<n)
        return 0;
    else
        return 1;
}


int topo2(ALG &g)//求vl（即不影响施工进度的情况下最晚开始的时间）
{


    int i, j, k;
    arcnode *p;
    if(!topo1(g)) return 0;
    printf("%d\n", ve[n]);
    for(i=1; i<=n; i++)  //初始化最迟发生时间
        vl[i] = ve[n];


    while(!t.empty()) //按拓扑逆序求vl
    {
        j = t.top();
        t.pop();
        for(p = g.a[j].first; 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<=n; j++) //求关键活动
    {
        for(p = g.a[j].first; 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)
                {
                    printf("%d %d\n", a, b);
                    a = j;
                    b = k;
                }
            }
        }
    }
    printf("%d %d\n", a, b);
}


int main()
{
    ALG g;
    while(~scanf("%d %d", &n, &m))
    {
        memset(indegree, 0, sizeof(indegree));
        create(g);

        topo2(g);

    }
    return 0;
}


/***************************************************
User name: 
Result: Accepted
Take time: 52ms
Take Memory: 2972KB
Submit time: 2017-07-25 19:33:52
****************************************************/