最小生成树 2429: [HAOI2006]聪明的猴子

BZOJ 2429: [HAOI2006]聪明的猴子

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 877 Solved: 566
[Submit][Status][Discuss]

Description

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

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

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

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

Input

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

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

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

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

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

Output

包括一个整数，表示可以在这个地区的所有树冠上觅食的猴子数

Sample Input

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

Sample Output

HINT

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

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

/*
思路：这个题目还是比较明显的，因为最小生成树的性质满足最大边最小，“能所有树冠上觅食”，所以一定这个猴子一定可以经过最小生成树的最大边，因为没有其他更小的距离来把全部的树都走一遍了。
只需记录最小生成树的最大边，统计有多少只猴子能够跨越这个最大边即可。
*/

 1 /*kruskal算法过程：建完边之后，把所有边从小到大拍一个序，按这个顺序往最小生成树里面加边，用并查集维护关系。*/
 2 #include<algorithm>
 3 #include<cmath>
 4 #include<iostream>
 5 using namespace std;
 6 #include<cstdio>
 7 #define N 1005
 8 #define M 505
 9 int t,n,m,dist[M],father[N];
10 struct Zb{
11     int x,y;
12 }zb[N];
13 struct Edge{
14     int u,v;
15     double w;
16     bool operator<(Edge Q)
17     const{return w<Q.w;}
18 }edge[N*N];
19 double ans;
20 void input()
21 {
22     scanf("%d",&m);
23     for(int i=1;i<=m;++i)
24     {
25         scanf("%d",&dist[i]);
26     }
27     scanf("%d",&n);
28     for(int i=1;i<=n;++i)
29     {
30         scanf("%d%d",&zb[i].x,&zb[i].y);
31     }
32 }
33 void add_edge(int a,int b)
34 {
35     ++t;
36     edge[t].u=a;edge[t].v=b;
37     edge[t].w=sqrt((zb[a].x-zb[b].x)*(zb[a].x-zb[b].x)+(zb[a].y-zb[b].y)*(zb[a].y-zb[b].y));
38 }
39 void build_edge()
40 {
41     for(int i=1;i<=n;++i)
42       for(int j=i+1;j<=n;++j)
43       {
44           add_edge(i,j);
45       }
46 }
47 int find(int x)
48 {
49     return(father[x]==x?x:father[x]=find(father[x]));
50 }
51 void kruskal()
52 {
53     for(int i=1;i<=n;++i)  father[i]=i;
54     sort(edge+1,edge+t+1);
55     int num=0;
56     for(int l=1;l<=t;++l)
57     {
58         int x1=find(edge[l].u);
59         int y1=find(edge[l].v);
60         if(x1==y1) continue;
61         num++;
62         father[y1]=x1;
63         if(num==n-1)
64         {
65             ans=edge[l].w;
66             break;
67         }
68     }
69 }
70 int main()
71 {
72     input();
73     build_edge();
74     kruskal();
75     int ret=0;
76     for(int i=1;i<=m;++i)
77     {
78         if(dist[i]>=ans) ret++;
79     }
80     printf("%d",ret);
81     return 0;
82 }