普利姆算法

普利姆算法

一:介绍:

(一):原理介绍

Prim算法构造最小生成树过程如下图所示:

​ 1、首先从图中任选一个顶点加入树T中,此时最小生成树T中就只含有一个顶点

​ 2、然后选择与当前最小生成树T中顶点集合距离最近的顶点,并将该顶点和相应的边加入最小生成树T中,每次操作树T中的顶点数量和边都加一。

​ 3、执行1,2两步,当图中的所有顶点都并入到T中,T就是最小生成树了。此时T中必有n-1条边
普利姆算法

(二):算法思路介绍

​ 1、先初始化辅助数组,即从起始点到其他顶点的权值保存到数组中,并标志起始点

​ 2、循环执行以下操作(循环次数为图中顶点的数量减一,因为起始点不算)

​ (1)、定位数组中权值最小的权值,从而找到对应的下一个要加入的顶点。

​ (2)、在辅助数组中将新加入的顶点的权值置为0(自己到自己的距离为0)

​ (3)、跟新辅助数组,在新加入定点之后,集合T中的顶点到其他顶点的最小距离会变化,所以跟新新加入的顶点到其他顶点的距离,如果比原来的小,则跟新,否则不跟新。

二:代码

(一):结构

#ifndef DATASTRUCTURE_MST_H
#define DATASTRUCTURE_MST_H

#include <stdio.h>
#include "../Graph/Graph.h"
//#include "../Graph/Graph.c"

//生成树的类型定义
typedef struct {			
    VertexType adjvex;      //较早加入当前集合T的顶点(权值最小边的起始点)
    VRType lowcost;         //以该顶点为起始点的边的权值
}Edge;                      //辅助数组

typedef Edge closeEdge[MAX_VERTEX_NUM];     //创建一个辅助数组,表示已经加入到T中集合的顶点到其他顶点的最近距离

//在辅助数组中找到权值最小的边的数组下标
int Minimum(MGraph G,closeEdge close);

//普里母算法
void MinSpanTree_Prim(MGraph G,VertexType v);

#endif //DATASTRUCTURE_MST_H

(二):算法

#include "MST.h"

/**
 * 普利姆算法
 * @param G     无向网
 * @param v     开始顶点(随意)	初始为顶点1
 */
void MinSpanTree_Prim(MGraph G,VertexType v){
    int k = LocateVex(G,v);
    closeEdge close;		//辅助数组
    //初始化辅助数组,将与该起始点有关的所有边的信息:边的起始点和权值,初始化到辅助数组中
    //如初始顶点为1,那么有顶点1相连的有顶点2,3,4,则顶点1到2,3,4为其边上的权值,而与顶点5,6不相邻,所以为INFINITY
    //其中辅助数组的大小为图中顶点的数量,因此对应起始点与图中顶点的权值,即close[i].adjvex表示的是起始顶点,i对应的是顶点在
    //G中顶点数组的下标,所以可以通过i来找到顶点值
    for (int i = 0; i < G.vexnum; ++i) {
        if (i!=k){			
            close[i].adjvex = k;
            close[i].lowcost = G.arcs[k][i].adj;
        }
    }

    //因为起始点已经在最小生成树中了,所以需要将其辅助数组置0,或者因为从该点开始的,所以距离为0
    close[k].lowcost = 0;

    //选择下一个点,并更新辅助数组中的信息
    for (int i = 1; i < G.vexnum; ++i) {        //为什么从1开始,因为目前已经有一个在最小生成树中了
        k = Minimum(G,close);                   //找到权值最小的边所在数组的下标
        printf("v%d \t v%d\t",G.vexs[close[k].adjvex],G.vexs[k]);           //输出选择的路径
        close[k].lowcost = 0;                   //入树之后需要将权值设置为0

        //因为加入了一个顶点,那么辅助数组中的权值可能就会发生改变,所以,需要从新加入的顶点开始,比较其到其他邻接顶点的权值
        //如果小的就需要更新
        for (int j = 0; j < G.vexnum; ++j) {
            if (G.arcs[k][j].adj<close[j].lowcost){
                close[j].lowcost = G.arcs[k][j].adj;
                close[j].adjvex = k;

            }
        }
        printf("\n");
    }
}


//在辅助数组中找到权值最小的边的数组下标
int Minimum(MGraph G,closeEdge close){
    int min = INFINITY;			//为无向网,所以INFINTIY表示最大距离
    int index;
    for (int i = 0; i < G.vexnum; ++i) {
        if (close[i].lowcost>0&&close[i].lowcost<min){  //当权值为0的时候,说明顶点已经归到了最下生成树中了
            min = close[i].lowcost;
            index = i;
        }
    }
    return index;
}

/**
 * 确定顶点v在图中的位置
 * @param G
 * @param v		为顶点值
 * @return		如果在图的顶点数组中找到,则返回在数组中的下标,否则返回-1;
 */
int LocateVex(MGraph G,VertexType v){		
    for (int i = 0; i < G.vexnum; ++i) {
        if (v == G.vexs[i]){
            return i;
        }
    }
    return -1;
}

上一篇:HttpWebResponse 获取网页乱码


下一篇:pymysql 链接mysql