C++ SPFA算法解析

news/2024/11/18 14:29:28/文章来源:https://www.cnblogs.com/Mono-Awen/p/18374186

前言

将了解 C++ 求最短路中 SPFA 的算法

SPFA

SPFA的一些说明

SPFA：适用于权值有负值，且没有负圈的图的单源最短路径，论文中的复杂度O(kE)，k为每个节点进入Queue的次数，且k一般<=2，但此处的复杂度证明是有问题的，其实SPFA的最坏情况应该是O(VE).!

引例：

输入格式

给出一个有向图，请输出从某一点出发到所有点的最短路径长度
三个整数 n， m， s，分别表示点的个数、有向边的个数、出发点的编号。
接下来 m 行, 包含三个整数 u，v， w， u --> v， 长度为 w;

输出格式

输出一行 n 个整数，第 i 个表示 s 到第 i 个点的最短路径，若不能到达则输出 2³¹ - 1;

输入样例	输出样例
4 6 1	0 2 4 3
1 2 2
2 3 2
2 4 1
1 3 5
3 4 3
1 4 4

此题的输入 数据范围 为 70% 即可通过

我们需要一个数组来记录 一个点 到 某个点 的最短距离，我们可以定义一个 cost数组 来记录权值

写好输入数据

我们的最大值可定义为 2e6 + 9 (即2000009), 不能到达则输出 2³¹ - 1 (即2147483647）将它定义为 INF;

const int N = 2e6 + 9; // const常量，不可改变,数尽量大一些
const int INF = 2147483647;
int n, m, s, u, v, w;void AddEdge() {} // 连边函数signed main() {cin >> n >> m >> s;for (int i = 1; i <= m; ++ i) {cin >> u >> v >> w;AddEdge(u, v, w); // 待会进行的连边所使用的函数}
}

连边

我们需要定义一个 结构体 Edge 然后进行连边操作

int cnt = 0, head[N];struct Edge {int nxt, to, val; 
} edge[N];void AddEdge(int from, int to, int val) {cnt++; // 记录操作次数edge[cnt].nxt = head[from];edge[cnt].to = to;edge[cnt].val = val;head[from] = cnt;    
}

SPFA函数的编写

主程序 main（续写）中连边后执行 SPFA 函数

signed main() {……for (int i = 1; i <= m; ++ i) {cin >> u >> v >> w;AddEdge(u, v, w);}SPFA();
}

vis 记录节点是否 进/出 队列，cost 记录节点从 s 到 某个节点 的总共权值

int vis[N], cost[N];

在 SPFA 函数中，我们要进行 vis 数组的清零，需用到 memset 函数来执行，让 cost[n] 里的各数值为 INF，即作为无穷大的数，当不可到达时，可直接输出 cost[i] (i为各个点），可到达时，使用 min（） 函数可求最短路径，将 cost[s] (即起点）的数值为 0，因为它到它自己本身的距离为 0

void SPFA() {memset(vis, 0, sizeof(vis));for (int i = 1; i <= n; ++ i)cost[i] = INF;cost[s] = 0;
}

定义一个 队列 q，每次取队头执行 松弛操作，添加起点

void SPFA() {……cost[s] = 0;queue <int> q;q.push(s);vis[s] = 1；// 让起点 s 标记为 1， 代表已进入队列
}

当队列不为空时，进行循环

void SPFA() {……vis[s] = 1;while (!q.empty()) {int x = q.front(); // 获取队头数字q.pop(); // 弹出队头vis[x] = 0; // 队头已取出，此时 取出的队头 数字 的 vis[队头数字] 改为 0 代表 出队}
}

以下代码段为松弛操作（以便不懂的小伙伴）

if (cost[y] > cost[x] + edge[i].val) {cost[y] = cost[x] + edge[i].val;if (vis[y] == 0) {vis[y] = 1;q.push(y);}
}

队头已取出，此时 取出的队头数字 的 vis[队头数字] 改为 0 代表出队

void SPFA() {……while (!q.empty()) {……vis[x] = 0;for (int i = head[x]; i; i = edge[i].nxt) {int y = edge[i].to; // 代表为 i 时的 edge的 下一个指向if (cost[y] > cost[x] + edge[i].val) { // 如果 指向（未更改 原本存储） 的权值 大于 x 的权值 + 指向的权值cost[y] = cost[x] + edge[i].val; //进行 权值 的 更改if (vis[y] == 0) { //判断是否曾进入队列vis[y] = 1;q.push(y);}}}}
}

最后在主程序中进行权值的输出

signed main() {……SPFA();for (int i = 1; i <= n; ++ i ) {cout << cost[i];}
}

最终代码

#include <bits/stdc++.h>
using namespace std;
const int N = 2e6 + 9; // const常量，不可改变
const int INF = 2147483647;
int n, m, s, u, v, w;
int cnt = 0, head[N], vis[N], cost[N];struct Edge {int nxt, to, val;
} edge[N];
void AddEdge(int from, int to, int val) {cnt++;edge[cnt].nxt = head[from];edge[cnt].to = to;edge[cnt].val = val;head[from] = cnt;
}
void SPFA() {memset(vis, 0, sizeof(0));for (int i = 1; i <= n; ++ i)cost[i] = INF;cost[s] = 0;queue<int> q;q.push(s);vis[s] = 1;while(!q.empty()) {int x = q.front();q.pop();vis[x] = 0;for (int i = head[x]; i; i = edge[i].nxt) {int y = edge[i].to;if (cost[y] > cost[x] + edge[i].val) {cost[y] = cost[x] + edge[i].val;if (vis[y] == 0) {vis[y] = 1;q.push(y);}}}}
}
signed main() {cin >> n >> m >> s;for (int i = 1; i <= m; ++ i) {cin >> u >> v >> w;AddEdge(u, v, w);}SPFA();for (int i = 1; i <= n; ++ i) {cout << cost[i] << " ";}
}

SPFA的负环判断

洛谷负环判断模板题目：传送门

我将上述代码的变量 cnt 循环次数改为 tim，创建了一个 cnt数组, 用来储存 每个数的入队次数

存储的是入队次数而不是松弛次数

int cnt[N];

题目部分要求

若 u >= 0，则表示存在一条从 u 至υ边权为 w 的边，还存在一条从υ至u 边权为 w 的边。
若 u < 0，则只表示存在一条从 u 至υ边权为 w 的边。

主程序 main 中的部分代码段应改为

for (int i = 1; i <= m; ++ i) {cin >> u >> v >> w;AddEdge(u, v, w);if (w >= 0) AddEdge(v, u, w);
}

最终程序

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 3e6 + 9;
const int INF = 3e6 + 10;
int T;
int n, m, u, v, w, tim, head[N], vis[N], cost[N], cnt[N]; //cnt 换成 tim, 且创建 cnt 数组
struct Edge
{int nxt, to, val;
} edge[N];
// 清空函数: start
void Clear() {memset(edge, 0, sizeof(edge));tim = 0;memset(head, 0, sizeof(head));memset(vis, 0, sizeof(vis));memset(cost, INF, sizeof(cost));memset(cnt, 0, sizeof(cnt));
}
//清空函数: end
void AddEdge(int from, int to, int val) {tim++; //原本的cnt需改成tim，含义一致edge[tim].nxt = head[from];edge[tim].to = to;edge[tim].val = val;head[from] = tim;
}
bool SPFA() {memset(vis,0, sizeof(vis));for (int i = 1; i <= n; ++ i) cost[i] = INF;cost[1] = 0;queue<int> q;q.push(1); vis[1] = 1;while (!q.empty()) {int x = q.front();q.pop();vis[x] = false;for (int i = head[x]; i; i = edge[i].nxt) {int v = edge[i].to;if (cost[v] > cost[x] + edge[i].val) {cost[v] = cost[x] + edge[i].val;if (!vis[v]) {vis[v] = 1;q.push(v);cnt[v] ++; // 记录 edge[i].to 的次数if (cnt[v] >= n) return true; // 如果 入队次数 已经 大于等于 n}}}}return false;
}signed main() {ios::sync_with_stdio(false); cin.tie(0);cin >> T;while (T--) { //重复循环次数Clear(); //进行清空cin >> n >> m;for (int i = 1; i <= m; ++ i) {cin >> u >> v >> w;AddEdge(u, v, w);if (w >= 0) AddEdge(v, u, w); //更改的地方}cout << (SPFA() ? "YES" : "NO") << endl; //更改的地方}
}