最大异或和

https://zybuluo.com/ysner/note/1238161

题面

给定一个初始长度为(N)的非负整数序列({a})。
有(m)个操作，操作分为两种：

• 在序列末尾加一个数，(++N)

• 给出(l,r,x)，找出满足(lleq pleq r)的位置(p),最大化(a_pbigoplus a_{p+1}bigoplus a_{p+2}bigoplus...bigoplus a_{N}bigoplus x)。

• (nleq3*10^5,a_ileq10^7)

  解析

  维护区间最大异或和当然要召唤可持久化(Trie)树啦。
  它的功能是，在完成建树后，能在给定一段区间和一个数的情况下得出它们异或得到的最大值。

具体建立方法有点像把(Trie)树和主席树结合起来。
插入一个数（新建一棵树），枚举到某一位时，若这个数这一位是(p)，则只用递归处理(p)子树，(pbigoplus1)子树可以使用前面已有的对应子树。
对于询问，递归时若(x)这一位为(p)，则把(r)树和(l-1)树对应子树(能贡献答案的)大小作差，若差为(0)，说明区间中不存在该位为(pbigoplus1)的数，答案该位为(0)。依此类推。
（其实看代码最清楚）

至于这个奇怪的询问方式，设异或前缀和为(S)，则询问可视为(S_Nbigoplus S_{p-1})。

注意数组(rt,sum,t)等数组大小要开到(n*50)左右。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#define re register
#define il inline
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,b) for(re int i=a;i>=b;i--)
using namespace std;
const int mod=1e9+7,N=3e7+100;
int n,m,tot,sum[N],t[N][2],rt[N],tim;
il ll gi()
{
   re ll x=0,t=1;
   re char ch=getchar();
   while(ch!=‘-‘&&(ch<‘0‘||ch>‘9‘)) ch=getchar();
   if(ch==‘-‘) t=-1,ch=getchar();
   while(ch>=‘0‘&&ch<=‘9‘) x=x*10+ch-48,ch=getchar();
   return x*t;
}
il void Build(re int x,re int &y,re int w,re int d)
{
  sum[y=++tim]=sum[x]+1;
  if(d<0) return;
  re int tmp=(w>>d)&1;
  t[y][tmp^1]=t[x][tmp^1];
  Build(t[x][tmp],t[y][tmp],w,d-1);
}
il int Query(re int x,re int y,re int d)
{
  if(d<0) return 0;
  re int tmp=(w>>d)&1,p=sum[t[y][tmp^1]]-sum[t[x][tmp^1]];
  if(p>0) return (1<<d)+Query(t[x][tmp^1],t[y][tmp^1],d-1);
  else return Query(t[x][tmp],d-1);
}
int main()
{
  n=gi();m=gi();
  Build(rt[0],rt[1],25);++n;
  fp(i,2,n)
    {
      re int x=gi();tot^=x;
      Build(rt[i-1],rt[i],25);
    }
  while(m--)
    {
      re char s[5];scanf("%s",s);
      if(s[0]==‘A‘)
    {
      re int x=gi();tot^=x;
      Build(rt[n],rt[n+1],25);++n;
    }
      else
    {
      re int l=gi(),r=gi(),x=gi();
      printf("%dn",Query(rt[l-1],rt[r],tot^x,25));
    }
    }
  return 0;
}