构建可持久化Trie所需的空间和时间复杂度都是字符串的总长度的线性函数。
比例系数为字符集大小。

例题：AcWing 256. 最大异或和：
给定一个非负整数序列 a，初始长度为 N。

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

1、”A x”：添加操作，表示在序列末尾添加一个数 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 行，每行描述一个操作，格式如题面所述。

输出格式
每个询问操作输出一个整数，表示询问的答案。

每个答案占一行。

数据范围
N,M≤3∗105,0≤a[i]≤107。

输入样例：
5 5
2 6 4 3 6
A 1
Q 3 5 4
A 4
Q 5 7 0
Q 3 6 6
输出样例：
4
5
6

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<set>
#include<deque>
#include<map>
#include<vector>
#include<cmath>
#define ll long long
#define llu unsigned ll
using namespace std;
const ll lnf=0x3f3f3f3f3f3f3f3f;
const int inf = 0x3f3f3f3f;
const int maxn=601000;
const int mod=1e9+7;
int trie[maxn*24][2],latest[maxn*24];
int s[maxn],root[maxn],n,m,tot;
//递归方便统计 latest

void _insert(int i,int k,int p,int q)
{
    if(k<0)
    {
        latest[q]=i;
        return ;
    }
    int c=(s[i]>>k)&1;
    if(p) trie[q][c^1]=trie[p][c^1];
    trie[q][c]=++tot;
    _insert(i,k-1,trie[p][c],trie[q][c]);
    latest[q]=max(latest[trie[q][0]],latest[trie[q][1]]);
}


int ask(int now,int val,int k,int limit)
{
    if(k<0) return s[latest[now]]^val;
    int c=(val>>k)&1;
    if(latest[trie[now][c^1]]>=limit)
        return ask(trie[now][c^1],val,k-1,limit);
    else return ask(trie[now][c],val,k-1,limit);
}

int main(void)
{
    char pos[20];
    int x,y,z;
    scanf("%d%d",&n,&m);
    latest[0]=-1;
    root[0]=++tot;
    _insert(0,23,0,root[0]);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&z);
        s[i]=s[i-1]^z;
        root[i]=++tot;
        _insert(i,23,root[i-1],root[i]);
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%s",pos);
        if(pos[0]=='A')
        {
            scanf("%d",&z);
            n++;
            s[n]=s[n-1]^z;
            root[n]=++tot;
            _insert(n,23,root[n-1],root[n]);
        }
        else
        {
            scanf("%d%d%d",&x,&y,&z);
            printf("%d\n",ask(root[y-1],s[n]^z,23,x-1));
        }
    }
    return 0;
}