题意:平面上有一些点,之后还会在平面上插入一些点,还会询问某一个点到平面中最近的点的距离
K−Dtree K − D t r e e ,最近的点像 [SDOI2010]Hide [ S D O I 2010 ] H i d e and a n d Seek S e e k 这样求就好了
然后这题是带插入的,和平衡树一样,可能会被卡掉
然而 K−Dtree K − D t r e e 又不能旋转,所以只能像替罪羊树那样重构了
然后不会 K−Dtree K − D t r e e 的看这里,再看这里,这个太长了,感觉将就着看吧
#include<bits/stdc++.h>
#define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
#define go(u) for(register int i=fi[u],v=e[i].to;i;v=e[i=e[i].nx].to)
#define file(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
char ss[1<<17],*A=ss,*B=ss;
inline char gc(){return A==B&&(B=(A=ss)+fread(ss,1,1<<17,stdin),A==B)?-1:*A++;}
template<class T>inline void sd(T&x){
char c;T y=1;while(c=gc(),(c<48||57<c)&&c!=-1)if(c==45)y=-1;x=c-48;
while(c=gc(),47<c&&c<58)x=x*10+c-48;x*=y;
}
char sr[1<<21],z[20];int C=-1,Z;
inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
template<class T>inline void we(T x){
if(C>1<<20)Ot();if(x<0)sr[++C]=45,x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++C]=z[Z],--Z);sr[++C]='\n';
}
const int N=5e5+5,M=2*N,inf=-1u>>1;
const double alpha=.77;
typedef int arr[N];
int n,m,rt,D,ans,cb,cn,tot,bin[M];
struct po{
int d[2],mi[2],mx[2],l,r,sz;
inline int&operator[](int x){return d[x];}
inline void clr(){fp(i,0,1)mi[i]=mx[i]=d[i];l=r=0,sz=1;}
inline bool operator<(po b)const{return d[D]<b[D];}
}a[N],tr[M],t;
inline int dis(po a,po b){int s=0;fp(i,0,1)s+=abs(a[i]-b[i]);return s;}
struct KD_Tree{
#define lc tr[p].l
#define rc tr[p].r
inline void up(int p){
tr[p].sz=tr[lc].sz+tr[rc].sz+1;
fp(i,0,1)
tr[p].mi[i]=min(tr[p][i],min(tr[lc].mi[i],tr[rc].mi[i])),
tr[p].mx[i]=max(tr[p][i],max(tr[lc].mx[i],tr[rc].mx[i]));
}
inline int New(){return cb?bin[cb--]:++cn;}
int build(int L,int R,int k){
int mid=(L+R)>>1,p=New();D=k;
nth_element(a+L,a+mid,a+R+1);tr[p]=a[mid];
lc=L<mid?build(L,mid-1,k^1):0;
rc=mid<R?build(mid+1,R,k^1):0;
return up(p),p;
}
void cle(int&p){if(!p)return;bin[++cb]=p;cle(lc),cle(rc),a[++tot]=tr[p],a[tot].clr(),p=0;}
void chk(int&p,int k){
if(!p)return;
if(tr[lc].sz>tr[p].sz*alpha||tr[rc].sz>tr[p].sz*alpha)
return tot=0,cle(p),p=build(1,tot,k),void();
if(t[k]<tr[p][k])chk(lc,k^1);
else chk(rc,k^1);
}
void ins(int&p,int k){
if(!p)return tr[p=New()]=t,void();
if(t[k]<tr[p][k])ins(lc,k^1);
else ins(rc,k^1);up(p);
}
inline int h(int p){return!p?inf:(max(tr[p].mi[0]-t[0],0)+max(t[0]-tr[p].mx[0],0)+max(tr[p].mi[1]-t[1],0)+max(t[1]-tr[p].mx[1],0));}
void qmin(int p){
int dl=h(lc),dr=h(rc),d0=dis(tr[p],t);
cmin(ans,d0);
if(dl<dr){
if(dl<ans)qmin(lc);
if(dr<ans)qmin(rc);
}else{
if(dr<ans)qmin(rc);
if(dl<ans)qmin(lc);
}
}
}kd;
int main(){
#ifndef ONLINE_JUDGE
file("s");
#endif
fp(i,0,1)a[0].mx[i]=-inf,a[0].mi[i]=inf;tr[0]=a[0];
sd(n);sd(m);int op;
fp(i,1,n){fp(j,0,1)sd(a[i][j]);a[i].clr();}rt=kd.build(1,n,0);
while(m--){
sd(op);sd(t[0]),sd(t[1]);t.clr();
if(op&1)kd.ins(rt,0),kd.chk(rt,0);
else ans=inf,kd.qmin(rt),we(ans);
}
return Ot(),0;
}