【1】最短路径
最短路径？别乱想哈，其实就是字面意思，一个带边值的图中从某一个顶点到另外一个顶点的最短路径。
官方定义：对于内网图而言，最短路径是指两顶点之间经过的边上权值之和最小的路径。
并且我们称路径上的第一个顶点为源点，最后一个顶点为终点。
由于非内网图没有边上的权值，所谓的最短路径其实是指两顶点之间经过的边数最少的路径。
别废话了！整点实际的哈，你能很快计算出下图中由源点V0到终点V8的最短路径吗？
【2】迪杰斯特拉算法
迪杰斯特拉算法是按路径长度递增的次序产生最短路径的思路求解。
具体算法及其详细讲解如下：
阅读程序前，先要搞明白几个数组作用：
final[w]=1; 表示V0到Vw顶点已经有最短路径的结果
ShortPathTable[w]; 表示V0到Vw顶点的最短路径权值和
Pathmatirx[w]; 表示V0到Vw顶点的前驱顶点下标值
开始调用算法前，我们需要为案例图创建邻接矩阵图，如下所示：
（1）程序开始运行，final数组是为了标记V0到某顶点是否已经求得最短路径。
如果V0到Vw已经有结果，那么final[w]=1;
（2）第5~10行，是对数据初始化工作。 此时final数组均赋值为0，表示所有点均未求得最短路径。
D数组为 {0,1,5,65515,65535,65535,65535,65535,65535}。因为V0与V1和V2的边权值为1和5。
P数组全为0，表示目前没有路径。
（3）第11行，表示V0至V0路径为0。
第12行，表示V0点到V0点已经求得最短路径，因此final[0]=1。
此时final数组为 {1,0,0,0,0,0,0,0,0}，此时整个初始化工作完成。
（4）第13~33行为主循环，每次循环求得V0与一个顶点的最短路径。除去V0，因此索引从1开始。
（5）第15~23行，先令min为65535的极大值，通过w控制循环，与D[w]比较找到最小值为min=1，同时确定k=1。
（6）第24行，由k=1可知与V0最近的顶点是V1，并且由D[1]=1知道此时V0到V1的最短路径是1。
因此再将对应的final[1]设置为1。此时final数组为 {1,1,0,0,0,0,0,0,0}
（7）第25~32行是一循环，此循环甚为关键。
它的目的是在刚才已经找到V0与V1的最短路径基础之上，对V1与其它顶点的边进行计算，得到V0与它们的当前最短距离，如图所示：
因为min=1,所以D[2]=5不再是V0到V2的最短路径，现在D[2]=V0->V1->V2=min+3=4, D[3]=V0->V1->V3=min+7=8,
D[4]=V0->V1->V4=min+5=6，于是D数组当前值为 {0,1,4,8,6,65535,65535,65535,65535}。
而P[2]=1，P[3]=1，P[4]=1，其表示V0到V2，V3，V4点的最短路径它们的前驱均是V1。
此时P数组为 {0,0,1,1,1,0,0,0,0}。
（8）新一轮循环，此时i=2。第15~23行，对w循环，注意因为final[0]=1和final[1]=1，由第18行的！final[w]可知：
V0与V1并不参与最小值的获取。通过循环比较，找到最小值min=4，k=2。
（9）第24行，由k=2，表示已经求出V0到V2的最短路径，并且由D[2]=4知道最短路径距离为4。
因此将V2对应的final[2]设置为1，此时final数组为 {1,1,1,0,0,0,0,0,0}。
（10）第25~32行，在刚才已经找到V0与V2的最短路径的基础上，对V2与其它顶点的边进行计算，得到V0与它们的最短距离。
如图所示：
因为min=4，所以D[4]=6不再是V0到V4的最短距离，现在D[4]=V0->V2->V4=min+1=5，D[5]=V0->V2->V5=min+7=11。
因此D数组当前值为 {0,1,4,8,5,11,65535,65535,65535,65535}。
而原本P[4]=1，此时P[4]=2，P[5]=2，它表示的意思是V0到V4和V5的最短路径前驱均为V2。
此时P数组为 {0,0,1,1,2,2,0,0,0}。
（11）重新再开始一轮循环，此时i=3。第15~23行，通过对w循环比较找到最小值min=5，k=4。
（12）第24行，由k=4表示已经求出V0到V4的最短路径，并且由D[4]=5知道最短路径距离为5。
因此将V4对应的final[4]设置为1。此时final数组为 {1,1,1,0,1,0,0,0,0}。
（13）第25~32行，对V4与其它顶点的边值进行计算，得到V0与它们的当前最短距离，如图所示：
因为min=5，所以D[3]=8不再是V0到V3的最短路径，现在D[3]=V0->V4->V3=min+2=7，同理：
D[5]=V0->V4->V5=min+3=8，D[6]=V0->V4->V6=min+6=11，
D[7]=V0->V4->V7=min+9=14，因此，D数组当前值为 {0,1,4,7,5,8,11,14,65535}。
而原本P[3]=1，此时P[3]=4，原本P[5]=2，此时P[5]=4。
另外P[6]，P[7]=4，它表示V0到V3，V5，V6，V7点的最短路径它们的前驱是V4。
此时P数组值为 {0,0,1,4,2,4,4,4,0}。
（14）之后的循环完全类似。得到最终的结果，如图所示：
此时final数组为 {1,1,1,1,1,1,1,1,1}，它表示所有的顶点均完成了最短路径的查找工作。
此时D数组为 {0,1,4,7,5,8,10,12,16}，它表示V0到各个顶点的最短路径数，比如D[8]=1+3+1+2+3+2+4=16。
此时的P数组为 {0,0,1,4,2,4,3,6,7}：
这个P数组值可能难理解一些。比如P[8]=7，它表示要查看V0到V8的最短路径时，顶点V8的前驱是V7；
再由P[7]=6表示要查看V0到V7的最短路径时，顶点V7的前驱是V6，同理，P[6]=3表示V6的前驱是V3。
这样就可以得到：V0到V8最短路径为V8<-V7<-V6<-V3<-V4<-V2<-V1<-V0
即就是： V0->V1->V2->V4->V3->V6->V7->V8。
【3】算法实现
实现代码如下：
1 #include <iostream>
2 #include "SeqList.h"
3 #include "Stack.h"
4 #include <iomanip>
5 using namespace std;
6
7 #define INFINITY 65535
8
9 template<class NameType, class DistType>
10 class Graph
11 {
12 private:
13 SeqList<NameType> Vertices;
14 DistType **Edges;
15 int nVer, nEdges;
16
17 public:
18 Graph()
19 : Edges(NULL)
20 , nEdges(0)
21 , nVer(0)
22 {}
23 ~Graph()
24 {}
25
26 public:
27 int GetVer() const
28 {
29 return nVer;
30 }
31
32 istream & operator>>(istream &in)
33 {
34 int v, u, value;
35 int i, j;
36 NameType item;
37 cout << "请输入顶点的个数: " << endl;
38 in >> nVer;
39 cout << "请输入顶点的数据信息: " << endl;
40 for (i = 0; i < nVer; ++i)
41 {
42 in >> item;
43 Vertices.push_back(item); // 保存全部顶点
44 }
45 /////二维数组的创建并初始化
46 Edges = new DistType*[nVer]; // DistType *ar[10];
47 for (i = 0; i < nVer; ++i)
48 {
49 Edges[i] = new DistType[nVer];
50 for (j = 0; j < nVer; ++j)
51 {
52 Edges[i][j] = 0;
53 }
54 }
55 cout << "请输入边的个数: " << endl;
56 in >> nEdges;
57 cout << "请输入边的信息：" << endl;
58 for (i = 0; i < nEdges; ++i)
59 {
60 in >> v >> u >> value;
61 Edges[v][u] = value;
62 Edges[u][v] = value;
63 }
64 return in;
65 }
66 ostream & operator<<(ostream &out) const
67 {
68 int i, j;
69 out << "顶点信息 " << endl;
70 for (i = 1; i <= nVer; ++i)
71 {
72 out << Vertices[i] << setw(5);
73 }
74 out << endl;
75 out << "矩阵信息：" << endl;
76 out << setw(10);
77 for (i = 1; i <= nVer; ++i)
78 {
79 out << Vertices[i] << setw(5);
80 }
81 out << endl;
82 for (i = 0; i < nVer; ++i)
83 {
84 out << Vertices[i+1] << setw(5);
85 for (j = 0; j < nVer; ++j)
86 {
87 if (0 == Edges[i][j] && i != j)
88 Edges[i][j] = INFINITY;
89 cout << Edges[i][j] << setw(5);
90 }
91 out << endl;
92 }
93 out << endl;
94
95 return out;
96 }
97 // 迪杰斯特拉算法实现
98 void ShortestPath_Dijkstra(int v0, int* final, int*p, int *D)
99 {
100 int v, w, k, min;
101 // 初始化数据
102 for (v = 0; v < nVer; ++v)
103 {
104 final[v] = 0; // 全部顶点初始化为未知对短路径状态
105 D[v] = Edges[v0][v]; //将与V0点有连线的顶点加上权值
106 p[v] = 0; // 初始化路径数组p为0
107 }
108 D[v0] = 0; // V0至V0路径为0
109 final[v0] = 1; // final[W]=1表示V0至V0不需要求路径
110 // 开始主循环，每次求得V0到某个V顶点的最短路径
111 for (v = 1; v < nVer; ++v)
112 {
113 min = INFINITY; // 当前所知离V0顶点最近距离
114 for (w = 0; w < nVer; ++w) // 寻找离V0最近的顶点
115 {
116 if (!final[w] && D[w] < min)
117 {
118 min = D[w]; // w顶点离V0顶点更近
119 k = w;
120 }
121 }
122
123 final[k] = 1; // 将目前找到的最近的顶点置为1
124 for (w = 0; w < nVer; ++w) // 修正当前最短路径距离
125 {
126 // 如果经过V顶点的路径比现在这条路径的长度短的话
127 if (!final[w] && (min + Edges[k][w] < D[w]))
128 {
129 // 说明找到了最短的路径，修改D[w] 和 p[w]
130 D[w] = min + Edges[k][w]; // 修改当前路径长度
131 p[w] = k;
132 }
133 }
134 }
135 }
136 };
137
138 template<class NameType, class DistType>
139 istream & operator>>(istream &in, Graph<NameType,DistType> &g)
140 {
141 g >> in;
142 return in;
143 }
144
145 template<class NameType, class DistType>
146 ostream & operator<<(ostream &out, const Graph<NameType,DistType> &g)
147 {
148 g << out;
149 return out;
150 }
151
152 void main()
153 {
154 Graph<char, int> myg;
155 cin >> myg;
156 cout << "打印所有输入信息：" << endl;
157 cout << myg << endl;
158 cout << "求最短路径....." << endl;
159 int numVer = myg.GetVer();
160 int* pFinal = new int[numVer];
161 int* pPathmatirx = new int[numVer];
162 int* pShortPath = new int[numVer];
163 myg.ShortestPath_Dijkstra(0, pFinal, pPathmatirx, pShortPath);
164 cout << "打印各顶点最短路径标记数组值：" << " ";
165 for (int i = 0; i < numVer; ++i)
166 {
167 cout << pFinal[i] << " ";
168 }
169 cout << endl;
170 cout << "打印最短路径数组值：" << " ";
171 for (int i = 0; i < numVer; ++i)
172 {
173 cout << pShortPath[i] << " ";
174 }
175 cout << endl;
176 cout << "打印最短路径前驱数组值：" << " ";
177 for (int i = 0; i < numVer; ++i)
178 {
179 cout << pPathmatirx[i] << " ";
180 }
181 cout << endl;
182 cout << "打印V0到各个顶点最短路径值以及路径信息：" << endl;
183 SeqStack<int> sQ;
184 for (int i = 1; i < numVer; ++i)
185 {
186 cout << "V0~V" << i << ": " << pShortPath[i] << endl;
187
188 sQ.Push(pPathmatirx[i]);
189 int n = 0;
190 while (sQ.GetTop(n) && n != 0)
191 {
192 sQ.Push(pPathmatirx[n]);
193 }
194
195 while (!sQ.IsEmpty())
196 {
197 int m = 0;
198 sQ.Pop(m);
199 cout << "V" << m << "->";
200 }
201 cout << "V" << i << endl;
202 }
203 delete []pFinal;
204 delete []pPathmatirx;
205 delete []pShortPath;
206 pFinal = NULL;
207 pPathmatirx = NULL;
208 pShortPath = NULL;
209 }
210 // 备注：
211 // 最短路径迪杰斯特拉算法实现
212 // 整理于2013-12-04
213 // 测试输入程序为：
214 /*
215 请输入顶点的个数:
216 9
217 请输入顶点的数据信息:
218 A B C D E F G H I
219 请输入边的个数:
220 16
221 请输入边的信息：
222 0 1 1
223 0 2 5
224 1 2 3
225 1 3 7
226 1 4 5
227 2 4 1
228 2 5 7
229 3 4 2
230 3 6 3
231 4 5 3
232 4 6 6
233 4 7 9
234 5 7 5
235 6 7 2
236 6 8 7
237 7 8 4
238 打印所有输入信息：
239 顶点信息
240 A B C D E F G H I
241 矩阵信息：
242 A B C D E F G H I
243 A 0 1 5655356553565535655356553565535
244 B 1 0 3 7 565535655356553565535
245 C 5 3 065535 1 7655356553565535
246 D65535 765535 0 265535 36553565535
247 E65535 5 1 2 0 3 6 965535
248 F6553565535 765535 3 065535 565535
249 G655356553565535 3 665535 0 2 7
250 H65535655356553565535 9 5 2 0 4
251 I655356553565535655356553565535 7 4 0
252
253
254 求最短路径.....
255 打印各顶点最短路径标记数组值： 1 1 1 1 1 1 1 1 1
256 打印最短路径数组值： 0 1 4 7 5 8 10 12 16
257 打印最短路径前驱数组值： 0 0 1 4 2 4 3 6 7
258 打印V0到各个顶点最短路径值以及路径信息：
259 V0~V1: 1
260 V0->V1
261 V0~V2: 4
262 V0->V1->V2
263 V0~V3: 7
264 V0->V1->V2->V4->V3
265 V0~V4: 5
266 V0->V1->V2->V4
267 V0~V5: 8
268 V0->V1->V2->V4->V5
269 V0~V6: 10
270 V0->V1->V2->V4->V3->V6
271 V0~V7: 12
272 V0->V1->V2->V4->V3->V6->V7
273 V0~V8: 16
274 V0->V1->V2->V4->V3->V6->V7->V8
275 */
关于实现代码中的SeqList.h文件和Stack.h文件从随笔《顺序表》和《栈》中查找拷贝一份即可。调试环境为VS2010。
Good Good Study, Day Day Up.
顺序 选择 循环 总结
转载请注明本文地址： 数据结构图之三（最短路径--迪杰斯特拉算法）