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
Solution
二分答案,然後把距離小於答案的都合併起來,再算並查集裏有幾個根,如果>=m就增加答案,否則減少答案。
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
#define db double
using namespace std;
int n,m,s,u,v;
db ans,Max;
int f[1005],a[1005],b[1005];
db dis(int x1,int y1,int x2,int y2)
{
return sqrt((db)((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)));
}
int get(int x)
{
if(f[x]==x) return x;else return f[x]=get(f[x]);
}
int main()
{
cin>>n>>m;
for(int i=1;i<=n;i++) scanf("%d%d",&a[i],&b[i]);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++) Max=max(Max,dis(a[i],b[i],a[j],b[j]));
db l=0,r=Max;
for(int w=1;w<=50;w++)
{
db mid=(l+r)/2.0;
s=0;
for(int i=1;i<=n;i++) f[i]=i;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(i!=j&&dis(a[i],b[i],a[j],b[j])<mid)
{
u=get(i);
v=get(j);
f[v]=u;
}
for(int i=1;i<=n;i++) s=s+(f[i]==i);
if(s>=m) l=mid; else r=mid;
}
printf("%0.2f",l);
return 0;
}