1821: [JSOI2010]Group 部落划分 Group
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
0 0
0 1
1 1
1 0
Sample Output
1.00
HINT
Source
贪心处理
按照两点之间距离从小到大排序
每次合并两个相距最小的 并且 n-- 知道减到k个部落
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
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) {;}
作者:乌鸦坐飞机
出处:http://www.cnblogs.com/whistle13326/
新的风暴已经出现
怎么能够停止不前
穿越时空 竭尽全力
我会来到你身边
微笑面对危险
梦想成真不会遥远
鼓起勇气 坚定向前
奇迹一定会出现