1. 程式人生 > >BZOJ 3261 最大異或和(可持久化Trie)

BZOJ 3261 最大異或和(可持久化Trie)

找到 比較 != oid NPU esp sca void con

Description

給定一個非負整數序列{a},初始長度為N。 有M個操作,有以下兩種操作類型: 1、Ax:添加操作,表示在序列末尾添加一個數x,序列的長度N+1。 2、Qlrx:詢問操作,你需要找到一個位置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保證沒有修改操作。
0<=a[i]<=10^7。

Sample Output

4
5
6
解題思路:   假如說就問你一個數在一群數中的最大異或,你就想到了01字典樹 (詳見HDU4825) 那要是這道題不就是把一個序列的後綴異或和做為一個數,問你區間最大異或 最大異或好辦,就是在字典樹上貪心查找,而區間異或,就是在一個只存了l~r這個區間的字典樹上貪心
這就是可持久化字典樹了(不要問我這是什麽,詳解網上遍地都是) 每個節點建一個新版本,詢問r和l-1版本的差集貪心即可。 當然插入數的時候要以前綴的形式插入,詢問時利用兩次異或值為0的性質即可。 可持久化數據結構耗費很大空間
 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 using namespace std;
 5 struct pnt{
 6     int child[2];
 7     int num;
 8 };
 9 struct KCJTrie{
10     pnt tr[17000000
]; 11 int siz; 12 void insert(int &spc,int rt,int x) 13 { 14 spc=++siz; 15 int root=spc; 16 for(int i=30;~i;i--) 17 { 18 int tmp=((x&(1<<i))!=0); 19 tr[root]=tr[rt]; 20 tr[root].num++; 21 rt=tr[rt].child[tmp]; 22 tr[root].child[tmp]=++siz; 23 root=tr[root].child[tmp]; 24 } 25 tr[root].num=tr[rt].num+1; 26 return ; 27 } 28 int query(int spc1,int spc2,int v) 29 { 30 int ans=0; 31 int rt1=spc1,rt2=spc2; 32 for(int i=30;~i;i--) 33 { 34 int tmp=((v&(1<<i))!=0); 35 if(tr[tr[rt1].child[tmp^1]].num-tr[tr[rt2].child[tmp^1]].num) 36 { 37 ans|=(1<<i); 38 rt1=tr[rt1].child[tmp^1]; 39 rt2=tr[rt2].child[tmp^1]; 40 }else{ 41 rt1=tr[rt1].child[tmp]; 42 rt2=tr[rt2].child[tmp]; 43 } 44 } 45 return ans; 46 } 47 }T; 48 int a[1000000]; 49 int b[1000000]; 50 int s[1000000]; 51 int n,m; 52 char cmd[5]; 53 int main() 54 { 55 scanf("%d%d",&n,&m); 56 n++;//看hzwer學長的博客,的確在第一位放0比較好處理,就不用特判l==1的情況了 57 for(int i=2;i<=n;i++) 58 scanf("%d",&a[i]); 59 for(int i=1;i<=n;i++) 60 { 61 b[i]=b[i-1]^a[i]; 62 T.insert(s[i],s[i-1],b[i]); 63 } 64 for(int i=1;i<=m;i++) 65 { 66 scanf("%s",cmd); 67 if(cmd[0]==A) 68 { 69 n++; 70 scanf("%d",&a[n]); 71 b[n]=b[n-1]^a[n]; 72 T.insert(s[n],s[n-1],b[n]); 73 }else{ 74 int l,r,x; 75 scanf("%d%d%d",&l,&r,&x); 76 printf("%d\n",T.query(s[r],s[l-1],x^b[n])); 77 } 78 } 79 return 0; 80 }

BZOJ 3261 最大異或和(可持久化Trie)