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;
}
         
         
              
posted @ 2011-04-22 10:30  CoderZhuang  阅读(140)  评论(0编辑  收藏  举报