• 单源无权最短路径


    单源无权最短路径是指对一个无权图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. }
  • 相关阅读:
    你知道this.$options吗?(Vue)
    数据结构和算法——线性结构
    VTK——angleWidget的3D转换
    【Java List与数组】List<T>数组和数组List<T>的区别(124)
    【JY】求?减隔震元件的滞回面积~
    低代码平台简介(10家国产化低代码平台详细介绍)
    【面试题】如何理解 前端设计模式-测策略模式?
    vue(十二)——vue3新特性之Teleport
    Python实现IP代理批量采集, 并检测代理是否可用 保存一下<嫖一波免费代理>
    国内什么牌子的ipad手写笔好用?开学性价比触控笔
  • 原文地址:https://blog.csdn.net/thdwx/article/details/127424732