BZOJ4571 [Scoi2016]美味

来源：互联网发布：码字软件大纲编辑：程序博客网时间：2024/04/28 05:04

区间最大异或值的经典做法是可持久化trie，但是加上一个数的话可持久化trie就变得非常的扯淡

考虑可持久化trie其实可以等价为一颗上限为2^k-1的主席树，在trie上确定一位其实相当于将答案的区间缩小的一半，也就是在主席树上向下走一层

当所有数加上x之后，我们在主席树上走的时候就不能直接调用siz[son[x][0]]，但是因为所有数都被加了，所以我们其实要查询的是左子树的区间向前窜x位之后的区间有没有数，这样的话每次走的时候在主席树上重新查[l-x,mid-x]或者[mid+1-x,r-x]来判断应该往哪边走即可，复杂度多了个log，但是n=2*10^5，不虚

然后就过了

时间复杂度O(m log^2 n)

#include<iostream>#include<cstdlib>#include<cstdio>#include<cstring>#include<cmath>#include<ctime>#include<algorithm>#include<iomanip>#include<vector>#include<stack>#include<queue>#include<map>#include<set>#include<bitset>using namespace std;#define MAXN 200010#define MAXM 6000010#define ll long long#define INF 1000000000#define MOD 1000000007#define eps 1e-8int n,m;int a[MAXN];int son[MAXM][2],siz[MAXM];int rt[MAXN],tot;int mx=(1<<19)-1;void change(int &x,int xx,int l,int r,int p){    x=++tot;    memcpy(son[x],son[xx],sizeof(son[x]));    siz[x]=siz[xx]+1;    int mid=l+r>>1;    if(l==r){        return ;    }    if(p<=mid){        change(son[x][0],son[xx][0],l,mid,p);    }else{        change(son[x][1],son[xx][1],mid+1,r,p);    }}int X,XX,D;int ask2(int x,int xx,int y,int z,int l,int r){    if(y==l&&z==r){        return siz[xx]-siz[x];    }    int mid=y+z>>1;    if(r<=mid){        return ask2(son[x][0],son[xx][0],y,mid,l,r);    }else if(l>mid){        return ask2(son[x][1],son[xx][1],mid+1,z,l,r);    }else{        return ask2(son[x][0],son[xx][0],y,mid,l,mid)+ask2(son[x][1],son[xx][1],mid+1,z,mid+1,r);    }}int ask(int l,int r,int av,int y){    if(l==r){        return av^l;    }    int mid=l+r>>1;    D--;    if(av&(1<<D)){        if(ask2(X,XX,0,mx,max(0,l-y),max(0,mid-y))){            return ask(l,mid,av,y);        }else{            return ask(mid+1,r,av,y);        }    }else{        if(ask2(X,XX,0,mx,max(0,mid+1-y),max(0,r-y))){            return ask(mid+1,r,av,y);        }else{            return ask(l,mid,av,y);        }    }     }int main(){    int i,x,y,l,r;    scanf("%d%d",&n,&m);    for(i=1;i<=n;i++){        scanf("%d",&a[i]);        change(rt[i],rt[i-1],0,mx,a[i]);    }         for(i=1;i<=m;i++){        scanf("%d%d%d%d",&x,&y,&l,&r);        X=rt[l-1];        XX=rt[r];        D=19;        printf("%d\n",ask(0,mx,x,y));    }    return 0;} /*4 41 2 3 41 4 1 42 3 2 33 2 3 34 1 2 4 */

0 0