- 【题目地址】
题意见题面。
首先,考虑异或最大可以按照二进制位重高位到低位贪心,所以我们对于一个询问的 b i b_i bi,假如当前答案要确定的这一位上 b i b_i bi为 0 0 0,那么我们只需要权值在 ( ⋯ 10000 ) 2 ∼ ( ⋯ 11111 ) 2 (\cdots10000)_2\sim(\cdots11111)_2 (⋯10000)2∼(⋯11111)2的菜品(也就是这一位是 1 1 1的范围,省略号表示前面已经确定了的),不断的修改左右范围,一位一位确定即可,复杂度为 l o g v logv logv的。
由于询问还要加上 x i x_i xi,所以我们反过来思考,那么真正的查询区间就是原来的 L ∼ R L\sim R L∼R变成了 L − x ∼ R − x L-x\sim R-x L−x∼R−x(注意和0取max)。
所以每一位我们只需查询在当前的可选取的菜品内有没有 a i a_i ai在这个范围内的,这个可以用主席树在 l o g n logn logn的时间内实现。
那么总复杂度就是 O ( n l o g n + m l o g 2 n ) O(nlogn+mlog^2n) O(nlogn+mlog2n),那么这个题就解决了。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int M=2e5+10,Log=20;
const int UP=2e5+1;
int sum[M*Log],ls[M*Log],rs[M*Log],tot,root[M];
int n,m;
void update(int pre,int &o,int l,int r,int p){
o=++tot;ls[o]=ls[pre];rs[o]=rs[pre];sum[o]=sum[pre];
if(l==r){++sum[o];return;}
int mid=l+r>>1;
if(p<=mid) update(ls[pre],ls[o],l,mid,p);
else update(rs[pre],rs[o],mid+1,r,p);
sum[o]=sum[ls[o]]+sum[rs[o]];
}
int query(int o,int l,int r,int L,int R){
if(L<=l&&r<=R) return sum[o];
int mid=l+r>>1;
if(L>mid) return query(rs[o],mid+1,r,L,R);
if(R<=mid) return query(ls[o],l,mid,L,R);
return query(ls[o],l,mid,L,R)+query(rs[o],mid+1,r,L,R);
}
int query_area(int l,int r,int vl,int vr){
vr=min(vr,UP);vl=min(vl,UP);
return query(root[r],0,UP,vl,vr)-query(root[l-1],0,UP,vl,vr);
}
int b,x,l,r;
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%d",&x);
update(root[i-1],root[i],0,UP,x);
}
while(m--){
scanf("%d%d%d%d",&b,&x,&l,&r);
int ans=0;
for(int i=19,low=0,up=(1<<20)-1;i>=0;i--){
if((b>>i)&1){
up^=(1<<i);
if(!query_area(l,r,max(0,low-x),max(0,up-x)))up^=(1<<i),low^=(1<<i);
else ans|=(1<<i);
}else{
low^=(1<<i);
if(!query_area(l,r,max(0,low-x),max(0,up-x)))low^=(1<<i),up^=(1<<i);
else ans|=(1<<i);
}
}
printf("%d\n",ans);
}
return 0;
}