这道题可以用kruskal做,不过我觉得二分答案复杂度第一些?
二分一个答案mid,检查的时候遍历所有点对,如果i,j的距离比mid小,那么他们一定在一个联通块内
此时让fa[find(i)]=ffind(j)
注意因为find(i)会变,所以要每次重新算,不能在前面先算find(i),因为这个wa了无数次
二分时可以二分平方保证精度最后再开根
代码如下:
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define maxn 1005
using namespace std;
int n,k,fa[maxn],x[maxn],y[maxn];
double ans;
const double eps=1e-4;
inline int rd(){
int x=0,f=1;char c=' ';
while(c<'0' || c>'9') {if(c=='-')f=-1;c=getchar();}
while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
return x*f;
}
inline int find(int x){
if(x==fa[x]) return x;
return fa[x]=find(fa[x]);
}
inline bool check(double mid){
for(int i=1;i<=n;i++) fa[i]=i;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if((y[i]-y[j])*(y[i]-y[j])+(x[i]-x[j])*(x[i]-x[j])<=mid)
fa[find(i)]=find(j);//这里写错wa了满屏
int num=0;
for(int i=1;i<=n;i++)
if(find(i)==i) num++;
if(num<k) return false;
return true;
}
int main(){
n=rd(); k=rd();
double mx=0.0,my=0.0;
for(int i=1;i<=n;i++)
x[i]=rd(),y[i]=rd(),
mx=mx>x[i]?mx:x[i],my=my>y[i]?my:y[i];//求上界
double l=0.0,r=mx*mx+my*my;
while(r-l>eps){//二分平方保证精度
double mid=(l+r)/2.0;
if(check(mid)) l=mid;
else r=mid;
}
printf("%.2lf",sqrt(l));
return 0;
}