【贪心】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=1000+5;
int n,d;
rec inter[MAXN];

int calculate()
{
    sort(inter,inter+n);
    int ans=1;
    int before=0;
    for (int i=1;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=0;
    while (scanf("%d%d",&n,&d))
    {
        if (n==d && d==0) break;
        t++;
        int f=1;
        for (int i=0;i<n;i++)
        {
            double x,y;
            cin>>x>>y;
            if (d<y)
            {
                f=0;    
            }
            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<<-1<<endl;
            
    }
    return 0;
}