P3293 [SCOI2016]美味

解析

有点懊恼。
见到异或后无脑上trie，想到了和正解差不多的值域限制，但在trie上复杂度总是不对…
生活在踹树上。

格局打开！
线段树是解决异或问题的优秀工具。
把值域限制从trie上扒下来，用主席树支持区间值域查询即可。
复杂度 $O(m{\log }^{2}V)$。

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned ll
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define ok debug("OK\n")


inline ll read() {
  ll x(0),f(1);char c=getchar();
  while(!isdigit(c)) {if(c=='-') f=-1;c=getchar();}
  while(isdigit(c)) {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
  return x*f;
}

const int N=2e5+100;
//const int mod=998244353;

bool mem1;

inline ll ksm(ll x,ll k,int mod){
  ll res=1;
  while(k){
    if(k&1) res=res*x%mod;
    x=x*x%mod;
    k>>=1;
  }
  return res;
}

int n,m;

int mx;
struct node{
  int ls,rs,siz;
}tr[N*30];
int rt[N],tot;
#define mid ((l+r)>>1)
inline int copy(int k){
  tr[++tot]=tr[k];return tot;
}
inline void upd(int &k,int l,int r,int p){
  k=copy(k);
  if(l==r){
    tr[k].siz++;return;
  }
  if(p<=mid) upd(tr[k].ls,l,mid,p);
  else upd(tr[k].rs,mid+1,r,p);
  tr[k].siz++;
}
int ask(int k,int l,int r,int x,int y){
  if(!k) return 0;
  if(x<=l&&r<=y) return tr[k].siz;
  int res(0);
  if(x<=mid) res+=ask(tr[k].ls,l,mid,x,y);
  if(y>mid) res+=ask(tr[k].rs,mid+1,r,x,y);
  return res;
}
#undef mid
int a[N],mi[50];
void init(){
  mi[0]=1;
  for(int i=1;i<=18;i++) mi[i]=mi[i-1]<<1;
  for(int i=1;i<=n;i++){
    a[i]=read();
    mx=max(mx,a[i]);
  }
  for(int i=1;i<=n;i++){
    rt[i]=rt[i-1];
    upd(rt[i],0,mx,a[i]);
  }
}
void work(){
  int b=read(),x=read(),p=read(),q=read(),ans(0),cur(0);
  for(int k=18;k>=0;k--){
    bool f=(b&mi[k])?1:0;    
    int goal=cur|(mi[k]*(!f));    
    int l=max(0,goal-x),r=min(mx,goal+mi[k]-x-1);
    //printf("  k=%d (%d %d) goal=%d cur=%d\n",k,l,r,goal,cur);
    if(l<=r&&ask(rt[q],0,mx,l,r)-ask(rt[p-1],0,mx,l,r)){
      ans|=mi[k];
      cur=goal;
    }
    else cur=goal^mi[k];     
  }
  printf("%d\n",ans);
}

bool mem2;
signed main(){
  #ifndef ONLINE_JUDGE
  freopen("a.in","r",stdin);
  freopen("a.out","w",stdout);
  #endif
  //debug("mem=%.4lf\n",abs(&mem2-&mem1)/1024./1024.);
  n=read();m=read();
  init();
  while(m--) work();
  return 0;
}
/*

*/