https://www.lydsy.com/JudgeOnline/problem.php?id=2648
如果是带插入的话就要写动态生成节点了,类似主席树的写法。
下题是求曼哈顿距离,曼哈顿距离可以维护每个子树中所有节点的x的最大最小值,y的最大最小值,然后query的时候计算左右子树的非法值,非法值是指在这个子树的超平面(也就是x,y范围外的和)外的和,左子树大还是右子树大,就先去哪里,如果非法值已经超过了当前记录的ans的话,就不去那个子树找了。
这个复杂度我感觉还是很玄学,
#include<bits/stdc++.h>
using namespace std;
const int maxl=5e5+10;
int n,m,idx,cnt,root,ans;
struct point
{
int x[2];
bool operator < (const point &b)const
{
return x[idx]<b.x[idx];
}
}a[maxl];
struct kdtree
{
int ls,rs,mi[2],mx[2];
point p;
}tree[maxl<<2];
inline void push_up(int k)
{
int s;
if(tree[k].ls)
{
s=tree[k].ls;
for(int i=0;i<2;i++)
{
tree[k].mx[i]=max(tree[k].mx[i],tree[s].mx[i]);
tree[k].mi[i]=min(tree[k].mi[i],tree[s].mi[i]);
}
}
if(tree[k].rs)
{
s=tree[k].rs;
for(int i=0;i<2;i++)
{
tree[k].mx[i]=max(tree[k].mx[i],tree[s].mx[i]);
tree[k].mi[i]=min(tree[k].mi[i],tree[s].mi[i]);
}
}
}
inline void build(int &k,int l,int r,int d)
{
idx=d;k=++cnt;
int mid=(l+r)>>1;
nth_element(a+l,a+mid,a+r+1);
tree[k].p=a[mid];
for(int i=0;i<2;i++)
tree[k].mx[i]=tree[k].mi[i]=a[mid].x[i];
if(l<mid)
build(tree[k].ls,l,mid-1,d^1);
if(mid<r)
build(tree[k].rs,mid+1,r,d^1);
push_up(k);
}
inline void insert(int &k,point pt,int d)
{
if(!k)
{
k=++cnt;
for(int i=0;i<2;i++)
tree[k].mx[i]=tree[k].mi[i]=pt.x[i];
tree[k].p=pt;
return;
}
idx=d;
if(pt<tree[k].p) insert(tree[k].ls,pt,d^1);
else insert(tree[k].rs,pt,d^1);
push_up(k);
}
inline void prework()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
scanf("%d%d",&a[i].x[0],&a[i].x[1]);
build(root,1,n,0);
}
inline int dist(point a,point b)
{
return abs(a.x[0]-b.x[0])+abs(a.x[1]-b.x[1]);
}
inline int calc(int k,point p)
{
int ret=0;
for(int i=0;i<2;i++)
{
ret+=max(0,tree[k].mi[i]-p.x[i]);
ret+=max(0,p.x[i]-tree[k].mx[i]);
}
return ret;
}
inline void query(int k,point pt,int d)
{
if(!k) return;
ans=min(ans,dist(pt,tree[k].p));
int lv=2e9,rv=2e9;
if(tree[k].ls) lv=calc(tree[k].ls,pt);
if(tree[k].rs) rv=calc(tree[k].rs,pt);
if(lv<rv)
{
if(lv<ans) query(tree[k].ls,pt,d^1);
if(rv<ans) query(tree[k].rs,pt,d^1);
}
else
{
if(rv<ans) query(tree[k].rs,pt,d^1);
if(lv<ans) query(tree[k].ls,pt,d^1);
}
}
inline void mainwork()
{
int op,x,y;point d;
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&op,&x,&y);
d.x[0]=x;d.x[1]=y;
if(op==1)
insert(root,d,0);
else
{
ans=2e9;
query(root,d,0);
printf("%d\n",ans);
}
}
}
int main()
{
prework();
mainwork();
return 0;
}
但是只能假装是sqrt(n)....