Floyd算法是一种利用动态规划思想寻找多源点之间最短路径的算法。该算法时间复杂度为O(n^3),对于稠密图来说效果比执行[V]次的Dijkstra算法以及SPFA算法要好,其可用于边权为负的图。代码编写简单
DP思路
f[k][i][j]表示i和j之间可以通过编号为1…k的节点的最短路径。
有2种转义方式:
1.f[k][i][j]由f[k-1][i][j]转移而来,表示i到j不经过点k
2.f[k][i][j]由f[k-1][i][k]+f[k-1][k][j]转移而来,表示经过k这个点
那么状态转移方程式为:
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]=min(f[k−1][i][j],f[k−1][i][k]+f[k−1][k][j])
状态转移方程的优化
实际上看别人代码会发现状态转移方程式是
f
[
i
]
[
j
]
=
f
[
i
]
[
k
]
+
f
[
k
]
[
j
]
f[i][j] = f[i][k] + f[k][j]
f[i][j]=f[i][k]+f[k][j]
明显进行了状态压缩
此时的转移方程实际有4种情况
f
[
k
]
[
i
]
[
j
]
=
f
[
k
−
1
]
[
i
]
[
k
]
+
f
[
k
−
1
]
[
k
]
[
j
]
f[k][i][j] = f[k-1][i][k] + f[k-1][k][j]
f[k][i][j]=f[k−1][i][k]+f[k−1][k][j]
f
[
k
]
[
i
]
[
j
]
=
f
[
k
]
[
i
]
[
k
]
+
f
[
k
−
1
]
[
k
]
[
j
]
f[k][i][j] = f[k][i][k] + f[k-1][k][j]
f[k][i][j]=f[k][i][k]+f[k−1][k][j]
f
[
k
]
[
i
]
[
j
]
=
f
[
k
−
1
]
[
i
]
[
k
]
+
f
[
k
]
[
k
]
[
j
]
f[k][i][j] = f[k-1][i][k] + f[k][k][j]
f[k][i][j]=f[k−1][i][k]+f[k][k][j]
f
[
k
]
[
i
]
[
j
]
=
f
[
k
]
[
i
]
[
k
]
+
f
[
k
]
[
k
]
[
j
]
f[k][i][j] = f[k][i][k] + f[k][k][j]
f[k][i][j]=f[k][i][k]+f[k][k][j]
仔细分析后会发现f[k][i][k] = f[k-1][i][k], 因为k不可能是f[k][i][k]的中间节点(无环)
代码模板
void SPFA(){
for(int k=0; k<n; k++){
for(int i=0; i<n; i++{
for(int j=0; j<n; j++){
dp[i][j] = min(dp[i][j], dp[i][k]+dp[k][j]);
}
}
}
}