南阳ACM-喷水装置（二）

题目：

喷水装置（二）

时间限制：3000 ms  |  内存限制：65535 KB

难度：4

描述

有一块草坪，横向长w,纵向长为h,在它的橫向中心线上不同位置处装有n(n<=10000)个点状的喷水装置，每个喷水装置i喷水的效果是让以它为中心半径为Ri的圆都被润湿。请在给出的喷水装置中选择尽量少的喷水装置，把整个草坪全部润湿。

输入

第一行输入一个正整数N表示共有n次测试数据。
每一组测试数据的第一行有三个整数n,w,h，n表示共有n个喷水装置，w表示草坪的横向长度，h表示草坪的纵向长度。
随后的n行，都有两个整数xi和ri,xi表示第i个喷水装置的的横坐标（最左边为0），ri表示该喷水装置能覆盖的圆的半径。

输出

每组测试数据输出一个正整数，表示共需要多少个喷水装置，每个输出单独占一行。
如果不存在一种能够把整个草坪湿润的方案，请输出0。

样例输入

22 8 61 14 52 10 64 56 5

样例输出

1
2

个人理解：需要注意的就是圆圈里必须套个方框这个圆才是可以用的 先选大圈 贪心算法

AC情况：

代码 C语言：

# include <stdio.h>  # include <math.h>  # define N 1000  typedef double elem;  # define max(a,b) (a)>(b)?(a):(b)  void  KuaiPai(elem A[][N],int left,int right);//快速排序 自身实现  void change(elem *A,elem *B)//交换函数  {     elem C=*A;     *A=*B;     *B=C;  }  double PS[2][N],d,hw,h,MAX,sum;  int T,num,W,H,i,j,X,R;  int main(){      //freopen("SASA.txt","r",stdin);      scanf("%d",&T);      while(T--)      {          scanf("%d%d%d",&num,&W,&H);          h=H*0.5;          hw=h*h;          for(j=i=0;i<num;i++){              scanf("%d%d",&X,&R);              if(R>h){//只对符合条件的装置计数              d=sqrt(R*R-hw);                PS[0][j]=X-d;                PS[1][j++]=X+d;//符合条件的装置的个数=j               }          }          sum=num=0;//每次初始化为0          KuaiPai(PS,0,j-1);//对符合条件的装置 下标[0,j-1]          while(sum<W)          {              MAX=0;//用MAX来记录最长的更新的距离              for( i=0;i<j&&PS[0][i]<=sum;i++) //每一次遍历找最长的更新距离，并且起点至少在sum以内                  MAX=max(MAX,(PS[1][i]-sum));//筛选最大值               if(MAX){//如果最大值不是 0  说明没有断层                 num++;//记数                 sum+=MAX;               }               else break;          }          printf("%d\n",MAX?num:0);//输出      }      return 0;  }  void  KuaiPai(elem A[][N],int left,int right)  {      int i=left,j=right;      elem temp=A[0][left];//记录一下      if(left>=right)  return;//如果重合      while(i!=j)      {          while(A[0][j]>=temp && i<j)j--;//从右向左找          while(A[0][i]<=temp && i<j)i++;//从左向右找          if(i<j)          {              change(&A[0][i],&A[0][j]);//交换              change(&A[1][i],&A[1][j]);          }  }      change(&A[0][left],&A[0][i]);//换回来      change(&A[1][left],&A[1][i]);      KuaiPai(A,left,i-1);//左递归      KuaiPai(A,i+1,right);//右递归  }