算法思想
每次找到离源点最近的一个顶点,然后以该顶点为中心进行拓展,最终得到源点到其余各点的最短路径
在开始时,有必要将代码中的state数组说明一下:
struct state{
int predecessor;//当前节点的上一个节点
int length;//从源节点到当前节点的距离
enum{
permanent,tentative
}label;//枚举类型,表示当前节点state的状态,
//permanent表示state里面的数据是最终确定的,不会发生变化
//tentative表示state里面的数据是暂时的,还会发生变化
}state[MAX_NODES];
现在假如想找节点0到节点3的最短路径
开始时,我们将 节点0 标记为 永久的,在代码上显示为 state[0].label = permanent。然后我们依次检查每一个与 节点0 相邻的节点,并且用它们与 节点0 的距离重新计算源节点到他的最短距离
然后找到一个离源节点最近的节点k(即state[k].length最小)。根据节点k,对与节点k相连的节点进行松弛操作(松弛操作即改变节点的)。
之后我们将 节点k 标记为 永久的。我们依次检查每一个与 节点k 相邻的节点,并且用它们与 节点k 的距离重新计算源节点到他的最短距离(即更新state[x].length)
代码如下
#define MAX_NODES 1024
#define INFINITY 1000000000
int n,dist[MAX_NODES][MAX_NODES];
//n表示结点个数
//dist[i][j]:表示从节点i到j的距离
void shortest_path(int s,int t,int path[]){
//s:表示终点结点的编号
//t:表示起点结点的编号
//path[]:存放最后结果(最短路径)
struct state{
int predecessor;
int length;
enum{
permanent,tentative
}label;
}state[MAX_NODES];
int i,k,min;
struct state *p;
//算法第一步:初始化state数组
for(p = &state[0];p<&state[n];p++){//初始化state数组
p->predecessor = -1;
p->length = INFINITY;
p->label = tentative;
}
state[t].length = 0;
state[t].label = permanent;//设置源节点state对应的值
k = t;//初始时,离源节点最近的是不是他自己呢?
do{
//算法第二步:k是离源节点最近的节点(即state[k].length最小)。
//根据节点k,对与节点k相连的节点进行松弛操作。
//例如存在一条从k到v的边,那么可以通过将边k-v添加到尾部来拓展一条从t(源节点)到v的路径,
//这条路径的长度是state[k].length + dist[k][v]。
//如果这个值比目前已知的state[v].length的值小,我们可以用新值来代替它
for(i = 0;i<n;i++){
if(dist[k][i] != 0 && state[i].label == tentative){
if(state[k].length + dist[k][i] < state[i].length){
state[i].predecessor = k;
state[i].length = state[k].length + dist[k][i];
}
}
}
//算法第三步:找到那个离源节点最近的节点k,
//这里通过打擂台的方式。而且能上擂台的节点必须是tentative的,你说对不对?
k = 0;min = INFINITY;
for(i = 0;i<n;i++)
if(state[i].label == tentative && state[i].length < min){
min = state[i].length;
k = i;
}
state[k].label = permanent;
}while(k!=s);
i = 0;k = s;
do{
path[i++] = k;
k = state[k].predecessor;
}while(k>=0);
}