POJ - 1328 Radar Installation 安装雷达 贪心+区间选点 重庆一中高2018级竞赛班第八次测试 2016.8.7 Problem 3

【问题描述】
在一个坐标系中，x轴代表海岸线，x轴以上是大海。在海上有n个岛屿。现在打算在x轴上建一些雷达，一部雷达的覆盖范围为半径为 d 的圆。请你计算在x轴上最少建多少部雷达，才能把全部的岛屿覆盖起来，如果不能覆盖，输出-1。

【输入格式】
输入包含若干组数据，每组数据的第一行为整数 n 和 d ，它们的意义如题目描述。接下来的 n 行，每行为两个整数 x,y，表示一个岛屿的坐标。输入以0 0结束。

【输出格式】
每组数据输出一行，如果建立的雷达能覆盖所有岛屿，输出最少需要安装的雷达数目，否则输出-1。

【输入样例】

3 2
1 2
-3 1
2 1

1 2
0 2

0 0

【输出样例】

Case 1: 2
Case 2: 1

【数据范围】
1<=n<=1000 -10^9<=x<1=0^9 0< d,y <=10^9

思路：很容易想到贪心来做，每个岛屿在X轴上可安装雷达的区域是[x[i]-sqrt(d*d-y[i]*y[i]),x[i]+sqrt(d*d-y[i]*y[i])]，于是问题变成给定n个区间，选择最少的点，使得每个区间至少有一个点！但是我没有另外用一个数组来存储这些区间，所以可能贪心时排序等出现了一些问题。
注：sqrt最好是实数double互相运算，不然可能会有误差。

代码如下：

/*
    Name: radar.cpp
    Copyright: Twitter & Instagram @stevebieberjr
    Author: @stevebieberjr
    Date: 07/08/16 20:00
*/
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;

int n;
double d;

struct island
{
    double x,y;
}a[10005];

struct qujian
{
    double z,y;
}b[10005];

bool cmp(qujian i,qujian j)
{
    if(i.y==j.y) return i.z<j.z;
    return i.y<j.y;
}

int main()
{
    freopen("radar.in","r",stdin);
    freopen("radar.out","w",stdout);

    int T=0;
    while(1)
    {
        scanf("%d%lf",&n,&d);
        T++;
        double maxy=0;
        if(n==0 && d==0) break;

        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf",&a[i].x,&a[i].y);
            b[i].z=a[i].x-sqrt(d*d-a[i].y*a[i].y);
            b[i].y=a[i].x+sqrt(d*d-a[i].y*a[i].y);
            maxy=max(maxy,a[i].y);
        }

        if(maxy>d)
        {
            printf("Case %d: -1\n",T);
            continue;
        }

        sort(b+1,b+1+n,cmp);

        int cnt=1;
        double last=b[1].y;
        for(int i=2;i<=n;i++)
        {
            if(b[i].z>last)
            {
                last=b[i].y;
                cnt++;
            }
        }
        printf("Case %d: %d\n",T,cnt);
    }
    return 0;
}