SPFA算法

目录

应用场景：

问题描述

思路分析：

实现代码：

---

应用场景：

SPFA算法应用于解决单源权值为负数且无负环的最短路问题，对于一般的正值权重问题也可以使用spfa算法来代替dijkstra算法。

核心思想是通过队列优化bellman-ford算法，在bellman-ford算法中需要遍历所有的边，但是对于所有的边，如果前一个点的边到源点的距离没有更改，那么此点对于后面的点的最短路问题是没有帮助的。简而言之，spfa算法是我更新过谁，再拿谁来更新后面的点。

---

问题描述

给定一个 n 个点 m 条边的有向图，图中可能存在重边和自环， 边权可能为负数。

请你求出 1 号点到 n 号点的最短距离，如果无法从 1 号点走到 n 号点，则输出 impossible。

数据保证不存在负权回路。

输入格式

第一行包含整数 n 和 m。

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

输出格式

输出一个整数，表示 1 号点到 n 号点的最短距离。

如果路径不存在，则输出 impossible。

数据范围

1≤n,m≤105,
图中涉及边长绝对值均不超过 10000。

输入样例：

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

输出样例：

2

---

思路分析：

1.spfa用队列存储 简化bellman-ford算法中遍历所有边的过程 spfa只会将前一个已经更新的点用于更新下面的边
2.队列中存储的为已经更新过的点 每次将此类点入队出队 再用此点去更新下面的点到源点的距离
3.spfa算法与dijkstra的区别在于1.st数组的不同含义 2.dijkstra基于贪心的思想 每个点只会使用一次 在spfa中每个点可能多次入队 不是单调的 在spfa算法中 st数组是表示当前节点是否在队列中

---

实现代码：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N = 100010;
int dist[N],w[N];
int idx,h[N],e[N],ne[N];
bool st[N];
int n,m;

//邻接表存储的插入操作 带有权值边的插入
void add(int a,int b,int c){
    e[idx] = b;
    w[idx] = c;
    ne[idx] = h[a];
    h[a] = idx++;
}

//spfa的核心代码
//spfa用队列存储 简化bellman-ford算法中遍历所有边的过程 spfa只会将前一个已经更新的点用于更新下面的边
//队列中存储的为已经更新过的点 每次将此类点入队出队 再用此点去更新下面的点到源点的距离
//spfa算法与dijkstra的区别在于1.st数组的不同含义 2.dijkstra基于贪心的思想 每个点只会使用一次
//在spfa中每个点可能多次入队 不是单调的
//在spfa算法中 st数组是表示当前节点是否在队列中
void spfa(){
    //初始化dist数组
    memset(dist,0x3f,sizeof(dist));
    dist[1] = 0;
    //定义队列
    queue<int> q;
    q.push(1);
    //当对列不为空时
    while(q.size()){
        //取出节点
        auto t = q.front();
        q.pop();
        //从队列中取出来之后该节点st被标记为false,代表之后该节点如果发生更新可再次入队
        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];
                //当前已经加入队列的结点，无需再次加入队列，即便发生了更新也只用更新数值即可，重复添加降低效率
                if(!st[j]){
                    q.push(j);
                    st[j] = true;
                }
            }
        }
    }
}
int main(){
    //初始化
    memset(h,-1,sizeof(h));
    //输入操作
    cin>>n>>m;
    for(int i=0;i<m;i++){
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }
    spfa();
    //spfa算法只会更新能走到源点的点 此时判断条件为0x3f3f3f3f
    if(dist[n]==0x3f3f3f3f) cout<<"impossible"<<endl;
    else cout<<dist[n]<<endl;
    return 0;
}