二维恶心强制在线卡空间数点（主席树+KD树）

为啥要主席树加KD树？
因为这个恶心出题人强制在线卡空间卡时间。
主席树处理询问$3$
KD树处理$=1,2$
勉强卡过时限。

#include<bits/stdc++.h>
#define maxn 800005
#define Ct const
#define rep(i,j,k) for(int i=(j),LIM=(k);i<=LIM;i++)
using namespace std;

char cb[1<<16],*cs=cb,*ct=cb;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<16,stdin),cs==ct)?0:*cs++)
void read(int &res){
   char ch;bool f =0;
	for(;!isdigit(ch=getc());) if(ch=='-') f=1;
	for(res=ch-'0';isdigit(ch=getc());res=res*10+ch-'0');
	(f) && (res=-res);
}

int tot,rtj,rtp,D,c[maxn];
struct node{
   
	int a[2],sz,sm,v,mn[2],mx[2],ch[2],d;
	void upd(node *P){
   
		sm=v,sz=1;
		rep(i,0,1) mn[i]=mx[i]=a[i];
		rep(j,0,1) if(ch[j]){
   
			rep(i,0,1)
				mn[i]=min(mn[i],P[ch[j]].mn[i]),
				mx[i]=max(mx[i],P[ch[j]].mx[i]);
			sm+=P[ch[j]].sm,sz+=P[ch[j]].sz;
		}
	}
}P[maxn];
#define upd() upd(P)
bool cmp(Ct int &u,Ct int &v){
    return P[u].a[D]<P[v].a[D]; }
void Build(int &u,int l,int r,int nD){
   
	if(l>r) return;
	int m=l+r>>1;D=nD,nth_element(c+l,c+m,c+r+1,cmp);u=c[m];
	Build(P[u].ch[0],l