当前位置:   article > 正文

单源无权最短路径_无权最短路径-广度优先搜索算法

无权最短路径-广度优先搜索算法

单源无权最短路径是指对一个无权图G=(V,E),给定一个顶点s作为输入(单源),从该顶点出发分别到各个顶点的最短无权路径。最短无权路径可以看成是最短赋权路径的特殊情况,只是所有的边的权值都为1.

以下图为例:

选择v_{3}为输入顶点,从v_{3}开始计算每一个顶点的最短无权路径。

显然,从v_{3}v_{3}的距离为0,把这个信息作个标记:

接下来找到与v_{3}距离为1的顶点,这些顶点可以通过考察v_{3}的邻接顶点找到:

 

接下来找到所有与v_{3}距离为3的顶点,这些点可以通过考察与v_{3}距离为1的顶点的邻接顶点找到:

 接下来找到与v_{3}距离为3的顶点:

 

到此就得到了所有顶点的最短无权路径,这种搜索一个图的方法称为广度优先搜索,这种方法按照层处理顶点,首先处理输入顶点的邻接顶点,再处理邻接顶点的邻接顶点,很像树的层序遍历。

代码实现:

  1. void Unweighted(g* p) {
  2. int currdis;
  3. for (currdis = 0; currdis < 7; currdis++) {//当前处理的顶点的距离
  4. for (int i = 0; i < 7; i++) {//遍历找到未知的对应距离元素
  5. if (p->v[i]->known == 0 && p->v[i]->d == currdis) {//如果元素未知且距离相符就处理
  6. p->v[i]->known = 1;//将该元素标记为已知
  7. l* tmp = p->v[i]->next;
  8. while (tmp != NULL) {//处理其邻接顶点
  9. if (p->v[tmp->val]->d == inf) {//如果邻接顶点的距离没有被更新才进行更新
  10. p->v[tmp->val]->d = currdis + 1;
  11. p->v[tmp->val]->path = i;//路径中的上一个顶点
  12. }
  13. tmp = tmp->next;
  14. }
  15. }
  16. }
  17. }
  18. }

 这个算法包含了两个for循环,所以算法的时间复杂度为O(|V|^2),由于在循环过程中,不断地考察了一些已知的顶点和一些已经处理过的未知顶点,所以算法的效率较低。我们可以参考拓扑排序的改进来对这个算法进行改进,我们希望在处理当前定点后处理那些未知的并且距离已经赋值的顶点,而不对已知的或是未被处理过的顶点进行考察,所以可以使用一个队列来存放将要被处理的顶点。首先将输入顶点入队,当队列不为空时,队头元素出队,并对队头元素的邻接顶点进行处理,如果邻接顶点是未知的,则改变它的距离并入队,重复这些工作直到队列为空。以上图为例,下面用表格来说明这个过程:

v初始状态
Knowndispath
v_{1}0inf-1
v_{2}0inf-1
v_{3}00-1
v_{4}0inf-1
v_{5}0inf-1
v_{6}0inf-1
v_{7}0inf-1
Queuev_{3}
vv_{3}出队
Knowndispath
v_{1}01

  v_{3}

v_{2}0inf-1
v_{3}10-1
v_{4}0inf-1
v_{5}0inf-1
v_{6}01

 v_{3}

v_{7}0inf-1
Queuev_{1}v_{6}

 

vv_{1}出队
Knowndispath
v_{1}11v_{3}
v_{2}02v_{1}
v_{3}10-1
v_{4}02v_{1}
v_{5}0inf-1
v_{6}01v_{3}
v_{7}0inf-1
Queuev_{6}v_{2}v_{4}
vv_{6}出队
Knowndispath
v_{1}11v_{3}
v_{2}02v_{1}
v_{3}10-1
v_{4}02v_{1}
v_{5}0inf-1
v_{6}11v_{3}
v_{7}0inf-1
Queuev_{2}v_{4}
vv_{2}出队
Knowndispath
v_{1}11v_{3}
v_{2}12v_{1}
v_{3}10-1
v_{4}02v_{1}
v_{5}03v_{2}
v_{6}11v_{3}
v_{7}0inf-1
Queuev_{4}v_{5}
vv_{4}出队
Knowndispath
v_{1}11v_{3}
v_{2}12v_{1}
v_{3}10-1
v_{4}12v_{1}
v_{5}03v_{2}
v_{6}11v_{3}
v_{7}03v_{4}
Queuev_{5}v_{7}
vv_{5}出队
Knowndispath
v_{1}11v_{3}
v_{2}02v_{1}
v_{3}10-1
v_{4}02v_{1}
v_{5}03v_{2}
v_{6}01v_{3}
v_{7}03v_{4}
Queuev_{7}
vv_{7}出队
Knowndispath
v_{1}11v_{3}
v_{2}02v_{1}
v_{3}10-1
v_{4}02v_{1}
v_{5}03v_{2}
v_{6}01v_{3}
v_{7}03v_{4}
Queueempty

 代码实现:

  1. void UnweightedQ(g* p) {
  2. q* pq = CreatQueue();
  3. enqueue(pq, 2);//首先将输入顶点入队
  4. while (!isempty(pq)) {//队列不为空就出队
  5. int n = dequeue(pq);
  6. p->v[n]->known = 1;
  7. l* tmp = p->v[n]->next;
  8. while (tmp != NULL) {//对出队顶点的邻接顶点进行考察,如果没有被处理过,就进行处理并入队
  9. if (p->v[tmp->val]->d == inf) {
  10. p->v[tmp->val]->d = p->v[n]->d + 1;
  11. p->v[tmp->val]->path = n;
  12. enqueue(pq, tmp->val);
  13. }
  14. tmp = tmp->next;
  15. }
  16. }
  17. }

 图的定义和测试代码:

  1. #define inf 99999999//定义无穷大来表示未到达的顶点的距离
  2. typedef struct list {//邻接表,用来存放该顶点的邻接顶点的编号
  3. int val;
  4. struct list* next;
  5. }l;
  6. typedef struct table {//用来存放各个顶点的信息
  7. int known;//顶点是否已知
  8. int d;//顶点到输入顶点的距离
  9. int path;//顶点在路径中的上一个顶点
  10. l* next;//指向邻接表的指针
  11. }t;
  12. typedef struct graph {//图,存放指向七个顶点信息的指针
  13. t* v[7];
  14. }g;
  15. g* CreatGraph() {//用暴力方法建立一个图
  16. g* pg = (g*)malloc(sizeof(g));
  17. for (int i = 0; i < 7; i++) {//在程序中,下标为0处代表编号为1的顶点
  18. t* p = (t*)malloc(sizeof(l));
  19. p->next = NULL;
  20. p->d = inf;
  21. p->known = 0;
  22. p->path = -1;
  23. pg->v[i] = p;
  24. }
  25. l* p = (l*)malloc(sizeof(l));
  26. p->val = 1;
  27. p->next = pg->v[0]->next;
  28. pg->v[0]->next = p;
  29. p = (l*)malloc(sizeof(l));
  30. p->val = 3;
  31. p->next = pg->v[0]->next;
  32. pg->v[0]->next = p;
  33. p = (l*)malloc(sizeof(l));
  34. p->val = 3;
  35. p->next = pg->v[1]->next;
  36. pg->v[1]->next = p;
  37. p = (l*)malloc(sizeof(l));
  38. p->val = 4;
  39. p->next = pg->v[1]->next;
  40. pg->v[1]->next = p;
  41. p = (l*)malloc(sizeof(l));
  42. p->val = 0;
  43. p->next = pg->v[2]->next;
  44. pg->v[2]->next = p;
  45. p = (l*)malloc(sizeof(l));
  46. p->val = 5;
  47. p->next = pg->v[2]->next;
  48. pg->v[2]->next = p;
  49. p = (l*)malloc(sizeof(l));
  50. p->val = 2;
  51. p->next = pg->v[3]->next;
  52. pg->v[3]->next = p;
  53. p = (l*)malloc(sizeof(l));
  54. p->val = 4;
  55. p->next = pg->v[3]->next;
  56. pg->v[3]->next = p;
  57. p = (l*)malloc(sizeof(l));
  58. p->val = 5;
  59. p->next = pg->v[3]->next;
  60. pg->v[3]->next = p;
  61. p = (l*)malloc(sizeof(l));
  62. p->val = 6;
  63. p->next = pg->v[3]->next;
  64. pg->v[3]->next = p;
  65. p = (l*)malloc(sizeof(l));
  66. p->val = 6;
  67. p->next = pg->v[4]->next;
  68. pg->v[4]->next = p;
  69. p = (l*)malloc(sizeof(l));
  70. p->val = 5;
  71. p->next = pg->v[6]->next;
  72. pg->v[6]->next = p;
  73. return pg;
  74. }
  75. //队列ADT
  76. typedef struct queue {
  77. int arr[7];
  78. int rear;
  79. int front;
  80. int num;
  81. }q;
  82. q* CreatQueue() {
  83. q* p = (q*)malloc(sizeof(q));
  84. p->front = 0;
  85. p->rear = -1;
  86. p->num = 0;
  87. return p;
  88. }
  89. void enqueue(q* p, int i) {
  90. p->arr[++p->rear] = i;
  91. p->num++;
  92. }
  93. int dequeue(q* p) {
  94. p->num--;
  95. return p->arr[p->front++];
  96. }
  97. int isempty(q* p) {
  98. return p->num == 0;
  99. }
  100. //打印图中的信息
  101. void print(g* p) {
  102. for (int i = 0; i < 7; i++) {
  103. printf("v%d path = v%d know = %d distance = %d\n", i + 1, p->v[i]->path + 1, p->v[i]->known, p->v[i]->d);
  104. }
  105. }
  106. //测试代码
  107. int main() {
  108. g* pg = CreatGraph();
  109. pg->v[2]->d = 0;//选择编号3的顶点为输入顶点
  110. //Unweighted(pg);
  111. UnweightedQ(pg);
  112. print(pg);
  113. return 0;
  114. }
声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/笔触狂放9/article/detail/152137
推荐阅读
相关标签
  

闽ICP备14008679号