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

来源

《算法艺术与信息学竞赛》

代码：

#include<iostream>
#include<cmath>
#include<algorithm>
#include<stdio.h>
#include<fstream>
using namespace std;
struct radius
{
  int a;
  int b;
}rad[10001];
bool cmp(radius x,radius y)
{
  if(x.a<y.a) return true;
  return false;
}
int main()
{
  int t;
  int n,w,h;  //n表示有多少喷水装置，w表示长，h表示宽 
  int a,b;
//  ifstream cin("test.txt");
  cin>>t;
  while(t--)
  {
    cin>>n>>w>>h;
    for(int i=0;i<n;i++)
    {
      cin>>a>>b;
      double s=(double)pow(b,2)-pow(h/2,2);
      if(s<0)
      {
        rad[i].a=0;
        rad[i].b=0;
      }
      else
      {
        rad[i].a=a-sqrt(s);
        rad[i].b=a+sqrt(s);   
      }
    }
    sort(rad,rad+n,cmp);
    for(int i=0;i<n;i++)
      cout<<rad[i].a<<" "<<rad[i].b<<endl;
    double sum=0;
    double max;
    int flag=1;
    int count=0;
    while(sum<w)
    {
      max=0;
      for(int i=0;i<n&&rad[i].a<=sum;i++)
      {
        if(rad[i].b-sum>max) max=rad[i].b-sum;
      } 
     
      if(max==0)
      {
        flag=0;
        break;   
      }
      else
      {
        sum+=max;count++;
      }
     
    }
    if(flag) cout<<count<<endl;
    else cout<<0<<endl;
  }
  return 0;
}