【GDOI2018模拟9.21】数列

Description

Input

Output

Sample Input

5 5 0
3 1 2 9 4
1 5
4 4
1 4
1 3
2 5

Sample Output

3
4
15
11
3

Solution

这题坑死我了
设输入的数组为a[]
一个区间l,r异或，相当于将异或前缀后（即将a[]前缀异或）a[r]^a[l]的值
把a二进制分开来加入trie
我现在有一个trie，在加入a[i]的时候，可以求出此时a[i]和i之前哪一个异或得到的结果最大
但这显然无法解决这题
考虑分块
设c[i][j]表示从第i个块到第j个块包含的数的答案
那么这个区间就不是从0开始了
因为trie只是维护有没有这个数字
所以把trie变成可持久化trie，类似主席树的前缀和模式，就可以维护了
那么每次询问的时候，先把c加上去，再把两边独立出来的再算一次答案，统计一下
细节很多（调了我两个小时）

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define N 20100
#define M 150
using namespace std;
int n,m,cty,a[N],ans=0,e[32],tot,c[M+100][M+100];
struct node{
    int l,r,d;
}t[N*100]; 
void ins(int x,int v)
{
    t[v].d++;
    fd(i,29,0)
    {
        if(x&e[i])
        {
            t[++tot]=t[t[v].r],t[v].r=tot;
            v=t[v].r;t[v].d++;
        }
        else
        {
            t[++tot]=t[t[v].l],t[v].l=tot;
            v=t[v].l;t[v].d++;
        }
    }
}
int get(int v1,int v2,int x)
{
    int an=0;
    fd(i,29,0)
    {
        if(an+e[i+1]<ans) return 0;
        if(x&e[i])
        {
            if(t[t[v2].l].d-t[t[v1].l].d) v1=t[v1].l,v2=t[v2].l,an+=e[i];
            else v1=t[v1].r,v2=t[v2].r;
        }
        else
        {
            if(t[t[v2].r].d-t[t[v1].r].d) v1=t[v1].r,v2=t[v2].r,an+=e[i];
            else v1=t[v1].l,v2=t[v2].l;
        }
    }
    return an;
}
int main()
{
    freopen("sequence.in","r",stdin);
    freopen("sequence.out","w",stdout);
    int ac;scanf("%d%d%d",&n,&ac,&cty);tot=n+1;
    e[0]=1;fo(i,1,31) e[i]=e[i-1]*2;ins(0,1);t[2]=t[1];
    fo(i,1,n) 
    {
        scanf("%d",&a[i]),a[i]^=a[i-1],t[i+1]=t[i],ins(a[i],i+1);
    }
    int m=n/M;
    fo(j,0,m)
    {
        fo(i,j*M+1,n)
        {
            c[j][i/M]=max(c[j][i/M],get(j*M,i,a[i]));
        }
    }
    fo(l,1,m) fo(i,0,m-l) c[i][i+l]=max(c[i][i+l],max(c[i+1][i+l],c[i][i+l-1]));
    for(;ac;ac--)
    {
        int x,y;scanf("%d%d",&x,&y);
        int l=(x+cty*ans)%n+1,r=(y+cty*ans)%n+1;
        if(l>r) swap(l,r);
        ans=c[l/M+1][r/M-1];
        fo(i,max((r)/M*M,l-1),r)
        ans=max(ans,get(l-1,i,a[i]));
        fo(i,l-1,min((l/M+1)*M,r))
        ans=max(ans,get(i,r+1,a[i]));
        printf("%d\n",ans);
    }
}