LibreOJ 2534 「CQOI2018」异或序列

莫队呀莫队。
此题双倍经验 CF617E。 设 \(xorsum_i = a_1 \oplus a_2 \oplus \dots \oplus a_i\)
那么根据异或运算的性质,求 \(a_l \oplus a_{l+1} \oplus \dots \oplus a_r\) 就相当于 \(xorsum_r \oplus xorsum_{l-1}\)

所以问题转化为求给定的 \([l,r]\) 区间内有多少个数对 \((x,y)\),满足 \(xorsum_y \oplus xorsum_{x-1} = k\)

那么,我们可以考虑在莫队处理 \(a_i\) 对答案产生的贡献的时候,设 \(cnt_x = \forall j \in [l,i]\)\(xorsum_j = x\) 的个数。
那么 \(a_i\) 对答案的贡献就是 \(cnt_{xorsum_i \oplus k}\)

删除同理。

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int N = 1e5;
const int M = 1e5;
const int A = 131071;
int n,m,k,f[N + 10],block,pos[N + 10],cnt[A + 10],cur;
int ans[M + 10];
struct s_query
{
int l,r,id;
inline bool operator<(const s_query &a) const
{
return pos[l] < pos[a.l] || (pos[l] == pos[a.l] && r < a.r);
}
} qry[M + 10];
int main()
{
scanf("%d%d%d",&n,&m,&k);
block = pow(n,2.0 / 3);
for(register int i = 1;i <= n;++i)
scanf("%d",f + i),pos[i] = (i - 1) / block + 1,f[i] ^= f[i - 1];
cnt[0] = 1;
for(register int i = 1;i <= m;++i)
scanf("%d%d",&qry[i].l,&qry[i].r),qry[i].id = i;
sort(qry + 1,qry + m + 1);
for(register int i = 1,l = 1,r = 0;i <= m;++i)
{
while(r < qry[i].r)
cur += cnt[f[++r] ^ k],++cnt[f[r]];
while(r > qry[i].r)
--cnt[f[r]],cur -= cnt[f[r--] ^ k];
while(l < qry[i].l)
--cnt[f[l - 1]],cur -= cnt[f[l - 1] ^ k],++l;
while(l > qry[i].l)
--l,cur += cnt[f[l - 1] ^ k],++cnt[f[l - 1]];
ans[qry[i].id] = cur;
}
for(register int i = 1;i <= m;++i)
printf("%d\n",ans[i]);
}