poj1328 Radar Installation

第一篇博客~

题意：有一条海岸线，还上有n个岛屿，雷达的探测范围是以d为半径的圆，当岛屿与雷达之间的距离小于等于d时则岛屿可以被探测到，在海岸线上建立尽可能少的雷达，并且每一个岛屿都必须能够被探测到。如果某个岛屿永远无法被探测到的情况输出-1。

贪心算法，先把每个岛屿对应圆心在x轴上的范围也就是求出每个岛屿能够被探测到的雷达建设范围，记录范围的左端与右端（要用double型），之后对每个岛屿的雷达建设范围的左端进行从小到大排序（排序时注意浮点型qsort排序），排序完成后就可以看成找尽可能少的点使所有的区段都能包含至少一个这样的点。

AC代码：

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<math.h>
#include<algorithm>
using namespace std;
struct node
{
    double left;
    double right;
};
struct node coor[1005];
int cmp(const void *a,const void *b)
{
    return (*(node*)a).left>=(*(node*)b).left?1:-1;/*这里注意浮点型的排序返回不能用a-b形式的*/
}
int main()
{
    int n,d,i,x,y,t=0,num,fail;
    double j,std;
    while(scanf("%d%d",&n,&d)!=EOF&&(n||d))
    {
        fail=0;
        t++;
        num=1;
        for(i=0;i<n;i++)
        {
            scanf("%d%d",&x,&y);
            if(y>d)/*当某一个岛屿的纵坐标超出探测范围，则这个岛屿无论如何都无法被探测到*/
            fail=1;
            else
            {
                /*计算雷达建设范围*/
                j=sqrt((double)(d*d-y*y));
                coor[i].left=x-j;
                coor[i].right=x+j;
            }
        }
        if(fail)
        {
            printf("Case %d: -1\n",t);
        }
        else
        {
          qsort(coor,n,sizeof(node),cmp);
          std=coor[0].right;/*以最左端的一个雷达建设范围的右端为参考值*/
          for(i=1;i<n;i++)
          {
              if(coor[i].left>std)/*当某一个雷达建设范围的左端超出参考值，说明这个范围对应的
              岛屿与前面的岛屿不能公用一个雷达，num++，并且更新参考值为该雷达建设范围的右端*/
              {
                  num++;
                  std=coor[i].right;
              }
              else if(coor[i].right<std)/*当某一雷达建设范围的右端小于参考值时说明这个雷达建设
              范围是套在参考值对应的雷达建设范围内的，为保证此雷达建设范围对应的岛屿可被探测，
              参考值应更新为此雷达建设范围的右端*/
              {
                  std=coor[i].right;
              }
          }
         printf("Case %d: %d\n", t, num);
        }
    }
    return 0;
}