本题有用到动态规划的思想,现对floyd()内 d[i][j] = min(d[i][j],d[i][k] + d[k][j]) 一行作出阐释
假设存在1~n共n个节点,d[k][i][j]是大小为n*n的矩阵,表示为经过k节点情况下i到j节点的最短路径长度,当i和j间是断路,对应值是无穷大。
经过min操作后得到d[k][i][j]的结果只会存在两种情况:
d[k-1][i][j] 即i到j的最短距离不经过k点
d[k-1][i][k]+d[k-1][k][j] 即i到j最短路经过k点
综上不难看出,d[k]只可能和d[k-1]有关,所以原式:d[k][i][j]=min(d[k−1][i][j],d[k−1][i][k]+d[k−1][k][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 d[N][N];//邻接矩阵
int n,m,k,a,b,w;
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 >> k;
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--)//录入图
{
scanf("%d%d%d",&a,&b,&w);
d[a][b] = min(d[a][b],w);//因为存在重边所以保留最小的
}
floyd();
while(k--)
{
scanf("%d%d",&a,&b);
if(d[a][b] > INF/2) puts("impossible");//设置 /2 的主要原因就是存在负权边,具体解释可参照笔记 853
else cout << d[a][b] << endl;
}
return 0;
}