数据结构与算法——图的最短路径之迪杰斯特拉（Dijkstra）算法

图的最短路径之迪杰斯特拉（Dijkstra）算法

看图分析：

1. 迪杰斯特拉算法的思想：是从起点开始，根据边的权值信息逐步推出到各个顶点的最小值，然后利用推导出来的最小权值顶点不断迭代更新到顶点的最小路径；

2. 结合上图来看下推导过程：
            a. 默认v0到v1、v2 ...v8各个顶点的最小路径为不可能的最大值记缓存信息为buffer[numVertex]（即未知）
            b. 从v0开始，遍历v0的所有边权值，取最小权值为1，为此时的最小路径；
            c. 当获取到v1的最小路径时，我们可以利用v1，推导出到顶点v2 、v3、v4此时已知条件下的最小路径并更新buffer：
                  v0---->v3 :  v0 -->v1--->v3    1 + 7 =8 
                  v0---->v4:   v0 -->v1--->v4     1 + 5 = 6
                  v0---->v2:   需要比较v0-->v2 与 v0-->v1--->v2两个路径，取最小值 ，（此处取v0-->v1--->v2：1 +3 = 4）  
            d.   经过第三步，我们在buffer中已经保存了v0 到 v1 、v2、v3、v4的已知条件下的最小路径，由于v0-->v1已经是已确定的最小路径，我们不用再处理它，现在对v0到 v2、v3、v4 选取最近路径值，显然是v2 路径为 v0-->v1--->v2

    e. 重复c步骤，检测更新与v2相连的顶点v1，v4、v5的最小路径，并更新buffer

    f. 重复d，从buffer中筛选最小路径。。。一直循环执行c、d步骤，直至numVertex-1 循环结束，找出最短路径

 

实现伪代码：

#include <iostream>
using namespace std;

#define MAXVER 9
#define INFINITY 65535

typedef char VertexType;
typedef int EdgeType;

typedef struct
{
    VertexType vexs[MAXVER];
    EdgeType arc[MAXVER][MAXVER];
    int numVertexes, numEdges;
}MGraph;

typedef int Pathmatirx[MAXVER];
typedef int ShortPathTable[MAXVER];

void ShortestPath_Dijkstra(MGraph G,int v0,Pathmatirx *P,ShortPathTable *D) 
{
    int v, w, k, min;
    int final[MAXVER];
    for ( v = 0; v < G.numVertexes; v++)
    {
        final[v] = 0;
        (*D)[v] = G.arc[v0][v];//v0顶点设置权值
        (*P)[v] = 0;
    }

    (*D)[v0] = 0;
    final[v0] = 1;

    for (v = 1; v < G.numVertexes; v++) 
    {
        min = INFINITY;//设置初始值为一个不可能的值
        for ( w = 0; w < G.numVertexes; w++)
        {//遍历所有顶点
            /*D[w]是V0->V的最小路径长度
            *  1. v0开始时：D[w],保存了v0所有边的权值，第一次简单的找出最小权值的顶点
            */
            if (!final[w] && (*D)[w] < min) 
            {
                k = w;
                min = (*D)[w];
            }
        }

        final[k] = 1; //更新目前已找到的v0到 顶点k的最小权值状态

        for (w = 0; w < G.numVertexes; w++) 
        {
            //D[w]:保存的是当前已获得到某个顶点w的最小权值
            //由于我们获得了一个新的最小权值顶点，因此需要更新D[w]缓存的信息
            //对于未找到最小路径的顶点进行更新，这里利用了已得到的最小值和k顶点边数组行的权值
            //理解 min + G.arc[k][w] 为 k 的边尾顶点到v0的最小距离
            if (!final[w] && (min + G.arc[k][w] < (*D)[w])) 
            {
                (*D)[w] = min + G.arc[k][w];//对应上述分析的C步骤，检测更新相连顶点的此时最优值
                (*P)[w] = k;//我们是为了寻找路径，所以必须知道当前前驱顶点是谁
            }
        }

    }

}