【JSOI2010】【BZOJ1821】Group 部落划分 Group

Description
聪聪研究发现，荒岛野人总是过着群居的生活，但是，并不是整个荒岛上的所有野人都属于同一个部落，野人们总是拉帮结派形成属于自己的部落，不同的部落之间则经常发生争斗。只是，这一切都成为谜团了——聪聪根本就不知道部落究竟是如何分布的。 不过好消息是，聪聪得到了一份荒岛的地图。地图上标注了N个野人居住的地点（可以看作是平面上的坐标）。我们知道，同一个部落的野人总是生活在附近。我们把两个部落的距离，定义为部落中距离最近的那两个居住点的距离。聪聪还获得了一个有意义的信息——这些野人总共被分为了K个部落！这真是个好消息。聪聪希望从这些信息里挖掘出所有部落的详细信息。他正在尝试这样一种算法： 对于任意一种部落划分的方法，都能够求出两个部落之间的距离，聪聪希望求出一种部落划分的方法，使靠得最近的两个部落尽可能远离。例如，下面的左图表示了一个好的划分，而右图则不是。请你编程帮助聪聪解决这个难题。
Input
第一行包含两个整数N和K(1<=N<=1000,1<=n)，分别代表了野人居住点的数量和部落的数量。 接下来n行，每行包含两个正整数x,y，描述了一个居住点的坐标(0<=”x,” y<=”10000)。” <=”” div=”“>
Output
输出一行，为最优划分时，最近的两个部落的距离，精确到小数点后两位。
Sample Input
4 2
0 0
0 1
1 1
1 0

Sample Output
1.00

HINT

Source

JSOI2010第二轮Contest1
<=n)，分别代表了野人居住点的数量和部落的数量。

自己感受一下题面吧T_T
其实就是一个最小生成树.
之前无意间看过Oxer大爷Blog给这题写的Tag,看见题瞬间明白怎么做了(虽然并不是标算233
二分一下距离,小于这个距离的边加进去,并查集判联通.

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 1010
#define GET (ch>='0'&&ch<='9')
using namespace std;
int n,k;
double dis[MAXN][MAXN],x[MAXN],y[MAXN];
int f[MAXN];
bool use[MAXN];
double l=0,r=1000000,ans;
int T=50;
int find(int x)
{
    return f[x]==x?f[x]:f[x]=find(f[x]);
}
void Union(int a,int b)
{
    int x=find(a),y=find(b),t=rand()%2;
    if (x!=y)   {if (t) swap(x,y);f[x]=y;}
}
double sqr(double x)
{
    return x*x;
}
double Dis(double x1,double y1,double x2,double y2)
{
    return sqrt(sqr(x2-x1)+sqr(y2-y1));
}
int main()
{
    scanf("%d%d",&n,&k);
    for (int i=1;i<=n;i++)  scanf("%lf%lf",&x[i],&y[i]);
    for (int i=1;i<=n;i++)
        for (int j=i+1;j<=n;j++)    dis[i][j]=dis[j][i]=Dis(x[i],y[i],x[j],y[j]);
    while (T--)
    {
        double mid=(l+r)/2.0;int num=0;
        for (int i=1;i<=n;i++)  f[i]=i;
        memset(use,0,sizeof(use));
        for (int i=1;i<=n;i++)
            for (int j=i+1;j<=n;j++)
                if (dis[i][j]<mid)  Union(i,j);
        for (int i=1;i<=n;i++)  num+=!use[find(i)],use[find(i)]=1;
        if (num<k)  r=mid;
        else    l=mid;
    }
    printf("%.2lf\n",l);
}