(原创)最短路径-Dijkstra算法,以Til the Cows Come Home为例

六月ゝ 毕业季﹏ 提交于 2020-10-18 02:44:10

 

 (1)首先先解释一下单源最短路径:

    1)容易的解释:指定一个点(源点)到其余各个顶点的最短路径,也叫做单源最短路径

    2)官方解释:给定一个带权有向图G=(V,E),其中每条边的权是一个实数。另外,还给定V中的一个顶点,称为源。现在要计算从源到其他所有各顶点的最短路径长度。这里的长度就是指路上各边权之和。这个问题通常称为单源最短路径问题。

 

2)解释一下Dijkstra算法:

 

 

例如求A点到B、C、D、E、F顶点的最短路径;

 

 

 

我们可以先这样设想:

1)先把所有的点到另一个点的长度全部初始化为无穷大本身到本身则初始化为0再输入数值;

第一:将所有点到点初始化为无穷大

 

 

 

代码大致如下:

1 const int INF = 0x3f3f3f3f;  //为无穷大;
 2 int G[2000][2000];
 3 int N  ;//N为点的个数;
 4 
 5         for( i = 1 ;i <= N ;i++)
 6         {
 7             for( j = 1 ;j <= N ;j++)
 8             {
 9                 G[i][j] = INF;
10             }
11         }

 

也可用下面这种:

1 #include<string.h>
2 const int INF = 0x3f3f3f3f;
3 int N; //点的个数;
4 int G[2000][2000];
5 memset(G,INF,sizeof(G));

 

 

第二:将本身到本身初始化为0 因为本身到本身的距离就为0

 

 

1 for(int i = 1 ;i <= N; i++)
2 {
3    G[i][i] = 0;    //本身到本身距离为0;
4 }

 

第三:输入数据:

 

 

 

 

 

对应下面的表格:

 

 

代码实现大概如下:

 

1 int M;        //M为边的个数;
 2 int x  , y;  // x ,y 为同一边的两个点;
 3 int D;    //D 为题目给的边的权值;
 4         for( i = 1; i <= M ;i++)
 5         {
 6             cin>>x>>y>>D;
 7            
 8                 if(G[x][y]>D)
 9                 {
10                     G[x][y] = D;
11                     G[y][x] = D;
12                 }
13            
14         }
15
16

 

17 //上面我们已经把点到其他点的距离初始化为无穷大,把本身到本身初始化为0;
   //那么上面这个循环的化,我们可以把题目给的两点之间的权值输入;

18 //题目给的权值一定为正值,所以比0大,故本身到本身的距离还是维持为0;
   //而小于无穷大,所以用上面那个循环可把边的权值输入;

   2)我们还需用一个一维数组d来存储A顶点到其余各个顶点的初始路程

 

 

下面我们来模拟一下:

这是输入最初数据的表格:

 

模拟:

 

 

1)既然是从A点到其余各个顶点的最短路径,那就先找一个离A号顶点最近的顶点。  AB的距离最小

 

      

       所以d[2]的值就已经从估计值变成了确定值,目前离A顶点最近的是B顶点。

 

2既然选了B点,接下来看B顶点有哪些出边呢,B->EB->C两条边。

 

 

先讨论B->E能否让A顶点到E顶点的路程变短;

A->B->E :d[2]+e[2][5]表示从ABE的距离;d[2]AB顶点的路程,e[2][5]表示BE的距离;  12+7 = 19;

A->E   :d[5] = 16;

我们发现 d[2]+e[2][5]=19d[5] = 16;    所以  d[5]  < d[2]+e[2][5];

所以AE目前的最短距离为16

     d    

 

 

因此d[5]更新为 16;(这个过程有个专业术语叫做:松弛);

d

 

 

 再判断B->C

A->B->C:d[2]+e[2][3] 表示ABC的距离;d[2]AB顶点的路程,e[2][3]表示BC 的距离;12+10 = 22

A->C:A没有直接到C ,所以为无穷大;

22<无穷大

所以AC目前的最短距离为22d[3] = 22;

 

 

这样从AB的当前最短路径就更新完了;

3)找除B外离A最近的点:由图可知是F点;

 

 

 

 

 由图可知F的出边有两条:F->E, F->D

先讨论 F->E能否让AE的距离变短;

A->F->E d[6]+e[6][5] = 14+ 9 = 23;

A->E :d[5] = 16;

所以d[5] < d[6]+e[6][5];

所以当前AE的最短距离为16  d[5] = 16;

 

 

 

再讨论E->D的距离能否让AD的距离变短

A->E->D  :   d[5]+e[5][4] = 16+2=18;

A->D   :   无穷大;

所以AD的当前最短距离更新为18 d[4] = 18;

 

 

 

4)找除BF外离A最近的点E

 

 

 

 E的出边有CD,看E->CE->D能否使AC的距离变短,使AD的距离变短

先讨论E->D

A->E->D d[5]+e[5][4] = 16+2 = 18

之前更新后的AD d[4] = 22;

所以AD 当前最短变为18 d[4] = 18;

 

 

 讨论E->C,看是否让AC变短

A->E->C d[5]+e[5][3] = 16+6 =22;

原更新的AC22

所以当前AC最短为22  d[3] = 22;

 

 

5)找除B,E,F外离A最近的点,D

 

 

 

 D的出边只有D->C

D->C,能否让AC变短;

A->D->C  d[4]+e[4][3] = 18+5 = 23;

AC   22

 所以AC当前最短为22  d[3] = 22;

 

 

 

6)最后剩下C

全部找完,确定了A到各个点的最短路径;

下面以一道题为例

Bessie is out in the field and wants to get back to the barn to get as much sleep as possible before Farmer John wakes her for the morning milking. Bessie needs her beauty sleep, so she wants to get back as quickly as possible.


Farmer John's field has N (2 <= N <= 1000) landmarks in it, uniquely numbered 1..N. Landmark 1 is the barn; the apple tree grove in which Bessie stands all day is landmark N. Cows travel in the field using T (1 <= T <= 2000) bidirectional cow-trails of various lengths between the landmarks. Bessie is not confident of her navigation ability, so she always stays on a trail from its start to its end once she starts it.

Given the trails between the landmarks, determine the minimum distance Bessie must walk to get back to the barn. It is guaranteed that some such route exists.


Input

* Line 1: Two integers: T and N

* Lines 2..T+1: Each line describes a trail as three space-separated integers. The first two integers are the landmarks between which the trail travels. The third integer is the length of the trail, range 1..100.

Output

* Line 1: A single integer, the minimum distance that Bessie must travel to get from landmark N to landmark 1.

 这道题的题意就是:给定N个地点(编号从1N),T条道路,先输入两个正整数TN,接下来T行,每行三个整数,分别代表这条道路的起点和终点和长度,要求求出从N1的最短路径;此题用Dijkstra算法(先懂得Dijkstra算法才能知道这道题怎么做)

代码如下:  

解法一:未优化的Dijkstra ,时间复杂度为N²;

以上面的思想写就好了,先找出最小,不断更新“当前最小”,并注意标记哪些点已经被访问过了,就是上面被涂红的点;

 

 1 #include<iostream>
 2 #include<string.h>
 3 #include<stdio.h>
 4 using namespace std;
 5 
 6 
 7 
 8 const int INF = 0x3f3f3f3f;
 9 int N ,M;
10 int x , y;
11 int min1;
12 int D;
13 int flag ;
14 int vis[2005];          //用来记录哪些“最近的点”被访问过了;
15 int d[2005];           //用来记录“当前最短距离”;
16 int G[2005][2005];        //输入点与点间的距离;
17 int main()
18 {
19     while(scanf("%d%d",&M,&N)!=EOF)  //输入边的条数和点的个数;
20     {
21             for(int i = 1 ;i <= N ;i++)
22         {
23             for(int j = 1 ; j <= N ;j++)
24             {
25                 G[i][j] = INF;           //先将全部点到点的距离初始化为无穷大;
26             }
27         }
28     for(int i  = 1 ; i <= N ;i++)
29     {
30         G[i][i] = 0;                  //本身点到本身点距离初始化为0 ;
31     }
32     for(int i = 1 ; i <= M ;i++)
33     {
34         cin>>x>>y>>D;               //输入题目给的点和两点间的距离;
35         if(G[x][y]>D)    
36         {
37             G[x][y] = D;          //输入数据;注意这是有向图,所以要正反来一遍;有一些题目从A→B,和从B→A的距离不一样;此题中A→B和B→A一样;
38             G[y][x] = D;
39         }
40     }
41     memset(d,INF,sizeof(d));       //就是上面思路中的一位数组d ,用来记录“当前最短距离”;每次更新;
42     d[1] = 0 ;          //先初始化d[1] = 0 ;起点1 到1 的距离为0;
43     
44     for(int i = 1 ; i <= N; i++)
45     {
46         min1 = INF;                      //每次找该点与其他点的最短距离时,先将min1初始化为无穷大;每次循环找新的点到其余点时,就需要min1初始化为INF;
如找i = 1 时与其连通距离最短的点,找到后min1 = d[1]; 后循环后找 i = 2 时与其连通距离最短的点时,应先将min1 再次赋为INF,找到的d[2]才对;
47 for(int j = 1 ; j<= N; j++) 48 { 49 if(min1>d[j]) 50 { 5152 if(vis[j]==0) //这是判断没有”被访问过”的点;即还没有被找到的“最近距离”的点;即上面思路中还没被涂红的点; 53 { 54 min1 = d[j]; //找到最小; 55 flag = j; //记录该序号; 56 } 57 58 } 59 } 60 for(int k = 1 ; k <= N ;k++) 61 { 62 if(d[k]>d[flag]+G[flag][k]) //不断更新“最近的点”; 63 d[k] = d[flag]+G[flag][k]; 64 } 65 vis[flag] = 1; //标记该点被访问过了;即图中“涂红”点 66 67 } 68 cout<<d[N]<<endl; //直接输出最短距离; 69 for(int i = 1 ; i <= N; i++) //多组输入的话,则应将所有“涂红的点”恢复;重新开始新的数据; 70 { 71 vis[i] = 0; 72 } 73 } 74 75 76 }

 

解法二:

 

下面是优化了的写法,可以用优先队列取优化;每次取最小的时候可以用优先队列去取,时间复杂度优化到了NlogN;

 

 1 #include<iostream>
 2 #include<stdio.h>
 3 #include<queue>
 4 #include<string.h>
 5 using namespace std ;
 6 
 7 const int INF = 0x3f3f3f3f;
 8 int G[2000][2000];
 9 int d[2000];
10 
11 int i ,j;
12 
13 struct node{
14     int num;              //记录结点的序号;
15     int dis;            //记录结点的“最短距离”
16     friend bool operator<(node a ,node b)
17     {
18         return a.dis>b.dis;//运算符重载,优先队列原本是从大到小输出,现在改写后变成从小到大;
19     }
20 };
21 
22 int main()
23 {
24     int M , N;
25     int x,y,D;
26     
27     priority_queue<node>que;//将等下的“最短路径”入队;
28     
29     
30     while(scanf("%d%d",&M,&N)!=EOF)
31     {
32         for( i = 1 ;i <= N ;i++)
33         {
34             for( j = 1 ;j <= N ;j++)
35             {
36                 G[i][j] = INF;//先全部初始化为无穷大;
37             }
38         }
39         
40             for(int i  = 1 ; i <= N ;i++)
41         {
42             G[i][i] = 0;
43         }
44         for( i = 1; i <= M ;i++)
45         {
46             scanf("%d%d%d",&x,&y,&D); //输入数据;
47             
48                 if(G[x][y]>D)//如果比无穷大小就更新距离;即有数据的就有数据,无直接连通就是无穷大;
49                 {
50                     G[x][y] = D;
51                     G[y][x] = D;
52                 }
53             
54         }
55         memset(d,0x3f,sizeof(d));
56         d[1] = 0;
57         que.push({1,0}); //将起点入队,起点和起点的距离是为0,因为是本身;
58         while(!que.empty())
59         {
60             node tp = que.top(); //每次取出队列最前的数,即最小的数;
61 
62             que.pop();
63             
64             for(i = 1 ;i <= N ;i++)
65             {
66                 if(G[tp.num][i])//如果两点间有距离
67                 {
68                     if(d[i]>d[tp.num]+G[tp.num][i])//每次更新“当前最短距离”
69                     {
70                         d[i] = d[tp.num] + G[tp.num][i];
71                         que.push({i,d[i]});//入队;
72                     }
73                 }
74             }
75         }
76         
77         printf("%d\n",d[N]); //直接输出1到N的最短距离;
78         while(!que.empty())  //注意队列要清空;
79         {
80             que.pop();
81         }
82         
83         
84     }
85     return 0;
86 }

 

 

 

 

 

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!