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

题意:
给出一段序列，多次询问区间[l,r]内一段子区间所有权值异或和最大，询问该异或和是多少。
解析:
《看完题就傻系列》
这题好不好想啊..
首先最容易想到的就是预处理出所有的前缀异或和。
然后每一次相当于在[l-1,r]中挑两个数的异或和最大。
所求东西的形式容易想到可持久化trie。
不过怎么应用？这点真是卡死我= =!
要知道我们肯定不可能暴力O(m*n*log(signed int))。
所以一定有特殊的优化技巧！
丫的我承认我分块做的少=w=
的确有一种分块做法的存在，至于是否有其他方法，反正我是没想到，因为这个强制在线简直蛋疼。
首先我们将整个序列分成sqrt(n)块。然后我们可以预处理出每个块的左端点到其右边任意一个点的最大选两个异或和。
为什么这么做一会再说。
显然这是一个递推。
b[i][j]表示第i块的左端点到第j个的最大两数异或和。
显然b[i][j]=max(b[i][j-1],query(root[l-1],root[j],a[j]))。
这个复杂度最坏肯定不会超过O(sqrt(n)*n*log(signed int))
然后我们再考虑每一个询问。
对于一个询问，我们可以找到在这段区间内出现的第一个块中左端点的位置。
然后我们可以直接拿出来b[pos][r]这个值。
剩下的该左端点前的所有点到该询问的左端点的点的个数显然不超过sqrt(n).
所以对于每一个我们都可以暴力在可持久化trie上跑。
这部分的复杂度大概是O(m*sqrt(n)*log(signed int))
所以总复杂度大概卡卡时是可以接受的。

感觉这道题的大分类应该算是分块，毕竟可持久化trie只是提供了一个在O(log(signed int))复杂度内能解决求某一部分值的方案。

另外吐槽这题那个lastans可能加爆int，从而使得re。
代码:

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 12100
#define M 800100
using namespace std;
int n,m;
int Block;
int a[N],size;
int lson[M],rson[M],siz[M];
int root[N],b[110][N];
void insert(int y,int &x,int num,int dep)
{
    x=++size;
    siz[x]=siz[y]+1;
    if(dep==-1)return;
    lson[x]=lson[y],rson[x]=rson[y];
    if((num&(1<<dep)))insert(rson[y],rson[x],num,dep-1);
    else insert(lson[y],lson[x],num,dep-1);
}
int query(int y,int x,int num)
{
    int ans=0;
    for(int i=30;i>=0;i--)
    {
        int tmp=num&(1<<i);
        if(!tmp)
        {
            if(siz[rson[x]]-siz[rson[y]]!=0)
            {
                ans+=(1<<i);
                x=rson[x],y=rson[y];
            }else x=lson[x],y=lson[y];
        }else
        {
            if(siz[lson[x]]-siz[lson[y]]!=0)
            {
                ans+=(1<<i);
                x=lson[x],y=lson[y];
            }else x=rson[x],y=rson[y];
        }
    }
    return ans;
}
int main()
{
    scanf("%d%d",&n,&m);
    Block=(int)sqrt(n);
    int Size_Block=n/Block+(n%Block!=0);

    for(int i=1;i<=n;i++)scanf("%d",&a[i]),a[i]^=a[i-1];
    for(int i=1;i<=n;i++)insert(root[i-1],root[i],a[i],30);
    for(int i=1;i<=Block;i++)
    {       
        for(int j=(i-1)*Size_Block+1;j<=n;j++)
        {
            b[i][j]=max(b[i][j-1],query(root[(i-1)*Size_Block],root[j],a[j]));
            if(i==1)b[i][j]=max(b[i][j],a[j]); 
        }
    }
    int lastans=0;
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        x=(x+lastans%n)%n+1,y=(y+lastans%n)%n+1;
        if(x>y)swap(x,y);
        x--;
        int numblock=x/Block+(x%Block!=0);
        int l1=numblock*Size_Block+1;
        lastans=0;
        if(l1<=y)lastans=b[numblock+1][y];
        l1=min(l1,y);
        for(int j=x;j<l1;j++)
        {
            lastans=max(lastans,query(root[x],root[y],a[j]));
        }
        printf("%d\n",lastans);
    }
}