题目链接: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条边的边权就是答案。
1 #include <iostream> 2 #include <cstdio> 3 #include <algorithm> 4 #include <cstring> 5 #include <cmath> 6 #define rep(i,l,r) for(int i=l; i<=r; i++) 7 #define clr(x,y) memset(x,y,sizeof(x)) 8 #define travel(x) for(Edge *p=last[x]; p; p=p->pre) 9 using namespace std; 10 const int INF = 0x3f3f3f3f; 11 const int maxn = 1010; 12 inline int read(){ 13 int ans = 0, f = 1; 14 char c = getchar(); 15 for(; !isdigit(c); c = getchar()) 16 if (c == '-') f = -1; 17 for(; isdigit(c); c = getchar()) 18 ans = ans * 10 + c - '0'; 19 return ans * f; 20 } 21 struct Edge{ 22 int from,to; 23 double cost; 24 inline bool operator < (const Edge &_Tp) const{ 25 return cost < _Tp.cost; 26 } 27 }edge[500010]; 28 struct Point{ 29 int x,y; 30 }t[maxn]; 31 int n,k,cnt=0,tot=0,fa[maxn]; 32 inline double getdis(Point a,Point b){ 33 return sqrt((a.x-b.x) * (a.x-b.x) + (a.y-b.y) * (a.y-b.y)); 34 } 35 inline void addedge(int x,int y,double dis){ 36 edge[++cnt].from = x; edge[cnt].to = y; edge[cnt].cost = dis; 37 } 38 int getfa(int x){ 39 return fa[x] == x ? x : fa[x] = getfa(fa[x]); 40 } 41 int main(){ 42 n = read(); k = read(); 43 rep(i,1,n) t[i].x = read(), t[i].y = read(), fa[i] = i; 44 rep(i,1,n-1) rep(j,i+1,n) addedge(i,j,getdis(t[i],t[j])); 45 sort(edge+1,edge+cnt+1); 46 k = n - k; 47 rep(i,1,cnt){ 48 int a = getfa(edge[i].from), b = getfa(edge[i].to); 49 if (a == b) continue; 50 fa[a] = b; 51 if (++tot > k){ 52 printf("%.2lf ",edge[i].cost); 53 return 0; 54 } 55 } 56 return 0; 57 }