如果求一个连通图中以某个顶点为根的深度最小的生成树_图的应用(2)-最小路径问题...

上篇文章是求图的最小生成树，是为了把所有的点以最小的成本连接起来，组成一个网络。

GhostClock：图的应用-连通图的最小生成树(Prim算法和Kruskal算法)​zhuanlan.zhihu.com

这篇文章是求图的最短路径。

求图的最短路径，其实就是求起点到终点，所经过的路径值最小。

把上图存入邻接矩阵中：

用邻接矩阵保存图

相应的数据结构

#define OK 1
#define ERROR 0
#define TRUE 1
#define FALSE 0

#define MAXEDGE 20
#define MAXVEX 20
#define INFINITYC 65535

typedef int Status;
typedef struct
{
    int vexs[MAXVEX];
    int arc[MAXVEX][MAXVEX];
    int numVertexes, numEdges;
}MGraph;

/*用于存储最短路径下标的数组*/
typedef int Patharc[MAXVEX];
/*用于存储到各点最短路径权值的和*/
typedef int ShortPathTable[MAXVEX];

创建邻近矩阵

/*创建邻近矩阵*/
void CreateMGraph(MGraph *G)
{
    int i, j;
    
    G->numEdges=16;
    G->numVertexes=9;
    
    for (i = 0; i < G->numVertexes; i++)
    {
        G->vexs[i]=i;
    }
    
    for (i = 0; i < G->numVertexes; i++)
    {
        for ( j = 0; j < G->numVertexes; j++)
        {
            if (i==j)
                G->arc[i][j]=0;
            else
                G->arc[i][j] = G->arc[j][i] = INFINITYC;
        }
    }
    
    G->arc[0][1]=1;
    G->arc[0][2]=5;
    G->arc[1][2]=3;
    G->arc[1][3]=7;
    G->arc[1][4]=5;
    
    G->arc[2][4]=1;
    G->arc[2][5]=7;
    G->arc[3][4]=2;
    G->arc[3][6]=3;
    G->arc[4][5]=3;
    
    G->arc[4][6]=6;
    G->arc[4][7]=9;
    G->arc[5][7]=5;
    G->arc[6][7]=2;
    G->arc[6][8]=7;
    
    G->arc[7][8]=4;
    
    
    for(i = 0; i < G->numVertexes; i++)
    {
        for(j = i; j < G->numVertexes; j++)
        {
            G->arc[j][i] =G->arc[i][j];
        }
    }
    
}

Dijkstra(迪杰斯特拉)算法

https://www.wikiwand.com/zh-hans/%E6%88%B4%E5%85%8B%E6%96%AF%E7%89%B9%E6%8B%89%E7%AE%97%E6%B3%95​www.wikiwand.com

该算法有三个数组：final数组，D数组，p数组。

final数组：表示V0到某个顶点Vw是否已经求得最短路径的标记。如果V0到Vw已经有结果，则final[w]=1。

D数组：表示V0到某个顶点Vw的路径。

p数组：当前顶点的前驱顶点的下标。

/*求得网图中2点间最短路径
 Dijkstra 算法
 G: 网图;
 v0: V0开始的顶点;
 p[v]: 前驱顶点下标;
 D[v]: 表示从V0到V的最短路径长度和;
 */
void ShortestPath_Dijkstra(MGraph G, int v0, Patharc *P, ShortPathTable *D)
{
    int v,w,k,min;
    k = 0;
    /*final[w] = 1 表示求得顶点V0~Vw的最短路径*/
    int final[MAXVEX];
    
    /*1.初始化数组*/
    for(v=0; v<G.numVertexes; v++)
    {
        //全部顶点初始化为未知最短路径状态0
        final[v] = 0;
        //将与V0 点有连线的顶点最短路径值;
        (*D)[v] = G.arc[v0][v];
        //初始化路径数组p = 0;
        (*P)[v] = 0;
    }
    
    //V0到V0的路径为0
    (*D)[v0] = 0;
    //V0到V0 是没有路径的.
    final[v0] = 1;
    //v0到V0是没有路径的
    (*P)[v0] = -1;
    
    //2. 开始主循环,每次求得V0到某个顶点的最短路径
    for(v=1; v<G.numVertexes; v++)
    {
        
        //当前所知距离V0顶点最近的距离
        min=INFINITYC;
        /*3.寻找离V0最近的顶点*/
        for(w=0; w<G.numVertexes; w++)
        {
            if(!final[w] && (*D)[w]<min)
            {
                k=w;
                //w顶点距离V0顶点更近
                min = (*D)[w];
            }
        }
        
        //将目前找到最近的顶点置为1;
        final[k] = 1;
        
        /*4.把刚刚找到v0到v1最短路径的基础上,对于v1 与 其他顶点的边进行计算,得到v0与它们的当前最短距离;*/
        for(w=0; w<G.numVertexes; w++)
        {
            //如果经过v顶点的路径比现在这条路径长度短,则更新
            if(!final[w] && (min + G.arc[k][w]<(*D)[w]))
            {
                //找到更短路径, 则修改D[W],P[W]
                //修改当前路径的长度
                (*D)[w] = min + G.arc[k][w];
                (*P)[w]=k;
            }
        }
    }
}

使用和打印：

printf("最短路径-Dijkstra算法n");
int i,j,v0;
MGraph G;
Patharc P;
ShortPathTable D;
v0=0;
CreateMGraph(&G);
ShortestPath_Dijkstra(G, v0, &P, &D);
printf("最短路径路线:n");
for(i=1;i<G.numVertexes;++i)
{
    printf("v%d -> v%d : ",v0,i);
    j=i;
    while(P[j]!=-1)
    {
        printf("%d ",P[j]);
        j=P[j];
    }
    printf("n");
}
printf("n最短路径权值和n");
for(i=1;i<G.numVertexes;++i)
    printf("v%d -> v%d : %d n",G.vexs[0],G.vexs[i],D[i]);
printf("n");

Dijkstra算法打印

Floyd(弗洛伊德)算法

https://www.wikiwand.com/zh-hans/Floyd-Warshall%E7%AE%97%E6%B3%95​www.wikiwand.com

公式:

Floyd算法公式

实现

/*Floyd算法，求网图G中各顶点v到其余顶点w的最短路径P[v][w]及带权长度D[v][w]。
 Patharc 和 ShortPathTable 都是二维数组;
 */
void ShortestPath_Floyd(MGraph G, Patharc *P, ShortPathTable *D)
{
    int v,w,k;
    
    /* 1. 初始化D与P 矩阵*/
    for(v=0; v<G.numVertexes; ++v)
    {
        for(w=0; w<G.numVertexes; ++w)
        {
            /* D[v][w]值即为对应点间的权值 */
            (*D)[v][w]=G.arc[v][w];
             /* 初始化P P[v][w] = w*/
            (*P)[v][w]=w;
        }
    }
    
    //2.K表示经过的中转顶点
    for(k=0; k<G.numVertexes; ++k)
    {
        for(v=0; v<G.numVertexes; ++v)
        {
            for(w=0; w<G.numVertexes; ++w)
            {
                /*如果经过下标为k顶点路径比原两点间路径更短 */
                if ((*D)[v][w]>(*D)[v][k]+(*D)[k][w])
                {
                    /* 将当前两点间权值设为更小的一个 */
                    (*D)[v][w]=(*D)[v][k]+(*D)[k][w];
                    /* 路径设置为经过下标为k的顶点 */
                    (*P)[v][w]=(*P)[v][k];
                }
            }
        }
    }
}

使用和打印：

printf("Hello,最短路径弗洛伊德Floyd算法");
int v,w,k;
MGraph G;
Patharc P;
ShortPathTable D; /* 求某点到其余各点的最短路径 */
CreateMGraph(&G);
ShortestPath_Floyd(G,&P,&D);
//打印所有可能的顶点之间的最短路径以及路线值
printf("各顶点间最短路径如下:n");
for(v=0; v<G.numVertexes; ++v)
{
    for(w=v+1; w<G.numVertexes; w++)
    {
        printf("v%d-v%d weight: %d ",v,w,D[v][w]);
        //获得第一个路径顶点下标
        k=P[v][w];
        //打印源点
        printf(" path: %d",v);
        //如果路径顶点下标不是终点
        while(k!=w)
        {
            //打印路径顶点
            printf(" -> %d",k);
            //获得下一个路径顶点下标
            k=P[k][w];
        }
        //打印终点
        printf(" -> %dn",w);
    }
    printf("n");
}
//打印最终变换后的最短路径D数组
printf("最短路径D数组n");
for(v=0; v<G.numVertexes; ++v)
{
    for(w=0; w<G.numVertexes; ++w)
    {
        printf("%dt",D[v][w]);
    }
    printf("n");
}
//打印最终变换后的最短路径P数组
printf("最短路径P数组n");
for(v=0; v<G.numVertexes; ++v)
{
    for(w=0; w<G.numVertexes; ++w)
    {
        printf("%d ",P[v][w]);
    }
    printf("n");
}

Floyd算法打印