洛谷P2504 [HAOI2006]聪明的猴子

#include<bits/stdc++.h>
using namespace std;
const int maxn=5050;
const int maxm=5000005;
int n,m,k,ans;
double Max;
int monkey[maxn];
typedef pair<int,int> P;
P point[maxn];
#define Dis(x,y) (sqrt((x.first-y.first)*(x.first-y.first)+(x.second-y.second)*(x.second-y.second)))
struct node{
    int cnt,fa;
}f[maxn];
int find(int x){return f[x].fa==x?x:f[x].fa=find(f[x].fa);}
void Union(int x,int y)
{
    x=find(x),y=find(y);if(x==y) return;
    if(f[x].cnt<=f[y].cnt) {f[x].fa=y;f[y].cnt+=f[x].cnt;}
    else {f[y].fa=x;f[x].cnt+=f[y].cnt;} 
}
struct edge{
    int x,y;double dis;
}E[maxm];
inline bool cmp(const edge &a,const edge &b) {return a.dis<b.dis;}
int main()
{
    scanf("%d",&k);
    for(int i=1;i<=k;i++) scanf("%d",&monkey[i]);
    sort(monkey+1,monkey+1+k);
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d%d",&point[i].first,&point[i].second);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=i;j++)
        {
            if(i==j) continue;
            E[++m].x=i;E[m].y=j;E[m].dis=Dis(point[i],point[j]);
        }
    for(int i=1;i<=n;i++) f[i].fa=i,f[i].cnt=1;
    sort(E+1,E+1+m,cmp);
    for(int i=1;i<=m;i++)
    {
        if(f[find(1)].cnt==n) break;
        if(find(E[i].x)!=find(E[i].y)) Union(E[i].x,E[i].y),Max=max(Max,E[i].dis);
    }
    for(int i=1;i<=k;i++) if(monkey[i]>=Max) ++ans;
    printf("%d",ans);
    return 0;
}