SDUTOJ 3362 - 数据结构实验之图论六：村村通公路

Problem Description

当前农村公路建设正如火如荼的展开，某乡镇政府决定实现村村通公路，工程师现有各个村落之间的原始道路统计数据表，表中列出了各村之间可以建设公路的若干条道路的成本，你的任务是根据给出的数据表，求使得每个村都有公路连通所需要的最低成本。

Input

连续多组数据输入，每组数据包括村落数目N(N <= 1000)和可供选择的道路数目M(M <= 3000)，随后M行对应M条道路，每行给出3个正整数，分别是该条道路直接连通的两个村庄的编号和修建该道路的预算成本，村庄从1～N编号。

Output

输出使每个村庄都有公路连通所需要的最低成本，如果输入数据不能使所有村庄畅通，则输出-1，表示有些村庄之间没有路连通。

Sample Input

5 8
1 2 12
1 3 9
1 4 11
1 5 3
2 3 6
2 4 9
3 4 4
4 5 6

Sample Output

19

Hint

Source

xam


思路：
prim算法

注：
设顶点集A为最小生成树的点集；顶点集B为待处理点集
lowcost[i]用来保存，顶点集A到i顶点的最小代价。
lowcost数组初始化为INF（无穷大）
随着结点的不断并入，lowcost数组的值不断更新
为0代表该结点已并入A顶点集
为INF代表该结点与A顶点集无直接边相连

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

#define INF 0x3f3f3f3f // 必须定义一个接近无穷大的数

int vexnum, edgenum; // 顶点数和边数
int map[1005][1005]; // 邻接矩阵
int lowcost[1005];
int result = 0; // 累加器

int prim() // 普里姆算法
{
    int i, j, k, min;
    lowcost[1] = 0; // 顶点集A中并入1号顶点，即从1号顶点开始
    for (i = 2; i <= vexnum; i++) // 更新lowcost数组
    {
        lowcost[i] = map[i][1];
    }

    for (i = 2; i <= vexnum; i++) // 通过n-1次循环，并入n-1个点
    {
        min = INF;

        j = 1;
        k = 0;
        while(j <= vexnum) // 找到顶点集A的最小代价边
        {
            if (lowcost[j] != 0&& lowcost[j] < min) // 保证不成环且最小
            {
                min = lowcost[j];
                k = j;
            }
            j++;
        }
        if (min == INF) // 进行判断，若顶点集A不再有边，退出循环
        {
            break;
        }
        result += min; // 累加权值
        lowcost[k] = 0; // 将K号顶点并入顶点集A
        for (j = 1; j <= vexnum; j++) // 通过新并入的K号顶点，更新lowcost数组
        {
            if (lowcost[j] != 0&&lowcost[j] > map[j][k])
            {
                lowcost[j] = map[j][k];
            }
        }
    }
    for (i = 1; i <= vexnum; i++) // 判断是否所有边都并入顶点集A
    {
        if (lowcost[i] != 0) // 若有顶点为并入顶点集A，则图不连通
        {
            return -1;
        }
    }
    return result;
}

int main()
{
    int i,u,v,c;
    while(~scanf("%d %d",&vexnum,&edgenum))
    {
        result=0;
        memset(map,INF,sizeof(map));
        memset(lowcost,INF,sizeof(lowcost));
        for(i=1;i<=edgenum;i++)
        {
            scanf("%d %d %d",&u,&v,&c);
            map[u][v]=c;
            map[v][u]=c;
        }
        int value = prim();
        printf("%d\n", value);
    }
    return 0;
}