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 l o g n ∗ α ( n ) ) O(n^2log\ n*α(n)) O(n2log n∗α(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;
}