HDU1875 畅通工程再续【Kruskal算法+并查集】

畅通工程再续

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 35029    Accepted Submission(s): 11633

 

Problem Description

相信大家都听说一个“百岛湖”的地方吧，百岛湖的居民生活在不同的小岛中，当他们想去其他的小岛时都要通过划小船来实现。现在政府决定大力发展百岛湖，发展首先要解决的问题当然是交通问题，政府决定实现百岛湖的全畅通！经过考察小组RPRush对百岛湖的情况充分了解后，决定在符合条件的小岛间建上桥，所谓符合条件，就是2个小岛之间的距离不能小于10米，也不能大于1000米。当然，为了节省资金，只要求实现任意2个小岛之间有路通即可。其中桥的价格为 100元/米。

 

 

Input

输入包括多组数据。输入首先包括一个整数T(T <= 200)，代表有T组数据。
每组数据首先是一个整数C(C <= 100),代表小岛的个数，接下来是C组坐标，代表每个小岛的坐标，这些坐标都是 0 <= x, y <= 1000的整数。

 

 

Output

每组输入数据输出一行，代表建桥的最小花费，结果保留一位小数。如果无法实现工程以达到全部畅通，输出”oh!”.

 

 

Sample Input

2 2 10 10 20 20 3 1 1 2 2 1000 1000

 

 

Sample Output

1414.2 oh!

 

 

Author

8600

 

 

Source

2008浙大研究生复试热身赛（2）——全真模拟

 

问题链接：HDU1875 畅通工程再续

问题描述：（略）

问题分析：

  这是一个最小生成树的为问题，解决的算法有Kruskal（克鲁斯卡尔）算法和Prim（普里姆） 算法。

程序说明：

  本程序使用Kruskal算法实现。有关最小生成树的问题,使用克鲁斯卡尔算法更具有优势，只需要对所有的边进行排序后处理一遍即可。程序中使用了并查集，用来判定加入一条边后会不会产生循环。程序中，图采用边列表的方式存储，排序一下就好了。

  这个题给的是坐标，需要计算距离。

AC的C++语言程序如下：

/* HDU1875 畅通工程再续 */

#include <iostream>
#include <algorithm>
#include <stdio.h>
#include <math.h>

using namespace std;

const int N = 100;

int f[N * (N - 1) + 1], fcnt;

void UFInit(int n)
{
    for(int i = 1; i <=n; i++)
        f[i] = i;
    fcnt = n;
}

int Find(int a) {
    return a == f[a] ? a : f[a] = Find(f[a]);
}

bool Union(int a, int b)
{
    a = Find(a);
    b = Find(b);
    if (a != b) {
        f[a] = b;
        fcnt--;
        return true;
    } else
        return false;
}

struct Point {
    int x, y;
} p[N];

struct Edge {
    int u, v;
    double w;
} e[N * (N - 1)];

bool cmp(Edge& a, Edge& b)
{
    return a.w - b.w < 1e-6;
}

int main()
{
    int t, c;
    scanf("%d", &t);
    while(t--) {
        scanf("%d", &c);
        for(int i = 0; i < c; i++)
            scanf("%d%d", &p[i].x, &p[i].y);

        int cnt = 0;
        for(int i = 0; i < c; i++)
            for(int j = i + 1; j < c; j++) {
                e[cnt].u = i + 1;
                e[cnt].v = j + 1;
                e[cnt++].w = sqrt((p[i].x - p[j].x) * (p[i].x - p[j].x) + (p[i].y - p[j].y) * (p[i].y - p[j].y));
            }

        // Kruscal算法
       double ans = 0;
       sort(e, e + cnt, cmp);
       UFInit(c);
       for(int i = 0; i < cnt; i++) {
           if(e[i].w >= 10 && e[i].w <= 1000)
               if(Union(e[i].u, e[i].v))
                   ans += e[i].w;
       }

       if(fcnt == 1)
            printf("%.1f\n", ans * 100);
       else
           printf("oh!\n");
    }

    return 0;
}