1.首先,引入一个辅助向量D,它的每个分量 D
表示当前所找到的
从起始点
(即源点
)到其它每个顶点
的长度。
![](https://i-blog.csdnimg.cn/blog_migrate/ee97014a38fbe2d5b732042f7fd5a3ae.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
2.D的初始状态为:若从
到
有弧(即从
到
存在连接边),则D
为弧上的权值(即为从
到
的边的权值);否则置D
为∞。
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
![](https://i-blog.csdnimg.cn/blog_migrate/ee97014a38fbe2d5b732042f7fd5a3ae.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
![](https://i-blog.csdnimg.cn/blog_migrate/ee97014a38fbe2d5b732042f7fd5a3ae.png)
显然,长度为 D
= Min{ D |
∈V } 的路径就是从
出发到顶点
的长度最短的一条路径,此路径为(
)。
![](https://i-blog.csdnimg.cn/blog_migrate/0749be93b0bc176b937b591204d4a27c.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/9cd7bb597f950ae30308e775cb068f36.png)
![](https://i-blog.csdnimg.cn/blog_migrate/16403fce3150543cad8e7b02be895526.png)
3.那么,下一条长度次短的是哪一条呢?也就是找到从源点
到下一个顶点的最短路径长度所对应的顶点,且这条最短路径长度仅次于从源点
到顶点
的最短路径长度。
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/9cd7bb597f950ae30308e775cb068f36.png)
假设该次短路径的终点是
,则可想而知,这条路径要么是(
),或者是(
)。它的长度或者是从
到
的弧上的权值,或者是D
加上从
到
的弧上的权值。
![](https://i-blog.csdnimg.cn/blog_migrate/57d9503337e3a36d86b4f6ff62c2a82f.png)
![](https://i-blog.csdnimg.cn/blog_migrate/a446e504782bee6bbc6acc377827b8ab.png)
![](https://i-blog.csdnimg.cn/blog_migrate/adfd5adc06d72d4461bb6a828f94f4a9.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/57d9503337e3a36d86b4f6ff62c2a82f.png)
![](https://i-blog.csdnimg.cn/blog_migrate/0749be93b0bc176b937b591204d4a27c.png)
![](https://i-blog.csdnimg.cn/blog_migrate/cb20b2c886a85ea68a7ae7bac1f0b95f.png)
![](https://i-blog.csdnimg.cn/blog_migrate/57d9503337e3a36d86b4f6ff62c2a82f.png)
4.一般情况下,假设S为已求得的从源点
出发的最短路径长度的顶点的集合,则可证明:下一条次最短路径(设其终点为
)要么是弧(
),或者是从源点
出发的中间只经过S中的顶点而最后到达顶点
的路径。
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/b16ba7064043a6234e2b47137ca3c82d.png)
![](https://i-blog.csdnimg.cn/blog_migrate/7da24e8adc88ae58c91a8e9df713eb9b.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/b16ba7064043a6234e2b47137ca3c82d.png)
因此,下一条长度次短的的最短路径长度必是D
= Min{ D
|
∈V-S },其中D
要么是弧(
)上的权值,或者是D
(
∈S)和弧(
,
)上的权值之和。
![](https://i-blog.csdnimg.cn/blog_migrate/0749be93b0bc176b937b591204d4a27c.png)
![](https://i-blog.csdnimg.cn/blog_migrate/ee97014a38fbe2d5b732042f7fd5a3ae.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
![](https://i-blog.csdnimg.cn/blog_migrate/ee97014a38fbe2d5b732042f7fd5a3ae.png)
![](https://i-blog.csdnimg.cn/blog_migrate/4efff6509e15ac17232d4b333326e44b.png)
![](https://i-blog.csdnimg.cn/blog_migrate/4fc4951c54e6932383f45e54d14cef25.png)
![](https://i-blog.csdnimg.cn/blog_migrate/57d9503337e3a36d86b4f6ff62c2a82f.png)
![](https://i-blog.csdnimg.cn/blog_migrate/57d9503337e3a36d86b4f6ff62c2a82f.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
算法描述如下:
1)令arcs表示弧上的权值。若弧不存在,则置arcs为∞(在本程序中为MAXCOST)。S为已找到的从
出发的的终点的集合,初始状态为空集。那么,从
出发到图上其余各顶点
可能达到的长度的初值为D=arcs[Locate Vex(G,
)],
∈V;
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/6974ded72c840d4ac670f57ada8bc835.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
2)选择
,使得D
=Min{ D |
∈V-S } ;
![](https://i-blog.csdnimg.cn/blog_migrate/9cd7bb597f950ae30308e775cb068f36.png)
![](https://i-blog.csdnimg.cn/blog_migrate/0749be93b0bc176b937b591204d4a27c.png)
![](https://i-blog.csdnimg.cn/blog_migrate/218fb1cc8dec27a1854e4bc0e1d4ad10.png)
按路径长度
递增次序产生算法:
把顶点集合V分成两组:
(1)S:已求出的顶点的集合(初始时只含有源点V0)
(2)V-S=T:尚未确定的顶点集合
将T中顶点按递增的次序加入到S中,保证:
(1)从源点V0到S中其他各顶点的长度都不大于从V0到T中任何顶点的最短路径长度
(2)每个顶点对应一个距离值
S中顶点:从V0到此顶点的长度
T中顶点:从V0到此顶点的只包括S中顶点作中间顶点的最短路径长度
依据:可以证明V0到T中顶点Vk的,或是从V0到Vk的直接路径的权值;或是从V0经S中顶点到Vk的路径权值之和
(
反证法可证)
求最短路径步骤
算法步骤如下:
G={V,E}
1. 初始时令 S={V0},T=V-S={其余顶点},T中顶点对应的距离值
若存在<V0,Vi>,d(V0,Vi)为<V0,Vi>弧上的权值
若不存在<V0,Vi>,d(V0,Vi)为∞
2. 从T中选取一个与S中顶点有关联边且权值最小的顶点W,加入到S中
3. 对其余T中顶点的距离值进行修改:若加进W作中间顶点,从V0到Vi的距离值缩短,则修改此距离值
重复上述步骤2、3,直到S中包含所有顶点,即W=Vi为止
大学经典教材<<数据结构>>(C语言版 严蔚敏 吴为民 编著) 中该算法的实现
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
|
/*
测试数据 教科书 P189 G6 的邻接矩阵 其中 数字 1000000 代表无穷大
6
1000000 1000000 10 100000 30 100
1000000 1000000 5 1000000 1000000 1000000
1000000 1000000 1000000 50 1000000 1000000
1000000 1000000 1000000 1000000 1000000 10
1000000 1000000 1000000 20 1000000 60
1000000 1000000 1000000 1000000 1000000 1000000
结果:
D[0] D[1] D[2] D[3] D[4] D[5]
0 1000000 10 50 30 60
*/
#include <iostream>
#include <cstdio>
#define MAX 1000000
using
namespace
std;
int
arcs[10][10];
//邻接矩阵
int
D[10];
//保存最短路径长度
int
p[10][10];
//路径
int
final[10];
//若final[i] = 1则说明 顶点vi已在集合S中
int
n = 0;
//顶点个数
int
v0 = 0;
//源点
int
v,w;
void
ShortestPath_DIJ()
{
for
(v = 0; v < n; v++)
//循环 初始化
{
final[v] = 0; D[v] = arcs[v0][v];
for
(w = 0; w < n; w++) p[v][w] = 0;
//设空路径
if
(D[v] < MAX) {p[v][v0] = 1; p[v][v] = 1;}
}
D[v0] = 0; final[v0]=0;
//初始化 v0顶点属于集合S
//开始主循环 每次求得v0到某个顶点v的最短路径 并加v到集合S中
for
(
int
i = 1; i < n; i++)
{
int
min = MAX;
for
(w = 0; w < n; w++)
{
//我认为的核心过程--选点
if
(!final[w])
//如果w顶点在V-S中
{
//这个过程最终选出的点 应该是选出当前V-S中与S有关联边
//且权值最小的顶点 书上描述为 当前离V0最近的点
if
(D[w] < min) {v = w; min = D[w];}
}
}
final[v] = 1;
//选出该点后加入到合集S中
for
(w = 0; w < n; w++)
//更新当前最短路径和距离
{
/*在此循环中 v为当前刚选入集合S中的点
则以点V为中间点 考察 d0v+dvw 是否小于 D[w] 如果小于 则更新
比如加进点 3 则若要考察 D[5] 是否要更新 就 判断 d(v0-v3) + d(v3-v5) 的和是否小于D[5]
*/
if
(!final[w] && (min+arcs[v][w]<D[w]))
{
D[w] = min + arcs[v][w];
// p[w] = p[v];
p[w][w] = 1;
//p[w] = p[v] + [w]
}
}
}
}
int
main()
{
cin >> n;
for
(
int
i = 0; i < n; i++)
{
for
(
int
j = 0; j < n; j++)
{
cin >> arcs[i][j];
}
}
ShortestPath_DIJ();
for
(
int
i = 0; i < n; i++)
printf
(
"D[%d] = %d\n"
,i,D[i]);
return
0;
}
|