【贪心】POJ1328-Radar Installation

【思路】

以每一座岛屿为圆心，雷达范围为半径作圆，记录下与x轴的左右交点。如果与x轴没交点，则直接退出输出“-1”。以左交点为关键字进行排序，从左到右进行贪心。容易知道，离每一个雷达最远的那一座岛与雷达相距恰巧为半径的时候，可以得到最优解。假设上一个雷达与第before座岛相距为半径大小，对于当前的岛屿i：

如果before岛的右交点在i岛左交点的左侧，此时必然需要一个新的雷达，将当前before暂定为i，雷达数加一；

如果before岛的右交点在i岛右交点的右侧，为了让雷达举例尽可能地广，将雷达移至当前岛屿与x轴的交点上，将当前before暂定为i，雷达数不变。

 #include<iostream>
 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 double sqrt(double n);
 struct rec
 {
     double left,right;
     bool operator < (const rec& x) const
     {
         return left<x.left;
     }
 };

 const int MAXN=+;
 int n,d;
 rec inter[MAXN];

 int calculate()
 {
     sort(inter,inter+n);
     int ans=;
     int before=;
     for (int i=;i<n;i++)
     {
         if (inter[before].right<inter[i].left)
         {
             ans++;
             before=i;
         }
         else if (inter[before].right>=inter[i].right)
         {
             before=i;
         }
     }
     return ans;
 }

 int main()
 {
     int t=;
     while (scanf("%d%d",&n,&d))
     {
         if (n==d && d==) break;
         t++;
         int f=;
         for (int i=;i<n;i++)
         {
             double x,y;
             cin>>x>>y;
             if (d<y)
             {
                 f=;
             }
             else
             {
                 inter[i].left=x*1.0-sqrt(d*d-y*y);
                 inter[i].right=x*1.0+sqrt(d*d-y*y);
             }
         }
         cout<<"Case "<<t<<": ";
         if (f) cout<<calculate()<<endl;else cout<<-<<endl;

     }
     return ;
 }