POJ1328 Radar Installation 贪心

Problem Address：http://poj.org/problem?id=1328

【思路】

贪心的思路。

（1）把所有点按 x 的递增排序。同时计算出每个点可能的圆心的区间，记录其最左值和最右值。

（2）计算每个点到最后一个点的最小的右坐标。即用一个数组保存，表示从这个点到最后一个的点的所有最右值中的最小值。

（3）扫描一遍数组，如果存在某个点的 y 的绝对值大于 d，则输出-1。结束。

（3）从最左端开始，首先取第一个点的最小右值作为当前圆心，圆个数初始化为1。（之所以要用最小右值是因为存在位于该点右边但是需要的圆心位置却比该点的小的情况。）如果某个点的圆心区间包含了这个圆心，则继续扫描下去。若不包含，则取当前点的最小右值作为当前圆心，这个圆心必定是包含当前点且可以作为最右圆心的圆心；同时圆个数加一。之后继续扫描下去。直到扫描完毕，即输出答案。

【代码】

#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;

const int maxn = 1000;

struct point
{
	int x;
	int y;
}pt[maxn+5];

double leftcenter[maxn+5], rightcenter[maxn+5], rightmin[maxn+5];

bool cmp(const point &a, const point &b)
{
	if (a.x!=b.x) return a.x<b.x;
	else return a.y<b.y;
}

#define min(a,b) ((a)<=(b)?(a):(b))

int main()
{
	int n, d;
	double temp;
	int i;
	int ct;
	int index = 0;
	while(scanf("%d %d", &n, &d)!=EOF)
	{
		if (n==0 && d==0) break;
		index++;
		for (i=0; i<n; i++)
		{
			scanf("%d %d", &pt[i].x, &pt[i].y);
		}
		for (i=0; i<n; i++)
		{
			if (pt[i].y>d || pt[i].y<-d) break;//查看是否应该输出-1
		}
		if (i<n)
		{
			printf("Case %d: -1\n", index);
			continue;
		}
		sort(pt, pt+n, cmp);//排序
		rightmin[n] = (double)INT_MAX;//讲标号为n的最小右值设为无限大，方面计算
		for (i=n-1; i>=0; i--)从后往前扫描，多个值同时计算。也可分开扫描。
		{
			temp = sqrt((double)(d*d-pt[i].y*pt[i].y));
			leftcenter[i] = (double)pt[i].x - temp;//计算左值
			rightcenter[i] = (double)pt[i].x + temp;//计算右值
			rightmin[i] = min(rightcenter[i], rightmin[i+1]);//计算最小右值，为当前右值与其右点的最小右值的最小值
		}
		ct = 1;//圆个数初始化为1
		temp = rightmin[0];//当前圆心初始化为第一个点的最小右值
		for (i=0; i<n; i++)
		{
			if (temp<=rightcenter[i] && temp>=leftcenter[i])//如果区间包含当前圆心，则继续扫描
				continue;
			else
			{
				temp = rightmin[i];//否则取当前点的最小右值作为当前圆心
				ct++;//并使圆个数加一
			}
		}
		printf("Case %d: %d\n", index, ct);
	}
	return 0;
}