最短路总结3 BellmanFord

最短路系列链接
最短路总结1 最短路问题概述与朴素dijkstra
最短路总结2 堆优化dijkstra
最短路总结3 BellmanFord
最短路总结4 SPFA及应用
最短路总结5 floyd

1.bellman-ford算法

1. Bellman-Ford算法原理简单来说，就是使用全部的边对于起点到其他n-1个点的路径进行松弛,重复n-1次。算法复杂度为O(VE)。
2. Bellman-Ford算法有什么用呢？和Dijkstra算法作用相同，可以求出图中（网络中）单源最短路径。Dijkstra算法要求所有权值不能为负数，但Bellman-Ford算法可以。
3. 当存在负权回路时，SPFA失效，只能用bellman-ford；当限制最多经过k条边时，SPFA失效，只能用bellman-ford。
4. 区别于其他最短路算法，bellman-ford算法一般按边存储.

2.算法流程

1. 存储边edge,初始化各点到源点的距离
2. 迭代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;
}