bzoj [JSOI2010]Group 部落划分 Group【二分+并查集】

我是zz吗这么简单都写错……
一眼二分，然后判断的话是枚举点，然后计算这个点到已有联通块的最小距离，如果这个点到一些联通块的距离小于当前二分的val，则把这些联通块合并起来，这里用并查集维护，最后看这样得出的部落数是否大于k（多出来的可以直接合并）
有个非常小的优化就是不用double二分，直接把点两两之间的距离存起来排个序，直接在排序后数组里二分即可

#include<iostream>
#include<cstdio>
#include<cmath>
#include<vector>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1005;
int n,m,x
,y
,f
,tot;
double mn
,d[N*N];
int read()
{
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
r=r*10+p-48;
p=getchar();
}
return r*f;
}
double dis(double x1,double y1,double x2,double y2)
{
return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
int zhao(int x)
{
return x==f[x]?x:f[x]=zhao(f[x]);
}
bool ok(double va)
{
int col=0;
for(int i=1;i<=n;i++)
f[i]=i;
for(int i=1;i<=n;i++)
{
int s=0,w;
for(int j=1;j<=n;j++)
mn[j]=1e9;
for(int j=1;j<i;j++)
{
int nw=zhao(j);
mn[nw]=min(mn[nw],dis(x[i],y[i],x[j],y[j]));
}
for(int j=1;j<i;j++)
if(f[j]==j&&mn[j]<va)
s++,w=j;
if(s>1)
{
for(int j=1;j<i;j++)
if(f[j]==j&&mn[j]<va&&j!=w)
f[j]=w;
f[i]=w;
}
else if(s==1)
f[i]=w;
}
for(int i=1;i<=n;i++)
if(f[i]==i)
col++;
// printf("%.2lf %d\n",va,col);
// for(int i=1;i<=n;i++)
// cerr<<zhao(i)<<" ";
// cerr<<endl;
return col>=m;
}
int main()
{
n=read(),m=read();
for(int i=1;i<=n;i++)
{
x[i]=read(),y[i]=read();
for(int j=1;j<i;j++)
d[++tot]=dis(x[i],y[i],x[j],y[j]);
}
sort(d+1,d+1+tot);
int l=0,r=tot,ans;
while(l<=r)
{
int mid=(l+r)>>1;
if(ok(d[mid]))
l=mid+1,ans=mid;
else
r=mid-1;
}
printf("%.2lf\n",d[ans]);
return 0;
}