poj 1328 Radar Installation 贪心

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

X轴的上方代表海，下方代表陆地。海中有多个岛屿，X轴上可任意放置半径为d的雷达，求至少放置多少个雷达可覆盖所有岛屿，不能完成输出－1。

    以岛屿为圆心做半径为d的圆，可在X轴上产生两个交点，即为X轴上可覆盖该岛屿的区间，求出所有岛屿所对应的区间，将交集不为空的区间分为一组，这样你会发现答案其实就是求区间可分为多少组。

#include<stdio.h>

#include<stdlib.h>

#include<algorithm>

#include<iostream>

#include<math.h>

using namespace std;

struct point

{

       double x,y;

}p[1005];

struct in

{

      double l,r;

}q[1005];

int cmp(in a,in b)

{

    if(a.l-b.l<-1e-8)

    return 1;

    else

    return 0;

}

double dis(point a,point b)

{

       return sqrt(a.x*a.x+b.y*b.y);

}

int main()

{

    int i,j,ans,n,d,flag,cas=1;

    double x,a,b;

    while(scanf("%d%d",&n,&d),n+d)

    {

         flag=0;ans=1;

         for(i=0;i<n;i++)

         {

              scanf("%lf%lf",&p[i].x,&p[i].y);

              if(p[i].y>d)

              flag=1;

         }

         if(flag)

         {

              printf("Case %d: -1\n",cas++);

              continue;

         }

         for(i=0;i<n;i++)

         {

              q[i].l=p[i].x-sqrt(d*d-p[i].y*p[i].y);

              q[i].r=sqrt(d*d-p[i].y*p[i].y)+p[i].x;

         }           //求解每个区间。

         sort(q,q+n,cmp);   //每个区间以左端点大小，从小到大排序。

         a=q[0].l;b=q[0].r;

         for(i=1;i<n;i++)

         {

              if(q[i].l>b)

              {

                   ans++;

                   a=q[i].l;

                   b=q[i].r;

              }

              else if(q[i].l<=b)

              {

                   a=q[i].l;

                   if(q[i].r<b)

                   b=q[i].r;

              }

         }

         printf("Case %d: %d\n",cas++,ans);

         

    }

    return 0;

}