zoukankan      html  css  js  c++  java
  • vijos:P1053Easy sssp(spfa判环)

    描述

    输入数据给出一个有N(2 <= N <= 1,000)个节点,M(M <= 100,000)条边的带权有向图. 
    要求你写一个程序, 判断这个有向图中是否存在负权回路. 如果从一个点沿着某条路径出发, 又回到了自己, 而且所经过的边上的权和小于0, 就说这条路是一个负权回路.
    如果存在负权回路, 只输出一行-1;
    如果不存在负权回路, 再求出一个点S(1 <= S <= N)到每个点的最短路的长度. 约定: S到S的距离为0, 如果S与这个点不连通, 则输出NoPath.

    格式

    输入格式

    第一行: 点数N(2 <= N <= 1,000), 边数M(M <= 100,000), 源点S(1 <= S <= N);
    以下M行, 每行三个整数a, b, c表示点a, b(1 <= a, b <= N)之间连有一条边, 权值为c(-1,000,000 <= c <= 1,000,000)

    输出格式

    如果存在负权环, 只输出一行-1, 否则按以下格式输出
    共N行, 第i行描述S点到点i的最短路: 
    如果S与i不连通, 输出NoPath;
    如果i = S, 输出0;
    其他情况输出S到i的最短路的长度.

    输入
    6 8 1
    1 3 4
    1 2 6
    3 4 -7
    6 4 2
    2 4 5
    3 6 3
    4 5 1
    3 5 4
    输出
    0
    6
    4
    -3
    -2
    7
    思路:spfa算法,不过这道题存在陷阱。第一存在重边,解决方法用邻接矩阵表示图,保留权值最小的一个。第二图可能不连通,解决方法对每个结点进行spfa,为了优化时间,设置一个used数组标志已经遍历过的结点,遍历过的结点无需再进行spfa。最后加上若d[s]<0则存在负环。
    #include<cstdio>
    #include<cstring>
    #include<queue>
    using namespace std;
    const int MAXN=1005;
    const int INF=0x7f7f7f7f;
    int n,m,s;
    int contain;
    int mp[MAXN][MAXN];
    int d[MAXN];
    int cnt[MAXN];
    bool vis[MAXN],used[MAXN];
    void addedge(int u,int v,int w)
    {
        if(mp[u][v]>w)
            mp[u][v]=w;
    }
    
    bool spfa(int scource)
    {
        for(int i=1;i<=n;i++)
        {
            d[i]=INF;
            vis[i]=false;
        }
        
        queue<int> que;
        que.push(scource);
        vis[scource]=true;
        used[scource]=true;
        d[scource]=0;
        while(!que.empty())
        {
            int now = que.front();que.pop();
            vis[now]=false;
            for(int i=1;i<=n;i++)
            {
                if(now==i||mp[now][i]==INF)    continue;
                if(d[i]>d[now]+mp[now][i])
                {
                    d[i]=d[now]+mp[now][i];
                    if(d[scource]<0)    return true;
                    if(!vis[i])
                    {
                        used[i]=vis[i]=true;
                        que.push(i);
                        cnt[i]++;
                        if(cnt[i]>=n)    return true;
                    }
                }
                
            }
        }
        return false;
    }
    
    void solve()
    {
        for(int i=1;i<=n;i++)
        {
            if(!used[i]&&spfa(i))
            {
                printf("-1
    ");
                return ;
            }
        }
        if(spfa(s))
        {
            printf("-1
    ");
        }
        else
        {
            for(int i=1;i<=n;i++)
                if(d[i]==INF)    printf("NoPath
    ");
                else    printf("%d
    ",d[i]);
        }
    }
    
    int main()
    {
        scanf("%d%d%d",&n,&m,&s);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(i==j)    mp[i][j]=0;
                else    mp[i][j]=mp[j][i]=INF;
        for(int i=0;i<m;i++)
        {
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            addedge(u,v,w);
        }
        solve();
        
        return 0;
    }

     spfa+前向星可解决重边问题。d[source][source]<0则存在负环

    #include <cstdio>
    #include <cstring>
    #include <queue>
    using namespace std;
    const int MAXN=1005;
    const int INF=0x3f3f3f3f;
    struct Edge{
        int to,w,next;
    }es[100005];
    int head[MAXN],tot;
    int n,m,source;
    int d[MAXN][MAXN],vis[MAXN],used[MAXN];
    void addedge(int u,int v,int w)
    {
        es[tot].to=v;
        es[tot].w=w;
        es[tot].next=head[u];
        head[u]=tot++;
    }
    bool spfa(int s)
    {
        for(int i=1;i<=n;i++)
        {
            d[s][i]=INF;
            vis[i]=0;
        }
        queue<int> que;
        que.push(s);
        d[s][s]=0;
        vis[s]=1;
        used[s]=1;
        while(!que.empty())
        {
            int u=que.front();que.pop();
            vis[u]=0;
            for(int i=head[u];i!=-1;i=es[i].next)
            {
                Edge e=es[i];
                if(d[s][e.to]>d[s][u]+e.w)    
                {
                    d[s][e.to]=d[s][u]+e.w;
                    if(d[s][s]<0)    return true;
                    if(!vis[e.to])
                    {
                        que.push(e.to);
                        vis[e.to]=1;
                        used[e.to]=1;
                    }
                }
            }
        }
        return false;
    }
    void solve()
    {
        if(spfa(source))
        {
            printf("-1
    ");
            return ;
        }
        for(int i=1;i<=n;i++)
        {
            if(!used[i]&&spfa(i))
            {
                printf("-1
    ");
                return ;
            }
        }
        for(int i=1;i<=n;i++)
        {
            if(d[source][i]==INF)    printf("NoPath
    ");
            else    printf("%d
    ",d[source][i]);
        }
    }
    int main()
    {
        memset(head,-1,sizeof(head));
        scanf("%d%d%d",&n,&m,&source);
        for(int i=0;i<m;i++)
        {
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            addedge(u,v,w);
        }
        solve();
        
        return 0;
    }
  • 相关阅读:
    想算法无头绪,咋办?
    使用R语言分析股价波动
    cmd控制台 wrapper | OpenSCManager failed
    Ubuntu使用vi命令时,不能正常编辑文件,使用方向键时老是出现很多字母解决方案
    Ubuntu su命令 Authentication failure的解决办法
    Springboot2.x整合RabbitMQ
    RabbitMQ 笔记
    Ubuntu安装RabbitMQ
    使用xshell连接ubuntu
    Java 泛型
  • 原文地址:https://www.cnblogs.com/program-ccc/p/5339675.html
Copyright © 2011-2022 走看看