5904. 【NOIP2018模拟10.15】刺客信条（AC）

题目大意：

故事发生在1486 年的意大利，Ezio 原本只是一个文艺复兴时期的贵族，后来因为家族成员受到圣殿骑士的杀害，决心成为一名刺客。最终，凭借着他的努力和出众的天赋，成为了杰出的刺客大师。刺客组织在他的带领下，为被剥削的平民声张正义，赶跑了原本统治意大利的圣殿骑士首领-教皇亚历山大六世。在他的一生中，经历了无数次惊心动魄、扣人心弦的探险和刺杀。

这次的故事就是他暗杀一位作恶多端的红衣主教。红衣主教可以吸取他周围人的生命力量，而他的红衣教徒也拥有这个力量。红衣主教的家是一个x*y 的长方形房间，也就是说，他的家的四个角坐标分别为(0,0)(x,0)(0,y)(x,y)。教堂的门在（0，0） ,而红衣主教就在 （x,y）的卧室休息。他的家中还有n个守护着他的红衣教徒，站在(ai,bi)。Ezio想要趁主教休息时，从门进入潜入到他的卧室刺杀他，因为主教休息时会脱下红衣，这样吸取生命的力量就消失了。可是守卫他的红衣教徒依然很危险，离红衣教徒太近就会被吸取生命。因此，Ezio想知道，在能刺杀主教的前提，从门到他的卧室的路上，他最远和离他最近的红衣教徒保持多远的距离。注意：教徒都在房间里。

思路：

这题是以前做过的一道题的升级版，好像是去年的有次比赛，叫奶酪。你只需要二分一个半径，然后判断园的集合$(0,0)$到$(x,y)$之间的路径是否被切断。并查集，搜索都可以实现。(分成四种情况讨论)。

程序：

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#define LL long long
using namespace std;
const LL N=2005;
LL a[N],b[N],x,y,n;
double ans;
bool flag[N],hjy;

double sqr(double x){
	return x*x;
}

bool dis(LL x,LL y,double r){
	return sqr(a[x]-a[y])+sqr(b[x]-b[y])<4*r*r;
}

void dfs(LL x,double r,LL p){
	flag[x]=1;
	if ((a[x]-r<0||b[x]+r>y)&&p==1) {
		hjy=1;
		return;
	}
	if ((b[x]-r<0)&&p==2){
		hjy=1;
		return;
	}
	if (b[x]-r<0&&p==3){
		hjy=1;
		return;
	}
	if (a[x]-r<0&&p==4) {
		hjy=1;
		return;
	}
	for (LL i=1;i<=n;i++) 
	  if (dis(i,x,r)&&!flag[i]&&!hjy)
		dfs(i,r,p);
}

bool cheak(double r){
	hjy=0;
	memset(flag,0,sizeof(flag));
	for (LL i=1;i<=n;i++)
		if (!hjy&&a[i]+r>x&&!flag[i]) dfs(i,r,1); 	
	memset(flag,0,sizeof(flag));
	for (LL i=1;i<=n;i++)
		if (!hjy&&b[i]+r>y&&!flag[i]) dfs(i,r,2);
	memset(flag,0,sizeof(flag));
	for (LL i=1;i<=n;i++)
		if (!hjy&&a[i]-r<0&&!flag[i]) dfs(i,r,3); 	
	memset(flag,0,sizeof(flag));
	for (LL i=1;i<=n;i++)
		if (!hjy&&b[i]-r<0&&!flag[i]) dfs(i,r,4);
	if (hjy) return 1;
	return 0;
}

int main(){
	freopen("AC.in","r",stdin);
	freopen("AC.out","w",stdout);
	scanf("%lld%lld%lld",&x,&y,&n);
	for (LL i=1;i<=n;i++)
		scanf("%lld%lld",&a[i],&b[i]);
	double l=0,r=x;
	while (r>=l){
		double mid=(r+l)/2;
		if (cheak(mid)) r=mid-0.001;
				else {
					ans=mid;
					l=mid+0.001;
				}
	}
	printf("%.2lf",ans);
}