POJ1328 Radar Installation(贪心)

Radar Installation

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 77816		Accepted: 17409

Description

Assume the coasting is an infinite straight line. Land is in one side of coasting, sea in the other. Each small island is a point locating in the sea side. And any radar installation, locating on the coasting, can only cover d distance, so an island in the sea can be covered by a radius installation, if the distance between them is at most d. 

We use Cartesian coordinate system, defining the coasting is the x-axis. The sea side is above x-axis, and the land side below. Given the position of each island in the sea, and given the distance of the coverage of the radar installation, your task is to write a program to find the minimal number of radar installations to cover all the islands. Note that the position of an island is represented by its x-y coordinates. 
 
Figure A Sample Input of Radar Installations

Input

The input consists of several test cases. The first line of each case contains two integers n (1<=n<=1000) and d, where n is the number of islands in the sea and d is the distance of coverage of the radar installation. This is followed by n lines each containing two integers representing the coordinate of the position of each island. Then a blank line follows to separate the cases. 

The input is terminated by a line containing pair of zeros 

Output

For each test case output one line consisting of the test case number followed by the minimal number of radar installations needed. "-1" installation means no solution for that case.

Sample Input

3 2
1 2
-3 1
2 1

1 2
0 2

0 0

Sample Output

Case 1: 2
Case 2: 1

好题！真是一道好题啊！这一道题的意思就是说，在x轴上方是海域，下方和x轴都属于陆地。现在给定一些海岛坐标，并且给出雷达可以搜索的长度，要求雷达尽可能覆盖海岛，并且所用雷达数目最少。这道题，一开始想真的觉得好困惑，这咋整？后来慢慢形成思路，我把他们全部按照x轴从小到大排序，从最左边开始贪心，在雷达包含选定的海岛的情况下，尽可能放入更多的海岛。好，这是一个初始方案，接下来写起来就没这么简单咯。对了，雷达肯定是在x轴上的，这样照到的海域肯定会比在x轴下面要大！

问题1，选定第一个海岛之后，他在x轴上有两个交点，然后如果放入第二个海岛，这个交点的范围会一直缩小的，直到其他全部都不能放入这个交点范围内。而我一开始却只保留了右边的交点，因为我是以x轴从小到大排序的嘛，天真的以为只要下一个海岛左交点小于这个右交点就行了，但是发现我左边的交点也是变化的，也需要算出来才行。

问题2，排序的时候还要注意，如果x轴是相同的，那就按照y的大小从大到小排序，这样最上面那个点和x轴的两个交点的范围是相对于下面的点来说是比较窄的，可以减少一下范围的更新次数嘛。。。

问题3，我也是看到别人给的数据才想的，那就是如果真的这么奇葩他的海岛给的y点是小于0的，也就是在陆地上的怎么办，然后就是我一开始是用的m+n来结束while读入循环，然后发现他给我出了这样一组数据，3 -3，三个岛，雷达范围是-3？什么鬼？反正这样也会让我的程序结束，这样的数据我觉得是不合法的，但是他给出来了就改一改吧，毕竟acm本来就是个坑。

然后我自己还有各种各样的粗心大意的错误，TLE了两发，while循环那忘记要j++了。。。RE的两发，这个。。。我也不知道为什么。。。wa了一发，还是太粗心。。。然后，终于终于过了！

代码如下：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
int mark[1005];
int n, m, ppp;
struct node
{
	int x, y;
	double left,right; 
	bool operator < (const node &a)const
	{
		if(x != a.x)
			return x < a.x;
		else
			return y > a.y;
	}
}island[1005]; 

void func(int i)
{
	ppp++;
	mark[i] = 1;
	double LEFT = island[i].left, RIGHT = island[i].right;
	int j = i + 1;
	while(j < n)
	{	
		if(!mark[j] && ((island[j].left >= LEFT && island[j].left <= RIGHT)					 
					|| (island[j].right >= LEFT && island[j].right <= RIGHT)
					|| (island[j].left <= LEFT && island[j].right >= RIGHT)))
		{
			mark[j] = 1;
			LEFT = island[j].left > LEFT ? island[j].left : LEFT;
			RIGHT = island[j].right < RIGHT ? island[j].right : RIGHT;
		} 
		j++;
	} 
}

int main()
{ 
	int flag, k = 1, i, j; 
	int tx, ty;
	while(scanf("%d%d", &n, &m), n != 0 || m != 0)
	{
		memset(mark, 0, sizeof(mark));
		flag = 0;
		ppp = 0;
		for(i = 0, j = 0; i < n; i++)
		{
			scanf("%d%d", &tx, &ty);
			if(ty > m)
				flag = 1;
			if(ty >= 0 && !flag)
			{
				island[j].x = tx;
				island[j].y = ty;
				island[j].left = 1.0 * tx - sqrt(m * m - ty * ty);
				island[j].right = 1.0 * tx + sqrt(m * m - ty * ty);
				j++;
			}		
		}
		n = j;
		printf("Case %d: ", k++);
		if(flag || m <= 0)
			printf("-1\n"); 
		else
		{
			sort(island, island + n);
			for(i = 0; i < n; i++)
			{
				if(!mark[i])
					func(i); 
			}
			printf("%d\n", ppp);
		}
	}
	return 0;
}

看到代码才想起，就是我们判定下一个海岛是否在这个范围内，有三种情况，一种是，下一个海岛的左交点坐标在范围内，一种是右交点在范围内，还有一种！！！就是下一个海岛的左交点小于范围的左交点，右交点大于范围的右交点！相当于下一个海岛的左右范围比原本的范围还要大！这也是有交集的！！！脑子不够用我当时就是少了这种情况啊！！！好像就是因为这个wa了一发。

一年多了，重写了一次，发现还是不会。。。不过代码风格变得更清晰了

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<cmath>
#include<utility> 
#include<vector>
using namespace std;
#define eps 1e-7
typedef pair<double, double> pii;
vector <pii> v;
int n, d;

int main() {
	int Case = 1, x, y;
	while(cin >> n >> d) {
		
		bool flag = 1;
		v.clear();
		if(n == 0 && d == 0)
			break;
		if(d < 0)
			flag = 0;
		for(int i = 0; i < n; i++) {
			cin >> x >> y;
			
			if(y <= 0)
				continue;
			if(y > d)
				flag = 0;
				
			if(!flag)
				continue;
				
			double tmp = sqrt((double)d * d - (double)y * y);
			v.push_back(pii(x - tmp, x + tmp));
		}
		
		if(flag) {
			int ans = 0;
			sort(v.begin(), v.end());
			for(int i = 0; i < v.size(); i++) {
				double l = v[i].first;
				double r = v[i].second;
				int k = i + 1;
				while(k < v.size() && v[k].second >= l && v[k].first <= r) {
					l = max(v[k].first, l);
					r = min(v[k].second, r);
					k++; 
				} 
				i = k - 1;
				ans++;
			}
			printf("Case %d: %d\n", Case++, ans);
		} else {
			printf("Case %d: -1\n", Case++);
		}
	}
	return 0;
}