BZOJ3261 最大异或和 解题报告(可持久化Trie树)

本题链接：https://www.lydsy.com/JudgeOnline/problem.php?id=3261

题目描述

给定一个非负整数序列{a}，初始长度为N。

有M个操作，有以下两种操作类型：

1、Ax：添加操作，表示在序列末尾添加一个数x，序列的长度N+1。

2、Q l r x：询问操作，你需要找到一个位置p，满足l<=p<=r，使得：

a[p] xor a[p+1] xor ... xor a[N] xor x 最大，输出最大是多少。

输入

第一行包含两个整数 N  ，M，含义如问题描述所示。   
第二行包含 N个非负整数，表示初始的序列 A 。 
接下来 M行，每行描述一个操作，格式如题面所述。

输出

假设询问操作有 T个，则输出应该有 T行，每行一个整数表示询问的答案。

样例输入

5 5
2 6 4 3 6
A 1
Q 3 5 4
A 4
Q 5 7 0
Q 3 6 6
对于测试点 1-2，N,M<=5 。
对于测试点 3-7，N,M<=80000 。
对于测试点 8-10，N,M<=300000 。
其中测试点 1, 3, 5, 7, 9保证没有修改操作。
0<=a[i]<=10^7。

样例输出

4
5
6

 

题解：

在这道题目之前，我想读者需要一个前置知识便于理解。请看我的另一篇博客：https://www.cnblogs.com/xxzh/p/9178838.html

如果学会了这个知识的话，也就是知道了如何使用一颗01 Trie树

那么考虑这道题，很显然，现在我们需要维护一颗可持久化trie树。如果不清楚静态主席树原理的，请移步我的另一篇博客：https://www.cnblogs.com/xxzh/p/9158819.html

 

正式开始对这题讲解：

对于每次操作Q，它给的算式太鬼了，实际上应该化成max( (a[n] xor x)xor a[p-1] ) p∈[l,r];

 

对每个前置建立Trie，建立原理参考主席树博客。

在建树的过程中，很重要的一点就是我们每次插入的节点其实是原来a数组中的前缀。也就是说我们维护的是前缀和的前缀（有点绕读者仔细琢磨）

 

然后注意看算式，我们每次找的其实是p-1，注意是p-1而不是p。

也就是本来我们query查找的应该是l-1到r，由于是p-1我们现在查找的是l-2到r-1

当然，对于处理p-1我们也有方法，只需要在原来的数组a之前加一个0，这样的话就可以直接算l-1到r了

 

具体的操作还是贪心的从高位到低位建立Trie，查询的时候注意减一下

下面附上代码：

#include<bits/stdc++.h>
using namespace std;

const int maxn=6e5+;
int n,m,sz;
int t[maxn<<][],sum[maxn<<],b[maxn],q[maxn];
inline int read(){
    char ch=getchar();
    int s=,f=;
    while (!(ch>=''&&ch<='')) {if (ch=='-') f=-;ch=getchar();}
    while (ch>=''&&ch<='') {s=(s<<)+(s<<)+ch-'';ch=getchar();}
    return s*f;
}
int ins(int last,int val)
{
    int res,k;
    res=k=++sz;
    for (int i=;i>=;i--)
    {
        sum[k]=sum[last]+;
        t[k][]=t[last][];t[k][]=t[last][];
        bool d=val&(<<i);
        k=t[k][d]=++sz;last=t[last][d];
    }
    sum[k]=sum[last]+;
    return res;
}
int query(int k1,int k2,int val)
{
    int res=;
    for (int i=;i>=;i--)
    {
        bool d=val&(<<i);
        if (sum[t[k2][d^]]-sum[t[k1][d^]]>){
            res|=(<<i);
            k1=t[k1][d^];
            k2=t[k2][d^];
        }
        else k1=t[k1][d],k2=t[k2][d];
    }
    return res;
}
int main()
{
    n=read()+;m=read();
    q[]=ins(q[],b[]);
    for (int i=;i<=n;i++)
    {
        b[i]=b[i-]^read();//注意插入trie的是前缀和
        q[i]=ins(q[i-],b[i]);
    }
    for (int i=;i<=m;i++)
    {
        char ch=getchar();
        while (!(ch=='A'||ch=='Q')) ch=getchar();
        if (ch=='A') {
            ++n;
            b[n]=b[n-]^read();
            q[n]=ins(q[n-],b[n]);
        }
        else {
            int l=read(),r=read(),x=read();
            printf("%d\n",query(q[l-],q[r],x^b[n]));
        }
    }
    return ;
}