1821: [JSOI2010]Group 部落划分 Group

Posted 2020-10-05 心之所向 素履以往

                                                    1821: [JSOI2010]Group 部落划分 Group

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

HINT

 

Source

JSOI2010第二轮Contest1

 

贪心处理

按照两点之间距离从小到大排序

 

每次合并两个相距最小的 并且 n-- 知道减到k个部落

 

 1 #include <cmath>
 2 #include <cstdio>
 3 #include <cctype>
 4 #include <cstring>
 5 #include <algorithm>
 6 
 7 using namespace std;
 8 
 9 const int MAXN=1010;
10 
11 int n,k,cnt,tot;
12 
13 int x[MAXN],y[MAXN],fa[MAXN];
14 
15 struct node {
16     int x,y;
17     double dis;
18     bool operator < (const node&a) const {
19         return dis<a.dis;
20     }
21 };
22 node e[MAXN*MAXN];
23 
24 inline void read(int&x) {
25     int f=1;register char c=getchar();
26     for(x=0;!isdigit(c);c==‘-‘&&(f=-1),c=getchar());
27     for(;isdigit(c);x=x*10+c-48,c=getchar());
28     x=x*f;
29 }
30 
31 inline double c(int i,int j) {
32     return sqrt((double)(x[i]-x[j])*(x[i]-x[j])+(double)(y[i]-y[j])*(y[i]-y[j]));
33 }
34 
35 inline int find(int x) {
36     if(x==fa[x]) return x;
37     return fa[x]=find(fa[x]);
38 }
39 
40 int hh() {
41 //    freopen("people.in","r",stdin);
42 //    freopen("people.out","w",stdout);
43     read(n);read(k);
44     for(int i=1;i<=n;++i) read(x[i]),read(y[i]),fa[i]=i;
45     for(int i=1;i<=n;++i)
46       for(int j=i+1;j<=n;++j) {
47             e[++cnt].x=i;
48             e[cnt].y=j;
49             e[cnt].dis=c(i,j);
50         }
51     sort(e+1,e+1+cnt);
52     n=n-k+1;
53     for(int i=1;i<=cnt;++i) {
54         int xx=find(e[i].x);
55         int yy=find(e[i].y);
56         if(xx!=yy) {
57             fa[xx]=yy;
58             ++tot;
59             if(tot==n) {
60                 printf("%.2lf\n",e[i].dis);
61                 goto END;
62             }
63         }
64     }
65 END:
66     return 0;
67 }
68 
69 int sb=hh();
70 int main(int argc,char**argv) {;}

代码