1、题目
有向图G有n个点、m条边,图中可以有重边、负权重边、自环,可能有负权回路,计算1号点到n号点的最短距离。
2、思路
for i:0~n-1{
for j:0~m-1
dis[b] = min(dis[b], backup[a] + w);
}
-
外层:for循环n次
内层:for循环m次,循环所有边a—w--->b
每次:dis[b]=min(dis[b],dis[a]+w)
(松弛操作)
(外层循环k次表示的含义:从1号点,经过不超过k条边的最短路的距离;第n次循环有更新:图中一定有负环) -
循环完的结果:
所有边满足dis[b]<=dis[a]+w (三角不等式) -
有负权回路:
1到n点的最短路不一定存在。
3、复杂度
O(n*m) time
4、实现
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 510, M = 10010;//点数和边数最大值
int n, m, k;
int dis[N], backup[N];//dis[i]:点1到i最短距离
struct Edge{
int a, b, w;
}edges[M];
int bellman_ford(){
memset(dis, 0x3f, sizeof(dis));
dis[1] = 0;
for (int i = 0; i < k; i++){
//特殊处理:每次循环时,先把dis数组备份到backup中
memcpy(backup, dis, sizeof(dis));
for (int j = 0; j < m; j++){
int a = edges[j].a, b = edges[j].b, w = edges[j].w;
dis[b] = min(dis[b], backup[a] + w);
}
}
if (dis[n]>0x3f3f3f3f / 2) //写dis[n]==0x3f3f3f3f不行,因为可能有负权把0x3f3f3f3f变小了
return -1;
return dis[n];
}
int main(){
ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
//输入
cin >> n >> m >> k;
for (int i = 0; i < m; i++){
int a, b, w; cin >> a >> b >> w;
edges[i] = { a, b, w };
}
//计算bellman-ford
int t = bellman_ford();
if (t == -1)
cout << "impossible" << endl;
else
cout << t << endl;
return 0;
}