Dijkstra算法

http://ghj19850926.blog.163.com/blog/static/1859156020141115522903/

 Dijkstra算法又称为单源最短路径,所谓单源是在一个有向图中,从一个顶点出发,求该顶点至所有可到达顶点的最短路径问题。 
      设G=(V,E)是一个有向图,V表示顶点,E表示边。它的每一条边(i,j)属于E,都有一个非负权W(I,j),在G中指定一个结点v0,要求把从v0到G的每一个接vj(vj属于V)的最短有向路径找出来(或者指出不存在)。
      Dijstra算法是运用贪心的策略,从源点开始,不断地通过相联通的点找出到其他点的最短距离
基本思想是:
      设置一个顶点的集合s,并不断地扩充这个集合,一个顶点属于集合s当且仅当从源点到该点的路径已求出。开始时s中仅有源点,并且调整非s中点的最短路径长度,找当前最短路径点,将其加入到集合s,直到终点在s中。
基本步骤:
1、把所有结点分成两组:
      第一组:包括已经确定最短路径的结点;
      第二组:包括尚未确定最短路径的结点。
2、开始时,第一组只包含起点,第二组包含剩余的点;
3、用贪心的策略,按最短路径长度递增的顺序把第二组的结点加到第一组去,直到v0可达的所有结点都包含于第一组中。在这个过程中,不断更新最短路径,总保持从v0到第一组各结点的最短路径长度dist都不大于从v0到第二组任何结点的路径长度。
4、每个结点对应一个距离值,第一组结点对应的距离就是v0到此结点的最短路径长度,第二组结点对应的距离值就是v0由第一组结点到此结点的最短路径长度。
5、直到所有的顶点都扫描完毕(v0可达的所有结点都包含于第一组中),找到v0到其它各点的所有最短路径。

 动画演示:http://www.jcc.jx.cn/kejiandb/Dijkstra.swf

 如图:求0点到其他点的最短路径。

技术分享技术分享

(1)开始时,s1={v0},s2={v1,v2,v3,v4},v0到各点的最短路径是{0,10,&,30,100};
(2)在还未进入s1的顶点之中,最短路径为v1,因此s1={v0,v1},由于v1到v2有路径,因此v0到各点的最短路径更新为{0,10,60,30,100};
(3)在还未进入s1的顶点之中,最短路径为v3,因此s1={v0,v1,v3},由于v3到v2、v4有路径,因此v0到各点的最短路径更新为{0,10,50,30,90};
(4)在还未进入s1的顶点之中,最短路径为v2,因此s1={v0,v1,v3,v2},由于v2到v4有路径,因此v0到各点的最短路径更新为{0,10,50,30,60};
数据结构:
(1)用一个二维数组a[i..j,i..j]来存储各点之间的距离,10000表示无通路:
(2)用数组dist[i..j]表示最短路径;
(3)用集合s表示找到最短路径的结点。 

技术分享
 1技术分享//单源最短路(dijkstra算法)
 2技术分享program dijkstra;
 3技术分享const max=10000;    //表示无通路
 4技术分享type jihe=set of 0..100;  //顶点数
 5技术分享var
 6技术分享   a:array[0..100,0..100] of integer;//各点之间的距离
 7技术分享   dist:array[0..100] of integer; //最短路径
 8技术分享   i,j,k,m,n:integer;
 9技术分享   s:jihe;
10技术分享
11技术分享procedure init;
12技术分享begin
13技术分享   s:=[0];    //开始集合只包含源点
14技术分享   readln(n);
15技术分享   assign(input,‘dijs.in‘);reset(input);
16技术分享   for i:=0 to n do    //读入各点之间的距离
17技术分享      for j:=0 to n do
18技术分享         read(a[i,j]);
19技术分享   for i:=0 to n do    //源点到各点之间的直接距离
20技术分享      dist[i]:=a[0,i];
21技术分享end;
22技术分享
23技术分享procedure dijsk;
24技术分享begin
25技术分享   for i:=0 to n do
26技术分享   begin
27技术分享      m:=0;
28技术分享      dist[m]:=max;  //初始化源点到各点的最小距离为无穷
29技术分享      for j:=1 to n do
30技术分享         if(not(j in s))and(dist[m]>dist[j]) then  //找出当前的最小距离
31技术分享        m:=j;    //记录找到的顶点    
32技术分享         s:=s+[m]; //把顶点加入集合
33技术分享      for k:=0 to n do    //更新经过新节点到其它点之间的最小距离
34技术分享         if(not(k in s))and(dist[k]>dist[m]+a[m,k]) then
35技术分享            dist[k]:=dist[m]+a[m,k];
36技术分享   end;
37技术分享end;
38技术分享
39技术分享begin
40技术分享   init;
41技术分享   dijsk;
42技术分享   for i:=1 to n do
43技术分享      writeln(i,‘:‘,dist[i]);
44技术分享   close(input);
45技术分享end.
 
http://ghj19850926.blog.163.com/blog/static/185915602014112351326191/
 

 Dijkstra算法解决了有向图上带正权值的单源最短路径问题,其运行时间要比Bellman-Ford算法低,但适用范围比Bellman-Ford算法窄。

迪杰斯特拉提出的按路径长度递增次序来产生源点到各顶点的最短路径的算法思想是:对有n个顶点的有向连通网络G=(V, E),首先从V中取出源点u0放入最短路径顶点集合U中,这时的最短路径网络S=({u0}, {?}); 然后从u?U和v?V-U中找一条代价最小的边(u*, v*)加入到S中去,此时S=({u0, v*}, {(u0, v*)})。每往U中增加一个顶点,则要对V-U中的各顶点的权值进行一次修正。若加进v*作为中间顶点,使得从u0到其他属于V-U的顶点vi的路径不加v*时最短,则修改u0到vi的权值,即以(u0, v*)的权值加上(v*, vi )的权值来代替原(u0, vi )的权值,否则不修改u0到vi的权值。接着再从权值修正后的V-U中选择最短的边加入S中,如此反复,直到U=V为止。

上面的说明都很抽象,下面图解算法思想:

      原始图为:

     技术分享

寻找最短路径的过程如下:

技术分享

        对第一个图中的有向网络按以上算法思想处理,所求得的从源点F到其余顶点的最短路径的过程如图13.16所示。其中单圆圈表示U中的顶点,而双圆圈表示V-U中的顶点。连接U中两个顶点的有向边用实线表示,连接U和V-U中两个顶点的有向边用虚线表示。圆圈旁的数字为源点到该顶点当前的距离值。
        初始时,S中只有一个源点F,它到V-U中各顶点的路径如图13.16(a)所示;选择图13.16(a)中最小代价边(F, B),同时由于路径(F, A)大于(F, B, A)和(F, C)大于(F, B, C),进行相应调整可得到图13.16(b);选择图13.16(b)中的最小代价边(B, C),同时由于(F, B, A)大于(F, B, C, A),进行相应调整可得到图13.16(c);选择图13.16(c)中最小代价边(C, A)即可得到图13.16(d);选择图13.16(d)中最小代价边(F, D) 即可得到图13.16(e); 最后选择(F, E)即可得到图13.16( f )。

具体的程序实现如下:

 

[cpp] view plaincopy技术分享技术分享
 
  1. #include<stdio.h>  
  2. #define M 12//边数  
  3. #define N 6//顶点数  
  4. #define MAX 10000  
  5.   
  6. void Dijkstra(int v, int dist[][N],int D[N],int p[N],int s[N]) ;  
  7. int flag[N]={0};  
  8. int flag1=0;  
  9. int flag2=0;  
  10. typedef struct  
  11. {  
  12.     int startvex;  
  13.     int endvex;  
  14.     int length;  
  15. }edge;//边的结构体  
  16. edge T[M];  
  17. void main()  
  18. {  
  19.     int dist[N][N]={{0,6,MAX,8,MAX,MAX},//图的邻接矩阵  
  20.                     {18,0,7,MAX,MAX,10},  
  21.                     {9,MAX,0,15,MAX,MAX},  
  22.                     {MAX,MAX,12,0,MAX,MAX},  
  23.                     {MAX,MAX,4,MAX,0,MAX},  
  24.                     {24,5,MAX,25,MAX,0}};  
  25.     int D[N]={0};  
  26.     int p[N]={0};  
  27.     int s[N]={0};  
  28.     int num=0;  
  29.     Dijkstra(5,dist,D, p,s) ;  
  30. }  
  31.   
  32.   
  33.  void Dijkstra(int v, int dist[][N],int D[N],int p[N],int s[N])   
  34.  {  int i, j, k, v1, min, max=10000, pre;   /* Max中的值用以表示dist矩阵中的值? */  
  35.     v1=v;   
  36.     for( i=0; i<N; i++)              /* 各数组进行初始化 */  
  37.     {   D[i]=dist[v1][i];   
  38.         if( D[i] != MAX )  p[i]= v1+1;   
  39.         else p[i]=0;   
  40.         s[i]=0;   
  41.     }  
  42.   
  43.     s[v1]=1;                         /* 将源点送U */  
  44.       for( i=0; i<N-1; i++)   /* 求源点到其余顶点的最短距离 */  
  45.     {   min=10001;    /* min>max, 以保证值为?的顶点也能加入U */  
  46.         for( j=0; j<N-1; j++)  
  47.               if ( ( !s[j] )&&(D[j]<min) )       /* 找出到源点具有最短距离的边 */  
  48.                 {min=D[j];   
  49.                         k=j;   
  50.                     }  
  51.                 s[k]=1;  /* 将找到的顶点k送入U */     
  52.     for(j=0; j<N; j++)  
  53.      if ( (!s[j])&&(D[j]>D[k]+dist[k][j]) ) /* 调整V-U中各顶点的距离值 */  
  54.         {D[j]=D[k]+dist[k][j];   
  55.         p[j]=k+1;                   /* k是j的前趋 */  
  56.                 }  
  57.             }                               /*  所有顶点已扩充到U中 */  
  58.             for( i=0; i<N; i++)  
  59.             {  
  60.                 printf(" %d : %d ", D[i], i);  
  61.                 pre=p[i];   
  62.             while ((pre!=0)&&(pre!=v+1))  
  63.             {   printf ("<- %d ", pre-1);   
  64.                 pre=p[pre-1];   
  65.             }  
  66.             printf("<-%d \n", v);   
  67.         }  
  68. }         

结果显示如下:

 

技术分享

注:如果程序出错,可能是使用的开发平台版本不同,请点击如下链接: 解释说明

 

 

原文:http://blog.csdn.net/tengweitw/article/details/17510891

作者:nineheadedbird

郑重声明:本站内容如果来自互联网及其他传播媒体,其版权均属原媒体及文章作者所有。转载目的在于传递更多信息及用于网络分享,并不代表本站赞同其观点和对其真实性负责,也不构成任何其他建议。