【ybt高效进阶1-2-2】【luogu P1325】雷达装置

雷达装置

题目链接：ybt高效进阶1-2-2 / luogu P1325

题目大意

就是问你要至少以 x 轴某个点为圆心画多少个半径是 d 的圆，才能使所有给出的点都至少在其中一个圆内。

思路

我们先考虑对于一个点，如果要让它包含在一个园内，这个圆的圆心的范围。

通过勾股定理可以得出，对于坐标 $(x,y)$ 的点，圆心的 y 坐标范围是 $(x-\sqrt{d^2-(|y|{)}^2})\sim (x+\sqrt{d^2-(|y|{)}^2})$。
同时，也可以得出，如果 $|y|>d$，则没有解。

那我们其实可以贪心，先把坐标按圆心 y 坐标范围最大值从小到大排序。
（这里的贪心不能按照点的 y 坐标，因为不一定 y 坐标越大，圆的圆心最右位置就一定更大）
如果一定要用放新的圆，圆心就放在最后边。
那这样就让这个圆可以尽可能包含更多的点。

怎样就一定要放新的圆的呢？
就是最后的圆包含不了这个点，或者说圆心不在范围内。

代码

#include<cmath>
#include<cstdio>
#include<algorithm>

using namespace std;

struct zb {
	double x, y, l, r;
}a[1001];
int n, ans;
double last, d;

double abss(double now) {
	if (now < 0.0) return -now;
	return now;
}

bool cmp(zb x, zb y) {
	return x.r < y.r;
}

int main() {
	scanf("%d %lf", &n, &d);
	for (int i = 1; i <= n; i++) {
		scanf("%lf %lf", &a[i].x, &a[i].y);
		a[i].y = abss(a[i].y);
		if (a[i].y > d) {
			printf("-1");
			return 0;
		}
		a[i].l = a[i].x - sqrt(d * d - a[i].y * a[i].y);
		a[i].r = a[i].x + sqrt(d * d - a[i].y * a[i].y);
	}
	
	sort(a + 1, a + n + 1, cmp);
	
	last = a[1].r;
	ans++;
	for (int i = 2; i <= n; i++) {
		if (a[i].l > last) {
			last = a[i].r;
			ans++;
		}
	}
	
	printf("%d", ans);
	
	return 0;
}