迪杰斯特拉(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 | 微博@从流域到海域