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
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);
}