假如有一个中转站(中转顶点 temp ),先从顶点 v 到达顶点 temp ,再从顶点 temp 到达顶点 u ,则这段距离我们可以表示为 dis[v][temp] + dis[temp][u]。
如果这段距离比目前的直接到达方式更小的话,则更新最短距离dis[v][u],那么 temp 这个中转顶点可以是图中的所有顶点。
如果我们不断把所有顶点插入作为中转顶点,再不断更新最短距离,最后,得到的 dis[n][n] 自然是多源点之间的最短路径了。
我们由上述推导也可求得状态转移方程 dis[v][u]=min(dis[v][u],dis[v][temp]+dis[temp][u])。
这就是Floyd算法的思想与原理,不断插入中转顶点,利用动态规划思想来实现。
这里还要记录了最短路径的方案,针对有些要求最短路径方案的题,用 path[v][u] 的值来记录顶点v到顶点u的中转顶点,若为 -1 则表示无中转顶点。
//Floyd算法核心
void floyd(int n){
//顶点数目。
int i,j,k;//循环变量。
memset(dis,inf,sizeof(dis));
memset(path,-1,sizeof(path));
for(i=0;i<n;i++){
for(int j=0;j<n;j++){
//初始化dis。
dis[i][j]=graph[i][j];
}
}
for(k=0;k<n;k++){
//不断插点
for(i=0;i<n;i++){
for(j=0;j<n;j++){
//对图中所有点之间的最短距离进行更新。
if(dis[i][j]>dis[i][k]+dis[k][j]){
path[i][j]=k;//记录中转顶点。
dis[i][j]=dis[i][k]+dis[k][j];//更新。
//等价于 d[i][j] |= d[i][k] & d[k][j]
}
}
}
}
}