带权无向图求最短路径_迪杰斯特拉算法求最短路径(C语言详解版)

如今出行已经不需要再为找不着路而担心了,车上有车载导航,手机中有导航 App。只需要确定起点和终点,导航会自动规划出可行的距离最短的道路。这是最短路径在人们实际生活中最典型的应用。在一个网(有权图)中,求一个顶点到另一个顶点的最短路径的计算方式有两种:迪杰斯特拉(Dijkstra算法)和弗洛伊德(Floyd)算法。迪杰斯特拉算法计算的是有向网中的某个顶点到其余所有顶点的最短路径;弗洛伊德算法计算的是任意两顶点之间的最短路径。

最短路径算法既适用于有向网,也同样适用于无向网。本节将围绕有向网讲解迪杰斯特拉算法的具体实现。

迪杰斯特拉算法计算的是从网中一个顶点到其它顶点之间的最短路径问题。

9f56822d4d002fe42446a7cab01867f6.png图 1 带权有向图

如图 1 所示是一个有向网,在计算 V0 到其它所有顶点之间的最小路径时,迪杰斯特拉算法的计算方式为:从 V0 出发,由于可以直接到达 V2 和 V5 ,而其它顶点和 V0 之间没有弧的存在,所以之间的距离设定为无穷大,可以得到下面这个表格:

faffcb693624c99be06d4f745aaeba71.png 

从表格中可以看到,V0 到 V2 的距离最近,所以迪杰斯特拉算法设定 V0-V2 为 V0 到 V2 之间的最短路径,最短路径的权值和为 10。已经判断 V0-V2 是最短路径,所以以 V2 为起始点,判断 V2 到除了 V0 以外的其余各点之间的距离,如果对应的权值比前一张表格中记录的数值小,就说明网中有一条更短的路径,直接更新表格;反之表格中的数据不变。可以得到下面这个表格:

0cf01d10bc30d33055988ea843d0c711.png 

例如,表格中 V0 到 V3 的距离,发现当通过 V2 到达 V3 的距离比之前的 ∞ 要小,所以更新表格。更新之后,发现 V0-V4 的距离最近,设定 V0 到 V4 的最短路径的值为 30。之后从 V4 出发,判断到未确定最短路径的其它顶点的距离,继续更新表格:

561d6c4a67775181057c3651cefd0b30.png 

更新后确定从 V0 到 V3 的最短路径为 V0-V4-V3,权值为 50。然后从 V3 出发,继续判断:

de4fea46a97d0b679699c8d3c050bf32.png 

对于 V5 来说,通过 V0-V4-V3-V5 的路径要比之前的权值 90 还要小,所以更新表格,更新后可以看到,V0-V5 的距离此时最短,可以确定 V0 到 V5 的最短路径为 60。最后确定 V0-V1 的最短路径,由于从 V0 无法到达 V1 ,最终设定 V0 到 V1 的最短路径为 ∞(无穷大)。在确定了 V0 与其他所有顶点的最短路径后,迪杰斯特拉算法才算结束。

事例中借用了图 1 的有向网对迪杰斯特拉算法进行了讲解,实际上无向网中的最短路径问题也可以使用迪杰斯特拉算法解决,解决过程和上述过程完全一致。

迪杰斯特拉算法的代码实现

#include #define MAX_VERtEX_NUM 20                   //顶点的最大个数#define VRType int                          //表示弧的权值的类型#define VertexType int                      //图中顶点的数据类型#define INFINITY 65535typedef struct {    VertexType vexs[MAX_VERtEX_NUM];        //存储图中顶点数据    VRType arcs[MAX_VERtEX_NUM][MAX_VERtEX_NUM];                         //二维数组,记录顶点之间的关系    int vexnum,arcnum;                      //记录图的顶点数和弧(边)数}MGraph;typedef int PathMatrix[MAX_VERtEX_NUM];     //用于存储最短路径中经过的顶点的下标typedef int ShortPathTable[MAX_VERtEX_NUM]; //用于存储各个最短路径的权值和//根据顶点本身数据,判断出顶点在二维数组中的位置int LocateVex(MGraph * G,VertexType v){    int i=0;    //遍历一维数组,找到变量v    for (; ivexnum; i++) {        if (G->vexs[i]==v) {            break;        }    }    //如果找不到,输出提示语句,返回-1    if (i>G->vexnum) {        printf("no such vertex.\n");        return -1;    }    return i;}//构造有向网void CreateUDG(MGraph *G){    scanf("%d,%d",&(G->vexnum),&(G->arcnum));    for (int i=0; ivexnum; i++) {        scanf("%d",&(G->vexs[i]));    }    for (int i=0; ivexnum; i++) {        for (int j=0; jvexnum; j++) {            G->arcs[i][j]=INFINITY;        }    }    for (int i=0; iarcnum; i++) {        int v1,v2,w;        scanf("%d,%d,%d",&v1,&v2,&w);        int n=LocateVex(G, v1);        int m=LocateVex(G, v2);        if (m==-1 ||n==-1) {            printf("no this vertex\n");            return;        }        G->arcs[n][m]=w;    }}//迪杰斯特拉算法,v0表示有向网中起始点所在数组中的下标void ShortestPath_Dijkstra(MGraph G,int v0,PathMatrix *p,ShortPathTable *D){    int final[MAX_VERtEX_NUM];//用于存储各顶点是否已经确定最短路径的数组    //对各数组进行初始化    for (int v=0; v        final[v]=0;        (*D)[v]=G.arcs[v0][v];        (*p)[v]=0;    }    //由于以v0位下标的顶点为起始点,所以不用再判断    (*D)[v0]=0;    final[v0]=1;    int k = 0;    for (int i=0; i        int min=INFINITY;        //选择到各顶点权值最小的顶点,即为本次能确定最短路径的顶点        for (int w=0; w            if (!final[w]) {                if ((*D)[w]                    k=w;                    min=(*D)[w];                }            }        }        //设置该顶点的标志位为1,避免下次重复判断        final[k]=1;        //对v0到各顶点的权值进行更新        for (int w=0; w            if (!final[w]&&(min+G.arcs[k][w]                (*D)[w]=min+G.arcs[k][w];                (*p)[w]=k;//记录各个最短路径上存在的顶点            }        }    }}int main(){    MGraph G;    CreateUDG(&G);    PathMatrix P;    ShortPathTable D;    ShortestPath_Dijkstra(G, 0, &P, &D);    for (int i=1; i        printf("V%d - V%d的最短路径中的顶点有(逆序):",0,i);        printf(" V%d",i);        int j=i;        //由于每一段最短路径上都记录着经过的顶点,所以采用嵌套的方式输出即可得到各个最短路径上的所有顶点        while (P[j]!=0) {            printf(" V%d",P[j]);            j=P[j];        }        printf(" V0\n");    }    printf("源点到各顶点的最短路径长度为:\n");    for (int i=1; i        printf("V%d - V%d : %d \n",G.vexs[0],G.vexs[i],D[i]);    }    return 0;}

运行以上代码,计算如图 1 所示的有向网,运行结果为:

6,8
0
1
2
3
4
5
0,5,100
0,4,30
0,2,10
1,2,5
2,3,50
3,5,10
4,3,20
4,5,60
V0 - V1的最短路径中的顶点有:V1 V0
V0 - V2的最短路径中的顶点有:V2 V0
V0 - V3的最短路径中的顶点有:V3 V4 V0
V0 - V4的最短路径中的顶点有:V4 V0
V0 - V5的最短路径中的顶点有:V5 V3 V4 V0
源点到各顶点的最短路径长度为:
V0 - V1 : 65535
V0 - V2 : 10
V0 - V3 : 50
V0 - V4 : 30
V0 - V5 : 60

总之,迪杰斯特拉算法解决的是从网中的一个顶点到所有其它顶点之间的最短路径,算法整体的时间复杂度为 O(n2)。但是如果需要求任意两顶点之间的最短路径,使用迪杰斯特拉算法虽然最终虽然也能解决问题,但是大材小用,相比之下使用弗洛伊德算法解决此类问题会更合适。

弗洛伊德算法实现的具体过程下一节会作详细介绍。

  • 5
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,以下是使用C语言编写的程序实现: ```c #include <stdio.h> #include <stdlib.h> #include <string.h> #define INF 0x3f3f3f3f const int MAXN = 1005; int n, m; // n表示顶点个数,m表示边数 int s, t; // s表示源点,t表示终点 int graph[MAXN][MAXN]; // 邻接矩阵存 int dist[MAXN]; // 存储源点到各个点的最距离 int vis[MAXN]; // 标记是否已访问 // 斯特算法 void dijkstra(int s) { memset(dist, INF, sizeof(dist)); memset(vis, 0, sizeof(vis)); dist[s] = 0; for (int i = 1; i <= n; i++) { int k = -1, minn = INF; for (int j = 1; j <= n; j++) { if (!vis[j] && dist[j] < minn) { k = j; minn = dist[j]; } } if (k == -1) break; vis[k] = 1; for (int j = 1; j <= n; j++) { if (!vis[j] && dist[k] + graph[k][j] < dist[j]) { dist[j] = dist[k] + graph[k][j]; } } } } int main() { scanf("%d %d %d %d", &n, &m, &s, &t); memset(graph, INF, sizeof(graph)); for (int i = 1; i <= m; i++) { int u, v, w; scanf("%d %d %d", &u, &v, &w); graph[u][v] = graph[v][u] = w; } dijkstra(s); printf("%d\n", dist[t]); return 0; } ``` 程序中使用了邻接矩阵来存储的信息。其中,`graph[i][j]`表示点i到点j的边的权值(如果不存在则为INF)。 在`dijkstra`函数中,首先将`dist`数组初始化为INF,并将源点s的距离设为0。然后对于每个未访问的点,找到距离源点最近的点k,并将其标记为已访问。然后遍历以k为起点的所有边,如果用k到v的距离加上源点到k的距离比源点到v的距离更小,则更新源点到v的距离。 最终,`dist[t]`即为源点s到终点t的最路径长度。 需要注意的是,本程序中使用了一个`INF`常量表示正无穷,其值可以根据具体情况设置为一个足够大的数。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值