[HAOI2006]聪明的猴子

题目描述

在一个热带雨林中生存着一群猴子，它们以树上的果子为生。昨天下了一场大雨，现在雨过天晴，但整个雨林的地表还是被大水淹没着，部分植物的树冠露在水面上。猴子不会游泳，但跳跃能力比较强，它们仍然可以在露出水面的不同树冠上来回穿梭，以找到喜欢吃的果实。

现在，在这个地区露出水面的有N棵树，假设每棵树本身的直径都很小，可以忽略不计。我们在这块区域上建立直角坐标系，则每一棵树的位置由其所对应的坐标表示(任意两棵树的坐标都不相同)。

在这个地区住着的猴子有M个，下雨时，它们都躲到了茂密高大的树冠中，没有被大水冲走。由于各个猴子的年龄不同、身体素质不同，它们跳跃的能力不同。有的猴子跳跃的距离比较远(当然也可以跳到较近的树上)，而有些猴子跳跃的距离就比较近。这些猴子非常聪明，它们通过目测就可以准确地判断出自己能否跳到对面的树上。

【问题】现已知猴子的数量及每一个猴子的最大跳跃距离，还知道露出水面的每一棵树的坐标，你的任务是统计有多少个猴子可以在这个地区露出水面的所有树冠上觅食。

输入输出格式

输入格式：

输入文件monkey.in包括：

第1行为一个整数，表示猴子的个数M(2<=M<=500)；

第2行为M个整数，依次表示猴子的最大跳跃距离（每个整数值在1--1000之间）；

第3行为一个整数表示树的总棵数N(2<=N<=1000)；

第4行至第N+3行为N棵树的坐标（横纵坐标均为整数，范围为：-1000--1000）。

（同一行的整数间用空格分开）

输出格式：

输出文件monkey.out包括一个整数，表示可以在这个地区的所有树冠上觅食的猴子数。

输入输出样例

输入样例#1： 

4
 1 2 3 4
6
0 0
1 0
1 2
-1 -1
-2 0
2 2

输出样例#1： 

3

说明

【数据规模】

对于40%的数据，保证有2<=N <=100，1<=M<=100

对于全部的数据，保证有2<=N <= 1000，1<=M=500

感谢@charlie003 修正数据

分析：

最小生成树模板题？？？em。。。本题注意数组要开够，至少1000*1000，其次就是在模板的基础上加一个变量统计答案即可。

CODE：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int M=1000005;
int xcor[M],ycor[M],ma[M];
int fa[M],n,m;
int tot;
struct node{
    int u,v,w;
}a[M];
inline int get(){
    char c=getchar();
    int res=0,f=1;
    while (c<'0'||c>'9') {
        if (c=='-') f=-1;
        c=getchar();
    }
    while (c>='0'&&c<='9'){
        res=(res<<3)+(res<<1)+c-'0';
        c=getchar();
    }
    return res*f;
}
int findr(int x){
    if (fa[x]==x) return x;
    return fa[x]=findr(fa[x]);
}
void merge (int x,int y){
    int fx=findr(x);
    int fy=findr(y);
    if (fa[fx]!=fa[fy]) fa[fy]=fx;
    return ;
}
bool cmp(node x,node y){return x.w<y.w;}
int main(){
    m=get();
    for (int i=1;i<=m;i++) ma[i]=get();
    n=get();
    for (int i=1;i<=n;i++){
        xcor[i]=get(),ycor[i]=get();
        fa[i]=i;
    }
    for (int i=1;i<=n;i++){
        for (int j=i+1;j<=n;j++){
            int cost=(xcor[i]-xcor[j])*(xcor[i]-xcor[j])+(ycor[i]-ycor[j])*(ycor[i]-ycor[j]);
            a[++tot].u=i;
            a[tot].v=j;
            a[tot].w=cost;
        }
    }
    sort(ma+1,ma+m+1);
    sort(a+1,a+tot+1,cmp);
    //for (int i=1;i<=tot;i++) cout<<a[i].u<<" "<<a[i].v<<" "<<a[i].w<<endl;
    int ans=m,j=1;
    for (int i=1;i<=tot;i++){
        if (fa[findr(a[i].u)]!=fa[findr(a[i].v)]){
            merge(a[i].u,a[i].v);
            for (;ma[j]*ma[j]<a[i].w;j++) ans--;
            //cout<<i<<" "<<ans<<endl;
        }
    }
    if (ans<0) cout<<0<<endl;
    else cout<<ans<<endl;
    //system("pause");
    return 0;
}