迪杰斯特拉(Dijkstra)算法求图中最短路径

迪杰斯特拉(Dijkstra )算法


对于图G=(V,E),将图的顶点分为两组:
顶点集S:已求出的最短路径的顶点集合(初始为{v0});
顶点集V-S:尚未求出最短路径的顶点集合(初始为V-{v0} )。
算法按最短路径长度的递增顺序逐个将V-S的顶点加入S中,直到所有顶点均被加入S为止。

算法需借助辅助数组dist[N], dist[i]表示目前已经找到的、从开始点v0到终点vi的当前最短路径的长度。
dist各元素的初值:若从v0到vi存在弧,则dist[i]为弧上的权值;否则dist[i]为∞。

算法执行过程

(1)当前下一条长度最短的路径必为 ( v0, … , vk ),vk满足如下条件:
dist[k]=Min{dist[i] | vi∈V-S}
求得顶点vk的最短路径后,将vk加入顶点集S中。
修正:每加入一个新的顶点vk到顶点集S,则对V-S中剩余的各顶点,多了一个“中转”结点vk,从而可能多一条“中转”路径,新的中转路径可能小于原来的路径,所以需对V-S剩余的各顶点的最短路径长度dist[i]进行修正。

(2)如何修正?
原来v0到vi的最短路径长度为dist[i],加入vk后,以vk为中间顶点的“中转”路径长度为:dist[k] + wki,(wki为弧”<”vk, vi>上的权值),若“中转”路径长度小于原dist[i](即 dist[k] + wki < dist[i]),则将顶点vi的最短路径长度修正为“中转”路径的长度。
(3)重复上述步骤(1)(2),直到所有顶点均加入S为止。

另外,为了记录从v0出发到各顶点的最短“路径”(顶点序列),使用辅助数组path[ ], path[i]表示当前找到的从开始顶点v0到顶点vi的当前最短路径顶点序列。
path[N]各元素的初值:如果从v0到vi有弧存在,则path[i]为(v0, vi);否则path[i]为空。

//Dijkstra算法
#define INFINITY 32768
typedef unsigned int WeightType
typedef WeightType Adjtex
typedef Seqlist VertexSet;

ShortestPath_DJS(AdjMartrix g, int v0, WeightType dist[MAX_VERTEX_NUM],
VertexSet path[MAX_VERTEX_NUM]) {
/*path[i]中存放当前顶点i的最短路径,是一个线性表。dist[i]中存放顶点i的当前最短路径长度*/
    VertexSet s;     /*s为已找到最短路径的终点集合*/
    for(i = 0; i < g.vexnum; i++) {
        InitList(&path[i]);               /*初始化dist[i]和path[i]*/
        dist[i] = g.arcs.[v0][i].adj;     /*dist[i]初始放当前顶点i到v0的路径长度*/
        if(dist[i] < INFINITY) {          /*如果经过上一步权值存在(路径存在)*/
            AddTail(&path[i], g.vertex[v0]);     /*AddTail是表尾添加操作*/
            AddTail(&path[i], g.vertex[i]);      /*构造从v0到i的路径序列,存储在顺序表数组的第i行*/
        }
    }
    InitList(&s);
    AddTail(&s, g.vertex[v0]);                /*将v0看成第一个已找到最短路径的终点*/
    for(t = 1; t <= g.vexnum - 1; t++) {      /*求v0到其余n-1个顶点额最短路径(n=g.vexnum)*/
        min = INFINITY;
        for(i = 0; t <= g.vexnum; i++) {
            if(!Member(g.vertex[i], s) && dist[i] < min) {//该点不存在s集合中且该点到v0的路径存在
                k = i; min = dist[i];
            }
        }
        if(min == INFINITY) return;   //如果经过上步没有修改min的值,说明当前顶点到v0没有路径 回退
        AddTail(&s, g.vertex[k]);     //把当前顶点加到s集合中
        for(i = 0; i < g.vertex; i++) {
        //如果该顶点不在s集合中且路径存在且该点通过已经加进s集合的点建立的新路径长度小于原路径
            if(!Member(g.vertex[i], s) && g.arcs[k][i].adj != INFINITY 
                       && (dist[k] + g.arc[k][i].adj < dist[i]) {
                dist[i] = dist[k] + g.arcs[k][i].adj;    //修改dist[i]为即最新的最短路径长度
                path[i] = path[k];                       //更新路径
                AddTail(&path[i], g.vertex[i]);     /*path[i] = path[k]∪{vi}*/                          
            }
        }
    }
}


知乎:Solo | 微博@从流域到海域

    原文作者:Dijkstra算法
    原文地址: https://blog.csdn.net/solo95/article/details/51449533
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞