BZOJ 1821 JSOI2010 Group 部落划分 Group

Posted 2020-10-18 zhangenming

1821: [JSOI2010]Group 部落划分 Group

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 3064  Solved: 1449
[Submit][Status][Discuss]

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的时候，如果再对于每一条边，如果不是同一个联通块的话，就直接输出这条边的边权即可

 1 #include <bits/stdc++.h>
 2 #define ll long long
 3 using namespace std;
 4 inline int read(){
 5     int x=0;int f=1;char ch=getchar();
 6     while(!isdigit(ch)) {if(ch==‘-‘) f=-1;ch=getchar();}
 7     while(isdigit(ch)) {x=x*10+ch-‘0‘;ch=getchar();}
 8     return x*f;
 9 }
10 const int MAXN=1e6+10;
11 struct node{
12     int x,y;
13     double v;
14 }e[MAXN];
15 int f[MAXN],x[MAXN],y[MAXN],cnt;
16 inline bool mycmp(node n,node m){
17     return n.v<m.v;
18 }
19 inline int find(int x){
20     return x==f[x]?x:f[x]=find(f[x]);
21 }
22 int main(){
23     int n=read();int m=read();
24     for(int i=1;i<=n;i++){
25         f[i]=i;
26     }
27     for(int i=1;i<=n;i++){
28         x[i]=read();y[i]=read();
29     }
30     for(int i=1;i<=n;i++){
31         for(int j=i+1;j<=n;j++){
32             double dis=sqrt(1.0*(x[i]-x[j])*(x[i]-x[j])+1.0*(y[i]-y[j])*(y[i]-y[j]));
33             e[++cnt].x=i;e[cnt].y=j;e[cnt].v=dis;
34         }
35     }
36     sort(e+1,e+cnt+1,mycmp);
37     for(int i=1;i<=cnt;i++){        
38     int fx=find(e[i].x);int fy=find(e[i].y);
39     //cout<<fx<<‘ ‘<<fy<<endl;
40         if(fx!=fy){
41             if(n>m){
42                 n--;f[fx]=fy;
43             }
44             else{
45                printf("%.2f\n",e[i].v);return 0;
46             }
47         }
48     }
49     return 0;
50 }

View Code