BZOJ1821 [JSOI2010] Group 部落划分 Group

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1821

Description

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

Input

第一行包含两个整数N和K(1< = N < = 1000,1< K < = N)，分别代表了野人居住点的数量和部落的数量。

接下来N行，每行包含两个正整数x,y，描述了一个居住点的坐标(0 < =x, y < =10000)

Output

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

贪心，把边权小的边划分到部落内部，第n - k + 1条边的边权就是答案。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#define rep(i,l,r) for(int i=l; i<=r; i++)
#define clr(x,y) memset(x,y,sizeof(x))
#define travel(x) for(Edge *p=last[x]; p; p=p->pre)
using namespace std;
const int INF = 0x3f3f3f3f;
const int maxn = 1010;
inline int read(){
    int ans = 0, f = 1;
    char c = getchar();
    for(; !isdigit(c); c = getchar())
    if (c == '-') f = -1;
    for(; isdigit(c); c = getchar())
    ans = ans * 10 + c - '0';
    return ans * f;
}
struct Edge{
    int from,to;
    double cost;
    inline bool operator < (const Edge &_Tp) const{
        return cost < _Tp.cost;
    }
}edge[500010];
struct Point{
    int x,y;
}t[maxn];
int n,k,cnt=0,tot=0,fa[maxn];
inline double getdis(Point a,Point b){
    return sqrt((a.x-b.x) * (a.x-b.x) + (a.y-b.y) * (a.y-b.y));
}
inline void addedge(int x,int y,double dis){
    edge[++cnt].from = x; edge[cnt].to = y; edge[cnt].cost = dis;
}
int getfa(int x){
    return fa[x] == x ? x : fa[x] = getfa(fa[x]);
}
int main(){
    n = read(); k = read();
    rep(i,1,n) t[i].x = read(), t[i].y = read(), fa[i] = i;
    rep(i,1,n-1) rep(j,i+1,n) addedge(i,j,getdis(t[i],t[j]));
    sort(edge+1,edge+cnt+1);
    k = n - k;
    rep(i,1,cnt){
        int a = getfa(edge[i].from), b = getfa(edge[i].to);
        if (a == b) continue;
        fa[a] = b;
        if (++tot > k){
            printf("%.2lf\n",edge[i].cost);
            return 0;
        }
    }
    return 0;
}