问题描述:
给定一个 n n n个点 m m m 条边的有向图,图中可能存在重边和自环, 边权可能为负数。
再给定 k k k 个询问,每个询问包含两个整数 x x x和 y y y,表示查询从点x到点y的最短距离,如果路径不存在,则输出“impossible”。
数据保证图中不存在负权回路。
输入格式:
第一行包含整数 n n n , m m m, k k k。
接下来 m m m 行每行包含三个整数 x x x, y y y, z z z,表示存在一条从点 x x x 到点 y y y 的有向边,边长为 z z z。
接下来 k k k行,每行包含两个整数 x x x, y y y,表示询问点 x x x到点 y y y的最短距离。
输出格式:
共
k
k
k行,每行输出一个整数,表示询问的结果,若询问两点间不存在路径,则输出 impossible
。
数据范围
1
≤
n
≤
200
,
1≤n≤200 ,
1≤n≤200,
1
≤
k
≤
n
2
,
1≤k≤n^2,
1≤k≤n2,
1
≤
m
≤
20000
,
1≤m≤20000,
1≤m≤20000,
图中涉及边长绝对值均不超过
10000
10000
10000。
输入样例:
3 3 2
1 2 1
2 3 2
1 3 1
2 1
1 3
输出样例:
impossible
1
时间复杂度 O ( n 3 ) O(n^3) O(n3)
思路:
1、Floyd使用了动态规划的思想求解多源汇最短路问题。
2、转移方程式: f [ k ] [ i ] [ j ] = m i n f [ k − 1 ] [ i ] [ j ] , f [ k − 1 ] [ i ] [ k ] + f [ k − 1 ] [ k ] [ j ] ) f[k][i][j] = min{f[k−1][i][j],f[k−1][i][k]+f[k−1][k][j]}) f[k][i][j]=minf[k−1][i][j],f[k−1][i][k]+f[k−1][k][j])
也就是说,从 i i i到 j j j的,且只经过 1 ∼ k 1 ∼ k 1∼k的点的最短路分为两类,如果其不经过点k,那么就是 f [ k − 1 ] [ i ] [ j ] f[k−1][i][j] f[k−1][i][j],如果经过k,由于不存在负环,所以只经过k一次,所以就是 f [ k − 1 ] [ i ] [ k ] + f [ k − 1 ] [ k ] [ j ] f[k-1][i][k]+f[k-1][k][j] f[k−1][i][k]+f[k−1][k][j],两者取最小者即可。
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 210, INF = 1e9;
int n, m, Q;
int d[N][N];
void floyd() {
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]);
}
int main() {
cin >> n >> m >> Q;
// 初始化的时候,把自环直接取成长度为0,也就是自环直接忽略
for (int i = 1; i <= n; i++)
for (int j = 1; j <= n; j++)
if (i == j) d[i][j] = 0;
else d[i][j] = INF;
while (m--) {
int a, b, w;
cin >> a >> b >> w;
// 读边的时候,有平行边则只存长度最小的那个边
d[a][b] = min(d[a][b], w);
}
floyd();
while (Q--) {
int a, b;
cin >> a >> b;
// 由于有负权边,所以即使到达不了,仍然有可能更新为比INF更小的值
if (d[a][b] > INF / 2) cout << "impossible" << endl;
else cout << d[a][b] << endl;
}
return 0;
}