[BZOJ4571][SCOI2016]美味(贪心+主席树)

经典问题，按位贪心，每次需要知道的是”在这一位之前的位都以确定的情况下，能否找到这一位是0/1的数”，这就是在询问[L,R]内某个值域区间是否有数，主席树即可。

#include<cstdio>
#include<algorithm>
#define rep(i,l,r) for (int i=(l); i<=(r); i++)
using namespace std;

const int N=200010,M=300000;
int n,m,b,x,l,r,nd,rt[N],a[N],ls[N*20],rs[N*20],v[N*20];

void ins(int &x,int y,int L,int R,int pos){
    x=++nd; v[x]=v[y]+1; ls[x]=ls[y]; rs[x]=rs[y];
    if (L==R) return;
    int mid=(L+R)>>1;
    if (pos<=mid) ins(ls[x],ls[y],L,mid,pos);
        else ins(rs[x],rs[y],mid+1,R,pos);
}

int que(int x,int y,int L,int R,int l,int r){
    if (L==l && r==R) return v[y]-v[x];
    int mid=(L+R)>>1;
    if (r<=mid) return que(ls[x],ls[y],L,mid,l,r);
    else if (l>mid) return que(rs[x],rs[y],mid+1,R,l,r);
        else return que(ls[x],ls[y],L,mid,l,mid)+que(rs[x],rs[y],mid+1,R,mid+1,r);
}

int main(){
    freopen("bzoj4571.in","r",stdin);
    freopen("bzoj4571.out","w",stdout);
    scanf("%d%d",&n,&m);
    rep(i,1,n) scanf("%d",&a[i]),ins(rt[i],rt[i-1],0,M,a[i]);
    rep(i,1,m){
        scanf("%d%d%d%d",&b,&x,&l,&r); int a=0;
        for (int j=17; ~j; j--){
            if (!(b&(1<<j))) a|=1<<j;
            int L=max(a-x,0),R=(a|((1<<j)-1))-x;
            if (R<0 || !que(rt[l-1],rt[r],0,M,L,R)) a^=1<<j;
        }
        printf("%d\n",a^b);
    }
    return 0;
}