【算法学习】【贪心算法】poj 1328:Radar Installation

问题描述

简化一下，给出x轴上方的若干点，再在x轴上布下雷达，要同时包含住这些点，求雷达最小数目。

样例输入输出

首先明确问题的输入输出：
问题输入：x轴上方岛屿的数目，雷达的半径，岛屿的xy坐标
问题输出：雷达的最小数量

问题分析

乍一看会以雷达为中心，分析岛屿在不在雷达的范围里面。这样倒是将问题想复杂了，换个思路，以岛屿为中心，看雷达可以在哪些范围里面，毕竟题目要求的是雷达的数量。
以岛屿为原点，给出的雷达半径为半径画圆，有如下两种情况：

1. 与x轴没有交点，即d<y，此时的情况按照题意为-1
2. 与x轴有交点，不妨设交点的横坐标为left,right（就算只有一个交点也是一样的，无非left=right）

显而易见，雷达可以安放在left与right之间。那么问题从二维平面转换到了一维直线上，每个岛屿都对应着自己的left和right。这就类似于活动选择问题，每个活动有自己的开始时间和结束时间，在这个问题里最早结束的活动必然要选上。同理，不妨将这n个岛屿的right从左到右排序（即是从小到大排序），后续思路就跟活动选择问题差不多了。

贴代码：

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
int ans[100];
struct node
{
	int x, y;//岛屿的x,y坐标
	double left, right;//以该点为圆心，r为半径时，圆与x轴相交的两点left和right
}po[1000];
bool cmp(node a, node b)
{
	return a.right < b.right;
}
int radar(int n,int r)
{
	for (int i = 0; i < n; i++)
	{
		//如果雷达扫描半径较小不能覆盖岛屿
		if (r < po[i].y)
			return -1;
		else
		{
			po[i].left = po[i].x - sqrt(double(r * r - po[i].y *po[i].y));
			po[i].right=po[i].x+ sqrt(double(r * r - po[i].y * po[i].y));
		}
	}
	sort(po, po + n, cmp);
	//按照right从小到大排序，先选right排第一的圆
	//按照left,right在x轴的顺序，有些类似活动选择问题
	//将left看成一件活动的开始时间，right看成一件活动的结束时间
	//最早结束的活动必选
	int num = 1;//雷达数目
	double bound = po[0].right;
	for (int i = 1; i < n; i++)
	{
		//如果下一个圆的左端不在边界内，说明没有相交的部分
		//雷达数加1，同时边界向右扩
		//如果下一个圆的边界在边界内，说明二者可以共用一个雷达
		if (po[i].left > bound)
		{
			num++;
			bound = po[i].right;
		}
	}
	return num;
}
int main()
{
	int n, d;//n为岛屿的数量，d为雷达的半径
	cin >> n >> d;
	int testCase = 0;
	while (n&&d)
		//只要n,d都不为0
	{
		for (int i = 0; i < n; i++)
			cin >> po[i].x >> po[i].y;
		ans[testCase++] = radar(n, d);
		cin >> n >> d;
	}
	for (int i = 0; i < testCase; i++)
	{
		cout << "Case " <<i+1 << ": " << ans[i] << endl;
	}
	return 0;
}

ERROR汇总

1. 一开始compile error，检查发现是sqrt函数内要使用double或者float类型数。。。好吧我在vs上一直都不注意这些，来个强转好了
2. 后来presentation error，服了。真就不愧是poj的输入输出，每次都能在这种地方上报一堆错。注意到了case后面的空格，没看到冒号后面的空格。。。