[可持久化字典树 分块] BZOJ 2741 【FOTILE模拟赛】L

先异或取前缀和

最大异或值 可持久化字典树

然后么 区间问题 分块

用f[i][j] 表示 第i块中的数与i-j块的数的最大异或值 然后边角处理下

复杂度一个根号一个log

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
 
inline char nc()
{
    static char buf[100000],*p1=buf,*p2=buf;
    if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
    return *p1++;
}
 
inline void read(int &x)
{
    char c=nc(),b=1;
    for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
    for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}
 
const int N=13005;
 
int ncnt;
int root[N];
int ch[N*35][2],sum[N*35];
 
inline void insert(int &rt,int x,int t){
    rt=++ncnt; 
    int y; sum[y=rt]=sum[x]+1;
    for (int i=30;~i;i--)
    {
        if ((t>>i)&1)
            ch[y][0]=ch[x][0],y=(ch[y][1]=++ncnt),x=ch[x][1];
        else
            ch[y][1]=ch[x][1],y=(ch[y][0]=++ncnt),x=ch[x][0];
        sum[y]=sum[x]+1;
    }
}
 
inline int query(int l,int r,int t){
    int ret=0,tmp;
    for (int i=30;~i;i--)
    {
        tmp=(t>>i)&1;
        if (sum[ch[r][tmp^1]]-sum[ch[l][tmp^1]])
            ret+=(1<<i),r=ch[r][tmp^1],l=ch[l][tmp^1];
        else
            r=ch[r][tmp],l=ch[l][tmp];
    }
    return ret;
}
 
int n,Q,B;
int cnt,pos[N],st[505],ed[505];
int a[N];
 
int f[505][505];
 
inline void Pre()
{
    B=sqrt(n);
    for (int i=1;i<=n;i++) pos[i]=(i-1)/B+1;
    cnt=pos[n];
    for (int i=1;i<=cnt;i++)
        st[i]=B*(i-1)+1,ed[i]=B*i;
    ed[cnt]=min(ed[cnt],n);
    for (int i=1;i<=n;i++)
    {
        int p=pos[i];
        for (int j=p;j<=cnt;j++)
            f[p][j]=max(f[p][j],query(root[ed[p-1]],root[ed[j]],a[i]));
    }
}
 
int lastans=0;
 
inline void Solve(){
    int l,r,lpos,rpos,ret=0;
    read(l); read(r);
    l=((ll)l+lastans)%(n-1)+1; r=((ll)r+lastans)%(n-1)+1;
    if (l>r) swap(l,r);
    r++;
    lpos=pos[l],rpos=pos[r];
    if (lpos==rpos)
    {
        for (int i=l;i<=r;i++)
            ret=max(ret,query(root[r],root[i],a[i]));
        printf("%d\n",lastans=ret);
        return;
    }
    for (int i=l;i<=ed[lpos];i++)
        ret=max(ret,query(root[r],root[l-1],a[i]));
    for (int i=st[rpos];i<=r;i++)
        ret=max(ret,query(root[r],root[l-1],a[i]));
    for (int i=lpos+1;i<=rpos-1;i++)
        ret=max(ret,f[i][rpos-1]);
    printf("%d\n",lastans=ret);
}
 
int main()
{
    read(n); read(Q); n++;
    insert(root[1],root[0],0);
    for (int i=2;i<=n;i++)
    {
        read(a[i]); a[i]^=a[i-1];
        insert(root[i],root[i-1],a[i]);
    }
    Pre();
    while (Q--)
        Solve();
    return 0;
}