当前位置:   article > 正文

求最短路径的三种算法_最短路径算法

最短路径算法

目录

一.单源最短路

1.dijkstra算法及实现

2.spfa算法及实现

(1)spafa负环判断及实现

二.多源最短路

1.floyd算法及实现

一.单源最短路

1.dijkstra算法及实现
求源点到图中其余各顶点的最短路径
dfs效率慢,解决规模小,bfs只能边权为1的图
Dijkstra算法——迪杰斯塔拉算法(非负全图)
 基本思想:
 首先假定源点为u,顶点集合V被划分为两部分:集合 S 和 V-S.
 初始时S中仅含有源点u,其中S中的顶点到源点的最短路径已经确定。
集合S 和V - S中所包含的顶点到源点的最短路径的长度待定,
 称从源点出发只经过S中的点到达V - S中的点的路径为特殊路径,
并用dist[]记录当前每个顶点对应的最短特殊路径长度。
实现从1到n的最短路径
 输入:
3 3
1 2 5
2 3 5
3 1 2
输出:
2

  1. #include<iostream>
  2. #include<cstring>
  3. using namespace std;
  4. const int N = 1001;
  5. const int M = 10001;
  6. int head[N];
  7. int size1;
  8. int n, m;
  9. struct edge {
  10. int v, w, next;
  11. edge() {};
  12. edge(int _v, int _w, int _next) {
  13. v = _v;
  14. w = _w;
  15. next = _next;
  16. }
  17. }e[M*2];
  18. void init() {
  19. memset(head, -1, sizeof(head));
  20. size1 = 0;;
  21. }
  22. void insert(int u,int v,int w) {
  23. e[size1] = edge(v, w, head[u]);
  24. head[u] = size1++;
  25. }
  26. void insert2(int u, int v, int w) {
  27. insert(u, v, w);
  28. insert(v, u, w);
  29. }
  30. int dis[N];
  31. bool vis[N];
  32. void dijkstra(int u)
  33. {
  34. memset(vis, false, sizeof(vis));
  35. memset(dis, 0x3f, sizeof(dis)); //0x3f用以表示正无穷
  36. dis[u] = 0;
  37. int mind = 1000000000;
  38. for (int i = 0; i < n; i++) {
  39. int minj = -1;
  40. for (int j = 1; j <= n; j++) {
  41. if (!vis[j] && dis[j] < mind) {
  42. minj = j;
  43. mind = dis[j];
  44. }
  45. }
  46. if (minj == -1)
  47. return;
  48. vis[minj] = true;
  49. for (int j = head[minj];~j ; j = e[j].next) { ///~j等价于j!=-1;
  50. int v = e[j].v;
  51. int w = e[j].w;
  52. if (!vis[v] && dis[v] > dis[minj] + w)
  53. dis[v] = dis[minj] + w;
  54. }
  55. }
  56. }
  57. int main()
  58. {
  59. init();
  60. int u, v, w;
  61. cin >> n >> m;
  62. while (m--) {
  63. cin >> u >> v >> w;
  64. insert2(u, v, w);
  65. }
  66. dijkstra(1);
  67. cout << dis[n] << endl;
  68. return 0;
  69. }


2.SPFA算法——Shortest Path Faster Algorithm 
思路:
用数组dis记录每个结点的最短路径估计值,用邻接表或邻接矩阵来存储图G。
我们采取的方法是动态逼近法:设立一个先进先出的队列用来保存待优化的结点,
优化时每次取出队首结点u,
并且用u点当前的最短路径估计值对离开u点所指向的结点v进行松弛操作,
如果v点的最短路径估计值有所调整,且v点不在当前的队列中,就将v点放入队尾。
这样不断从队列中取出结点来进行松弛操作,直至队列空为止
优点:
(可以解决带负权的有向图并且在稀疏图优于dijsktra)
反例:
 1.出题特殊数据使SPFA算法慢,2.有负环不能解决
推荐后续优化的dijkstra算法
实现:

  1. #include<iostream>
  2. #include<cstring>
  3. #include<queue>
  4. using namespace std;
  5. const int N = 1001;
  6. const int M = 10001;
  7. int head[N];
  8. int size1;
  9. int n, m;
  10. struct edge {
  11. int v, w, next;
  12. edge() {};
  13. edge(int _v, int _w, int _next) {
  14. v = _v;
  15. w = _w;
  16. next = _next;
  17. }
  18. }e[M*2];
  19. void init() {
  20. memset(head, -1, sizeof(head));
  21. size1 = 0;;
  22. }
  23. void insert(int u,int v,int w) {
  24. e[size1] = edge(v, w, head[u]);
  25. head[u] = size1++;
  26. }
  27. void insert2(int u, int v, int w) {
  28. insert(u, v, w);
  29. insert(v, u, w);
  30. }
  31. int dis[N];
  32. bool vis[N];
  33. void spfa(int u) {
  34. memset(vis, false, sizeof(vis));
  35. vis[u] = true;
  36. memset(dis, 0x3f, sizeof(dis));
  37. dis[u] = 0;
  38. queue<int>q;
  39. q.push(u);
  40. while (!q.empty()) {
  41. u = q.front();
  42. q.pop();
  43. vis[u] = false;
  44. for (int j = head[u]; ~j; j = e[j].next) {
  45. int v = e[j].v;
  46. int w = e[j].w;
  47. if (dis[v] > dis[u] + w) {
  48. dis[v] = dis[u] + w;
  49. if (!vis[v]) {
  50. q.push(v);
  51. vis[v] = true;
  52. }
  53. }
  54. }
  55. }
  56. }
  57. int main()
  58. {
  59. init();
  60. int u, v, w;
  61. cin >> n >> m;
  62. while (m--) {
  63. cin >> u >> v >> w;
  64. insert2(u, v, w);
  65. }
  66. spfa(1);
  67. cout << dis[n] << endl;
  68. return 0;
  69. }


(1)spfa判断负环
在进行spfa用一个数组cnt标记每个顶点的入队次数,如果一个顶点的入队次数大于顶点总数
则表示该图含有负环

  1. #include<iostream>
  2. #include<cstring>
  3. #include<queue>
  4. using namespace std;
  5. const int N = 1001;
  6. const int M = 10001;
  7. int head[N];
  8. int size1;
  9. int n, m;
  10. struct edge {
  11. int v, w, next;
  12. edge() {};
  13. edge(int _v, int _w, int _next) {
  14. v = _v;
  15. w = _w;
  16. next = _next;
  17. }
  18. }e[M*2];
  19. void init() {
  20. memset(head, -1, sizeof(head));
  21. size1 = 0;;
  22. }
  23. void insert(int u,int v,int w) {
  24. e[size1] = edge(v, w, head[u]);
  25. head[u] = size1++;
  26. }
  27. void insert2(int u, int v, int w) {
  28. insert(u, v, w);
  29. insert(v, u, w);
  30. }
  31. int dis[N], in[N];
  32. bool vis[N];
  33. bool spfa(int u) {
  34. memset(vis, false, sizeof(vis));
  35. vis[u] = true;
  36. memset(dis, 0x3f, sizeof(dis));
  37. dis[u] = 0;
  38. memset(in, 0, sizeof(in));
  39. in[u] = u;
  40. queue<int>q;
  41. q.push(u);
  42. while (!q.empty()) {
  43. u = q.front();
  44. q.pop();
  45. vis[u] = false;
  46. for (int j = head[u]; ~j; j = e[j].next) {
  47. int v = e[j].v;
  48. int w = e[j].w;
  49. if (dis[v] > dis[u] + w)
  50. dis[v] = dis[u] + w;
  51. if (!vis[v]) {
  52. q.push(v);
  53. vis[v] = false;
  54. ++in[v];
  55. if (in[v] > n)
  56. return true;
  57. }
  58. }
  59. }
  60. return false;
  61. }
  62. int main()
  63. {
  64. init();
  65. int u, v, w;
  66. cin >> n >> m;
  67. while (m--) {
  68. cin >> u >> v >> w;
  69. insert2(u, v, w);
  70. }
  71. if (spfa(1))
  72. cout << "Yes" << endl;
  73. else
  74. cout << "No" << endl;
  75. return 0;
  76. }


二.多源最短路
1.floyd算法——(Floyd-Warshall algorithm),中文亦称弗洛伊德算法
利用动态规划的思想解决带权图中任意两个点之间最短路径算法
优势:代码简短,高效,可以解决带负权

反例:不能解决带负环
思路:
dp[0][i][j]=图
dp[k][i][j]=从i到j可能经过1...k最短路径,
min(dp[k][i][j],dp[k-1][i][k]+dp[k-1][k][j])
优化为二维:
min(dp[i][j],dp[i][k]+dp[k][j])
实现:
 

  1. #include<iostream>
  2. #include<cstring>
  3. using namespace std;
  4. const int N = 101;
  5. int g[N][N];
  6. void floyd(int n) {
  7. for (int k = 1; k <= n; k++) {
  8. for (int i = 1; i <= n; i++) {
  9. for (int j = 1; j <= n; j++) {
  10. g[i][j] = min(g[i][j], g[i][k] + g[j][k]);
  11. }
  12. }
  13. }
  14. }
  15. int main()
  16. {
  17. memset(g, 0x3f, sizeof(g));
  18. for (int i = 0; i < N; i++) {
  19. g[i][i] = 0;
  20. }
  21. int n, m;
  22. int u, v, w;
  23. cin >> n >> m;
  24. while (m--) {
  25. cin >> u >> v >> w;
  26. g[u][v] = g[v][u] = w;
  27. }
  28. floyd(n);
  29. for (int i = 1; i <= n; i++) {
  30. for (int j = 1; j <= n; j++) {
  31. cout << g[i][j] << ' ';
  32. }
  33. cout << endl;
  34. }
  35. return 0;
  36. }

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/你好赵伟/article/detail/710956
推荐阅读
相关标签
  

闽ICP备14008679号