BZOJ 2741 L (可持久化01Trie+分块)

题目大意：给你一个序列，共有$q$个询问，每次询问区间$[L,R]$内最大连续字段异或和，强制在线，$n<=12000,m<=5000$

有个细节没处理好$WA$了好久..还有一次$ans$没清零

先对序列建出可持久化$01Trie$

分块预处理出，任意两块所覆盖区域的最大$xor$和，枚举右侧块内的每个数，然后在$01Trie$里查找即可，预处理总时间$O(n \sqrt n)$

对于每次询问，中间部分的答案可以$O(1)$得到

边界情况，左侧不完整块内的每个数都要分别作为区间左端点和右端点在$01Trie$里查找最大值

右侧不完整块内的每个数作为区间右端点在01Trie里查找最大值即可，因为一旦出现它作为左端点的情况，它右面的数一定会作为右端点查询到它

时间$O(q\sqrt n+n\sqrt n)$

数据还挺强的

 #include <cmath>
 #include <queue>
 #include <vector>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define N1 12100
 #define N2 420100
 #define M1 120
 #define ll long long
 #define dd double
 #define uint unsigned int
 using namespace std;

 int gint()
 {
     int ret=,fh=;char c=getchar();
     while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
     while(c>=''&&c<=''){ret=ret*+c-'';c=getchar();}
     return ret*fh;
 }

 int n,m,sq,tq;
 uint bin[];

 struct Trie{
 int ch[N2][],num[N2],root[N1],tot;
 void build()
 {
     root[]=tot=;int x=;
     for(int i=;i>=;i--){
         ch[x][]=++tot;
         x=ch[x][],num[x]=;
     }
 }
 void insert(uint s,int rt1,int rt2,int w)
 {
     int x,y,p;
     y=root[rt1];
     x=root[rt2]=++tot;
     for(int i=;i>=;i--){
         p=(s&bin[i])?:;
         ch[x][p]=++tot;
         ch[x][p^]=ch[y][p^];
         num[ch[x][p]]=num[ch[y][p]]+w;
         x=ch[x][p],y=ch[y][p];
     }
 }
 uint query(uint s,int l,int r)
 {
     int x,y,p;uint ans=;
     x=root[r];
     y=l<?:root[l];
     for(int i=;i>=;i--){
         p=(s&bin[i])?:;
         if(num[ch[x][p^]]-num[ch[y][p^]]>){
             x=ch[x][p^],y=ch[y][p^];
             ans|=bin[i];
         }else if(num[ch[x][p]]-num[ch[y][p]]>){
             x=ch[x][p],y=ch[y][p];
         }else break;
     }return ans;
 }
 }T;

 uint a[N1],pa[N1];
 uint s1[M1][M1];
 int L[M1],R[M1];

 int main()
 {
     scanf("%d%d",&n,&m);
     for(int i=;i<=;i++)
         bin[i]=(<<i);
     int x,y;
     T.build();
     for(int i=;i<=n;i++)
     {
         a[i]=gint();
         pa[i]=pa[i-]^a[i];
         T.insert(pa[i],i-,i,);
     }
     sq=sqrt(n);tq=n/sq;
     for(int i=;i<=tq+;i++)
         L[i]=(i-)*sq+,R[i]=min(n,i*sq);
     for(int i=;i<=tq;i++)
     for(int j=i;j<=tq;j++)
     {
         s1[i][j]=s1[i][j-];
         for(int k=L[j];k<=R[j];k++)
             s1[i][j]=max(s1[i][j],T.query(pa[k],L[i]-,k-));
     }
     uint ans=;int px,py;
     for(int j=;j<=m;j++)
     {
         x=gint(),y=gint();
         x=(ans+x)%n+,y=(ans+y)%n+;
         if(x>y) swap(x,y);
         px=(x-)/sq+,py=(y-)/sq+;
         if(px!=py){
             ans=s1[px+][py-];
             for(int i=x;i<=R[px];i++)
                 ans=max(ans,T.query(pa[i],x-,i-));
             for(int i=x-;i<=R[px]-;i++)
                 ans=max(ans,T.query(pa[i],i,y));
             for(int i=L[py];i<=y;i++)
                 ans=max(ans,T.query(pa[i],x-,i-));
         }else{
             ans=;
             for(int i=x;i<=y;i++)
                 ans=max(ans,T.query(pa[i],x-,i-));
         }
         printf("%u\n",ans);
     }
     return ;
 }