【bzoj 4066】 简单题

题目

显然这就是让我们在二维上数个点

如果没有强制在线就随便做啦，扫描线+主席树应该是最好的选择

但是现在强制在线并且卡了树套树的空间，于是只能上\(kdt\)了

我们还是维护一下每个子树分割出来的矩形的四个坐标

查询的时候发现如果当前整个子树都被查询的矩形包含，我们就直接加入答案

如果完全没有交，我们就直接退出

否则我们判断一下是否要把当前子树的根加入答案，之后递归左右两边的子树就好了

在随机数据下拥有一个\(log\)的美妙复杂度，但是卡卡就变成根号了

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define re register
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
inline int read() {
	char c=getchar();int x=0;while(c<'0'||x>'9') c=getchar();
	while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
const double alph=0.7;
const int maxn=300005;
struct Point{int x[2],w;}p[maxn],a[maxn];
int l[maxn],r[maxn],sz[maxn],mx[maxn][2],mi[maxn][2],d[maxn],st[maxn];
int n,m,cnt,top,rt,lst,op,num;
inline int cmp(Point A,Point B) {return A.x[op]<B.x[op];}
inline int newnode() {
	if(top) return st[top--];
	return ++cnt;
}
inline void pushup(int k) {
	mi[k][0]=mx[k][0]=p[k].x[0];
	mi[k][1]=mx[k][1]=p[k].x[1];
	for(re int i=0;i<2;i++) {
		if(l[k]) mx[k][i]=max(mx[k][i],mx[l[k]][i]),
				 mi[k][i]=min(mi[k][i],mi[l[k]][i]);
		if(r[k]) mx[k][i]=max(mx[k][i],mx[r[k]][i]),
				 mi[k][i]=min(mi[k][i],mi[r[k]][i]);
	}
	sz[k]=sz[l[k]]+sz[r[k]]+1;
	d[k]=d[l[k]]+d[r[k]]+p[k].w;
}
int build(int x,int y,int o) {
	if(x>y) return 0;
	int mid=x+y>>1,k=newnode();
	op=o;std::nth_element(a+x,a+mid,a+y+1,cmp);p[k]=a[mid];
	l[k]=build(x,mid-1,o^1),r[k]=build(mid+1,y,o^1);
	pushup(k);return k;
}
void kill(int k) {
	if(l[k]) kill(l[k]);
	st[++top]=k;a[++num]=p[k];
	if(r[k]) kill(r[k]);
}
inline int check(int k,int o) {
	if(max(sz[l[k]],sz[r[k]])>sz[k]*alph)
		return num=0,kill(k),build(1,num,o);
	return k;
}
int ins(Point a,int k,int o) {
	if(!k) {
		k=newnode();p[k]=a;pushup(k);
		return k;
	}
	if(a.x[o]<=p[k].x[o]) l[k]=ins(a,l[k],o^1);
		else r[k]=ins(a,r[k],o^1);
	pushup(k);return check(k,o);
}
inline int pd(int x,int y,int xx,int yy,int k) {
	return (x<=mi[k][0]&&xx>=mx[k][0]&&y<=mi[k][1]&&yy>=mx[k][1]);
}
inline int chk(int x,int y,int xx,int yy,int k) {
	return (mi[k][0]>xx||mx[k][0]<x||mi[k][1]>yy||mx[k][1]<y);
}
int query(int x,int y,int xx,int yy,int k) {
	if(!k||chk(x,y,xx,yy,k)) return 0;
	if(pd(x,y,xx,yy,k)) return d[k];
	int now=0;
	if(x<=p[k].x[0]&&xx>=p[k].x[0]&&y<=p[k].x[1]&&yy>=p[k].x[1]) now+=p[k].w;
	if(pd(x,y,xx,yy,l[k])&&l[k]) now+=d[l[k]];else now+=query(x,y,xx,yy,l[k]);
	if(pd(x,y,xx,yy,r[k])&&r[k]) now+=d[r[k]];else now+=query(x,y,xx,yy,r[k]);
	return now;
}
int main() {
	n=read();
	int opt,x,y,xx,yy,val;
	while(1) {
		opt=read();
		if(opt==1) {
			x=read()^lst,y=read()^lst;val=read()^lst;
			Point a;a.x[0]=x,a.x[1]=y,a.w=val;
			rt=ins(a,rt,0);
		}
		if(opt==2) {
			x=read()^lst,y=read()^lst,xx=read()^lst,yy=read()^lst;
			lst=query(x,y,xx,yy,rt);
			printf("%d\n",lst);
		}
		if(opt==3) return 0;
	}
	return 0;
}