P2498 [SDOI2012] 拯救小云公主并查集二分

BestMonkey

已于 2024-05-16 15:45:12 修改

阅读量355

点赞数 3

分类专栏：题解文章标签：算法数据结构 c++ c语言

于 2024-05-16 15:44:26 首次发布

本文链接：https://blog.csdn.net/BestMonkey/article/details/138963613

版权

题解专栏收录该内容

64 篇文章 0 订阅

订阅专栏

[SDOI2012] 拯救小云公主

传送门

题目描述

英雄又即将踏上拯救公主的道路……

这次的拯救目标是——爱和正义的小云公主。

英雄来到boss的洞穴门口，他一下子就懵了，因为面前不只是一只boss，而是上千只boss。当英雄意识到自己还是等级1的时候，他明白这就是一个不可能完成的任务。

但他不死心，他在想，能不能避开boss去拯救公主呢，嘻嘻。

Boss的洞穴可以看成一个矩形，英雄在左下角（1,1），公主在右上角（row，line）。英雄为了避开boss，当然是离boss距离越远越好了，所以英雄决定找一条路径使到距离boss的最短距离最远。

Ps:英雄走的方向是任意的，但是不能走出矩形的范围。即英雄可以到达矩形范围内的任意一个点（没有必要是整点）

你可以帮帮他吗？

当英雄找到了美丽漂亮的小云公主，立刻就被boss包围了！！！英雄缓闭双眼，举手轻挥，白光一闪后使用了回城卷轴，回到了城堡，但只有小云公主回去了……因为英雄忘了进入回城的法阵了。

输入格式

第一行，输入三个整数，n表示boss的数目，row，line表示矩形的大小；

接下来n行，每行分别两个整数表示boss的位置坐标。

输出格式

输出一个小数，表示英雄的路径离boss的最远距离，精确到小数点后两位。

样例 #1

样例输入 #1

1 3 3
2 2

样例输出 #1

1.00

样例 #2

样例输入 #2

1 3 3
3 1

样例输出 #2

2.00

提示

数据范围：

20%数据，boss坐标范围小于等于50；

60%数据，n<=1500；

100%数据，n<=3000；

注明

$以上来自洛谷。$

好久没更新了，是更不动了？

这个 LaTeX 不太行啊？

解题思路

转换一下题意：

将每个 boss 看成是是以 $(x, y)$ 圆心作一个半径为r的圆，我们想要求的就是让这些圆尽可能大,并且不能影响我们从 $(1, 1)$ 到 $(n, m)$ ，即不能覆盖。

考虑边界 $(n, 1)$ 和 $(1, m)$ ，如果这两个点没有被覆盖，那我一定可以到达 $(n, m)$ 。

用并查集维护各圆间是否连通。

最后二分求 boss 最大的合法半径 r 的大小即可。

AC Code

#include<bits/stdc++.h>
using namespace std;
const int Maxn = 3000 + 5;
const double eps = 1e-8;
struct Project_Boss_Pos {
	double x, y;
} Boss_Pos[Maxn];
int Fa[Maxn];
inline int Get_Fa(int now) {
	return Fa[now] == now ? now : Fa[now] = Get_Fa(Fa[now]);
}
inline void Union(int x, int y) {
	int x_fa = Get_Fa(x), y_fa = Get_Fa(y);
	if (x_fa != y_fa)
		Fa[x_fa] = y_fa;
}
inline double Get_Dis(int i, int j) {
	return (Boss_Pos[i].x - Boss_Pos[j].x) * (Boss_Pos[i].x - Boss_Pos[j].x) + (Boss_Pos[i].y - Boss_Pos[j].y) * (Boss_Pos[i].y - Boss_Pos[j].y);
}
int n, row, line;
inline bool Check(double round_s_r) {
	for (register int i = 0; i <= n + 1; ++i)
		Fa[i] = i;
	for (register int i = 1; i <= n; ++i) {
		for (register int j = 1; j < i; ++j)
			if (Get_Dis(i, j) <= (4 * round_s_r * round_s_r))
				Union(i, j);
		if (Boss_Pos[i].x - round_s_r <= 1 || Boss_Pos[i].y + round_s_r >= line)
			Union(i, 0);
		if (Boss_Pos[i].x + round_s_r >= row || Boss_Pos[i].y - round_s_r <= 1)
			Union(i, n + 1);
	}
	return Get_Fa(0) != Get_Fa(n + 1);
}
double Answer = 0.000000000000000000000000;
signed main() {
	ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
	cin >> n >> row >> line;
	for (register int i = 1; i <= n; ++i)
		cin >> Boss_Pos[i].x >> Boss_Pos[i].y;
	double l = 0.000000000000000000000000, r = min(row, line), mid = 0.000000000000000000000000;
	while (fabs(l - r) > eps) {
		mid = (l + r) / 2;
		if (Check(mid)) {
			l = mid;
			Answer = mid;
		} else
			r = mid;
	}
	printf("%.2f\n", Answer);
	return 0;
}