poj1328解题报告（贪心、线段交集）

POJ 1328，题目链接http://poj.org/problem?id=1328

题意：

有一海岸线（x轴），一半是陆地（y<0）、一半是海（y>0），海上有一些小岛（用坐标点表示P1、P2...），现要在海岸线上建雷达（覆盖半径R）。给出所有小岛的位置，和雷达半径，求最少需要多少个雷达？

思路：

1. 知道小岛位置，和雷达半径，那么以小岛为圆心，雷达覆盖半径为半径画圆，可以求出小岛与x轴有0（雷达无法覆盖）、1（雷达只能在这个点上才能覆盖）、2个交点（雷达在两点之间都能覆盖该小岛）

2. 要求最少雷达多少个，即把雷达放在1中线段的交集内。

那么这就变成了线段交集问题。（贪心）

代码：

//404k 79ms
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>

typedef struct tagLINE{
	double left;
	double right;
}Line;
void sortLineBuf(Line *p, int num)
{
	Line temp;
	for (int i=0; i<num; ++i)
	{
		for (int j=i+1; j<num; ++j)
		{
			if (p[j].left < p[i].left){
				temp = p[i];
				p[i] = p[j];
				p[j] = temp;
			}
		}
	}
}
int main()
{
	int caseNum = 0;
	double tempPoint;
	Line tempLine;
	while (true)
	{
		int islandNum = 0, r = 0;
		scanf("%d%d", &islandNum, &r);
		if (islandNum == 0 && r == 0) break;

		double *p = (double*)malloc(sizeof(double) * islandNum * 2);
		double *pX = p;
		double *pY = p+islandNum;
		for(int i=0; i<islandNum; ++i){
			scanf("%lf%lf", &pX[i], &pY[i]);
		}
		//
		int rapar = 0;
		bool bImpossible = true;
		Line* pLine = (Line*)malloc(sizeof(Line) * islandNum);
		//1 转换为直线
		for(int i=0; i<islandNum; ++i){
			if (fabs(pY[i]) > r){
				bImpossible = false;
				rapar = -1;
				break;
			}
			tempPoint = sqrt(r*r - pY[i]*pY[i]);
			pLine[i].left = pX[i]-tempPoint;
			pLine[i].right = pX[i]+tempPoint;
		}
		if (bImpossible)
		{
			rapar = 1;
			//2
			sortLineBuf(pLine, islandNum);
			//3 求解线段交集
			tempLine = pLine[0];
			for (int i=1; i<islandNum; ++i)
			{
				if (pLine[i].left > tempLine.right)
				{
					++rapar;
					tempLine = pLine[i];
				}
				else if (pLine[i].right < tempLine.right)
				{
					tempLine = pLine[i];
				}
			}
		}
		printf("Case %d: %d\n", ++caseNum, rapar);
		free(p);
		free(pLine);
	}

	return 0;
}