NYOJ 喷水装置（二）

喷水装置（二）

时间限制：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<algorithm>
#include <math.h> 
using namespace std;
#define max(a,b) (a>b?a:b)
struct Ww{
  float xi,ri;
       };
bool com(const Ww &a,const Ww &b){
     if(a.xi<b.xi) return true;
     if(a.xi==b.xi&&a.ri>b.ri) return true;
     
     return false;
      
     }
float dis(float rr,float hh){
      if(rr*rr-hh*hh/4>=0)
    return sqrt(rr*rr-hh*hh/4);  
    else
    return 0;
      
      }


int main(){
int m;
freopen("12.txt","r",stdin);
cin>>m;
while(m--){
   int n;
   cin>>n;
   Ww a[n];
   float w,h;
   cin>>w;
   cin>>h;
   for(int i=0;i<n;i++)
   cin>>a[i].xi>>a[i].ri;   
   
 sort(a,a+n,com);
 
 
  
 float sum=0;
 int count=0;
 int first=0,rear=n-1;
 
 for(int i=0;i<=n-1;i++){
   
          if(rear==first&&(a[first].ri>=sqrt(w/2+h*h/4))) { count++; break;                   
                          } 
          if(rear==first&&a[first].ri<sqrt(w/2+h*h/4)) {break;}
          
         if(dis(a[first].ri,h)>=a[first].xi){
          if(dis(a[rear].ri,h)>=w-a[rear].xi)     {                             
           if(a[i].xi!=a[i+1].xi){           
          if(dis(a[i].ri,h)+dis(a[i+1].ri,h)>=a[i+1].xi-a[i].xi){
           if(first==i) count+=2;
           else
           count++;
           
           
           }}
     
           }else{
                 rear--;
                 }
           
           
           }else 
           {
                first++;
                }
                
                    }

     
 cout<<count<<endl;
           }

      
}

 

我后来重写的：

#include<stdio.h>
#include<string.h>
#include<math.h> 
#include<algorithm>
using namespace std;
struct ps
{
    double left; 
    double right; 
}w[10001];
bool comp(ps a,ps b) 
{
    if(a.left<b.left) return true;
    return false;
}
int main()
{
       freopen("12.txt","r",stdin);
    int ncases,n,i,width,high,x,r,count,flag;
    double length,sum,max;
    scanf("%d",&ncases);
    while(ncases--)
    {
        flag=1;sum=0;count=0;
        scanf("%d %d %d",&n,&width,&high);
        for(i=0;i<=n-1;i++) 
        {
            scanf("%d %d",&x,&r);
            length=(double)r*r-(double)high/2*high/2;
            if(length>=0) {length=sqrt(length);} 
            if(length<0) {length=0;} 
            w[i].left=x-length;
            w[i].right=x+length;
        }
        sort(w,w+n,comp);
        while(sum<width) 
        {
            max=0; 
            for(i=0;i<=n-1&&w[i].left<=sum;i++) 
            {
                if((w[i].right-sum)>max) 
                {
                    max=w[i].right-sum; 
                }
            }
            if(max==0) 
            {
                flag=0;
                break;
            }
            else
            {
                sum=sum+max; 
                count++;
            }
        }
        if(flag==1)
        {
            printf("%d\n",count);
        }
        else
        {
            printf("0\n");
        }
    }
    return 0;
}

 

最优程序：

 
#include<iostream>
#include<vector>
#include<algorithm>
#include<cmath>
using namespace std;
double r,x,y,width,height;
class Circle
{
public:
    Circle(double x,double r):x(x),r(r){}
    double Left() const{
        if(r*r-height*height/4<0) return x;return x-sqrt(r*r-(height*height/4));}
    double Right() const{if(r*r-height*height/4<0) return x;return x+sqrt(r*r-(height*height/4));}
    friend ostream & operator<<(ostream& out,const Circle& circle);
private:
    double x,r;
};
ostream & operator<<(ostream& out,const Circle& circle)
{
    out<<circle.x<<" "<<circle.r;
    return out;
}
struct CircleLess
{
    bool operator()(const Circle& c1,const Circle& c2)
    {
        return c1.Right()<c2.Right();
    }
};

int main()
{

    int m;
    cin>>m;
    while(m--)
    {
        int n;
        
        vector<Circle> Rs;
        cin>>n>>width>>height;
        Rs.reserve(n);
        
        
        while(n--)
        {
            cin>>x>>r;
            Rs.push_back(Circle(x,r));
        }
        sort(Rs.begin(),Rs.end(),CircleLess());
        vector<int> choosed;
        choosed.push_back(-1);
        double len=0;
        try{
            while(len<width)
            {
                int last=-1;
                for(vector<Circle>::size_type i=choosed.back()+1;i!=Rs.size();i++)
                {
                    if (Rs[i].Left()<=len)
                    {
                        //cout<<Rs[i].Left()<<endl;
                        last=i;
                    }
                }
                if (last==-1) throw -1;
                else
                {
                    choosed.push_back(last);
                    len=Rs[last].Right();
                }
            }    
            cout<<choosed.size()-1<<endl;
        
        }catch(...)
        {
            cout<<0<<endl;
        }
    }

}