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

本文主要是介绍迪杰斯特拉(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}*/                          }}}
}