[JSOI2010][BZOJ1821]Group 部落划分 Group

1821: [JSOI2010]Group 部落划分 Group

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 1477  Solved: 709
[Submit][Status][Discuss]

Description

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

Input

第一行包含两个整数N和K(1<=N<=1000,1<k<=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

 

分析：

一开始看到最大最小还以为是二分，后来发现直接暴力就可以搞过……

把所有的点对的距离都存下来，sort一下，然后从小到大往里加，知道加的点对数=n-k时，就不能继续加了。这时不能把edge[i+1]输出，而要继续往后扫，扫到第一个可以继续往里加的点对，把距离输出就好了。

 

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
int n,k,top,x,y,ans=0,a[1001][3],f[1001];
double dis;
struct node
{
       int x,y;
       double d;
};
node edge[500005];
double calc(int x1,int y1,int x2,int y2)
{
       return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
void insert(int xx,int yy,double dd)
{
     edge[++top].x=xx;
     edge[top].y=yy;
     edge[top].d=dd;
}   
bool cmp(node xx,node yy)
{
     return xx.d<yy.d;
}
int find(int x)
{
    return x==f[x]?x:find(f[x]);
}
int main()
{
    scanf("%d%d",&n,&k);
    for (int i=1;i<=n;i++)
        scanf("%d%d",&a[i][1],&a[i][2]);
    top=0;
    for (int i=1;i<=n;i++)
        for (int j=i+1;j<=n;j++)
        {
            dis=calc(a[i][1],a[i][2],a[j][1],a[j][2]);
            insert(i,j,dis);
        }
    sort(edge+1,edge+top+1,cmp);
    for (int i=1;i<=n;i++) f[i]=i;
    for (int i=1;i<=top;i++)
    {
        int p=find(edge[i].x);
        int q=find(edge[i].y);
        if (p!=q)
        {
                 f[p]=q;
                 ans++;
        }
        if (ans==n-k) 
                      for (int j=i+1;j<=top;j++)
                      {
                          int p=find(edge[j].x),q=find(edge[j].y);
                          if (p!=q) printf("%.2lf",edge[j].d);
                          break;
                      }
    }
    return 0;
}