【xsy1098】第k小 可持久化trie

题目大意：你要维护一个长度为$n$的序列，资瓷对整个序列$xor,and,or$一个数，以及区间第k小查询。

数据范围：$n≤50000$，所有数字$＜2^{31}$。

 

此题甚妙

我们不难想出没有位运算的区间第k大查询，直接可持久化trie就可以了。

考虑此题只有xor操作。

我们记一个$last$表示之前所有异或操作的数的异或和。

查询区间第$k$小时，随便查询下就行了（在比较时假装交换$trie$的左右子树）

 

考虑存在$or$和$and$操作，我们设对序列操作的数字为$x$，设$vis[k]$表示整个序列的第$k$位是否被全部归零过。

不难发现，$and$操作对第i位有意义，当且仅当数字$x$的第$i$位为$0$，（$or$情况类似等下会讲），在这种情况下，对于序列上的所有数，第$i$位都会归零，（last的第$i$位也会归零）

同理，$or$操作对第i位有意义，当且仅当数字$x$的第$i$位为$1$，对于序列上所有数，第$i$位都会变为1（$last$的第$i$位也会变为1）

实际上两个操作本质上是相同的，我们可以用$and$操作和$last|=x$去实现$or$操作。

不难发现，对于数列的每一位，归零操作至多只需要执行一次。

那么，对于每个需要执行归零操作的$or$或$and$操作，我们重建一次可持久化$trie$就可以了。

 

时间复杂度$O(n\ log^2 V+m\ log V)$其中V为数列中的最大值。

#include<bits/stdc++.h>
#define M 50005
#define INF 2147483647
using namespace std;

struct trie{int a[2],siz;}a[M*33]={0}; int use=0;
void add(int &x,int dep,int zhi){
    a[++use]=a[x]; a[x=use].siz++;
    if(dep==-1) return;
    bool k=zhi&(1<<dep);
    if(k) add(a[x].a[1],dep-1,zhi);
    else add(a[x].a[0],dep-1,zhi);
}    
int num[M]={0},n,m,zhi=0,root[M]={0},vis[32]={0},pls=INF;

int query(int x,int y,int dep,int K){
    if(dep==-1) return 0;
    bool k=zhi&(1<<dep);
    int now=a[a[y].a[k]].siz-a[a[x].a[k]].siz,res=1<<dep;
    if(K<=now) res=query(a[x].a[k],a[y].a[k],dep-1,K);
    else res+=query(a[x].a[k^1],a[y].a[k^1],dep-1,K-now);
    return res;
}

void build(){
    for(int i=1;i<=n;i++) num[i]=num[i]&pls;
    use=0;
    for(int i=1;i<=n;i++){
        root[i]=root[i-1];
        add(root[i],30,num[i]);
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",num+i);
    build();
    while(m--){
        char op[4]; int x,l,r,k; pls=INF;
        scanf("%s",op);
        if(op[0]=='X'){scanf("%d",&x); zhi^=x;}
        if(op[0]=='O'){
            scanf("%d",&x);
            for(int i=0;i<=30;i++)
            if(((1<<i)&x)){
                if(vis[i]==0) vis[i]=1,pls^=(1<<i);
                zhi|=(1<<i); 
            }
            if(pls==INF) continue;
            build();
        }
        if(op[1]=='n'){
            scanf("%d",&x);
            for(int i=0;i<=30;i++)
            if(((1<<i)&x)==0){
                if(vis[i]==0) vis[i]=1,pls^=(1<<i);
                zhi&=(INF-(1<<i)); 
            }
            if(pls==INF) continue;
            build();
        }
        if(op[1]=='s'){
            scanf("%d%d%d",&l,&r,&k);
            printf("%d\n",query(root[l-1],root[r],30,k));
        }
    }
}