漫画:图的 “多源” 最短路径

简介: 下面我们来看一看Floyd算法的详细步骤。

640.jpg640.jpg



—————  第二天  —————


640.jpg640.png640.jpg640.jpg




小灰的思路如下:

第一步,利用迪杰斯特拉算法的距离表,求出从顶点A出发,到其他各个顶点的最短距离:



640.jpg


第二步,继续使用迪杰斯特拉算法,求出从顶点B出发,到其他各个顶点的最短距离。

第三步,从顶点C出发,到各个顶点的最短距离。

第四步,从顶点D出发......

.......

就像这样,一直遍历到顶点G。

这个思路的时间复杂度是多少呢?

假如图中有n个顶点,如果不考虑堆优化,一次迪杰斯特拉算法的时间复杂度是O(n^2)。所以,把每一个顶点都计算一遍,总的时间复杂度是O(n^3)。

640.jpg640.jpg640.jpg



————————————


640.jpg640.jpg640.jpg640.jpg640.jpg640.jpg640.jpg640.jpg



举一个栗子:

640.png



上图的顶点A和顶点C没有直接相连的边,它们之间的直接距离是无穷大。

如果以B作为“中继顶点”,此时A到C的最短路径就是A-B-C,最短距离是3+2=5。

640.jpg


再举一个栗子:

640.png


上图的顶点A和顶点C直接相连,距离是6。但是存在一条“迂回”路径A-B-C,距离是3+2=5<6。

所以,经过中继顶点B,从A到C的最短距离可以是5。

640.jpg


下面我们来看一看Floyd算法的详细步骤。



1.要实现Floyd算法,首先需要构建带权图的邻接矩阵:


640.png

在邻接矩阵当中,每一个数字代表着从某个顶点到另一个顶点的直接距离,这个距离是没有涉及到任何中继顶点的。

2.此时假定只允许以顶点A作为中继顶点,那么各顶点之间的距离会变成什么样子呢?

B和C之间的距离原本是无穷大,此时以A为中继,距离缩短为AB距离+AC距离=

5+2=7。

更新对应矩阵元素(橙色区域代表顶点A到其他顶点的临时距离):

640.png


3.接下来以顶点A、B作为中继顶点,那么各顶点之间的距离会变成什么样子呢?

A和D之间的距离原本是无穷大,此时以B为中继,距离缩短为AB距离+BD距离=5+1=6。

A和E之间的距离原本是无穷大,此时以B为中继,距离缩短为AB距离+BE距离=5+6=11。

更新对应矩阵元素(橙色区域代表顶点B到其他顶点的临时距离):

640.png

4.接下来以顶点A、B、C作为中继顶点,那么各顶点之间的距离会变成什么样子呢?

A和F之间的距离原本是无穷大,此时以C为中继,距离缩短为AC距离+CF距离=2+8=10。

更新对应矩阵元素(橙色区域代表顶点C到其他顶点的临时距离):

640.png

.........

.........

以此类推,我们不断引入新的中继顶点,不断刷新矩阵中的临时距离。

最终,当所有顶点都可以作为中继顶点时,我们的距离矩阵更新如下:


640.png

此时,矩阵中每一个元素,都对应着某顶点到另一个顶点的最短距离。640.jpg640.jpg



为什么这么说呢?让我们回顾一下动态规划的两大要素:

问题的初始状态问题的状态转移方程式

对于寻找图的所有顶点之间距离的问题,初始状态就是顶点之间的直接距离,也就是邻接矩阵。

而问题的状态转移方程式又是什么呢?

假设新引入的中继顶点是n,那么:

顶点i 到 顶点j 的新距离 = Min(顶点i 到 顶点j 的旧距离,顶点i 到 顶点n 的距离+顶点n 到 顶点j 的距离)

640.jpg

final
static
int
 INF 
=
Integer
.
MAX_VALUE
;
public
static
void
 floyd
(
int
[][]
 matrix
){
//循环更新矩阵的值
for
(
int
 k
=
0
;
 k
<
matrix
.
length
;
 k
++){
for
(
int
 i
=
0
;
 i
<
matrix
.
length
;
 i
++){
for
(
int
 j
=
0
;
 j
<
matrix
.
length
;
 j
++){
if
(
matrix
[
i
][
k
]
==
 INF 
||
 matrix
[
k
][
j
]
==
 INF
)
{
continue
;
}
                matrix
[
i
][
j
]
=
Math
.
min
(
matrix
[
i
][
j
],
 matrix
[
i
][
k
]
+
 matrix
[
k
][
j
]);
}
}
}
// 打印floyd最短路径的结果
System
.
out
.
printf
(
"最短路径矩阵: \n"
);
for
(
int
 i 
=
0
;
 i 
<
 matrix
.
length
;
 i
++)
{
for
(
int
 j 
=
0
;
 j 
<
 matrix
.
length
;
 j
++)
System
.
out
.
printf
(
"%3d  "
,
 matrix
[
i
][
j
]);
System
.
out
.
printf
(
"\n"
);
}
}
public
static
void
 main
(
String
[]
 args
)
{
int
[][]
 matrix 
=
{
{
0
,
5
,
2
,
 INF
,
 INF
,
 INF
,
 INF
},
{
5
,
0
,
 INF
,
1
,
6
,
 INF
,
 INF
},
{
2
,
 INF
,
0
,
6
,
 INF
,
8
,
 INF
},
{
INF
,
1
,
6
,
0
,
1
,
2
,
 INF
},
{
INF
,
6
,
 INF
,
1
,
0
,
 INF
,
7
},
{
INF
,
 INF
,
8
,
2
,
 INF
,
0
,
3
},
{
INF
,
 INF
,
 INF
,
 INF
,
7
,
3
,
0
}
};
    floyd
(
matrix
);
}
640.jpg 640.jpg
相关文章
|
7月前
|
并行计算 数据可视化 算法
CMplot & rMVP | 全基因组曼哈顿图和QQ图轻松可视化!
`CMplot`和`rMVP`是R语言中的两个包,用于全基因组关联分析(GWAS)的数据可视化。`CMplot`专注于曼哈顿图和QQ图的绘制,支持多种图表类型,如常见的SNP密度图、环状曼哈顿图、矩阵图、单条染色体图和多重曼哈顿图等。`rMVP`不仅包含了`CMplot`的功能,还支持更复杂的GWAS方法,如线性/混合线性模型和基因组选择算法,优化了内存管理和计算效率,特别适合大规模数据集。此外,它还提供PCA图和柱状图。两者都提供了丰富的参数定制图表。
340 1
CMplot & rMVP | 全基因组曼哈顿图和QQ图轻松可视化!
|
人工智能 数据挖掘
这图怎么画 | 相关分析棒棒糖图
这图怎么画 | 相关分析棒棒糖图
131 0
|
程序员 数据库
如何快速画出一副漂亮的架构图
为什么要画好一幅架构图?一幅漂亮的架构图既是创作者的深度结构化思考和表达,对于读者来说也更加容易理解架构所要表达的意思。
30372 16
|
人工智能 数据可视化 数据挖掘
这图怎么画| 富集分析之双向柱状图
这图怎么画| 富集分析之双向柱状图
198 0
|
存储 索引 Cloud Native
【刷题日记】1791. 找出星型图的中心节点
昨天加班比较晚来不及刷每日一题,今天加班也回来晚了点,但是看到是一个简单的图,那么还是刷刷吧
火山图|给你geneList,帮我标到火山图上
火山图|给你geneList,帮我标到火山图上
211 0
|
机器学习/深度学习
离散数学_十章-图 ( 4 ):图的表示和图的同构
离散数学_十章-图 ( 4 ):图的表示和图的同构
297 0
|
机器学习/深度学习 传感器 算法
基于 Bowyer-Watson算法实现delaunay德劳内三角网络和Voronoi泰森多边形的建立附matlab代码
基于 Bowyer-Watson算法实现delaunay德劳内三角网络和Voronoi泰森多边形的建立附matlab代码
|
存储 人工智能 缓存
图是什么?地图?网络图?还是……
图是什么?地图?网络图?还是……
954 0
图是什么?地图?网络图?还是……
最少转机(无路径之和) 图的广度优先遍历
最少转机(无路径之和) 图的广度优先遍历
最少转机(无路径之和) 图的广度优先遍历