[BZOJ1328]Radar Installation（贪心）

传送门

---

题意：有n个海岛，分别分布在坐标系上，在x轴上放雷达，使得雷达可以覆盖到所有的海岛，雷达的辐射范围是一个圆有一个固定的辐射半径R，问最少需要打多少个雷达可以覆盖所有的海岛

对于每一个海岛，可以计算出x轴上的一段能管辖它的区间，计算方法如下：

用一个点的横坐标减去/加去$\sqrt{R^2-y^2}$就的得到了范围[l,r]，这样子[l,r]之间的点画圆都可以包括这个点（为什么，半径嘛…直角三角形斜边最长嘛）。
那么无解的情况就是R够不到这个点以及R是负数（smg…
然后先按照每个点的$l[i]$从小到大排序
那么对于每个区间，如果当前区间$i$的左端点$l[i]$大于最后一台雷达的坐标$pos$，那么也就是说最后一台无法够到它，那么就雷达，然后更新$pos=r[i]$。否则就够的到，让最后一台来管辖这个点，然后维护$pos=min(r[i],pos)$（因为实际上我们没有建这个雷达，所以我们并不能够用pos=r[i]来更新他）

---

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=1010;
struct node
{
    double x,y;
    double l,r;//直线上l到r的区域管辖他
}a[maxn];
int cmp(node a,node b){return a.l<b.l;}
int main()
{
    int n; double R; int ti=0;
    while(scanf("%d%lf",&n,&R)!=EOF)
    {
        if(n==0 && R==0) break; 
        ti++; bool flag=true;
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf",&a[i].x,&a[i].y);
            double d=sqrt(R*R-a[i].y*a[i].y);
            a[i].l=a[i].x-d,a[i].r=a[i].x+d; 
            if(a[i].y>R || a[i].y<-R) flag=false;
        }
        sort(a+1,a+n+1,cmp);
        double pos=-1e10; //已经安放的最后一台监控设备的坐标 
        int ans=0;
        for(int i=1;i<=n;i++)
        {
            if(a[i].l>pos)//如果最后一台无法检测到那么新建一台 
            {
                ans++;
                pos=a[i].r;
            }
            else pos=min(a[i].r,pos);//可以检测到不需新建，维护pos 
        }
        if(R<0 || !flag) printf("Case %d: -1\n",ti);
        else printf("Case %d: %d\n",ti,ans);
    }
    return 0;
}