bzoj 2741: 【FOTILE模拟赛】L

Description

FOTILE得到了一个长为N的序列A，为了拯救地球，他希望知道某些区间内的最大的连续XOR和。
即对于一个询问，你需要求出max(Ai xor Ai+1 xor Ai+2 ... xor Aj)，其中l<=i<=j<=r。
为了体现在线操作，对于一个询问(x,y)：
l = min ( ((x+lastans)
mod N)+1 , ((y+lastans) mod N)+1 ).

r = max ( ((x+lastans) mod N)+1 , ((y+lastans) mod N)+1 ).
其中lastans是上次询问的答案，一开始为0。

Input

第一行两个整数N和M。
第二行有N个正整数，其中第i个数为Ai，有多余空格。
后M行每行两个数x,y表示一对询问。

Output

共M行，第i行一个正整数表示第i个询问的结果。

Sample Input

3 3

1 4 3

0 1

0 1

4 3

Sample Output

5

7

7

HINT

HINT

N=12000，M=6000，x,y,Ai在signed longint范围内。

分块+可持久化trie
用f[i][j]表示第i块开头到位置j的最大ans
f[i][j]=max(f[i][j-1],ask(head[i],j-1,a[j]))
这里处理是nsqrt(n)logn
然后每次询问的时候把整块的一起处理
再枚举前面部分的位置，更新ans即可

#include<cmath>
#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
struct tree
{
int l,r;
int ll,rr;
}tr[5000001];
int root[600001];
long long px[32];
int tot;
inline int add(int p,int d,int dep)
{
if(dep==31)
return 0;
if(d>=px[dep])
{
tot++;
int pp=tot;
tr[pp].l=tr[p].l;
tr[pp].ll=tr[p].ll;
tr[pp].r=add(tr[p].r,d-px[dep],dep+1);
tr[pp].rr=tr[p].rr+1;
return pp;
}
else
{
tot++;
int pp=tot;
tr[pp].l=add(tr[p].l,d,dep+1);
tr[pp].ll=tr[p].ll+1;
tr[pp].r=tr[p].r;
tr[pp].rr=tr[p].rr;
return pp;
}
}
inline int ask(int l,int r,int x,int dep)
{
if(dep==31)
return 0;
if(x>=px[dep])
{
if(tr[r].ll-tr[l].ll>0)
return px[dep]+ask(tr[l].l,tr[r].l,x-px[dep],dep+1);
else
return ask(tr[l].r,tr[r].r,x-px[dep],dep+1);
}
else
{
if(tr[r].rr-tr[l].rr>0)
return px[dep]+ask(tr[l].r,tr[r].r,x,dep+1);
else
return ask(tr[l].l,tr[r].l,x,dep+1);
}
}
int a[100001];
int x[100001];
int belong[100001];
int head[1001];
int f[1001][20001];
int main()
{
//	 freopen("data.in","r",stdin);
//	 freopen("data.out","w",stdout);
long long n;
int m;
scanf("%lld%d",&n,&m);
int i,j;
px[0]=1;
for(i=1;i<=30;i++)
px[i]=px[i-1]*2;
for(i=0;i<=15;i++)
{
int t=px[i];
px[i]=px[30-i];
px[30-i]=t;
}
root[0]=1;
add(0,0,0);
for(i=1;i<=n;i++)
{
scanf("%lld",&x[i]);
a[i]=(a[i-1]^x[i]);
root[i]=tot+1;
add(root[i-1],a[i],0);
}
int nt=sqrt(n);
int nx=0;
for(i=1;i<=n;i++)
{
belong[i]=(i-1)/nt+1;
if(belong[i]!=belong[i-1])
{
nx++;
head[belong[i]]=i;
}
}
head[nx+1]=n+1;
for(i=1;i<=nx;i++)
{
f[i][head[i]]=x[head[i]];
for(j=head[i]+1;j<=n;j++)
{
int rt=head[i]-2;
if(rt<0)
rt=0;
else
rt=root[rt];
f[i][j]=max(f[i][j-1],ask(rt,root[j-1],a[j],0));
}
}
long long s,t;
long long last=0;
for(i=1;i<=m;i++)
{
scanf("%lld%lld",&s,&t);
int ss=(s+last)%n+(long long)1,tt=(t+last)%n+(long long)1;
s=min(ss,tt);
t=max(ss,tt);
/*	  if(s==t)
{
printf("%d\n",x[s]);
continue;
}*/
int nxt=belong[s];
if(head[nxt]!=s)
nxt++;
int ans=f[nxt][t];
for(j=s;j<=min(head[nxt]-1,(int)t);j++)
{
ans=max(ans,ask(root[s-1],root[t],a[j-1],0));
}
printf("%lld\n",ans);
last=ans;
}
return 0;
}