最短路径迪杰斯特拉算法(Dijkstra),用c语言实现

首先,迪杰斯特拉算法是用来解决单源最短路经问题的,主要是通过边的松弛来实现。

我们来看这个问题:

这个问题求得是从1号顶点到达所有其他顶点的最短距离,我们用邻接矩阵来存储这个图,如下:

我们用一个dis数组来存储从一号顶点到其他各个顶点的初始路径,如图

先找一个离一号顶点最近的顶点,通过dis我们知道最近的顶点是2号顶点,从第2个顶点有两条边是2-->3和2-->4。先通过2-->3这个边看能否是从1到3的路程变短,也就是比较dis[3]与dis[2]+map[2][3]的大小,很明显dis[3]=12,而dis[2]+map[2][3]=10,所以我们把dis[3]的值更新为10,这个过程就是我们所说的“松弛”,同样对于2-->4,dis[4]的初始值为无穷大,而dis[2]+dis[2][4]=4,所以我们把dis[4]的值松弛为4,经过这一个松弛后,dis变成了:

把已经找过的点进行标记,然后在剩下的3,4,5,6顶点钟找出离1号顶点最近的顶点,很明显最近的是4,然后根据上面的思路继续进行松弛,一直松弛到边完,这时dis变为

扫描二维码关注公众号,回复: 3438114 查看本文章

整个过程用完整的代码来表示如下:

 
  1. #include <stdio.h>

  2. #include <string.h>

  3. #include <string>

  4. #include <iostream>

  5. #include <stack>

  6. #include <queue>

  7. #include <vector>

  8. #include <algorithm>

  9. #define mem(a,b) memset(a,b,sizeof(a))

  10. using namespace std;

  11. const int inf=1<<29;

  12. int main()

  13. {

  14. int map[10][10],t1,t2,t3,min,u,n,m;

  15. int dis[10];

  16. int vis[10];

  17. scanf("%d%d",&n,&m);

  18. //初始化

  19. for(int i=1; i<=n; i++)

  20. for(int j=1; j<=n; j++)

  21. if(i==j)

  22. map[i][j]=0;

  23. else

  24. map[i][j]=inf;

  25. //读入边

  26. for(int i=1; i<=m; i++)

  27. {

  28. scanf("%d%d%d",&t1,&t2,&t3);

  29. map[t1][t2]=t3;

  30. }

  31. //初始化dis数组,表示1号顶点到其余各个顶点的最短路程

  32. for(int i=1; i<=n; i++)

  33. dis[i]=map[1][i];

  34. //初始化vis

  35. for(int i=1; i<=n; i++)

  36. vis[i]=0;

  37. vis[1]=1;//标记起始点1已经被访问过

  38. //迪杰斯特拉算法(Dijkstra)的核心内容

  39. for(int i=1; i<=n-1; i++) //因为松弛的是边数,所以是n-1

  40. {

  41. min=inf;

  42. for(int j=1; j<=n; j++)

  43. {

  44. if(vis[j]==0&&dis[j]<min)

  45. {

  46. min=dis[j];

  47. u=j;

  48. }

  49. }

  50. vis[u]=1;

  51. for(int v=1; v<=n; v++)

  52. if(map[u][v]<inf)

  53. if(dis[u]+map[u][v]<dis[v])

  54. dis[v]=dis[u]+map[u][v];

  55. }

  56. for(int i=1;i<=n;i++)

  57. printf("%d ",dis[i]);

  58. return 0;

  59. }

给出一组样例进行测试:

输入:

 
  1. 6 9

  2. 1 2 1

  3. 1 3 12

  4. 2 3 9

  5. 2 4 3

  6. 3 5 5

  7. 4 3 4

  8. 4 5 13

  9. 4 6 15

  10. 5 6 4


输出:

0 1 8 4 13 17

利用邻接表,我们可以把时间复杂度优化到O(M+N)logN,以下是用邻接表来优化这个算法的代码:

 
  1. #include <stdio.h>

  2. #include <string.h>

  3. #include <string>

  4. #include <iostream>

  5. #include <stack>

  6. #include <queue>

  7. #include <vector>

  8. #include <algorithm>

  9. #define mem(a,b) memset(a,b,sizeof(a))

  10. using namespace std;

  11. const int inf=0x3f3f3f3f;

  12. int u[10],v[10],w[10],first[10],next[10],dist[10],vis[10],n,m,k,minn;

  13. int main()

  14. {

  15. scanf("%d%d",&n,&m);

  16. //初始化

  17. for(int i=1; i<=n; i++)

  18. first[i]=-1;

  19. //读入边

  20. mem(dist,inf);

  21. for(int i=1; i<=m; i++)

  22. {

  23. scanf("%d%d%d",&u[i],&v[i],&w[i]);

  24. next[i]=first[u[i]];

  25. first[u[i]]=i;

  26. if(u[i]==1)

  27. dist[v[i]]=w[i];//初始化dis数组,表示1号顶点到其余各个顶点的最短路程

  28. }

  29. //初始化vis

  30. for(int i=1; i<=n; i++)

  31. vis[i]=0;

  32. vis[1]=1;//标记起始点1已经被访问过

  33. dist[1]=0;

  34. //迪杰斯特拉算法(Dijkstra)的核心内容

  35. for(int i=1; i<=n-1; i++)

  36. {

  37. minn=inf;

  38. for(int j=1; j<=n; j++)

  39. {

  40. if(vis[j]==0&&dist[j]<minn)

  41. {

  42. minn=dist[j];

  43. k=j;

  44. }

  45. }

  46. vis[k]=1;

  47. for(int l=first[k]; l!=-1; l=next[l])

  48. {

  49. if(w[l]<inf)

  50. if(dist[k]+w[l]<dist[v[l]])

  51. dist[v[l]]=dist[k]+w[l];

  52. }

  53. }

  54. for(int i=1; i<=n; i++)

  55. printf("%d ",dist[i]);

  56. return 0;

  57. }

--------------------- 本文来自 riba2534 的CSDN 博客 ,全文地址请点击:https://blog.csdn.net/riba2534/article/details/54563524?utm_source=copy

猜你喜欢

转载自blog.csdn.net/u014535666/article/details/82870854