Floyd算法需要两个矩阵,a[i][j]表示顶点i到顶点j的距离,b[i][j]表示由i到j经过了b[i][j]记录的值所表示的顶点。
对a初始化,若ij不相邻则值为∞,b[i][j]初始值为j的值。如果有N个顶点,那么对a,b两个矩阵进行N次更新:先以v1点作为中介点假如找到a[1][v1]+a[v1][6]<a[1][6],那么更新a[1][6]=a[1][v1]+a[v1][6],b[1][6]=v1,把全部如此更新后,继续以v2为中介点遍历......
Floyd 算法本质是dp思想,可以回溯故可以处理负权图。Dijkstra是贪心思想每次选取最短的边去松弛其它边是不能处理负权图的,因为负值原因可能松弛出来的边比最短边还短。
核心代码 (O(N^3)):
for(int k=1;k<=n;k++) //有N个顶点,把每个顶点都作为中介点遍历一遍更新最短路
{
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]); //最后得到的dp数组记录所有的点之间的最短距离
}
}
HDU2544 最短路
Problem Description
在每年的校赛里,所有进入决赛的同学都会获得一件很漂亮的t-shirt。但是每当我们的工作人员把上百件的衣服从商店运回到赛场的时候,却是非常累的!所以现在他们想要寻找最短的从商店到赛场的路线,你可以帮助他们吗?
Input
输入包括多组数据。每组数据第一行是两个整数N、M(N<=100,M<=10000),N表示成都的大街上有几个路口,标号为1的路口是商店所在地,标号为N的路口是赛场所在地,M则表示在成都有几条路。N=M=0表示输入结束。接下来M行,每行包括3个整数A,B,C(1<=A,B<=N,1<=C<=1000),表示在路口A与路口B之间有一条路,我们的工作人员需要C分钟的时间走过这条路。
输入保证至少存在1条商店到赛场的路线。
Output
对于每组输入,输出一行,表示工作人员从商店走到赛场的最短时间
Sample Input
2 1
1 2 3
3 3
1 2 5 2 3 5 3 1 2
0 0
Sample Output
3
2
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define inf 0x3f3f3f3f
int n, m, dp[101][101];
int main()
{
while (scanf("%d%d", &n, &m) != EOF)
{
memset(dp, inf, sizeof(dp));
if (n == 0 && m == 0)
break;
for (int i = 1; i <= m; i++)
{
int a, b, c;
cin >> a >> b >> c;
dp[a][b] = dp[b][a] = c; //可知应是无向图dp[a][b]=dp[b][a]
}
for (int i = 1; i <= n; i++) //对角线上到自己的距离是0
dp[i][i] = 0;
for (int k = 1; k <= n; k++) //数据量不大,可用Floyd遍历
{
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]);
}
}
cout << dp[1][n] << endl;
}
return 0;
}