HYSBZ 1821 Group 部落划分 Group

1821: [JSOI2010]Group 部落划分 Group

Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 1513 Solved: 720
[Submit][Status][Discuss]
Description

聪聪研究发现，荒岛野人总是过着群居的生活，但是，并不是整个荒岛上的所有野人都属于同一个部落，野人们总是拉帮结派形成属于自己的部落，不同的部落之间则经常发生争斗。只是，这一切都成为谜团了——聪聪根本就不知道部落究竟是如何分布的。 不过好消息是，聪聪得到了一份荒岛的地图。地图上标注了N个野人居住的地点（可以看作是平面上的坐标）。我们知道，同一个部落的野人总是生活在附近。我们把两个部落的距离，定义为部落中距离最近的那两个居住点的距离。聪聪还获得了一个有意义的信息——这些野人总共被分为了K个部落！这真是个好消息。聪聪希望从这些信息里挖掘出所有部落的详细信息。他正在尝试这样一种算法： 对于任意一种部落划分的方法，都能够求出两个部落之间的距离，聪聪希望求出一种部落划分的方法，使靠得最近的两个部落尽可能远离。 例如，下面的左图表示了一个好的划分，而右图则不是。请你编程帮助聪聪解决这个难题。
Input

第一行包含两个整数N和K(1<=N<=1000,1
Output

输出一行，为最优划分时，最近的两个部落的距离，精确到小数点后两位。
Sample Input

4 2

0 0

0 1

1 1

1 0

Sample Output

1.00

数据规模比较小，总共50w条边，一个Kruskal搞定

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
#define inf 0xffffff
#define maxn 1005

struct Point
{
    int x,y;
};

struct road
{
    int p,q;
    int x;
};

Point P[maxn];
road R[maxn*maxn/2];
int Father[maxn],Top;

int length(Point x,Point y)
{
    return ((x.x-y.x)*(x.x-y.x) + (x.y-y.y)*(x.y-y.y));
}

bool cmp(road a,road b)
{
    return a.x<b.x;
}

int Find(int x)
{
    return x == Father[x]? x : Father[x] = Find(Father[x]);
}

bool updata(int a,int b)
{
    int x = Find(a);
    int y = Find(b);
    if(x == y) return false;
    Father[y] = x;
    return true;
}

double Kruskal(int n,int k)
{
    int flag = 0,s = 0;
    for(int i = 0; i < n; i++)
        Father[i] = i;
    for(int i =  0; i < Top; i++)
        if(updata(R[i].p,R[i].q)){
            if(n>k) n--;
            else return sqrt(R[i].x);
        }
}

int main()
{
    int n ,k;
    while(~scanf("%d%d",&n,&k))
    {
        for(int i = 0; i < n; i++)
            scanf("%d%d",&P[i].x,&P[i].y);
        if(n == 1 || k == 1) {printf("0.00\n");continue;}
        Top = 0;
        for(int i = 0; i < n; i++)
            for(int j = i + 1; j < n; j++){
                R[Top].p = i;
                R[Top].q = j;
                R[Top++].x = length(P[i],P[j]);
            }
        sort(R,R+Top,cmp);
        printf("%.2f\n",Kruskal(n,k));
    }
    return 0;
}

/*
5 2
-1 -1
1 -1
-1 2
2 1
2 2

*/