南阳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。

样例输入

2
2 8 6
1 1
4 5
2 10 6
4 5
6 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);//右递归  
}