喷水装置（二）

喷水装置（二）

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

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. }