P4047 [JSOI2010]部落划分
斯德哥尔摩
2018-08-02 09:07:57
[P4047 [JSOI2010]部落划分](https://www.luogu.org/problemnew/show/P4047)
本蒟蒻表示没有好的法子,直接 二分答案+并查集 水过。。。
注意到这句话:靠得最近的两个部落尽可能远离。
不就是二分嘛!
二分最远距离的平方,并查集维护部落。
还有,距离最好不要上来就开根号,不然会丢失精度的。
注:图片上有彩蛋$JSOI$哦。
附代码:
```cpp
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#define MAXN 1010
#define eps (1e-6)
using namespace std;
int n,m,fa[MAXN];
double l=0,r=0,mid;
struct Point{
double x,y;
}a[MAXN];
inline int read(){
int date=0,w=1;char c=0;
while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
return date*w;
}
inline double dis(int i,int j){
return ((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));
}
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
inline void uniun(int x,int y){x=find(x);y=find(y);if(x!=y)fa[y]=x;}
int check(double x){
int s=0;
for(int i=1;i<=n;i++)fa[i]=i;
for(int i=1;i<=n;i++)
for(int j=1;j<i;j++)
if(dis(i,j)<=mid)uniun(j,i);
for(int i=1;i<=n;i++)if(find(i)==i)s++;
return s;
}
void work(){
double ans;
while(l<=r){
mid=(l+r)/2.0;
int s=check(mid);
if(s==m)ans=mid;
if(s>=m)l=mid+eps;
else r=mid-eps;
}
printf("%.2lf\n",sqrt(mid));
}
void init(){
n=read();m=read();
for(int i=1;i<=n;i++){
a[i].x=read();a[i].y=read();
for(int j=1;j<i;j++)r=max(r,dis(i,j));
}
}
int main(){
init();
work();
return 0;
}
```