BZOJ 3261 最大异或和(算竞进阶习题)

可持久化Trie

需要知道一个异或的特点,和前缀和差不多
a[p] xor a[p+1] xor....xor a[n] xor x = a[p-1] xor a[n] xor x

所以我们把a[1...n]的异或和预处理出来,用s[i]表示,用一个可持久化Trie维护
问题就转化成s[n] xor x = val,求一个序列与val异或的最大值

第i个Trie的root对应维护s[1..i],这样我们在查询值的时候为了保证在[...r-1]之类,只要查询r-1及之前的版本
为了保证在[l-1...]内,我们还需要一个数组latest维护末尾节点是最近一次第几个s的末尾,查询的时候跳过小于l-1的版本即可

这题洛谷上卡常。。根据题目给的范围,24位二进制足够了,没必要开32位,节约时间。。还有各种inline,快读,尽量用上
为了保证能够查询[1..r],我们还需要给空树插入一个0,保证覆盖整个区间(类似主席树)

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
inline int lowbit(int x){ return x & (-x); }
inline int read(){
    int X = 0, w = 0; char ch = 0;
    while(!isdigit(ch)) { w |= ch == '-'; ch = getchar(); }
    while(isdigit(ch)) X = (X << 3) + (X << 1) + (ch ^ 48), ch = getchar();
    return w ? -X : X;
}
inline int gcd(int a, int b){ return a % b ? gcd(b, a % b) : b; }
inline int lcm(int a, int b){ return a / gcd(a, b) * b; }
template<typename T>
inline T max(T x, T y, T z){ return max(max(x, y), z); }
template<typename T>
inline T min(T x, T y, T z){ return min(min(x, y), z); }
template<typename A, typename B, typename C>
inline A fpow(A x, B p, C lyd){
    A ans = 1;
    for(; p; p >>= 1, x = 1LL * x * x % lyd)if(p & 1)ans = 1LL * x * ans % lyd;
    return ans;
}
const int N = 600005;
int n, m, trie[N*24][2], latest[N*24], s[N], root[N], tot;

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

inline int query(int cur, int i, int val, int lmt){
    if(i < 0) return s[latest[cur]] ^ val;
    int t = (val >> i) & 1;
    if(latest[trie[cur][t^1]] >= lmt) return query(trie[cur][t^1], i - 1, val, lmt);
    return query(trie[cur][t], i - 1, val, lmt);
}

int main(){
    memset(latest, -1, sizeof latest);
    n = read(), m = read();
    root[0] = ++tot;
    insert(0, 23, root[0], 0);
    for(int i = 1; i <= n; i ++){
        s[i] = s[i - 1] ^ read();
        root[i] = ++tot;
        insert(i, 23, root[i], root[i - 1]);
    }
    while(m --){
        char opt[2]; scanf("%s", opt);
        if(opt[0] == 'A'){
            int x = read();
            root[++n] = ++tot, s[n] = s[n - 1] ^ x;
            insert(n, 23, root[n], root[n - 1]);
        }
        else if(opt[0] == 'Q'){
            int l = read(), r = read(), x = read();
            printf("%d
", query(root[r - 1], 23, s[n] ^ x, l - 1));
        }
    }
    return 0;
}
原文地址:https://www.cnblogs.com/onionQAQ/p/10555320.html