BZOJ2741 【FOTILE模拟赛】L 可持久化Trie+分块

题意：给定长度为n的序列，m次询问区间最大连续子段异或和

n<=12000, m<=6000

Sol：

首先搞前缀异或和，问题变成询问区间选两个数最大异或和

然后考虑暴力 设f[i][j]表示区间i - j 答案 f[i][j]=max(f[i+1][j],calc(a[i],i+1,j)) calc(x,l,r)表示x与[l,r]中数的最大异或值

然后分块大法好，f[i][j]表示i块到j块的，询问时整块直接得到，两边的暴力一下

复杂度O(n^1.5 logn + m*n^0.5 log n)

Code:

#include<bits/stdc++.h>
#define debug(x) cout<<#x<<"="<<x<<endl
using namespace std;
const int N = 12009 , M = 309;
const int Log = 30;

int a[N],rt[N],bel[N],f[M][M];
struct Blocks{int l,r;}blo[M];
int n,m,siz,cnt,last,lim;
struct Trie
{
	int tot;
	int son[N*Log*6][2],cnt[N*Log*6];
	int insert(int x,int num)
	{
		int tmp,y;tmp=y=++tot;
		for(int i=Log;i>=0;i--)
		{
			son[y][0]=son[x][0];son[y][1]=son[x][1];
			cnt[y]=cnt[x]+1;
			int wh=(num>>i)&1;
			son[y][wh]=++tot;
			x=son[x][wh];y=son[y][wh];
		}cnt[y]=cnt[x]+1;
		return tmp;
	}
	int query(int l,int r,int num)
	{
		int res=0;
		for(int i=Log;i>=0;i--)
		{
			int wh=(num>>i)&1;
			if(cnt[son[r][wh^1]]-cnt[son[l][wh^1]]) res|=(1<<i),r=son[r][wh^1],l=son[l][wh^1];
			else r=son[r][wh],l=son[l][wh];
		}
		return res;
	}
}uuz;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

int query(int l,int r)
{
	int res=0;
	if(bel[l]+1>=bel[r])
	{
		for(int i=l;i<=r;i++) res=max(res,uuz.query(rt[l-1],rt[r],a[i]));
		return res;
	}
	for(int i=l;bel[i]==bel[l];i++) res=max(res,uuz.query(rt[l-1],rt[r],a[i]));
	for(int i=r;bel[i]==bel[r];i--) res=max(res,uuz.query(rt[l-1],rt[r],a[i]));
	return max(res,f[bel[l]+1][bel[r]-1]);
}

int main()
{
	n=read()+1;m=read();lim=n-1;
	for(int i=2;i<=n;i++) a[i]=read()^a[i-1];
	for(int i=1;i<=n;i++) rt[i]=uuz.insert(rt[i-1],a[i]);
	siz=sqrt(n*1.0);cnt=(n-1)/siz+1;
	for(int i=1;i<=n;i++) bel[i]=(i-1)/siz+1;
	for(int i=1;i<=n;i++) blo[bel[i]].r=i;
	for(int i=n;i>=1;i--) blo[bel[i]].l=i;
	for(int i=1;i<=cnt;i++)
	{
		int L=blo[i].l,R=blo[i].r;
		for(int j=L;j<=R;j++) f[i][i]=max(f[i][i],uuz.query(rt[L-1],rt[R],a[j]));
	}
	for(int len=2;len<=cnt;len++)
		for(int l=1;l+len-1<=cnt;l++)
		{
			int r=l+len-1,L=blo[l].l,R=blo[r].r;
			for(int i=blo[l].l;i<=blo[l].r;i++) f[l][r]=max(f[l][r],uuz.query(rt[L-1],rt[R],a[i]));
			f[l][r]=max(f[l][r],f[l+1][r]);
		}

	while(m--)
	{
		int x=(1ll*last+read())%lim+1,y=(1ll*last+read())%lim+1;
		int l=min(x,y),r=max(x,y);
		printf("%d\n",last=query(l,r+1));
	}
	return 0;
}