问题 I: qwb VS 去污棒
时间限制: 2 Sec 内存限制: 256 MB提交: 74 解决: 26
[ 提交][ 状态][ 讨论版]
题目描述
qwb表白学姐失败后,郁郁寡欢,整天坐在太阳底下赏月。在外人看来,他每天自言自语,其实他在和自己的影子“去污棒”聊天。
去污棒和qwb互相出题考验对方,去污棒问了qwb这样一个问题:
现已知一个有n个正整数的序列a[1],a[2]...a[n],接下来有m个操作
suf[t]表示从t开始的后缀的异或和,即 suf[t]=a[t] xor a[t+1] xor ...xor a[len],len为序列长度。
去污棒和qwb互相出题考验对方,去污棒问了qwb这样一个问题:
现已知一个有n个正整数的序列a[1],a[2]...a[n],接下来有m个操作
操作一共有两种:
1.在序列末尾添加一个数x。
2.查询suf[p] xor x的最大值,其中xor是异或 ,l<=p<=r,suf[t]表示从t开始的后缀的异或和,即 suf[t]=a[t] xor a[t+1] xor ...xor a[len],len为序列长度。
输入
第一行一个整数T(<=5),表示一共有T组数据。
每组数据第一行两个整数n(<=200000),m(<=200000),意义如上所述。
随后一行有n个数,表示初始序列。随后m行,每行表示一个操作。
操作有两种,1: x 表示在末尾添加一个x,2: l r x表示查询suf[p] xor x的最大值,其中l<= p <= r,
所有数及x不超过224 且保证 所有操作合法。
输出
每组测试数据的第一行输出"Case x:",x为数据组数的标号,从1开始。
接下来,对每个操作2输出一行答案。
样例输入
1
5 5
1 2 3 4 5
2 1 3 4
1 10
1 7
2 4 4 5
2 1 5 19
样例输出
Case 1:
6
9
31
题目链接:I题
如果用过可持久化Trie就可以发现这题解法已经非常明显了,但是他的序列是会变动的,那么由于题目是后缀异或和的性质,可以倒着做,用可持久化Trie维护最后的后缀异或值序列,然后考虑倒着做的时候每一次是减掉末尾的数x的影响,记后缀异或值序列为{后缀和},那么有:$max(x⊕({后缀和}⊕y))$与$max((x⊕y)⊕({后缀和}))$是等价的,那么我们记录一下倒着做的时候去掉的x的影响即用一个变量维护一下x的倒着的前缀和即可,然后由于是倒着维护的序列且是倒着离线,因此查询区间为$[r+1,l]$,答案也是倒着输出,感动~学的数据结构比赛的时候终于派上用场了……
代码:
#include <bits/stdc++.h>
using namespace std;
const int N = 400010;
struct Trie
{
int nxt[2];
int cnt;
};
Trie L[N * 27];
int tot, root[N];
int suf[N], arr[N];
int Ans[N];
struct info
{
int ops;
int l, r, x;
};
info P[N];
void init()
{
memset(L, 0, sizeof(L));
tot = 0;
}
void update(int &cur, int ori, int step, int n, int v)
{
cur = ++tot;
L[cur] = L[ori];
L[cur].cnt += v;
if (step < 0)
return ;
int t = (n >> step) & 1;
update(L[cur].nxt[t], L[ori].nxt[t], step - 1, n, v);
}
int Find(int S, int E, int step, int n)
{
if (step < 0)
return 0;
int t = (n >> step) & 1;
if (L[L[E].nxt[t ^ 1]].cnt - L[L[S].nxt[t ^ 1]].cnt > 0)
return (1 << step) + Find(L[S].nxt[t ^ 1], L[E].nxt[t ^ 1], step - 1, n);
else
return Find(L[S].nxt[t], L[E].nxt[t], step - 1, n);
}
int main(void)
{
//printf("%d\n",31^19);
int tcase, n, m, i;
scanf("%d", &tcase);
for (int q = 1; q <= tcase; ++q)
{
init();
scanf("%d%d", &n, &m);
for (i = 1; i <= n; ++i)
scanf("%d", &arr[i]);
for (i = 1; i <= m; ++i)
{
scanf("%d", &P[i].ops);
if (P[i].ops == 1)
{
scanf("%d", &P[i].x);
arr[++n] = P[i].x;
}
else if (P[i].ops == 2)
{
scanf("%d%d%d", &P[i].l, &P[i].r, &P[i].x);
}
}
suf[n + 1] = 0;
for (i = n; i >= 1; --i)
{
suf[i] = suf[i + 1] ^ arr[i];
update(root[i], root[i + 1], 25, suf[i], 1);
}
int last = 0;
printf("Case %d:\n", q);
int sz = 0;
for (i = m; i >= 1; --i)
{
if (P[i].ops == 1)
{
last ^= P[i].x;
}
else if (P[i].ops == 2)
{
int one = P[i].x ^ last;
Ans[++sz] = Find(root[P[i].r + 1], root[P[i].l], 25, one);
}
}
for (i = sz; i >= 1; --i)
printf("%d\n", Ans[i]);
}
}