迪杰斯特拉算法2

【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为0107         }108         D[v0] = 0;    // V0至V0路径为0109         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; // 将目前找到的最近的顶点置为1124             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-04213 // 测试输入程序为：214 /*215 请输入顶点的个数:216 9217 请输入顶点的数据信息:218 A B C D E F G H I219 请输入边的个数:220 16221 请输入边的信息：222 0 1 1223 0 2 5224 1 2 3225 1 3 7226 1 4 5227 2 4 1228 2 5 7229 3 4 2230 3 6 3231 4 5 3232 4 6 6233 4 7 9234 5 7 5235 6 7 2236 6 8 7237 7 8 4238 打印所有输入信息：239 顶点信息240 A    B    C    D    E    F    G    H    I241 矩阵信息：242 A    B    C    D    E    F    G    H    I243 A    0    1    5655356553565535655356553565535244 B    1    0    3    7    565535655356553565535245 C    5    3    065535    1    7655356553565535246 D65535    765535    0    265535    36553565535247 E65535    5    1    2    0    3    6    965535248 F6553565535    765535    3    065535    565535249 G655356553565535    3    665535    0    2    7250 H65535655356553565535    9    5    2    0    4251 I655356553565535655356553565535    7    4    0252 253 254 求最短路径.....255 打印各顶点最短路径标记数组值： 1 1 1 1 1 1 1 1 1256 打印最短路径数组值： 0 1 4 7 5 8 10 12 16257 打印最短路径前驱数组值： 0 0 1 4 2 4 3 6 7258 打印V0到各个顶点最短路径值以及路径信息：259 V0~V1: 1260 V0->V1261 V0~V2: 4262 V0->V1->V2263 V0~V3: 7264 V0->V1->V2->V4->V3265 V0~V4: 5266 V0->V1->V2->V4267 V0~V5: 8268 V0->V1->V2->V4->V5269 V0~V6: 10270 V0->V1->V2->V4->V3->V6271 V0~V7: 12272 V0->V1->V2->V4->V3->V6->V7273 V0~V8: 16274 V0->V1->V2->V4->V3->V6->V7->V8275  */

关于实现代码中的SeqList.h文件和Stack.h文件从随笔《顺序表》和《栈》中查找拷贝一份即可。调试环境为VS2010。

http://www.cnblogs.com/Braveliu/p/3443311.html 顺序表

http://www.cnblogs.com/Braveliu/p/3452693.html 栈