循环每一个点,求出据他最近的点的距离和最远的点的距离。用KD-tree来求。
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
inline 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;
}
int n,D,ans=inf,rt=0,mn,mx;
struct point{
int d[2];
int& operator[](int x){return d[x];}
friend bool operator<(point a,point b){return a[D]<b[D];}
friend int dis(point a,point b){return abs(a[0]-b[0])+abs(a[1]-b[1]);}
}P[N];
struct node{
point x;int lc,rc,mn[2],mx[2];
}tr[N];
inline void update(int p){
int l=tr[p].lc,r=tr[p].rc;
for(int i=0;i<2;++i){
if(l){
tr[p].mn[i]=min(tr[p].mn[i],tr[l].mn[i]);
tr[p].mx[i]=max(tr[p].mx[i],tr[l].mx[i]);
}if(r){
tr[p].mn[i]=min(tr[p].mn[i],tr[r].mn[i]);
tr[p].mx[i]=max(tr[p].mx[i],tr[r].mx[i]);
}
}
}
inline void build(int &p,int l,int r,int op){
int mid=l+r>>1;p=mid;D=op;
nth_element(P+l,P+mid,P+r+1);tr[p].x=P[mid];
for(int i=0;i<2;++i) tr[p].mx[i]=tr[p].mn[i]=tr[p].x[i];
if(l<mid) build(tr[p].lc,l,mid-1,op^1);
if(r>mid) build(tr[p].rc,mid+1,r,op^1);
update(p);
}
inline int calcmin(int p,point x){
if(!p) return inf;int res=0;
for(int i=0;i<2;++i){
if(x[i]<tr[p].mn[i]) res+=tr[p].mn[i]-x[i];
if(x[i]>tr[p].mx[i]) res+=x[i]-tr[p].mx[i];
}return res;
}
inline int calcmax(int p,point x){
if(!p) return 0;int res=0;
for(int i=0;i<2;++i)
res+=max(abs(tr[p].mn[i]-x[i]),abs(tr[p].mx[i]-x[i]));return res;
}
inline void qmin(int p,point x){
int tmp=dis(tr[p].x,x);
if(tmp) mn=min(mn,tmp);
int dl=calcmin(tr[p].lc,x),dr=calcmin(tr[p].rc,x);
if(dl<dr){
if(dl<mn) qmin(tr[p].lc,x);
if(dr<mn) qmin(tr[p].rc,x);
}else{
if(dr<mn) qmin(tr[p].rc,x);
if(dl<mn) qmin(tr[p].lc,x);
}
}
inline void qmax(int p,point x){
int tmp=dis(tr[p].x,x);
if(tmp) mx=max(mx,tmp);
int dl=calcmax(tr[p].lc,x),dr=calcmax(tr[p].rc,x);
if(dl>dr){
if(dl>mx) qmax(tr[p].lc,x);
if(dr>mx) qmax(tr[p].rc,x);
}else{
if(dr>mx) qmax(tr[p].rc,x);
if(dl>mx) qmax(tr[p].lc,x);
}
}
int main(){
// freopen("a.in","r",stdin);
n=read();
for(int i=1;i<=n;++i)
for(int k=0;k<2;++k) P[i][k]=read();
build(rt,1,n,0);
for(int i=1;i<=n;++i){
mn=inf;qmin(rt,P[i]);mx=0;qmax(rt,P[i]);
ans=min(ans,mx-mn);
}printf("%d\n",ans);
return 0;
}