最小生成树---普里姆算法

【普里姆算法】
本算法以无向图为例，图的存储方式采用邻接矩阵
1)以第一个顶点开始本算法，将第一个顶点(i = 0)加入最小生成树集low_cost[i] = 0，0代表已经加入生成树集;
2)遍历第i行邻接矩阵，更新与顶点i有边的顶点的权值(该权值集的初始化值为与第一个顶点有边连接的顶点的权值)
3)遍历这个low_cost权值集合，找到最小权值并输出

为了找到与当前顶点有边的前一个顶点，算法用Vertex_Information[]进行了记录

我们先为左图构造邻接矩阵，如右图所示

算法从A点开始生成最小生成树，代码如下：

/*********************************************************************Author:tmwdate:2017-10-17********************************************************************************************/#include <stdio.h>#include <stdlib.h>#define MAX_VERTEX 100#define inf 65535  //用65535代表无穷大/*************************建立无向图邻接矩阵**************************/typedef struct Matrix_Graph{    char vertex[MAX_VERTEX];//存储顶点信息    int edge[MAX_VERTEX][MAX_VERTEX];//存储边信息    int vertex_number,edge_number;//存储图中的顶点数和边数}Matrix_Graph;void create_non_direction_matrix_Graph( Matrix_Graph *G ){    int i,j,w,k;    printf("请输入顶点数和边数：\n");    scanf("%d %d",&G->vertex_number,&G->edge_number);    printf("请输入无向图顶点信息(如ABCD...)：\n");    char ch;    while( ( ch = getchar() ) != '\n' );    for( i = 0 ; i < G->vertex_number ; i++ )        scanf("%c",&G->vertex[i]);    //初始化边信息    for( i = 0 ; i < G->vertex_number ; i++ )        for( j = 0 ; j < G->vertex_number ; j++ )        {            if( i == j )                G->edge[i][j] = 0;            else                G->edge[i][j] = inf;        }    printf("请输入无向图中相连的两个顶点下标值(Vi,Vj)，以及权值---即边信息\n");    for( k = 0 ; k < G->edge_number ; k++ )    {        scanf("%d %d %d",&i,&j,&w);        G->edge[i][j] = w;        G->edge[j][i] = G->edge[i][j];//由于是无向图，故是对称阵    }    //打印邻接矩阵    for( i = 0 ; i < G->vertex_number ; i++ )    {        for( j = 0 ; j < G->vertex_number ; j++ )            printf("%d\t",G->edge[i][j]);        printf("\n");    }}void Min_Create_Tree_Prim( Matrix_Graph *G ){    int i,j,min;    int mark = 0;    int low_cost[MAX_VERTEX];    int Vertex_Information[MAX_VERTEX];//便于找到权值最小边的前一个顶点    low_cost[0] = 0;     //low_cost记录已经用于生成最小生成树的顶点，若为0，则标记为用过                        //这里从第一个顶点开始生成最小生成树，因此首先标记它为0    for( i = 1 ; i < G->vertex_number ; i++ )//因为第一个顶点已经纳入生成树，因此从第二个顶点开始遍历    {        low_cost[i] = G->edge[0][i];//将与第一个顶点有边相连的顶点权值放入low_cost数组        Vertex_Information[i] = 0; //初始化相关顶点的下标都从第一个顶点开始    }    for( i = 1 ; i < G->vertex_number ; i++ )    {        min = inf;        j = 1;        while( j < G->vertex_number )//找出与当前顶点权值最小的顶点        {            if( low_cost[j]!=0 && low_cost[j] < min )            {                min = low_cost[j];                mark = j;//将最小权值顶点做个标记            }            j++;        }        low_cost[mark] = 0;//标记新加入顶点已经被加入生成树        printf("[%d %d] ---",Vertex_Information[mark],mark);//打印顶点下标号        printf("%c %c\n",G->vertex[Vertex_Information[mark]],G->vertex[mark]);//打印找到的当前顶点的最小权值边        //找 与当前顶点和新加入顶点中 还相连的其他顶点        for( j = 1 ; j < G->vertex_number ; j++ )        {            if( low_cost[j] != 0 && low_cost[j] > G->edge[mark][j] )//更新 与当前顶点和新加入顶点间 其他与之相连顶点的权值            {                low_cost[j] = G->edge[mark][j];                Vertex_Information[j] = mark;//标记边的第一个顶点            }        }    }}int main(){    Matrix_Graph *G;    G = (Matrix_Graph*)malloc(sizeof(Matrix_Graph));    create_non_direction_matrix_Graph(G);    Min_Create_Tree_Prim(G);    return 0;}

【程序执行结果】

从代码结果可看出，本次构造的最小生成树路径为：[A B]、[A F]、[B I]、[I C]、[B G]、[G H]、[H E]、[H D]