spfa求最短路与判断负环

spfa

1.基本思路

spfa是对Bellman-Ford算法的改进，Bellman-Ford的基本思路为：

for n-1:
   for 所有边(a, b, w)边 //随便存，如定义一个结构体
      dist[b] = min[dist[b], dist[a] + w]

由上式知，并非每次min运算都会使dist[b]改变，只有当dist[a]更小了，dist[b]才能变小。
因此可以开一个队列(或者其他的数据结构)来存储变小的a，当队列不空时，取出队头t，用t更新其外邻居的dist并加入队列(若队列已存在该邻居则不加入，优化时间，且dist也已更新)。

2.spfa求最短路

给定一个 n 个点 m 条边的有向图，图中可能存在重边和自环， 边权可能为负数。
请你求出 1 号点到 n 号点的最短距离，如果无法从 1 号点走到 n 号点，则输出 impossible。
数据保证不存在负权回路。
输入格式
第一行包含整数 n 和 m。
接下来 m 行每行包含三个整数 x,y,z，表示存在一条从点 x 到点 y 的有向边，边长为 z。
输出格式
输出一个整数，表示 1 号点到 n 号点的最短距离。
如果路径不存在，则输出 impossible。
数据范围
1≤n,m≤10⁵,
图中涉及边长绝对值均不超过 10000。
输入样例：

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

输出样例：

2

思路
有负权边，则考虑用Bellman-Ford或spfa，用spfa最好的时间复杂度为$O(m)$，最差为$O(nm)$与Bellman-Ford一致。该题的n、m范围较大，考虑用spfa。
代码

// spfa求最短路
#include <iostream>
#include <cstring>
#include <algorithm>
#include <queue>

using namespace std;

typedef pair<int,int> PII;

const int N = 1e5 + 10;

int n, m;
int h[N], e[N], ne[N], w[N], idx;
int dist[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 ++;
}

int spfa()
{
    memset(dist, 0x3f, sizeof dist);
    dist[1] = 0;
    
    queue<int> q;
    q.push(1);
    st[1] = true;
    
    while(q.size()) // 当队列不空
    {
        int t = q.front(); // 取出队头
        q.pop();
        
        st[t] = false; // 取出则不在队列中
        
        // 用t更新其外邻居j的dist并加入队列
        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];
                if(!st[j]) // 已在队列中就不用加入
                {
                    q.push(j);
                    st[j] = true;
                }
            }
        }
    }
    
    if(dist[n] == 0x3f3f3f3f) return -1; // 无最短路，则更新不到
    return dist[n];
}

int main()
{
    cin >> n >> m;
    
    memset(h, -1, sizeof h);
    
    while(m --)
    {
        int a, b, c;
        cin >> a >> b >> c;
        add(a, b, c);
    }
    
    int t = spfa();
    
    if(t == -1) puts("impossible");
    else printf("%d\n", t);
    
    return 0;
    
}

3.spfa判断负环

给定一个 n 个点 m 条边的有向图，图中可能存在重边和自环，边权可能为负数。
请你判断图中是否存在负权回路。
输入格式
第一行包含整数 n 和 m。
接下来 m 行每行包含三个整数 x,y,z，表示存在一条从点 x 到点 y 的有向边，边长为 z。
输出格式
如果图中存在负权回路，则输出 Yes，否则输出 No。
数据范围
1≤n≤2000,
1≤m≤10000,
图中涉及边长绝对值均不超过 10000。
输入样例：

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

输出样例：

Yes

思路
判断负环一般用spfa，在上一题的基础上加入计数数组cnt[]即可。
当dist[j] > dist[t] + w[i]，即t可更新其外邻居j的dist[j]，则：

dist[j] = dist[t] + w[i];
cnt[j] = cnt[t] + 1;

若cnt[j] >= n，则说明存在一条弧数为n的最短路，其中n+1个点必有两点是同一点，即存在环，且该环必为负环。
而该负环未必是顶点1可到达的，因此将所有点加入队列。
代码

#include <iostream>
#include <cstring>
#include <algorithm>
#include <queue>

using namespace std;

const int N = 1e5 + 10;

int n, m;
int h[N], w[N], e[N], ne[N], 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 ++;
}

int spfa()
{
    queue<int> q;
    
    for(int i = 1; i <= n; i ++) // 将所有点加入队列
    {
        st[i] = true;
        q.push(i);
    }
    
    while(q.size())
    {
        int 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])
                {
                    q.push(j);
                    st[j] = true;
                }
            }
        }
    }
    
    return false;
}

int main()
{
    cin >> n >> m;
    
    memset(h, -1, sizeof h);
    
    while(m --)
    {
        int a, b, c;
        cin >> a >> b >> c;
        add(a, b, c);
    }
    
    if(spfa()) puts("Yes");
    else puts("No");
    
    return 0;
}