POJ 1328 Radar Installation 贪心 难度:1

http://poj.org/problem?id=1328

思路:

1.肯定y大于d的情况下答案为-1,其他时候必定有非负整数解

2.x,y同时考虑是较为麻烦的,想办法消掉y,用d^2-y^2获得圆心允许范围,问题转化为在许多圆心允许范围内取尽可能少的点,也即在许多线段上取尽可能少的点,使得所有线段上都有点被取到

3.从左往右考虑,完全在左边的线段肯定要取点,如果这个点在当前线段上已经取了,明显就可以忽略当前线段,明显在线段上的最优点是右端点

 #include <iostream>
 #include <cmath>
 using namespace std;

 typedef pair<double,double> p;
 p pos[];
 p temp[];
 int n,d;

 int MERGE(int s,int e){
     int mid=(e+s)/;
     int i=s,j=mid,k=s;
     while(i<mid&&j<e){
         if(pos[i].first>pos[j].first||(pos[i].first==pos[j].first&&pos[i].second>pos[j].second)){
             temp[k].first=pos[j].first;
             temp[k].second=pos[j].second;
             k++;
             j++;
         }
         else {
                 temp[k].first=pos[i].first;
                 temp[k].second=pos[i].second;
                 k++;
                 i++;
         }
     }
     int kk=k;
     for(;i<mid;i++){
         pos[kk].first=pos[i].first;
         pos[kk].second=pos[i].second;
         kk++;
     }
     for(int ii=s;ii<k;ii++){
         pos[ii].first=temp[ii].first;
         pos[ii].second=temp[ii].second;
     }
     return ;
 }
 int MERGE_SORT1(int s,int e){
     if(e-s<){
         return ;
     }
     int mid=(s+e)/;
     MERGE_SORT1(s,mid);
     MERGE_SORT1(mid,e);
     MERGE(s,e);
     return ;
 }

 int solve(){
     double x,y;
     double sq;
     bool flag=false;
     for(int i=;i<n;i++){
         cin>>x>>y;
         if(y>d){
             flag=true;
             continue;
         }
         if(y<){
             i--;
             n--;
             continue;
         }
         sq=sqrt((double)(d*d-y*y));
         pos[i].first=x-sq;
         pos[i].second=x+sq;
     }
     if(flag)return ;
     MERGE_SORT1(,n);
     int ans=;
     double mx;
     mx=pos[].second;
     for(int i=;i<n;i++){
         if(mx<pos[i].first){
             ans++;
             mx=pos[i].second;
         }
         else if(mx>pos[i].second){
             mx=pos[i].second;
         }
     }
     return ans;
 }

 int main(){
     int ci=;
     while(ci++){
         cin>>n>>d;
         if(n<=)break;
         int ans=solve();
         if(ans){
             cout<<"Case "<<ci-<<":"<<" "<<ans<<endl;
         }
         else{
             cout<<"Case "<<ci-<<":"<<" "<<-<<endl;
          }
     }
     return ;
 }