「CQOI 2018」异或序列

forever_dreams

于 2019-10-19 14:44:54 发布

阅读量154

点赞数

分类专栏： --------数据结构-------- # 莫队文章标签：莫队

本文链接：https://blog.csdn.net/forever_dreams/article/details/102625473

版权

--------数据结构-------- 同时被 2 个专栏收录

149 篇文章 1 订阅

订阅专栏

莫队

5 篇文章 1 订阅

订阅专栏

传送门

problem

已知一个长度为 $n$ 的整数数列 $a_1,a_2,...,a_n$ ，给定查询参数 $l, r$ ，问在 $a_l,a_{l+1},...,a_r$ 区间内，有多少子序列满足异或和等于 $k$ 。也就是说，对于所有的 $x, y (l \leq x \leq y \leq r)$ ，能够满足 $a_x \oplus a_{x+1} \oplus \cdots \oplus a_y = k$ 的 $x, y$ 有多少组。

数据范围： $1 ≤ n, m ≤ 10^5$ ， $0 ≤ k, a_i ≤ 10^5$ 。

solution

这题其实挺简单的。

我们先求一个前缀异或和，把问题转化成求 $[l - 1, r]$ 中有多少对异或和为 $k$ 的点。

注意到 $k$ 是全局的，也即，每次询问不是完全独立的，我们也许可以从上一次询问得到一些信息来优化复杂度。

这便启示我们考虑莫队算法了。

维护一个桶（记作 $c n t$ ），插入一个数 $x$ 时，我们在 $k\oplus x$ 位置 $+ 1$ ，同时把 $cnt_x$ 加进答案里。删除反着做即可。

时间复杂度 $O(n\sqrt n)$ 。

code

#include<bits/stdc++.h>
#define N 100005
using namespace std;
int n,m,k,S,Xor[N],in[N],ans[N],cnt[N];
struct Query{int l,r,id;}Q[N];
bool operator<(const Query &p,const Query &q){
	return (in[p.l]==in[q.l])?p.r<q.r:in[p.l]<in[q.l];
}
void Del(int pos,int &now)  {now-=cnt[Xor[pos]],cnt[k^Xor[pos]]--;}
void Add(int pos,int &now)  {now+=cnt[Xor[pos]],cnt[k^Xor[pos]]++;}
int main(){
	scanf("%d%d%d",&n,&m,&k),S=sqrt(n);
	for(int i=1,x;i<=n;++i)  scanf("%d",&x),Xor[i]=Xor[i-1]^x;
	for(int i=1;i<=m;++i)  scanf("%d%d",&Q[i].l,&Q[i].r),Q[i].l--,Q[i].id=i;
	for(int i=1;i<=n;++i)  in[i]=(i-1)/S+1;
	sort(Q+1,Q+m+1);
	int L=1,R=0,now=0;
	for(int i=1;i<=m;++i){
		while(L<Q[i].l)  Del(L++,now);
		while(L>Q[i].l)  Add(--L,now);
		while(R<Q[i].r)  Add(++R,now);
		while(R>Q[i].r)  Del(R--,now);
		ans[Q[i].id]=now;
	}
	for(int i=1;i<=m;++i)  printf("%d\n",ans[i]);
	return 0;
}