sdut AOE网上的关键路径(spfa+前向星)

http://acm.sdut.edu.cn/sdutoj/showproblem.php?pid=2498&cid=1304

题目描述

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

输入

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

输出

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

示例输入

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

示例输出

18
1 2
2 5
5 7
7 9

题目分析:
由题意可以知道这是要求从起点s到终点e的最长路径,因为有10000个点,有50000条边,用spfa进行最短路,但是有一个问题就是要求路径的字典序最小。
如果正向建图的话,
比如下图:如果现在终点是6,现在2和3都能使6的距离达到最大且值相同。我们处理的时候会选2,但还是2这条路径却不是最优的,反而3是最优的。
所以我们逆向见图,求一个最短路然后倒着输出就好了。


   解决方式:

     倒序建图,当松弛时(u,v),遇到相同的情况,尽量使u变的更小,那么最终得到就是最小的字典序。

     对于求最长路径,将dis设为-INF,dis[s] = 0 

代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#define inf 0x3f3f3f3f
#define maxx 200001
using namespace std;
struct node
{
    int x,y,c,next;
} eg[maxx];
int n,m,flag,tt,pre[20002],ru[20002],ch[20002],dis[20002],v[20002],head[20002];
void init()
{
    tt=0;
    flag=0;
    memset(ru,0,sizeof(ru));
    memset(ch,0,sizeof(ch));
    memset(head,-1,sizeof(head));
}
void add(int xx,int yy,int zz)
{
    eg[tt].x=xx;
    eg[tt].y=yy;
    eg[tt].c=zz;
    eg[tt].next=head[xx];
    head[xx]=tt++;
}
void SPFA(int s,int e)
{
    int ff;
    v[s]=1;
    queue<int>q;
    while(!q.empty())
        q.pop();
    for(int i=1; i<=n; i++)
    {
        dis[i]=-inf;
        pre[i]=-1;
        v[i]=0;
    }
    dis[s]=0;
    q.push(s);
    while(!q.empty())
    {
        ff=q.front();
        q.pop();
        v[ff]=0;
        for(int i=head[ff]; i!=-1; i=eg[i].next)
        {
            int vv=eg[i].y;
            int w=eg[i].c;
            if(dis[ff]+w>dis[vv])
            {
                pre[vv]=ff;//
                dis[vv]=dis[ff]+w;
                if(v[vv]==0)
                {
                    q.push(vv);
                    v[vv]=1;
                }
            }
            else if(dis[ff]+w==dis[vv]&&pre[vv]>ff)
            {
                pre[vv]=ff;
                if(v[vv]==0)
                {
                    v[vv]=1;
                    q.push(vv);
                }
            }
        }
    }
    cout<<dis[e]<<endl;
    int T=e;
    while(T!=s)
    {
        printf("%d %d\n",T,pre[T]);
        T=pre[T];
    }
}
int main()
{
    int xx,yy,zz,RU,CH;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init();
        for(int i=0; i<m; i++)
        {
            scanf("%d%d%d",&xx,&yy,&zz);
            add(yy,xx,zz);
            ru[xx]++;
            ch[yy]++;
        }
        for(int i=1; i<=n; i++)
        {
            if(ru[i]==0)
                RU=i;
            if(ch[i]==0)
                CH=i;
        }
        SPFA(RU,CH);
    }
    return 0;
}

 

posted @ 2014-11-30 12:01  人艰不拆_zmc  阅读(641)  评论(0编辑  收藏  举报