bzoj3261 最大异或和

$题目传送门bzoj3261$

---

题目

$3261: 最大异或和$

$Time Limit: 10 Sec$
$Memory Limit: 512 MB$

$Description$

给定一个非负整数序列${a}$，初始长度为$N$。
有$M$个操作，有以下两种操作类型：
$1$、$A$ $x$：添加操作，表示在序列末尾添加一个数$x$，序列的长度$N+1$。
$2$、$Q$ $l$ $r$ $x$：询问操作，你需要找到一个位置$p$，满足$l\leq p\leq 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

$Sample Output$

4
5
6

对于测试点 $1-2，N，M\leq 5$。
对于测试点 $3-7，N，M\leq 80000$。
对于测试点 $8-10，N，M\leq 300000$。
其中测试点 $1，3，5，7，9$ 保证没有修改操作。

$0\leq a[i]\leq 10^7$。

---

题解

看到$xor$，本能反应就是按位$xor$。
首先用最裸的做法，用前缀和处理一下。
设$s[i]$ $=$ $a[1]$ $xor$ $a[2]$ $xor$ $...$ $xor$ $a[i]$
则$a[p]$ $xor$ $a[p+1]$ $xor$ $...$ $xor$ $a[N]$ $xor$ $x$
$=$ $s[p-1]$ $xor$ $s[N]$ $xor$ $x$
添加$O(1)$，询问$O(N)$，时间$O(MN)$，空间$O(N)$。
观察$N$的取值范围，可以想到用$Trie$。
构造一棵可持久化$Trie$，询问时询问$s[l-1]$到$s[r-1]$与$s[N]$ $xor$ $x$。
因为高位的$1$比低位的$1$加起来还要大，所以高位的$1$先取。
添加$O(logN)$，询问$O(logN)$，时间$O(MlogN)$，空间$O(MlogN)$。

---

总结

看到异或，本能反应就是位运算的方法。
例如：按位$xor$，数位$DP$，$Trie$，可持久化$Trie$……
这样的关于异或的题目，一般情况下，都可以用$O(logN)$的时间解决……
简单地说，就是先想$Trie$。做不了，那我也没办法了……

---

标程

#include <bits/stdc++.h>
using namespace std;
const int N = 6e5 + 50;
struct tree{
    int c[2], cnt;
}t[N * 25];
int top, s[25], a[N], b[N], r[N];
int update(int y, int z)
{
    int ret = ++top, x = ret;
    for (int i = 23; i >= 0; i--)
    {
        t[x] = t[y]; t[x].cnt++;
        bool tmp = (z >> i) & 1;
        y = t[y].c[tmp];
        x = t[x].c[tmp] = ++top;
    }
    t[x] = t[y]; t[x].cnt++;
    return ret;
}
int query(int x, int y, int z)
{
    int ret = 0;
    for (int i = 23; i >= 0; i--)
    {
        bool tmp = (z >> i) & 1;
        if (t[t[x].c[tmp ^ 1]].cnt != t[t[y].c[tmp ^ 1]].cnt)
            {ret += s[i]; x = t[x].c[tmp ^ 1]; y = t[y].c[tmp ^ 1];}
        else {x = t[x].c[tmp]; y = t[y].c[tmp];}
    }
    return ret;
}
int main()
{
    int n, m; s[0] = 1;
    for (int i = 1; i <= 23; i++)
        s[i] = s[i - 1] << 1;
    scanf("%d%d", &n, &m); n++;
    r[1] = update(r[0], 0);
    for (int i = 2; i <= n; i++)
    {
        scanf("%d", &a[i]);
        b[i] = b[i - 1] ^ a[i];
        r[i] = update(r[i - 1], b[i]);
    }
    while (m--)
    {
        char c[10]; scanf("%s", c);
        if (c[0] == 'A')
        {
            scanf("%d", &a[++n]);
            b[n] = b[n - 1] ^ a[n];
            r[n] = update(r[n - 1], b[n]);
        }
        else
        {
            int x, y, z; scanf("%d%d%d", &x, &y, &z);
            printf("%d\n", query(r[x - 1], r[y], b[n] ^ z));
        }
    }
    return 0;
}

---