洛谷P2504 [HAOI2006]聪明的猴子

 #include<bits/stdc++.h>
 using namespace std;
 const int maxn=;
 const int maxm=;
 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=;i<=k;i++) scanf("%d",&monkey[i]);
     sort(monkey+,monkey++k);
     scanf("%d",&n);
     for(int i=;i<=n;i++) scanf("%d%d",&point[i].first,&point[i].second);
     for(int i=;i<=n;i++)
         for(int j=;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=;i<=n;i++) f[i].fa=i,f[i].cnt=;
     sort(E+,E++m,cmp);
     for(int i=;i<=m;i++)
     {
         if(f[find()].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=;i<=k;i++) if(monkey[i]>=Max) ++ans;
     printf("%d",ans);
     return ;
 }