【图论】Floyd算法求任意两点间最短路

一、前言

对于最短路问题，我们已经介绍了很多的算法来求解。其中每一种算法都有其不可替代性，分别适用于不同的情况。而在这篇博客里我们需要讲解的是求任意两点最短路的一个很简单的算法：Floyd算法。

二、算法概念介绍

Floyd算法其实是用到了动态规划的方法去解决图论问题，对于确定的起点与终点，我们可以通过状态的转移由之前求得的已知最短路来求得未知的最短路。如果要搞懂Floyd算法，你需要对于动态规划的知识点有所了解，否则可能会对其核心思想----状态转移无法理解。由于动态规划算法的特殊性，由代码来理解会更加轻松一些，因此笔者会在代码中写下详细的注释来帮助大家理解Floyd算法。

三、算法的代码实现

例题链接：Acwing Floyd求最短路

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

再给定 k 个询问，每个询问包含两个整数 x 和 y，表示查询从点 x 到点 y 的最短距离，如果路径不存在，则输出 impossible。

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

输入格式
第一行包含三个整数 n,m,k。

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

接下来 k 行，每行包含两个整数 x,y，表示询问点 x 到点 y 的最短距离。

输出格式
共 k 行，每行输出一个整数，表示询问的结果，若询问两点间不存在路径，则输出 impossible。

数据范围
1≤n≤200,
1≤k≤n2
1≤m≤20000,
图中涉及边长绝对值均不超过 10000。

输入样例：
3 3 2
1 2 1
2 3 2
1 3 1
2 1
1 3
输出样例：
impossible
1

这一道题就是一道Floyd算法的板子题，搞懂了这道题就可以说理解了Floyd算法。

#include<iostream>
#include<iomanip>
#include<cstdio>
#include<string>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#include<ctime>
#include<cstring>
#include<list>
#define ll long long
#define ull unsigned long long
#define INF 0x3f3f3f3f
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef  pair<int, int> PII;
const int N = 1e3 + 7;
int n, m, k;
int dp[N][N];   //动态规划数组，dp[i][j]代表从i到j的最短距离

void floyd()    //Floyd算法实现函数
{
	for (int k = 1; k <= n; k++)   //枚举中间的点，即一定经过k点
		for (int i = 1; i <= n; i++)    //枚举起点
			for (int j = 1; j <= n; j++)   //枚举终点
				dp[i][j] = min(dp[i][j], dp[i][k] + dp[k][j]);   //动态规划状态转移方程
}

void solve()
{
	cin >> n >> m >> k;
	for(int i=1;i<=n;i++)    //初始化操作
		for (int j = 1; j <= n; j++)
		{
			if (i == j)dp[i][j] = 0;   //自己到自己的距离为0
			else dp[i][j] = INF;  //把每两个不同点的距离初始化为无限大
		} 
	for (int i = 0; i < m; i++)   //输入边的操作
	{
		int x, y, z;
		cin >> x >> y >> z;
		dp[x][y] = min(dp[x][y], z);   //如果两点不止一条边，保留最短的边
	}
	
	floyd();   //实现算法
	
	while (k--)
	{
		int x, y;
		cin >> x >> y;
		if (dp[x][y] > INF / 2)    //这里在下面解释
			cout << "impossible" << endl;
		else
			cout << dp[x][y] << endl;
	}
}

int main()
{
	//std::ios::sync_with_stdio(false);
	//cin.tie(0), cout.tie(0);
	solve();
	return 0;
}

解释为什么不是">INF"而是“>INF/2”：由于题目中可能会有负数边存在，而在方程“dp[i][j] = min(dp[i][j], dp[i][k] + dp[k][j]); ”中我们可以发现，如果边为负数居然会把INF更新为INF-x，但是这个值实际上并不是我们能到达的距离。根据题目中数据范围，我们可以发现，我们虽然不能写成”>INF"，但是“>INF/2”是可以判断是否能到达终点的。

作者：Avalon Demerzel，喜欢我的博客就点个赞吧，更多图论与数据结构知识点请见作者专栏《图论与数据结构》