bzoj 3261 （可持久化trie树）

3261: 最大异或和

Time Limit: 10 Sec   Memory Limit: 512 MB
Submit: 1422   Solved: 612
[ Submit][ Status][ Discuss]

Description

     

给定一个非负整数序列 {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 最大，输出最大是多少。  

Input

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

Output

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

Sample Input

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保证没有修改操作。
对于 100% 的数据， 0<=a[i]<=10^7。

Sample Output

4
5
6

HINT

对于      100%  的数据，     0<=a[i]<=10^7  。

Source

[ Submit][ Status][ Discuss]

题解：可持久化trie树

这个题只要转化一下就是可持久化trie树的裸题。

我们需要维护的是一段后缀的值，但是我们是不断加入新点的，所以直接维护后缀不好做。

因为前缀是不变的所以我们维护前缀，因为要求[l,r]的异或值，可以通过求[1,r]的异或前缀和^[1,l-1]的异或前缀和得到。

所以我们对于每个位置的前缀和建立trie树，因为要^x最大，所以我们可以先得到[1,n]异或和^x的01串， 然后在要求区间的trie树走相反的边，即01串中为1，如果连向0的方向有节点就向0的方向走。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define N 20000003
#define bit 24
using namespace std;
int sz,n,m,a[N],root[N];
int size[N],ch[N][3];
void insert(int i,int x)
{
	int pre=root[i-1];
	root[i]=++sz;
	int now=root[i];
	for (int i=bit-1;i>=0;i--){
		int t=(x>>i)&1;
		size[now]=size[pre]+1;
		ch[now][t^1]=ch[pre][t^1];
		ch[now][t]=++sz;
		now=ch[now][t];
		pre=ch[pre][t];
	}
	size[now]=size[pre]+1;
}
int find(int i,int j,int x)
{
	int ans=0;
	int l=root[i]; int r=root[j];
	for (int i=bit-1;i>=0;i--) {
		int t=((x>>i)&1)^1;
		if (size[ch[r][t]]-size[ch[l][t]]==0) t^=1;
	    else ans+=1<<i;
	    l=ch[l][t];
	    r=ch[r][t];
	}
	return ans;
}
int main()
{
	freopen("a.in","r",stdin);
	freopen("my.out","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++) scanf("%d",&a[i]);
	insert(1,0);
	for (int i=1;i<=n;i++){
	 a[i]^=a[i-1]; 
	 insert(i+1,a[i]);
    }
    for (int i=1;i<=m;i++){
    	char s[10]; int x,y,v;
    	scanf("%s%d",s,&x);
    	if (s[0]=='A') {
		 a[++n]=a[n-1]^x;
		 insert(n+1,a[n]);
	    }
    	else {
    		scanf("%d%d",&y,&v);
			printf("%d\n",find(x-1,y,v^a[n]));
		}
	}
}