BZOJ2429 [HAOI2006] 聪明的猴子

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

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

下课时水一水放松身心……其实这个Kruskal是乱写的，应该加到n-1条边就break，然而也能AC就懒得改了。

1 #include <iostream> 2 #include <cstdio> 3 #include <algorithm> 4 #include <cstring> 5 #include <cmath> 6 #define rep(i,l,r) for(int i=l; i<=r; i++) 7 #define clr(x,y) memset(x,y,sizeof(x)) 8 #define travel(x) for(Edge *p=last[x]; p; p=p->pre) 9 using namespace std; 10 const int INF = 0x3f3f3f3f; 11 const int maxn = 1010; 12 inline int read(){ 13 int ans = 0, f = 1; 14 char c = getchar(); 15 for(; !isdigit(c); c = getchar()) 16 if (c == ‘-‘) f = -1; 17 for(; isdigit(c); c = getchar()) 18 ans = ans * 10 + c - ‘0‘; 19 return ans * f; 20 } 21 struct Edge{ 22 int from,to; 23 double cost; 24 inline bool operator < (const Edge &_Tp) const{ 25 return cost < _Tp.cost; 26 } 27 }edge[500010]; 28 struct Point{ 29 int x,y; 30 }t[maxn]; 31 int n,m,cnt=0,h[maxn],fa[maxn]; 32 inline double getdis(Point a,Point b){ 33 return sqrt((a.x-b.x) * (a.x-b.x) + (a.y-b.y) * (a.y-b.y)); 34 } 35 inline void addedge(int x,int y,double dis){ 36 edge[++cnt].from = x; edge[cnt].to = y; edge[cnt].cost = dis; 37 } 38 int getfa(int x){ 39 return fa[x] == x ? x : fa[x] = getfa(fa[x]); 40 } 41 int main(){ 42 m = read(); 43 rep(i,1,m) h[i] = read(); 44 n = read(); 45 rep(i,1,n) t[i].x = read(), t[i].y = read(); 46 rep(i,1,n-1) rep(j,i+1,n) addedge(i,j,getdis(t[i],t[j])); 47 sort(edge+1,edge+cnt+1); 48 rep(i,1,n) fa[i] = i; 49 double mx = -INF; 50 rep(i,1,cnt){ 51 int a = getfa(edge[i].from); 52 int b = getfa(edge[i].to); 53 if (a != b){ 54 fa[a] = b; 55 mx = max(mx,edge[i].cost); 56 } 57 } 58 int ans = 0; rep(i,1,m) if (h[i] >= mx) ans++; 59 printf("%d\n",ans); 60 return 0; 61 }