bzoj1821: [JSOI2010]Group 部落划分 Group

题面在这里

题意：

平面上n个点要划分为k个部落，定义部落间的距离为部落中最近的两点距离。
现在要求部落间最近的距离最远。求这个距离。

做法：

显然可以二分答案。
二分这个最近的距离，于是任意两点距离比它还小的都要在一个部落里。
然后用并查集数一下连通块个数即可。
时间复杂度$O(n^2\times{log})$

代码：

/*************************************************************
    Problem: bzoj 1821 [JSOI2010]Group 部落划分 Group
    User: fengyuan
    Language: C++
    Result: Accepted
    Time: 1676 ms
    Memory: 1312 kb
    Submit_Time: 2018-01-05 16:28:29
*************************************************************/

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cctype>
#include<cstdlib>
#define x first
#define y second
#define sqr(x) ((x)*(x))
#define PII pair<int, int>
using namespace std;
typedef long long LL;

const int N = 1010;
int n, K;
PII a[N];
int fa[N];

inline int getFa(int v) { return fa[v] == v ? v : fa[v] = getFa(fa[v]); }
inline double dis(PII a, PII b) { return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y)); }
inline bool check(double x)
{
    int sum = n; for(int i = 1; i <= n; i ++) fa[i] = i;
    for(int i = 1; i <= n; i ++)
        for(int j = i+1; j <= n; j ++) if(dis(a[i], a[j]) <= x) {
            int fx = getFa(i), fy = getFa(j);
            if(fx != fy) fa[fx] = fy, sum --;
        }
    return sum >= K;
}
int main()
{
    scanf("%d%d", &n, &K);
    for(int i = 1; i <= n; i ++) scanf("%d%d", &a[i].x, &a[i].y);
    double l = 0, r = 1e6, mid;
    for(int i = 0; i < 70; i ++) {
        mid = (l+r)/2;
        if(check(mid)) l = mid+0.01; else r = mid-0.01;
    }
    mid = (l+r)/2;
    printf("%.2lf\n", mid);
    return 0;
}