图之Dijkstra算法 图之Dijkstra算法 最短路径之Dijkstra算法 =========================

========================================================================

       这主要是用来找起点与终点最短路径的方法,搜索是以目标节点全遍历的方式进行搜索,一步步确定每个节点的最短路径才终止.

       操作过程中分两集合:  确定最短路径顶点集合U 与 未确定最短路径顶点集合V.  利用U里面的有向向量链分别搜索V里面的每个节点,形成的最短链对应的那个节点就可以确定最短路径了,其就可以加入集合U,逐步搜索整个集合V,直到V中数量减小到零,全部确定.

======================================================================== 

   先上基本搜索步骤图对算法理解.转自:图之Dijkstra算法   

   图之Dijkstra算法

     Dijkstra算法是一种求单源最短路算法,即从一个点开始到所有其他点的最短路。其步骤如下: (原文链接有参考代码.)


1

2
 

3

4
 

5

6
 


    既然理解了搜索步骤,那我们看看其中的思想.  参见: 最短路径之Dijkstra算法详细讲解

最短路径之Dijkstra算法详细讲解

1  最短路径算法

在日常生活中,我们如果需要常常往返A地区和B地区之间,我们最希望知道的可能是从A地区到B地区间的众多路径中,那一条路径的路途最短。最短路径问题是图论研究中的一个经典算法问题, 旨在寻找图(由结点和路径组成的)中两结点之间的最短路径。 算法具体的形式包括:

(1)确定起点的最短路径问题:即已知起始结点,求最短路径的问题。

(2)确定终点的最短路径问题:与确定起点的问题相反,该问题是已知终结结点,求最短路径的问题。在无向图中该问题与确定起点的问题完全等同,在有向图中该问题等同于把所有路径方向反转的确定起点的问题。

(3)确定起点终点的最短路径问题:即已知起点和终点,求两结点之间的最短路径。

(4)全局最短路径问题:求图中所有的最短路径。

用于解决最短路径问题的算法被称做“最短路径算法”, 有时被简称作“路径算法”。 最常用的路径算法有:Dijkstra算法、A*算法、Bellman-Ford算法、Floyd-Warshall算法、Johnson算法。 

2  Dijkstra算法

2.1  Dijkstra算法

  Dijkstra算法是典型最短路算法,用于计算一个节点到其他所有节点的最短路径。主要特点是以起始点为中心向外层层扩展,直到扩展到终点为止。Dijkstra算法能得出最短路径的最优解,但由于它遍历计算的节点很多,所以效率低

Dijkstra算法是很有代表性的最短路算法,在很多专业课程中都作为基本内容有详细的介绍,如数据结构,图论,运筹学等等。 

2.2  Dijkstra算法思想

Dijkstra算法思想为:

           设G=(V,E)是一个带权有向图,把图中顶点集合V分成两组,

         第一组为已求出最短路径的顶点集合(用S表示,初始时S中只有一个源点,以后每求得一条最短路径 , 就将 加入到集合S中,直到全部顶点都加入到S中,算法就结束了),

          第二组为其余未确定最短路径的顶点集合(用U表示),按最短路径长度的递增次序依次把第二组的顶点加入S中。在加入的过程中,总保持从源点v到S中各顶点的最短路径长度不大于从源点v到U中任何顶点的最短路径长度。

            此外,每个顶点对应一个距离,S中的顶点的距离就是从v到此顶点的最短路径长度,U中的顶点的距离,是从v到此顶点只包括S中的顶点为中间顶点的当前最短路径长度。

2.3  Dijkstra算法具体步骤  

1)初始时,S只包含源点,即S=,v的距离为0。U包含除v外的其他顶点,U中顶点u距离为边上的权(若v与u有边)或 )(若u不是v的出边邻接点)。

(2)从U中选取一个距离v最小的顶点k,把k,加入S中(该选定的距离就是v到k的最短路径长度)。

(3)以k为新考虑的中间点,修改U中各顶点的距离;若从源点v到顶点u(u U)的距离(经过顶点k)比原来距离(不经过顶点k)短,则修改顶点u的距离值,修改后的距离值的顶点k的距离加上边上的权。

(4)重复步骤(2)和(3)直到所有顶点都包含在S中。

2.4  Dijkstra算法举例说明

如下图,设A为源点,求A到其他各顶点(B、C、D、E、F)的最短路径。线上所标注为相邻线段之间的距离,即权值。(注:此图为随意所画,其相邻顶点间的距离与图中的目视长度不能一一对等)

图一:Dijkstra无向图

 最短路径之Dijkstra算法详细讲解 - 绿岩 - 永远的绿岩

算法执行步骤如下表:【注:图片要是看不到请到“相册--日志相册”中,名为“Dijkstra算法过程”的图就是了】



=========================  转自最短路径之Dijkstra算法 =========================

dijkstra算法两个应用题:

HDOJ 1874 畅通工程续,现有解法:www.wutianqi.com/?p=1894
HDOJ 2544 最短路,现有解法:www.wutianqi.com/?p=1892

参考:http://hi.baidu.com/zealot886/item/c8a499ee5795bcddeb34c950

           数据结构(C语言版)

          http://www.cnblogs.com/biyeymyhjob/archive/2012/07/31/2615833.html

推荐几篇搜索算法相关的非常好的博文:

一、A*搜索算法

一(续)、A*,Dijkstra,BFS算法性能比较及A*算法的应用

二、Dijkstra 算法初探          (Dijkstra算法系列4篇文章)

二(续)、彻底理解Dijkstra算法

二(再续)、Dijkstra 算法+fibonacci堆的逐步c实现

二(三续)、Dijkstra 算法+Heap堆的完整c实现源码



==============

代码

 

[cpp]  view plain  copy
  1. #include <stdio.h>  
  2. #include <malloc.h>  
  3. #define VERTEXNUM 6  
  4.   
  5. //存放最短路径的边元素  
  6. typedef struct edge{  
  7.         int vertex;  
  8.         int value;  
  9.         struct edge* next;  
  10. }st_edge;  
  11.   
  12. void createGraph(int (*edge)[VERTEXNUM], int start, int end, int value);  
  13. void displayGraph(int (*edge)[VERTEXNUM]);  
  14. void displayPath(st_edge** path, int startVertex,int* shortestPath);  
  15. void dijkstra(int (*edge)[VERTEXNUM], st_edge*** path, int** shortestPath, int startVertex, int* vertexStatusArr);  
  16. int getDistance(int value, int startVertex, int start, int* shortestPath);  
  17. void createPath(st_edge **path, int startVertex, int start, int end, int edgeValue);  
  18.   
  19. int main(void){  
  20.         //动态创建存放边的二维数组  
  21.         int (*edge)[VERTEXNUM] = (int (*)[VERTEXNUM])malloc(sizeof(int)*VERTEXNUM*VERTEXNUM);  
  22.         int i,j;  
  23.         for(i=0;i<VERTEXNUM;i++){  
  24.                 for(j=0;j<VERTEXNUM;j++){  
  25.                         edge[i][j] = 0;  
  26.                 }  
  27.         }  
  28.         //存放顶点的遍历状态,0:未遍历,1:已遍历  
  29.         int* vertexStatusArr = (int*)malloc(sizeof(int)*VERTEXNUM);  
  30.         for(i=0;i<VERTEXNUM;i++){  
  31.                 vertexStatusArr[i] = 0;  
  32.         }  
  33.   
  34.         printf("after init:\n");  
  35.         displayGraph(edge);  
  36.         //创建图  
  37.         createGraph(edge,0,1,6);  
  38.         createGraph(edge,0,3,5);  
  39.         createGraph(edge,0,2,1);  
  40.         createGraph(edge,1,2,5);  
  41.         createGraph(edge,1,4,3);  
  42.         createGraph(edge,2,4,6);  
  43.         createGraph(edge,2,3,5);  
  44.         createGraph(edge,2,5,4);  
  45.         createGraph(edge,3,5,2);  
  46.         createGraph(edge,4,5,6);  
  47.   
  48.         printf("after create:\n");  
  49.         displayGraph(edge);  
  50.     //最短路径  
  51.         /*存储的结构如下: 
  52.             path[0]:edge0->NULL 
  53.             path[1]:edge1->NULL 
  54.             path[2]:edge1->edge2->NULL 
  55.             path[3]:edge1->edge2->edge3->NULL 
  56.             path[4]:edge4->NULL 
  57.             从顶点0到0的最短路径:从0出发直接到0 
  58.             从顶点0到1的最短路径:从0出发直接到1 
  59.             从顶点0到2的最短路径:从0出发到1,从1出发到2 
  60.             ...... 
  61.         */  
  62.     st_edge** path = NULL;  
  63.     //存储最短路径的权值  
  64.         /* 
  65.         shortestPath[0] = 0; 
  66.         shortestPath[1] = 8; 
  67.         shortestPath[2] = 12; 
  68.         从顶点0到0的路径是0 
  69.         从顶点0到1的路径是8 
  70.         从顶点0到2的路径是12 
  71.         */  
  72.     int* shortestPath = NULL;  
  73.     //从顶点0开始寻找最短路径  
  74.     int startVertex = 0;  
  75.     //最短路径  
  76.     dijkstra(edge, &path, &shortestPath, startVertex, vertexStatusArr);  
  77.     printf("the path is:\n");  
  78.     displayPath(path,startVertex,shortestPath);  
  79.   
  80.         free(edge);  
  81.         free(path);  
  82.         return 0;  
  83. }  
  84. //创建图  
  85. void createGraph(int (*edge)[VERTEXNUM], int start, int end, int value){  
  86.         edge[start][end] = value;  
  87.         edge[end][start] = value;  
  88. }  
  89. //打印存储的图  
  90. void displayGraph(int (*edge)[VERTEXNUM]){  
  91.         int i,j;  
  92.         for(i=0;i<VERTEXNUM;i++){  
  93.                 for(j=0;j<VERTEXNUM;j++){  
  94.                         printf("%d ",edge[i][j]);  
  95.                 }  
  96.                 printf("\n");  
  97.         }  
  98. }  
  99. //打印最短路径  
  100. void displayPath(st_edge** path, int startVertex,int* shortestPath){  
  101.         int i;  
  102.         st_edge* p;  
  103.         for(i=0;i<VERTEXNUM;i++){  
  104.                 printf("Path from %d to %d:",startVertex,i);  
  105.                 p = *(path+i);  
  106.                 while(p != NULL){  
  107.                         printf("%d(%d) ",p->vertex,p->value);  
  108.                         p = p->next;  
  109.                 }  
  110.                 printf("\n");  
  111.         printf("the count is:%d\n",shortestPath[i]);  
  112.         }  
  113. }  
  114. //最短路径  
  115. void dijkstra(int (*edge)[VERTEXNUM], st_edge*** path, int** shortestPath, int startVertex, int* vertexStatusArr){  
  116.     //初始化最短路径  
  117.     *path = (st_edge**)malloc(sizeof(st_edge*)*VERTEXNUM);  
  118.         int i,j;  
  119.     for(i=0;i<VERTEXNUM;i++){  
  120.         if(i == startVertex){  
  121.             st_edge* e = (st_edge*)malloc(sizeof(st_edge));  
  122.             e->vertex = startVertex;  
  123.             e->value = 0;  
  124.             e->next = NULL;  
  125.             (*path)[i] = e;  
  126.         }else{  
  127.             (*path)[i] = NULL;  
  128.         }  
  129.     }  
  130.     //初始化最短路径的权值  
  131.     *shortestPath = (int *)malloc(sizeof(int)*VERTEXNUM);  
  132.     for(i=0;i<VERTEXNUM;i++){  
  133.         if(i == startVertex){  
  134.             (*shortestPath)[i] = 0;  
  135.         }else{  
  136.             (*shortestPath)[i] = -1;  
  137.         }  
  138.     }  
  139.     //从顶点0开始,则顶点0就是已访问的  
  140.     vertexStatusArr[startVertex] = 1;  
  141.   
  142.     int shortest, distance,start, end, edgeValue, vNum = 1;  
  143.         //如果还顶点还没有访问完  
  144.         while(vNum < VERTEXNUM){  
  145.                 shortest = 9999;  
  146.                 for(i=0;i<VERTEXNUM;i++){  
  147.                         //选择已经访问过的点  
  148.                         if(vertexStatusArr[i] == 1){  
  149.                                 for(j=0;j<VERTEXNUM;j++){  
  150.                                         //选择一个没有访问过的点  
  151.                                         if(vertexStatusArr[j] == 0){  
  152.                                                 //选出一条value最小的边  
  153.                                                 if(edge[i][j] != 0 && (distance = getDistance(edge[i][j], startVertex, i,  *shortestPath)) < shortest){  
  154.                                                         shortest = distance;  
  155.                                                         edgeValue = edge[i][j];  
  156.                                                         start = i;  
  157.                                                         end = j;  
  158.                                                 }  
  159.                                         }  
  160.                                 }  
  161.                         }  
  162.                 }  
  163.                 vNum++;  
  164.             //将点设置为访问过  
  165.             vertexStatusArr[end] = 1;  
  166.             //保存最短路径权值  
  167.             (*shortestPath)[end] = shortest;  
  168.             //保存最短路径  
  169.             createPath(*path, startVertex, start, end, edgeValue);  
  170.         }  
  171. }  
  172.   
  173. //返回从startVertex到新的顶点的距离  
  174. int getDistance(int value, int startVertex, int start, int* shortestPath){  
  175.     if(start == startVertex){  
  176.         return value;  
  177.     }else{  
  178.         return shortestPath[start] + value;  
  179.     }  
  180. }  
  181.   
  182. //保存最短路径  
  183. void createPath(st_edge **path, int startVertex, int start, int end, int edgeValue){  
  184.     if(start == startVertex){  
  185.         st_edge* newEdge = (st_edge*)malloc(sizeof(st_edge));  
  186.         newEdge->vertex = end;  
  187.         newEdge->value = edgeValue;  
  188.         newEdge->next = NULL;  
  189.   
  190.         st_edge** p = path + end;  
  191.         while((*p) != NULL){  
  192.             p = &((*p)->next);  
  193.         }  
  194.         *p = newEdge;  
  195.     }else{  
  196.         st_edge** pCopySrc = path + start;  
  197.         st_edge** pCopyDes = path + end;  
  198.         st_edge* newEdge = NULL;  
  199.         while((*pCopySrc) != NULL){  
  200.             newEdge = (st_edge*)malloc(sizeof(st_edge));  
  201.             *newEdge = **pCopySrc;  
  202.             newEdge->next = NULL;  
  203.             *pCopyDes = newEdge;  
  204.             pCopySrc = &((*pCopySrc)->next);  
  205.             pCopyDes = &((*pCopyDes)->next);  
  206.         }  
  207.         newEdge = (st_edge*)malloc(sizeof(st_edge));  
  208.         newEdge->vertex = end;  
  209.         newEdge->value = edgeValue;  
  210.         newEdge->next = NULL;  
  211.         *pCopyDes = newEdge;  
  212.     }  
  213. }  

猜你喜欢

转载自blog.csdn.net/wangjingqi930330/article/details/78515346