迪科斯彻算法

这个算法是通过为每个顶点 v 保留目前为止所找到的从s到v的最短路径来工作的。初始时，原点 s 的路径长度值被赋为 0 （d[s] = 0），若存在能直接到达的边（s,m），则把d[m]设为w（s,m）,同时把所有其他（s不能直接到达的）顶点的路径长度设为无穷大，即表示我们不知道任何通向这些顶点的路径（对于 V 中所有顶点 v 除 s 和上述 m 外 d[v]= ∞）。当算法退出时，d[v] 中存储的便是从 s 到 v 的最短路径，或者如果路径不存在的话是无穷大。 Dijkstra 算法的基础操作是边的拓展：如果存在一条从 u 到 v 的边，那么从 s 到 v 的最短路径可以通过将边（u, v）添加到尾部来拓展一条从 s 到 u 的路径。这条路径的长度是 d[u] + w(u, v)。如果这个值比目前已知的 d[v] 的值要小，我们可以用新值来替代当前 d[v] 中的值。拓展边的操作一直运行到所有的 d[v] 都代表从 s 到 v 最短路径的花费。这个算法经过组织因而当 d[u] 达到它最终的值的时候每条边（u,v）都只被拓展一次。

算法维护两个顶点集 S 和 Q。集合 S 保留了我们已知的所有 d[v] 的值已经是最短路径的值顶点，而集合 Q 则保留其他所有顶点。集合S初始状态为空，而后每一步都有一个顶点从 Q 移动到 S。这个被选择的顶点是 Q 中拥有最小的 d[u] 值的顶点。当一个顶点 u 从 Q 中转移到了 S 中，算法对每条外接边 (u, v) 进行拓展。

迪科斯彻算法使用了广度优先搜索解决非负权有向图的单源最短路径问题，算法最终得到一个最短路径树。

算法实现：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>


#define MAX_VEX_NUM 30
#define Infinity 30000
#define YES 1
#define NO 0


typedef struct Graph * G;


struct Graph{


  int vexnum;
  int arcnum;
  int arcs[MAX_VEX_NUM][MAX_VEX_NUM];


};


void create_graph(G g);
void dij(G g,int d[],int p[],int v0);
int find_min(G g,int final[],int d[]);
void print_graph(G g,int p[]);


void main(){
  G g = (G) malloc(sizeof(struct Graph));
  memset(g,0,sizeof(struct Graph));
  create_graph(g);
  int d[MAX_VEX_NUM];
  int p[MAX_VEX_NUM];
  dij(g,d,p,0); 
  print_graph(g,p);
}


void create_graph(G g){
  printf("input vexnum and arcnum:\n");
  scanf("%d %d",&g->vexnum,&g->arcnum);
  int i = 0;
  int j = 0;
  for(i=0;i<MAX_VEX_NUM;i++){
     for(j=0;j<MAX_VEX_NUM;j++){
       g->arcs[i][j] = Infinity;
     }
  }
  
  printf("now, init the graph\n");
  
  int start = 0;
  int end = 0;
  int weight = 0;
  for(i=0;i<g->arcnum;i++){
    printf("input edge info:start end weight:\n");
    scanf("%d %d %d",&start,&end,&weight);
    g->arcs[start][end] = weight;
  }
}


void dij(G g,int d[],int p[],int v0){
  int final[MAX_VEX_NUM];
  int i = 0;
  for(i=0;i<MAX_VEX_NUM;i++){
    final[i] = NO;
    d[i] = Infinity;
    p[i] = -1;
  }
  d[v0] = 0;
  int j = 0;
  for(j=0;j<g->vexnum;j++){
     int v = find_min(g,final,d);
     if(v==0&&final[v]==YES){
       break;
     }
     final[v] = YES;
     int s = 0;  
     for(s=0;s<g->vexnum;s++){
        if((g->arcs[v][s]+d[v])<d[s]){
          d[s] = g->arcs[v][s] + d[v];
          p[s] = v;
        }
     }
  }
}




int find_min(G g,int final[],int d[]){
  int v = 0;
  int i = 0;
  for(i=0;i<g->vexnum;i++){
    if(final[i]==NO){
      if(final[v]==YES) v = i;
      if(d[v]>d[i]) v = i;
    }
  }
  return v;
}


void print_graph(G g,int p[]){
  int i = 0;
  int tmp = -1;
  printf("print_graph-dij-start\n");
  for(i=0;i<g->vexnum;i++){
    printf("%d",i);
    tmp = i;
    while(p[tmp]!=-1){
      printf("<-%d",p[tmp]);
      tmp = p[tmp];
    }
    printf("\n");  
  }
}