BZOJ 1941 [Sdoi2010] Hide and Seek

Description

小猪iPig在PKU刚上完了无聊的猪性代数课，天资聪慧的iPig被这门对他来说无比简单的课弄得非常寂寞，为了消除寂寞感，他决定和他的好朋友giPi（鸡皮）玩一个更加寂寞的游戏---捉迷藏。 但是，他们觉得，玩普通的捉迷藏没什么意思，还是不够寂寞，于是，他们决定玩寂寞无比的螃蟹版捉迷藏，顾名思义，就是说他们在玩游戏的时候只能沿水平或垂直方向走。一番寂寞的剪刀石头布后，他们决定iPig去捉giPi。由于他们都很熟悉PKU的地形了，所以giPi只会躲在PKU内n个隐秘地点，显然iPig也只会在那n个地点内找giPi。游戏一开始，他们选定一个地点，iPig保持不动，然后giPi用30秒的时间逃离现场（显然，giPi不会呆在原地）。然后iPig会随机地去找giPi，直到找到为止。由于iPig很懒，所以他到总是走最短的路径，而且，他选择起始点不是随便选的，他想找一个地点，使得该地点到最远的地点和最近的地点的距离差最小。iPig现在想知道这个距离差最小是多少。 由于iPig现在手上没有电脑，所以不能编程解决这个如此简单的问题，所以他马上打了个电话，要求你帮他解决这个问题。iPig告诉了你PKU的n个隐秘地点的坐标，请你编程求出iPig的问题。

Input

第一行输入一个整数N 第2~N+1行，每行两个整数X，Y，表示第i个地点的坐标

Output

一个整数，为距离差的最小值。

Sample Input

4
0 0
1 0
0 1
1 1

Sample Output

1

HINT

对于30%的数据,N<=1000 对于100%的数据，N<=500000,0<=X,Y<=10^8 保证数据没有重点保证N>=2

Source

SDOI2010 第二轮Day 1

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

KDtree~

对于每个点，求出它的最大距离-最小距离，再更新答案。

注意update里面要先把c[k].mn和c[k].mx赋值，因为这个调了好久呢。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

const int inf=999999999;

int n,ans,x,y,kkz,maxx,minn,root;

struct node{
	int a[2],mn[2],mx[2],l,r;
}a[500001],c[500001],tmp;

int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0' && ch<='9') {x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

bool operator < (node u,node v)
{
	return u.a[kkz]<v.a[kkz];
}

int dis(node u,node v)
{
	return abs(u.a[0]-v.a[0])+abs(u.a[1]-v.a[1]);
}

void update(int u)
{
	node l=c[c[u].l],r=c[c[u].r];
	for(int k=0;k<2;k++)
	{
		c[u].mx[k]=c[u].mn[k]=c[u].a[k];
		if(c[u].l) c[u].mn[k]=min(c[u].mn[k],l.mn[k]),c[u].mx[k]=max(c[u].mx[k],l.mx[k]);
		if(c[u].r) c[u].mn[k]=min(c[u].mn[k],r.mn[k]),c[u].mx[k]=max(c[u].mx[k],r.mx[k]);
	}
}

int build(int l,int r,bool num)
{
	kkz=num;
	int mid=l+r>>1;
	nth_element(a+l,a+mid,a+r+1);
	c[mid]=a[mid];
	if(l<mid) c[mid].l=build(l,mid-1,num^1);
	if(r>mid) c[mid].r=build(mid+1,r,num^1);
	update(mid);
	return mid;
}

int cal(int u)
{
	int now=0;
	for(int k=0;k<2;k++) now+=max(0,c[u].mn[k]-tmp.a[k]);
	for(int k=0;k<2;k++) now+=max(0,tmp.a[k]-c[u].mx[k]);
	return now;
}

void querymin(int u)
{
	int dl=inf,dr=inf;
	if(c[u].a[0]!=tmp.a[0] || c[u].a[1]!=tmp.a[1]) minn=min(minn,dis(c[u],tmp));
	if(c[u].l) dl=cal(c[u].l);
	if(c[u].r) dr=cal(c[u].r);
	if(dl<dr)
	{
		if(dl<minn) querymin(c[u].l);
		if(dr<minn) querymin(c[u].r);
	}
	else
	{
		if(dr<minn) querymin(c[u].r);
		if(dl<minn) querymin(c[u].l);
	}
}

int cal1(int u)
{
	int now=0;
	for(int k=0;k<2;k++) now+=max(abs(c[u].mx[k]-tmp.a[k]),abs(c[u].mn[k]-tmp.a[k]));
	return now;
}

void querymax(int u)
{
	int dl=0,dr=0;
	maxx=max(maxx,dis(c[u],tmp));
	if(c[u].l) dl=cal1(c[u].l);
	if(c[u].r) dr=cal1(c[u].r);
	if(dl>dr)
	{
		if(dl>maxx) querymax(c[u].l);
		if(dr>maxx) querymax(c[u].r);
	}
	else
	{
		if(dr>maxx) querymax(c[u].r);
		if(dl>maxx) querymax(c[u].l);
	}
}

int calmin()
{
	minn=inf;querymin(root);
	return minn;
}

int calmax()
{
	maxx=0;querymax(root);
	return maxx;
}

int main()
{
	n=read();ans=inf;
	for(int i=1;i<=n;i++) a[i].a[0]=read(),a[i].a[1]=read();
	root=build(1,n,0);
	for(int i=1;i<=n;i++) tmp=a[i],ans=min(ans,calmax()-calmin());
	printf("%d\n",ans);
}