AcWing 854. Floyd求最短路

最新推荐文章于 2022-03-03 16:10:01 发布

铁头娃撞碎南墙

最新推荐文章于 2022-03-03 16:10:01 发布

阅读量115

点赞数

分类专栏：简单算法文章标签：动态规划 Floyd 多源最短路

本文链接：https://blog.csdn.net/m0_46680603/article/details/119683870

版权

简单算法专栏收录该内容

53 篇文章 8 订阅

订阅专栏

本文介绍了如何使用Floyd-Warshall算法求解有向图中从一个点到另一个点的最短路径，包括算法原理、时间复杂度分析以及C++代码实现。针对给定的AcWing竞赛题目，展示了如何处理负权边和查询操作。

摘要由CSDN通过智能技术生成

题目来源：AcWing 854. 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 \leq n \leq 200,$
$1≤k≤n^2$
$1 \leq m \leq 20000,$
图中涉及边长绝对值均不超过 $10000$ 。

输入样例：

输出样例：

impossible
1

二、Floyd最短路

算法时间复杂度 $O(n^3)$ ，其基本思路是动态规划 $d [k, i, j]$ ，也就是说 $i$ 这个点，逐一经过 $1$ ~ $k$ 点到达 $j$ 这个点的最短距离；然后通过 $d [k, i, j] = m i n (d [k - 1, i, j], d [k - 1, i, k] + d [k - 1, k, j]);$ 然后第一维度可以优化掉。

for (int k = 1; k <= n; k++)
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= n; j++)
			d[i][j]=min(d[i][j], d[i][k] + d[k][j]);

三、代码

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;

const int N = 210, INF= 1e9;
int g[N][N], n, m, s;

/*
	解题思路：动态规划的思想
    假设节点序号是从1到n。
    假设f[0][i][j]是一个n*n的矩阵，第i行第j列代表从i到j的权值，如果i到j有边，那么其值就为边ij的权值。
    如果没有边，那么其值就为无穷大。

    f[k][i][j]代表（k的取值范围是从1到n），在考虑了从1到k的节点作为中间经过的节点时，从i到j的最短路径的长度。

    比如，f[1][i][j]就代表：在考虑了1节点作为中间经过的节点时，从i到j的最短路径的长度。
    分析可知，f[1][i][j]的值无非就是两种情况，而现在需要分析的路径也无非两种情况，i=>j，i=>1=>j：
    【情况1】f[0][i][j]：i=>j这种路径的长度 < i=>1=>j这种路径的长度
    【情况2】f[0][i][1]+f[0][1][j]：i=>1=>j这种路径的长度 < i=>j这种路径的长度
    
    形式化说明如下：
    f[k][i][j]可以从两种情况转移而来：
    【情况1】从f[k−1][i][j]转移而来，表示i到j的最短路径不经过k这个节点
    【情况2】从f[k−1][i][k]+f[k−1][k][j]转移而来，表示i到j的最短路径经过k这个节点

    总结就是：f[k][i][j]=min(f[k−1][i][j],f[k−1][i][k]+f[k−1][k][j])
    从总结上来看，发现f[k]只可能与f[k−1]有关。
*/

int main()
{
	scanf("%d%d%d", &n, &m, &s);
	
	// 初始化工作，注意对角线是0
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= n; j++)
			if (i == j) g[i][j] = 0;
			else g[i][j] = INF;

	for (int i = 0; i < m; i++)
	{
		int a, b, c;
		scanf("%d%d%d", &a, &b, &c);
		g[a][b] = min(c, g[a][b]);
	}

	for (int k = 1; k <= n; k++)
		for (int i = 1; i <= n; i++)
			for (int j = 1; j <= n; j++)
				g[i][j] = min(g[i][j], g[i][k] + g[k][j]);

	while (s--)
	{
		int a, b;
		scanf("%d%d", &a, &b);
		if (g[a][b] > INF / 2) puts("impossible");
		else printf("%d\n", g[a][b]);
	}

	return 0;
}