普利姆-求最小生成树

一、普利姆算法的基本思想：

 从连通网络 N = { V, E }中的某一顶点 u0 出发, 选择与它关联的具有最小权值的边 ( u0, v ), 将其顶点加入到生成树顶点集合U中。

 以后每一步存一个顶点在 U 中,而另一个顶点不在 U 中的各条边中选择权值最小的边(u, v), 把它的顶点加入到集合 U 中。如此继续下去, 直到网络中的所有顶点都加入到生成树顶点集合 U 中为止。

采用邻接矩阵作为图的存储表示。

二、普利姆算法模板：

首先，普利姆算法采用的存储结构是邻接矩阵；所以要定义一个二维数组 map如： map [100 ][100]；

从某一顶点 u 出发找到与其相连边的最小权值时；需要将与其相连的边的权值存到一个数组里；所以定义数组 s 如:s[100];

然后从s数组中寻找与该顶点相连边的最小权值；找到后，标记该顶点u，表示u已加入集合中，然后从新的顶点继续开始，直到所有点都放入集合中（被标记后就代表放入了集合）

需要注意的是map数组刚开始要初始化 min(如下面代码中min=0x3f3f3f3f) 的值；这样在 prim函数中与 max 进行比较的时候就不会出错啦

如果有 n个结点；由于只用到 n*n的矩阵就行了，所以，初始化的时候只对前 n*n 初始化就行啦；

如hdu-1875-畅通工程再续（就可以作为普利姆的模板）

#include <stdio.h>
#include <string.h>
#include <math.h>
#define MAX 1000
#define N 110
#define INF 0x3f3f3f3f
int n;
int x[N],y[N];//保存每个岛的坐标 
double map[N][N];//保存岛与岛之间的距离 

double dis(int i , int j)
{
	return sqrt( 1.0*((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j])) );
}

void input()
{
    int i,j;
    double len;
    scanf("%d",&n);
    for(i=1;i<=n;i++)       //map数组初始化； 
     for(j=1;j<=n;j++)
       map[i][j]=INF; 
    for(i=1; i<=n; i++)
    scanf("%d%d",&x[i],&y[i]);//小岛i的坐标 
    for(i=1; i<=n; i++)
    for(j=1 ;j<=n; j++)
    {
      len=dis(i,j);
      if(10<=len&&len<=1000)
		map[i][j]=map[j][i]=len;//求出岛屿i j之间的距离 
	}
	return ;
}

void prim()
{
    double dis[N],min,sum;
    int vis[N];
    int v,i,j,k;
	memset(vis, 0, sizeof(vis)); //初始化岛屿全为未在点集合内 
    for(i=1; i<=n; i++)
	dis[i]=map[1][i];//记录i点到集合1的距离 
    dis[1]=0;//本身到本身的距离为0
	vis[1] = 1; //标记1已在集合内 
    for(v=1; v<n; v++)  //还要纳入n-1个点
    {
        min=INF; 
		k=1;
        for(i=1; i<=n; i++)
        if(!vis[i] && dis[i] < min)//如果i点未在集合内  并且小于最小值 
        {
			min = dis[i];//记录最小权值 
			k = i;//记录编号 
		}
		if(min == INF)
		{
			printf("oh!\n");
			return ;
	     }
        vis[k] = 1;//标记为已在集合内 
        for(i=1; i<=n; i++)//更新其余点到集合1的距离 
        if(!vis[i] && dis[i] > map[k][i])//如果i未在集合1内并且 当前保存的i到集合的距离大于新更新的值 
        dis[i]=map[k][i];//更新最小值 
    }
    sum = 0;
	for(i = 2; i <= n; ++i)
	sum += dis[i] * 100.0;
	printf("%.1lf\n", sum);
    return ;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        input();
        prim();
    }
    return 0;
}