BZOJ 1821 部落划分（二分+并查集）

Description

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

Input

第一行包含两个整数N和K(1< = N < = 1000,1< K < = N)，分别代表了野人居住点的数量和部落的数量。
接下来N行，每行包含两个正整数x,y，描述了一个居住点的坐标(0 < =x, y < =10000)
Output

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

Sample Input

4 2

0 0

0 1

1 1

1 0

Sample Output

1.00

思路：预处理出野人两两之间的距离，二分部落之间的最小距离，用并查集维护集合。如果两个野人之间的距离比当前答案小，则将其分在同一个集合中。最后统计一共分了几个部落，如果分的部落数量大于等于k，那再让答案变大，最后注意一下精度问题。

题解：

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int maxn=1005;
struct cc{
    double x,y;
}num[maxn];
double a[maxn][maxn];
int f[maxn];
int find(int w)
{
    if(f[w]!=w)
    {
        f[w]=find(f[w]);
    }
    return f[w];
}
int n,k;
bool check(double mid)
{
    for(int i=1;i<=n;i++)
    {
        f[i]=i;
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(a[i][j]<=mid)
            {
                if(find(i)!=find(j))
                {
                    f[find(i)]=find(j);
                }
            } 
        }
    }
    int tot=0;
    for(int i=1;i<=n;i++)
    {
        if(f[i]==i)
        {
            tot++;
        }
    }
    if(tot>=k)
    {
        return 1;
    }
    else
    {
        return 0;
    }
}
int main()
{

    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%lf",&num[i].x,&num[i].y);
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            a[i][j]=sqrt((num[i].x-num[j].x)*(num[i].x-num[j].x)+(num[i].y-num[j].y)*(num[i].y-num[j].y));
        }
    }
    double l=0,r=1000001;
    double ans=0;
    while(r-l>=0.0001)
    {
        double mid=(l+r)/2;
        if(check(mid))
        {
            l=mid+0.001;
            ans=max(ans,mid);
        }
        else
        {
            r=mid-0.001;
        }
    }
    printf("%.2lf",ans);
    return 0;
}