hdu 1879 最小生成树 prim算法实现

http://acm.hdu.edu.cn/showproblem.php?pid=1879

这题好诡异，用G++直接超时，c++也就400ms过。。。。

说一下最小生成树的prim算法

按照我的理解，它的主要思想就是说：任意取一个点（通常是第一个），用一个dis[i]来存储终点为 i 的到生成树中任意一点的最小的距离， 在已经加入生成树的节点外找一个点，这个点到生成树的距离要是最小，（不论是到生成树的哪一个节点）。这样直到找到最小的权值，就是一个贪心思想，每次都是最小的和肯定是最小的。。

hdu 1879　prim

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define M 1000
#define INF 100000000
int n;
int map[M][M];
int dis[M];
int used[M];
int sum;
void mst()
{
    for(int i = 1;i <= n;i++)
    {
        dis[i] = map[1][i];  //最短的距离初始化成其他点到1号点的距离
    }
    used[1] = 1;//标记1号点加入生成树
    for(int i = 2;i <= n;i++)  //需要n-1条边
    {
        int min = INF;
        int j;
        for(int k = 1;k <= n;k++)
        {
            if(!used[k]&&dis[k]<min) //找出边的权值最小且不在生成树中的点
            {
                min = dis[k];
                j = k;
            }
        }
        used[j] = 1;  //标记点加入生成树
        //还可以在此处加个判断 if(min==INF)说明根本就没有最小生成树，之间肯定有不连通的地方 hdu1863中可以用到。
        sum += min;
        for(int k = 1;k <= n;k++)　　
        {
            if(!used[k]&&map[j][k]<dis[k])  //更新dis数组（也就是最小的权值）
                dis[k] = map[j][k];
        }
    }
}
int main()
{
    while(scanf("%d",&n)==1 && n)
    {
        sum = 0;
        memset(used,0,sizeof(used));
        //for(int i = 0;i < n;i++)
        //for(int j = 0;j < n;j++)    最小生成树prim算法 应该是要把两点之间的距离都变为正无穷，如果两点之间没有连通的话就是INF
        //map[i][j] = INF;             但是这个给了n*(n-1)/2条边 所以每两个点之间都会被附上权值，不初始化也行
        for(int i = 1;i <= n*(n-1)/2;i++)
        {
            int a,b,c,d;
            scanf("%d %d %d %d",&a,&b,&c,&d);
            if(d==1)
                map[a][b] = map[b][a] = 0;//如果已经建好的话权值为零
            else
                map[a][b] = map[b][a] = c;
        }
        mst();
        printf("%d\n",sum);
    }
    return 0;
}

这里在附上一个别人写的prim算法，不针对这题，但是觉得注释蛮详细的也好理解。

#include <stdio.h>
#include <stdlib.h>
#define MAX 100
#define MAXCOST 0x7fffffff

int graph[MAX][MAX];

int Prim(int graph[][MAX], int n)
{
	/* lowcost[i]记录以i为终点的边的最小权值，当lowcost[i]=0时表示终点i加入生成树 */
	int lowcost[MAX];

	/* mst[i]记录对应lowcost[i]的起点，当mst[i]=0时表示起点i加入生成树 */
	int mst[MAX];

	int i, j, min, minid, sum = 0;

	/* 默认选择1号节点加入生成树，从2号节点开始初始化 */
	for (i = 2; i <= n; i++)
	{
		/* 最短距离初始化为其他节点到1号节点的距离 */
		lowcost[i] = graph[1][i];

		/* 标记所有节点的起点皆为默认的1号节点 */
		mst[i] = 1;
	}

	/* 标记1号节点加入生成树 */
	mst[1] = 0;

	/* n个节点至少需要n-1条边构成最小生成树 */
	for (i = 2; i <= n; i++)
	{
		min = MAXCOST;
		minid = 0;

		/* 找满足条件的最小权值边的节点minid */
		for (j = 2; j <= n; j++)
		{
			/* 边权值较小且不在生成树中 */
			if (lowcost[j] < min && lowcost[j] != 0)
			{
				min = lowcost[j];
				minid = j;
			}
		}
		/* 输出生成树边的信息:起点，终点，权值 */
		printf("%c - %c : %d\n", mst[minid] + 'A' - 1, minid + 'A' - 1, min);

		/* 累加权值 */
		sum += min;

		/* 标记节点minid加入生成树 */
		lowcost[minid] = 0;

		/* 更新当前节点minid到其他节点的权值 */
		for (j = 2; j <= n; j++)
		{
			/* 发现更小的权值 */
			if (graph[minid][j] < lowcost[j])
			{
				/* 更新权值信息 */
				lowcost[j] = graph[minid][j];

				/* 更新最小权值边的起点 */
				mst[j] = minid;
			}
		}
	}
	/* 返回最小权值和 */
	return sum;
}

int main()
{
	int i, j, k, m, n;
	int x, y, cost;
	char chx, chy;

	/* 读取节点和边的数目 */
	scanf("%d%d", &m, &n);
	getchar();

	/* 初始化图，所有节点间距离为无穷大 */
	for (i = 1; i <= m; i++)
	{
		for (j = 1; j <= m; j++)
		{
			graph[i][j] = MAXCOST;
		}
	}

	/* 读取边信息 */
	for (k = 0; k < n; k++)
	{
		scanf("%c %c %d", &chx, &chy, &cost);
		getchar();
		i = chx - 'A' + 1;
		j = chy - 'A' + 1;
		graph[i][j] = cost;
		graph[j][i] = cost;
	}

	/* 求解最小生成树 */
	cost = Prim(graph, m);

	/* 输出最小权值和 */
	printf("Total:%d\n", cost);
}