前言:在求任意两点间的最短路问题中,图一般较为稠密,使用Floyd算法可以在O(N ^ 3)的时间实现。当然也可以把每个点作为起点,求解N次单源最短路径问题,但较为复杂。这里介绍Floyd算法以及使用Floyd算法打印路径和解决选址问题
算法分析
假设用d[k, i, j]表示“经过若干编号不超过k的结点”从i 到 j的最短路长度,这个问题可以划分为两个子问题,经过编号不超过 k - 1 的结点从i到j或者从i 先到k ,再到j,所以有:
d[k, i, j] = min(d[k - 1, i, j], d[k - 1, i, k] + d[k - 1, k, j])
因为一般解决稠密图,所以使用邻接矩阵。
Floyd算法的本质是动态规划,k表示我们所划分的阶段,所以在循环的最外层,但我们写的时候已经将其优化,像揹包问题等动态规划中那样,使用“滚动数组”优化。
所以才有: d[i][j] = min(d[i][j] , d[i][k] + d[k][j]);
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
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 210, INF = 1e9;
int n, m, q;
int d[N][N];
void flody()
{
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;
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, c;
cin >> a >> b >> c;
d[a][b] = min(d[a][b], c);
}
flody();
while(q -- )
{
int x, y;
cin >> x >> y;
int t = d[x][y];
if(t > INF / 2) puts("impossible");
else cout << t << endl;
}
return 0;
}
Floyd打印最短路径
path[i][j]数组存储 i ~ j这条路上 j 前面那个点,假设为k,i ~ j 拆开就是i~ k 和 k ~ j,可以求出k前面的点,递归将i ~ j上经过的所有点全部打印出来。
样例:
4 5
1 2 1
1 3 4
2 3 1
3 4 1
1 4 6
1 4
输出:
1 -> 4的最短路径为:1 2 3 4
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 210, INF = 1e6;
int n, m;
int d[N][N],path[N][N];
void init();
void floyd();
void printpath(int a, int b);
int main()
{
memset(d, 0x3f, sizeof d);
cout << "请输入点、边的数量:" << endl;
cin >> n >> m;
init();
cout << "请输入边的信息:" << endl;
for(int i = 0; i < m; i ++ )
{
int a, b, c;
cin >> a >> b >> c;
d[a][b] = min(d[a][b], c);
}
floyd();
cout << "请输入查询路径:" << endl;
int x, y;
cin >> x >> y;
printf("%d -> %d的最短路径为:", x, y);
printpath(x, y);
cout << endl;
return 0;
}
void init()
{
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;
path[i][j] = i;
}
}
void floyd()
{
for(int k = 1; k <= n; k ++ )
for(int i = 1; i <= n; i ++ )
for(int j = 1; j <= n ; j++ )
{
if(d[i][j] > d[i][k] + d[k][j])
{
d[i][j] = d[i][k] + d[k][j]; // 从i ~ j 的最短路长度
path[i][j] = k;
}
}
}
void printpath(int a, int b)
{
if(a == b)
{
cout << a << " ";
return ;
}
int k = path[a][b];
printpath(a, k);
cout << b << " ";
}
选址问题
比如建造一个加油站,使所有村庄到加油站的距离最近,问加油站应该建在哪里。
通过Floyd计算出每个点到所有点的最短距离总和,找到总和最小点的就是加油站的地址。
样例可以自己模拟一个:
第一行输入点数n和边数m
然后m行依次输入每个点的信息
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 210, INF = 1e6;
int n, m;
int d[N][N];
void floyd();
int main()
{
cin >> n >> m;
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, c;
cin >> a >> b >> c;
d[a][b] = min(d[a][b], c);
}
floyd();
int ans = INF; //存路经总和
int location = 0; //加油站位置
//确定加油站的位置
for(int i = 1; i <= n; i ++ )
{
int sum = 0; //用来更新每个点到其他所有点的距离之和
for(int j = 1; j <= n; j ++ )
sum += d[i][j]; //将此点到所有村庄的路程加起来
//求出路径之和加起来最短的
if(sum < ans)
{
ans = sum;
location = i; //i即为医院所在位置
}
}
cout << "位置及最短路程总和为:";
cout << location << " " << ans;
return 0;
}
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]);
}