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

Description

故事发生在1486 年的意大利，Ezio 原本只是一个文艺复兴时期的贵族，后来因为家族成员受到圣殿骑士的杀害，决心成为一名刺客。最终，凭借着他的努力和出众的天赋，成为了杰出的刺客大师。刺客组织在他的带领下，为被剥削的平民声张正义，赶跑了原本统治意大利的圣殿骑士首领-教皇亚历山大六世。在他的一生中，经历了无数次惊心动魄、扣人心弦的探险和刺杀。
这次的故事就是他暗杀一位作恶多端的红衣主教。红衣主教可以吸取他周围人的生命力量，而他的红衣教徒也拥有这个力量。红衣主教的家是一个x*y 的长方形房间，也就是说，他的家的四个角坐标分别为(0,0)(x,0)(0,y)(x,y)。教堂的门在（0，0） ,而红衣主教就在 （x,y）的卧室休息。他的家中还有n个守护着他的红衣教徒，站在(ai,bi)。Ezio想要趁主教休息时，从门进入潜入到他的卧室刺杀他，因为主教休息时会脱下红衣，这样吸取生命的力量就消失了。可是守卫他的红衣教徒依然很危险，离红衣教徒太近就会被吸取生命。因此，Ezio想知道，在能刺杀主教的前提，从门到他的卧室的路上，他最远和离他最近的红衣教徒保持多远的距离。注意：教徒都在房间里。

Input

第一行三个整数x,y,n。之后n行，每行两个整数ai，bi ,意义见题目描述。

Output

一行一个数D，表示Ezio能保持的最大距离，保留两位小数。

Sample Input

10 20 2
3 3
6 14

Sample Output

3.00

样例说明
贴着墙走

Data Constraint

数据范围
对 10%的数据n<=10,
对 30%的数据n<=100
对 100%的数据n<=2000
保证输入合法，x,y属于[1,10^6].

Solution

• 先二分答案，转化成判断性问题。

• 把能连通的点用并查集并起来（下右边界看成一点，上左边界看成一点），

• 若最后边界并在一起了则二分上界下调，否则下界上调。

• 注意一下精度问题即可。

• 理论复杂度 $O(n^{2}logn\ast \alpha (n))$ ，实际上中途退出情况很多，跑得很快。

Code

#include<cstdio>
#include<algorithm>
#include<cctype>
using namespace std;
typedef long long LL;
const int N=2005;
struct data
{
	int x,y;
}a[N];
int n,m,k;
double ans;
int f[N],size[N];
inline int read()
{
	int X=0,w=0; char ch=0;
	while(!isdigit(ch)) w|=ch=='-',ch=getchar();
	while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
	return w?-X:X;
}
inline bool cmp(data x,data y)
{
	return x.x<y.x || x.x==y.x && x.y<y.y;
}
inline int max(int x,int y)
{
	return x>y?x:y;
}
inline int get(int x)
{
	return f[x]==x?x:f[x]=get(f[x]);
}
inline void merge(int x,int y)
{
	int xx=get(x),yy=get(y);
	if(xx^yy)
	{
		if(size[xx]<size[yy]) swap(xx,yy);
		f[yy]=xx;
		size[xx]+=size[yy];
	}
}
inline bool check(double x)
{
	for(int i=1;i<=k+4;i++) size[f[i]=i]=1;
	merge(k+1,k+2);
	merge(k+3,k+4);
	for(int i=1;i<=k;i++)
	{
		if(a[i].y<x) merge(i,k+1);
		if(m-a[i].y<x) merge(i,k+3);
		if(a[i].x<x) merge(i,k+4);
		if(n-a[i].x<x) merge(i,k+2);
	}
	double dis=x*x*4;
	for(int i=1;i<=k;i++)
		for(int j=i+1;j<=k;j++)
		{
			LL num=((LL)a[j].x-a[i].x)*(a[j].x-a[i].x);
			if(num>=dis) break;
			num+=((LL)a[j].y-a[i].y)*(a[j].y-a[i].y);
			if(num<dis) merge(i,j);
			if(get(k+2)==get(k+3)) return false;
		}
	return true;
}
int main()
{
	freopen("AC.in","r",stdin);
	freopen("AC.out","w",stdout);
	n=read(),m=read(),k=read();
	for(int i=1;i<=k;i++) a[i].x=read(),a[i].y=read();
	sort(a+1,a+1+k,cmp);
	if(a[1].x+a[1].y==0) return 0&puts("0.00");
	int l=1,r=max(n,m)*300;
	while(l<=r)
	{
		int mid=l+r>>1;
		if(check(mid*0.003))
		{
			l=mid+1;
			ans=mid*0.003;
		}else r=mid-1;
	}
	printf("%.2lf",ans);
	return 0;
}