bzoj 2741: 【FOTILE模拟赛】L 分块+可持久化trie

题意

给出一个序列A,每次询问一段区间[l,r]内最大的连续异或和。强制在线。
n<=12000,m<=6000,A[i]<=maxint

分析

可以先预处理出前缀异或和,这样就变成了在区间内选两个数使其异或和尽量大。这个问题显然要用到可持久化trie。
我没有往分块上面想,但只要想到分块这题就很水了。
可以预处理出mx[i,j]表示第i块的开头到j这一段的答案,那么每次询问的时候找到第一个作为块开头的位置,从这里开始往后的答案都知道了,剩下部分只要暴力即可。
注意强制在线的时候可能会炸int。

代码

#include
#include
#include
#include
#include
#include
using namespace std;

typedef long long LL;

const int N=12005;
const int M=205;

int n,m,block,bin[35],a[N],root[N],pos[N],sta[M],end[M],mx[M][N],sz;
struct tree{int l,r,s;}t[N*35];

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;
}

void ins(int &d,int p,int x,int y)
{
    d=++sz;t[d]=t[p];t[d].s++;
    if (x<0) return;
    if ((y&bin[x])) ins(t[d].r,t[p].r,x-1,y);
    else ins(t[d].l,t[p].l,x-1,y);
}

int query(int d,int p,int x,int y)
{
    if (x<0) return 0;
    if ((y&bin[x]))
        if (t[t[d].l].s-t[t[p].l].s) return query(t[d].l,t[p].l,x-1,y)+bin[x];
        else return query(t[d].r,t[p].r,x-1,y);
    else
        if (t[t[d].r].s-t[t[p].r].s) return query(t[d].r,t[p].r,x-1,y)+bin[x];
        else return query(t[d].l,t[p].l,x-1,y);
}

int solve(int l,int r)
{
    int ans=0;
    if (pos[l]==pos[r])
    {
        for (int i=l;i30,a[i]));
        return ans;
    }
    int w;
    for (int i=l;i<=r;i++)
        if (i==sta[pos[i]])
        {
            ans=mx[pos[i]][r];w=i-1;
            break;
        }
    for (int i=l;i<=w;i++) ans=max(ans,query(root[r],root[i],30,a[i]));
    return ans;
}

int main()
{
    n=read()+1;m=read();
    bin[0]=1;block=sqrt(n);
    for (int i=1;i<=30;i++) bin[i]=bin[i-1]*2;
    for (int i=2;i<=n;i++) a[i]=a[i-1]^read();
    for (int i=1;i<=n;i++)  
    {
        ins(root[i],root[i-1],30,a[i]);
        pos[i]=(i+block-1)/block;
        if (!sta[pos[i]]) sta[pos[i]]=i;
        end[pos[i]]=i;
    }
    for (int i=1;i<=pos[n];i++)
        for (int j=sta[i]+1;j<=n;j++)
            mx[i][j]=max(mx[i][j-1],query(root[j-1],root[sta[i]-1],30,a[j]));
    int lastans=0;
    while (m--)
    {
        int x=read(),y=read();
        int l=min(((LL)x+(LL)lastans)%(n-1)+1,((LL)y+(LL)lastans)%(n-1)+1);
        int r=max(((LL)x+(LL)lastans)%(n-1)+1,((LL)y+(LL)lastans)%(n-1)+1)+1;
        printf("%d\n",lastans=solve(l,r));
    }
    return 0;
}

你可能感兴趣的:(可持久化字典树,分块)