算法学习系列(三十九):spfa算法判断负环

引言

关于这个算法学习系列,确实大部分的知识点基本上都讲完了,就是讲的非常的宏观,比如说图论的知识,基础的是讲完了,但是更高阶的当然还没写,也还没学,并且目前处于刷题量还不是很大的情况下,还是不要盲目的去学习更难的知识,所以以后的算法学习系列,打算以零碎知识点的方式来写文章,毕竟也写不了像之前那么多的内容了,也太难了。言归正传,今天学习的是如何使用 s p f a spfa spfa 算法判断一个有向图中是否存在负环。


一、题目描述

概括:给一个有向图,存在重边和自环,判断该图中是否存在负环。
给定一个 n 个点 m 条边的有向图,图中可能存在重边和自环,边权可能为负数。 给定一个 n 个点 m 条边的有向图,图中可能存在重边和自环, 边权可能为负数。 给定一个n个点m条边的有向图,图中可能存在重边和自环,边权可能为负数。

请你判断图中是否存在负权回路。 请你判断图中是否存在负权回路。 请你判断图中是否存在负权回路。

输入格式 输入格式 输入格式
第一行包含整数 n 和 m 。 第一行包含整数 n 和 m。 第一行包含整数nm

接下来 m 行每行包含三个整数 x , y , z ,表示存在一条从点 x 到点 y 的有向边,边长为 z 。 接下来 m 行每行包含三个整数 x,y,z,表示存在一条从点 x 到点 y 的有向边,边长为 z。 接下来m行每行包含三个整数x,y,z,表示存在一条从点x到点y的有向边,边长为z

输出格式 输出格式 输出格式
如果图中存在负权回路,则输出 Y e s ,否则输出 N o 。 如果图中存在负权回路,则输出 Yes,否则输出 No。 如果图中存在负权回路,则输出Yes,否则输出No

数据范围 数据范围 数据范围
1 ≤ n ≤ 2000 , 1 ≤ m ≤ 10000 , 图中涉及边长绝对值均不超过 10000 。 1≤n≤2000,1≤m≤10000,图中涉及边长绝对值均不超过 10000。 1n2000,1m10000,图中涉及边长绝对值均不超过10000

输入样例: 输入样例: 输入样例:

3 3
1 2 -1
2 3 4
3 1 -4

输出样例: 输出样例: 输出样例:

Yes

二、解题思路

思路:我们可以维护一个 c n t [ i ] cnt[i] cnt[i] 数组,表示从某个点走到 i i i 号点的边数,使得 c n t [ i ] = c n t [ t ] + 1 cnt[i] = cnt[t] + 1 cnt[i]=cnt[t]+1 只要 c n t [ i ] ≥ n cnt[i] \geq n cnt[i]n,即存在负环,因为总共有 n n n 个点,也就是最多存在 n − 1 n-1 n1 条边。又因为 s p f a spfa spfa 算法求的是最短路,所以只有负环才会被识别出来。并且要注意:我们的起始点要从所有的点开始,因为1号点可能到不了负环,负自环就到不了。


三、示例代码

spfa算法参考博客:spfa算法

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

const int N = 2010, M = 1e4+10;

int n, m;
int h[N], e[M], w[M], ne[M], idx;
int dist[N], cnt[N];
bool st[N];

void add(int a, int b, int c)
{
    e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx++;
}

bool spfa()
{
    memset(dist, 0x3f, sizeof dist);
    dist[1] = 0;

    queue<int> q;
    for(int i = 1; i <= n; ++i)
    {
        st[i] = true;
        q.push(i);
    }
    while(q.size())
    {
        auto t = q.front(); q.pop();
        st[t] = false;
        
        for(int i = h[t]; i != -1; i = ne[i])
        {
            int j = e[i];
            if(dist[j] > dist[t] + w[i])
            {
                dist[j] = dist[t] + w[i];
                cnt[j] = cnt[t] + 1;
                if(cnt[j] >= n) return true;
                if(!st[j])
                {
                    st[j] = true;
                    q.push(j);
                }
            }
        }
    }
    
    return false;
}

int main()
{
    ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    memset(h, -1, sizeof h);
    
    
    cin >> n >> m;
    while(m--)
    {
        int a, b, c;
        cin >> a >> b >> c;
        add(a,b,c);
    }
    
    
    if(spfa()) puts("Yes");
    else puts("No");
    
    return 0;
}
  • 25
    点赞
  • 15
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

lijiachang030718

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值