bzoj2648 SJY摆棋子

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)....