贪心 POJ 1328 Radar Installation

 /*
贪心
(转载)题意：有一条海岸线，在海岸线上方是大海，海中有一些岛屿，
这些岛的位置已知，海岸线上有雷达，雷达的覆盖半径知道，
问最少需要多少个雷达覆盖所有的岛屿。
(错误)思路：我开始是想从最左边的点雷达能探测的到的最右的位置出发，判断右边其余的点是否与该点距离小于d
是，岛屿数-1；不是，雷达数+1，继续。。。
(正确)思路：每个岛屿的座标已知，以雷达半径为半径画圆，与x轴有两个交点。
也就是说，若要覆盖该岛，雷达的位置范围是这两个交点。因此转化为覆盖区间的问题。
参考代码：http://www.cnblogs.com/kuangbin/archive/2011/07/30/2121838.html
*/
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>
using namespace std;

const int MAXN = 1e3 + 10;
const int INF = 0x3f3f3f3f;
struct NODE
{
int x, y;
double l, r;
}node[MAXN];

bool cmp(NODE a, NODE b)
{
return a.l < b.l;
}

bool ok(int n, int d)
{
if (d < 0)    return false;
for (int i=1; i<=n; ++i)
{
if (node[i].y > d)    return false;
}

return true;
}

void work(int n, int d)
{
int cnt = 1;
double now = node[1].r;
for (int i=2; i<=n; ++i)
{
if (now > node[i].r)    now = node[i].r;
if (now < node[i].l)
{
now = node[i].r;
cnt++;
}
}
printf ("%d\n", cnt);
}

int main(void)        //POJ 1328 Radar Installation
{
//freopen ("I.in", "r", stdin);

int n, d;
int cnt = 0;
while (~scanf ("%d%d", &n, &d) && n && d)
{
for (int i=1; i<=n; ++i)
{
scanf ("%d%d", &node[i].x, &node[i].y);
node[i].l = (double)node[i].x - sqrt ((double)d * d - node[i].y * node[i].y);
node[i].r = (double)node[i].x + sqrt ((double)d * d - node[i].y * node[i].y);
}
sort (node+1, node+1+n, cmp);
printf ("Case %d: ", ++cnt);
if (!ok (n, d))
{
printf ("%d\n", -1);    continue;
}
work (n, d);
}

return 0;
}

/*
void work(int n, int d)
{
int i = 1;
int j = 1;
int next = 1;
int num = 0;
int cnt = 0;
while (num < n)
{
double res = node[i].x + sqrt (d * d - node[i].y * node[i].y);
cnt++;    num++;
int flag = 0;
fors (j=i+1; j<=n; ++j)
{
if (pow (node[j].x - res, 2) + pow (node[j].y, 2) <= d * d)
{
num++;    next = j;    flag = 1;
}
}
if (flag)    i = next + 1;
else    i++;
}
printf ("%d\n", cnt);
}
*/