喷水装置（二）

喷水装置（二）

时间限制：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

样例输出

12

01.#include <stdio.h>

02.#include <algorithm>

03.#include <math.h>

04.using namespace std;

05.struct use

06.{

07.double left;

08.double right;

09.}a[10005];

10.int cmp(use l,use r)

11.{

12.if(l.left<r.left) return 1;

13.if(l.left==r.left&&l.right>r.right) return 1;

14.else return 0;

15.}

16.int main()

17.{

18.int i,n,w,h,x,r,max,t;

19.scanf("%d",&t);

20.while(t--)

21.{

22.scanf("%d %d %d",&n,&w,&h);

23.for(i=0;i<n;i++)

24.{

25.scanf("%d %d",&x,&r);

26.if(r>h/2)

27.{

28.a[i].left=x-sqrt(r*r-h*h/4);

29.a[i].right=x+sqrt(r*r-h*h/4);//如果圆能和矩形有交点，算出左右交点 

30.}

31.else

32.{

33.a[i].left=x;

34.a[i].right=x; //没有交点则把这个圆变成点，有效洒水距离为0 

35.}

36.}

37.sort(a,a+n,cmp);

38.int ans=0,sum=0,flag=1;

39.while(sum<w)

40.{

41.max=0;

42.for(i=0;i<n&&a[i].left<=sum;i++)

43.{

44.if(a[i].right-sum>max)

45.max=a[i].right-sum;

46.}//每次都找出符合 a[i].left<=sum（即能和上一个找出的圆的右边接在一起）的圆中，有效洒水长度最长的圆，计算出新的洒水距离 

47.if(max==0)

48.{

49.flag=0;

50.break;

51.}//如果新找出的最优圆有效洒水长度为0并且sum仍没有达到w（草坪宽度），那往后就不可能完成覆盖全部的草坪，退出循环，flag=0 

52.else

53.{

54.sum=sum+max;//每次都加上最长有效洒水距离，即算出总的洒水有效长度

55.ans++;

56.} 

57.}

58.if(flag) printf("%d\n",ans);

59.else printf("0\n");

60.}

61.return 0;

62.}