最小生成树
什么是最小生成树,就是一个图中,将所有顶点连结并且其权重总和最小。
当然前提是图是一个连通图
1. Prim算法
1.1核心思路
从一个顶点出发,不断贪婪已知可以连接的最小权值边。
- 一个adjvex数组不断更新相连,且权值更小的顶点。
- adjvex中是起点,adjuer[k] -> k
- 一个lowcost数组不断更新可选的最小权值,通过0标记已加入最小生成树。
- lowcost 中存储 adjuex[k] --> k 的权值,新顶点指向k的权值更小,就会更新。
思路:
1、定义2个数组;adjvex 用来保存相关顶点下标;lowcost 保存顶点之间的权值
2、初始化2个数组,从v0 开始寻找最小生成树,默认v0 是最小生成树上第一个顶点
3、循环lowcost 数组,根据权值,找到顶点k;
4、更新lowcost 数组
5、循环所有顶点,找到与顶点k有关系的顶点,并更新lowcost 数组与 adjvex 数组;
注意更新lowcost 数组与adjvex 数组的条件
1、与顶点k之间有连接
2、当前结点j 没有加入过最小生成树;
3、顶点k 与 当前定点j 之间的权值 小于顶点j与其他顶点k 之间的权值,则更新。简单说就是要比较值钱存储的值药效,则更新
代码实现
#include "stdio.h"
#include "stdlib.h"
#include "math.h"
#include "time.h"
#define OK 1
#define ERROR 0
#define TRUE 1
#define FALSE 0
#define MAXEDGE 20
#define MAXVEX 20
#define INFINITYC 65535
typedef int Status; /* Status是函数的类型,其值是函数结果状态代码,如OK等 */
typedef struct
{
int arc[MAXVEX][MAXVEX];
int numVertexes, numEdges;
}MGraph;
/*9.1 创建邻接矩阵*/
void CreateMGraph(MGraph *G)/* 构件图 */
{
int i, j;
/* printf("请输入边数和顶点数:"); */
G->numEdges=15;
G->numVertexes=9;
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]=10;
G->arc[0][5]=11;
G->arc[1][2]=18;
G->arc[1][8]=12;
G->arc[1][6]=16;
G->arc[2][8]=8;
G->arc[2][3]=22;
G->arc[3][8]=21;
G->arc[3][6]=24;
G->arc[3][7]=16;
G->arc[3][4]=20;
G->arc[4][7]=7;
G->arc[4][5]=26;
G->arc[5][6]=17;
G->arc[6][7]=19;
for(i = 0; i < G->numVertexes; i++)
{
for(j = i; j < G->numVertexes; j++)
{
G->arc[j][i] =G->arc[i][j];
}
}
}
/* Prim算法生成最小生成树 */
void MiniSpanTree_Prim(MGraph G)
{
int min, i, j, k;
int sum = 0;
/* 保存相关顶点下标 */
int adjvex[MAXVEX];
/* 保存相关顶点间边的权值 */
int lowcost[MAXVEX];
/* 初始化第一个权值为0,即v0加入生成树 */
/* lowcost的值为0,在这里就是此下标的顶点已经加入生成树 */
lowcost[0] = 0;
/* 初始化第一个顶点下标为0 */
adjvex[0] = 0;
//1. 初始化
for(i = 1; i < G.numVertexes; i++) /* 循环除下标为0外的全部顶点 */
{
lowcost[i] = G.arc[0][i]; /* 将v0顶点与之有边的权值存入数组 */
adjvex[i] = 0; /* 初始化都为v0的下标 */
}
//2. 循环除了下标为0以外的全部顶点, 找到lowcost数组中最小的顶点k
for(i = 1; i < G.numVertexes; i++)
{
/* 初始化最小权值为∞, */
/* 通常设置为不可能的大数字如32767、65535等 */
min = INFINITYC;
j = 1;k = 0;
while(j < G.numVertexes) /* 循环全部顶点 */
{
/* 如果权值不为0且权值小于min */
if(lowcost[j]!=0 && lowcost[j] < min)
{
/* 则让当前权值成为最小值,更新min */
min = lowcost[j];
/* 将当前最小值的下标存入k */
k = j;
}
j++;
}
/* 打印当前顶点边中权值最小的边 */
printf("(V%d, V%d)=%d\n", adjvex[k], k ,G.arc[adjvex[k]][k]);
sum+=G.arc[adjvex[k]][k];
/* 3.将当前顶点的权值设置为0,表示此顶点已经完成任务 */
lowcost[k] = 0;
/* 循环所有顶点,找到与顶点k 相连接的顶点
1. 与顶点k 之间连接;
2. 该结点没有被加入到生成树;
3. 顶点k 与 顶点j 之间的权值 < 顶点j 与其他顶点的权值,则更新lowcost 数组;
*/
for(j = 1; j < G.numVertexes; j++)
{
/* 如果下标为k顶点各边权值小于此前这些顶点未被加入生成树权值 */
if(lowcost[j]!=0 && G.arc[k][j] < lowcost[j])
{
/* 将较小的权值存入lowcost相应位置 */
lowcost[j] = G.arc[k][j];
/* 将下标为k的顶点存入adjvex */
adjvex[j] = k;
}
}
}
printf("sum = %d\n",sum);
}
int main(void)
{
printf("Hello,最小生成树_Prim算法\n");
MGraph G;
CreateMGraph(&G);
MiniSpanTree_Prim(G);
return 0;
}
Hello,最小生成树_Prim算法
(V0, V1)=10
(V0, V5)=11
(V1, V8)=12
(V8, V2)=8
(V1, V6)=16
(V6, V7)=19
(V7, V4)=7
(V7, V3)=16
sum = 99