最短路系列链接
最短路总结1 最短路问题概述与朴素dijkstra
最短路总结2 堆优化dijkstra
最短路总结3 BellmanFord
最短路总结4 SPFA及应用
最短路总结5 floyd
1.bellman-ford算法
- Bellman-Ford算法原理简单来说,就是使用全部的边对于起点到其他n-1个点的路径进行松弛,重复n-1次。算法复杂度为O(VE)。
- Bellman-Ford算法有什么用呢?和Dijkstra算法作用相同,可以求出图中(网络中)单源最短路径。Dijkstra算法要求所有权值不能为负数,但Bellman-Ford算法可以。
- 当存在负权回路时,SPFA失效,只能用bellman-ford;当限制最多经过k条边时,SPFA失效,只能用bellman-ford。
- 区别于其他最短路算法,bellman-ford算法一般按边存储.
2.算法流程
- 存储边edge,初始化各点到源点的距离
- 迭代n次
循环所有边(a,b,w),对所有点做松弛操作dist[b]=min(dist[b],dist[a]+w);
3.模板
int n, m; // n表示点数,m表示边数
int dist[N]; // dist[x]存储1到x的最短路距离
struct Edge // 边,a表示出点,b表示入点,w表示边的权重
{
int a, b, w;
}edges[M];
// 求1到n的最短路距离,如果无法从1走到n,则返回-1。
int bellman_ford()
{
memset(dist, 0x3f, sizeof dist);
dist[1] = 0;
// 如果第n次迭代仍然会松弛三角不等式,就说明存在一条长度是n+1的最短路径,由抽屉原理,路径中至少存在两个相同的点,说明图中存在负权回路。
for (int i = 0; i < n; i ++ )
{
for (int j = 0; j < m; j ++ )
{
int a = edges[j].a, b = edges[j].b, w = edges[j].w;
if (dist[b] > dist[a] + w)
dist[b] = dist[a] + w;
}
}
if (dist[n] > 0x3f3f3f3f / 2) return -1;
return dist[n];
}
4.板子题
给定一个n个点m条边的有向图,图中可能存在重边和自环, 边权可能为负数。
请你求出从1号点到n号点的最多经过k条边的最短距离,如果无法从1号点走到n号点,输出impossible。
注意:图中可能 存在负权回路 。
输入格式
第一行包含三个整数n,m,k。
接下来m行,每行包含三个整数x,y,z,表示存在一条从点x到点y的有向边,边长为z。
输出格式
输出一个整数,表示从1号点到n号点的最多经过k条边的最短距离。
如果不存在满足条件的路径,则输出“impossible”。
数据范围
1≤n,k≤500,
1≤m≤10000,
任意边长的绝对值不超过10000。
输入样例:
3 3 1
1 2 1
2 3 1
1 3 3
输出样例:
3
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=550,M=10010,inf=0x3f3f3f3f;
struct Edge{
int a,b,w;
}edges[M];
int n,m,k;
int dist[N];
int backup[N];
int ford() {
memset(dist,inf,sizeof(dist));
dist[1]=0;
for(int i=1;i<=k;i++) {
memcpy(backup,dist,sizeof(dist));
for(int j=1;j<=m;j++) {
int a=edges[j].a;
int b=edges[j].b;
int w=edges[j].w;
dist[b]=min(dist[b],backup[a]+w);
}
}
if(dist[n]>k*10000) return -1;
else return dist[n];
}
int main(){
cin >> n >> m >> k;
// memset(g,inf,sizeof(g));
for(int i=1;i<=m;) {
int x,y,z;
cin >> x >> y >> z;
if(x!=y) {
edges[i].a=x,edges[i].b=y,edges[i].w=z;
i++;
}
}
int t=ford();
// for(int i=1;i<=n;i++) cout<<dist[i]<<" ";
if(t==-1) cout<<"impossible";
else cout<<t;
return 0;
}