BZOJ4571：[SCOI2016]美味

浅谈主席树：https://www.cnblogs.com/AKMer/p/9956734.html

题目传送门：https://www.lydsy.com/JudgeOnline/problem.php?id=4571

直接从高位到低位贪心。我们可以一位一位确定还没与\(b\)异或起来的\(ans\)，最后再与\(b\)异或。假设前面的高位确定完之后，根据当前位贪心需要什么，可以在一段区间内找到我们需要的数字。如果存在，那么这一位就是贪心的结果，否则就是另一种。在一段区间内找值域在某一段区间的数字，正好可以用主席树搞定。

时间复杂度：\(O(nlognlogv)\)

空间复杂度：\(O(nlogv)\)

代码如下：

#include <cstdio>
#include <algorithm>
using namespace std;
const int maxn=2e5+5;
int n,m;
int rt[maxn];
int read() {
	int x=0,f=1;char ch=getchar();
	for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
	for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
	return x*f;
}
struct tree_node {
	int cnt,ls,rs;
};
struct chairman_tree {
	int tot;
	tree_node tree[maxn*18];
	void ins(int lst,int &now,int l,int r,int pos) {
		now=++tot;tree[now]=tree[lst];
		tree[now].cnt++;
		if(l==r)return;
		int mid=(l+r)>>1;
		if(pos<=mid)ins(tree[lst].ls,tree[now].ls,l,mid,pos);
		else ins(tree[lst].rs,tree[now].rs,mid+1,r,pos);
	}
	int query(int x,int y,int l,int r,int L,int R) {
		if(L<=l&&r<=R)return tree[y].cnt-tree[x].cnt;
		int mid=(l+r)>>1,res=0;
		if(L<=mid)res+=query(tree[x].ls,tree[y].ls,l,mid,L,R);
		if(R>mid)res+=query(tree[x].rs,tree[y].rs,mid+1,r,L,R);
		return res;
	}
}T;
int main() {
	n=read(),m=read();
	for(int i=1;i<=n;i++) {
		int x=read();
		T.ins(rt[i-1],rt[i],0,1e5,x);
	}
	for(int i=1;i<=m;i++) {
		int b=read(),x=read(),l=read(),r=read(),ans=0;
		for(int j=17;~j;j--) {
			int L,R,bit;
			if(b&(1<<j))L=ans,R=ans+(1<<j)-1,bit=0;
			else L=ans+(1<<j),R=ans+(1<<(j+1))-1,bit=1;
			L=max(L-x,0);R=min(100000,R-x);
			bool bo=T.query(rt[l-1],rt[r],0,1e5,L,R);
			if(!bo)bit^=1;ans|=bit<<j;
		}ans^=b;
		printf("%d\n",ans);
	}
	return 0;
}