P1053Easy sssp

描述

输入数据给出一个有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的最短路的长度.

样例输入

1
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

坑人的题目,结果连WA了18次9T_T),spfa判环(巨坑)
两次spfa第一次跑从s点与其他点的距离,第二次跑整个图,
判环。

#include<cstdio>
#include<queue>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define M 1010
#define N 100010
#define ll long long
using namespace std;
ll next[N],to[N],val[N],head[N],tot;
ll vis[N],dis[N],in[N],dis2[N];
int n,m,s;
int read() {
    int x=0,f=1;char c=getchar();
    while(c>'9'||c<'0') {if(c=='-')f=-1;c=getchar();}
    while(c<='9'&&c>='0') x=10*x+c-48,c=getchar();
    return x*f;
}
void add(int x,int y,int z) {
    to[++tot]=y;
    val[tot]=z;
    next[tot]=head[x];
    head[x]=tot;
}
void spfa(int s) {
    queue<int> q;
    q.push(s);
    dis[s]=0;vis[s]=1;in[s]++;
    while(!q.empty()) {
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u];i;i=next[i]) {
            int v=to[i];
            if(dis[v]>dis[u]+val[i]) {
                dis[v]=dis[u]+val[i];
                if(!vis[v]) {
                    q.push(v);
                    vis[v]=1;
                    in[v]++;
                    if(in[v]>n) {
                        printf("-1\n");
                        exit(0);
                    }
                }
            }
        }
    }
}
void spfa1(int s){
    queue<int>q;
    memset(in,0,sizeof(in));
    memset(vis,false,sizeof(vis));
    q.push(s);
    vis[s]=1;dis2[s]=0;in[s]++;
    while(!q.empty()) {
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u];i;i=next[i]) {
            int v=to[i];
            if(dis2[v]>dis2[u]+val[i]) {
                if(!vis[v]) {
                    q.push(v);
                    vis[v]=1;
                    in[v]++;
                    if(in[v]>n) {
                        printf("-1\n");
                        exit(0);
                    }
                }
            }
        }                   
    }
}
int main() {
    n=read();m=read();s=read();
    for(int i=1;i<=m;i++) {
        int x,y,z;
        x=read();y=read();z=read();
        add(x,y,z);
    }
    for(int i=0;i<=n;i++) dis[i]=dis2[i]=M*M;
    spfa(s);
    for(int i=1;i<=n;i++)
      if(dis[i]==dis[0])
        spfa1(i);
    for(int i=1;i<=n;i++) {
        if(dis[i]!=M*M) printf("%lld\n",dis[i]);
        if(dis[i]==M*M) printf("NoPath\n");
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值