Floyd算法的理解
应用场景
F l o y d Floyd Floyd算法经常用于求多源最短路问题,时间复杂度为 O ( n 3 ) O(n^3) O(n3),其中 n n n为图的顶点数。
算法原理
基于动态规划的思想,我们令
d
p
[
k
]
[
i
]
[
j
]
dp[k][i][j]
dp[k][i][j]表示,经过前
k
k
k个点
i
i
i到
j
j
j的最短路的长度。那么根据是否经过
k
k
k点,进行状态转移:
d
p
[
k
]
[
i
]
[
j
]
=
m
i
n
(
d
p
[
k
−
1
]
[
i
]
[
j
]
,
d
p
[
k
−
1
]
[
i
]
[
k
]
+
d
p
[
k
−
1
]
[
k
]
[
j
]
)
dp[k][i][j] = min(dp[k - 1][i][j], dp[k - 1][i][k] + dp[k - 1][k][j])
dp[k][i][j]=min(dp[k−1][i][j],dp[k−1][i][k]+dp[k−1][k][j])
d
p
dp
dp状态的初始化:
把
d
p
dp
dp数组初始化为图的邻接矩阵
时间复杂度
O
(
n
3
)
O(n^3)
O(n3)
空间复杂度
O
(
n
3
)
O(n^3)
O(n3)
优化
在枚举
d
p
dp
dp状态的时候,我们发现,每次用到的都是
k
−
1
k - 1
k−1的状态,也就是前一个点是否中转的状态,我们可以算到当前状态的时候把算出来的值直接覆盖掉前一个状态,即:
d
p
[
i
]
[
j
]
:
=
m
i
n
(
d
p
[
i
]
[
j
]
,
d
p
[
i
]
[
k
]
+
d
p
[
k
]
[
j
]
)
dp[i][j] := min(dp[i][j], dp[i][k] + dp[k][j])
dp[i][j]:=min(dp[i][j],dp[i][k]+dp[k][j])
因为
d
p
[
k
−
1
]
[
i
]
[
j
]
dp[k - 1][i][j]
dp[k−1][i][j]后面已经用不到了。
时间复杂度
O
(
n
3
)
O(n^3)
O(n3)
空间复杂度
O
(
n
2
)
O(n^2)
O(n2)
代码
memset(f, 127, sizeof f);
for(int i = 1; i <= n; i ++ ) f[i][i] = 0;
for(int i = 1; i <= m; i ++ ){
int x, y, z;
scanf("%d %d %d",&x, &y, &z);
f[x][y] = min(f[x][y], z);
}
for(int k = 1; k <= n; k ++ )
for(int i = 1; i <= n; i ++ )
for(int j = 1; j <= n; j ++ )
if(f[i][k] < 1 << 30 && f[k][j] < 1 << 30)
f[i][j] = min(f[i][j], f[i][k] + f[k][j]);