poj -1328 Radar Installation

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

海洋在海岸线上面，陆地在海岸线下面，然后有一些岛屿分布在海洋的各个地方，有一些雷达在海岸线上，覆盖半径为d，问至少需要多少雷达才能全部覆盖这些岛屿。

以岛屿坐标为圆心，d为半径的圆与x轴的两个交点之间的范围就是雷达在这个范围内任何一个地方都能覆盖这个岛屿，求出每一条线段，然后按起点排序，贪心求解。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <vector>
#include <cstring>
#include <string>
#include <algorithm>
#include <string>
#include <set>
#include <functional>
#include <numeric>
#include <sstream>
#include <stack>
#include <map>
#include <queue>

#define CL(arr, val)    memset(arr, val, sizeof(arr))

#define ll long long
#define inf 0x7f7f7f7f
#define lc l,m,rt<<1
#define rc m + 1,r,rt<<1|1
#define pi acos(-1.0)

#define L(x)    (x) << 1
#define R(x)    (x) << 1 | 1
#define MID(l, r)   (l + r) >> 1
#define Min(x, y)   (x) < (y) ? (x) : (y)
#define Max(x, y)   (x) < (y) ? (y) : (x)
#define E(x)        (1 << (x))
#define iabs(x)     (x) < 0 ? -(x) : (x)
#define OUT(x)  printf("%I64d
", x)
#define lowbit(x)   (x)&(-x)
#define Read()  freopen("a.txt", "r", stdin)
#define Write() freopen("b.txt", "w", stdout);
#define maxn 1000000000
#define N 1010
using namespace std;

struct point
{
    double x,y;
}p[1010];
bool cmp(const point &a,const point &b)
{
    return a.x<b.x;
}

int main()
{
  //Read();
   //Write()
   int n,d,a,b,j=1;
   while(~scanf("%d%d",&n,&d))
   {
       if(n==0&&d==0) break;
       bool flag=0;
       for(int i=0;i<n;i++)
       {
           scanf("%d%d",&a,&b);
           if(b>d||b<0) flag=1;
           p[i].x=a-sqrt(1.0*(d*d-b*b));
           p[i].y=a+sqrt(1.0*(d*d-b*b));
       }
       if(flag) {printf("Case %d: -1
",j++); continue;}
       sort(p,p+n,cmp);
       // for(int i=0;i<n;i++)
       // printf("%lf %lf
",p[i].x,p[i].y);
        double ans=p[0].y;
        int m=1;
       for(int i=1;i<n;i++)
       {
           if(p[i].y<ans) ans=p[i].y;  //这里需要注意，ans必须更新保证覆盖当前终点小的岛屿。
           else if(p[i].x>ans)
           {
               m++;
               ans=p[i].y;
           }
       }
       printf("Case %d: %d
",j++,m);
   }
   return 0;
}