6731. 【2020.06.17省选模拟】T3 拟赛

题目

---

正解

先求出前缀和。考虑分界点前缀和的取值。
可以发现它的取值只有$2^k$种。具体证明考虑前面若干段的异或和都在集合$B$中，所以它一定是$B$中若干个数异或起来。

思考一下如何判定是否有答案。
如果整个序列的异或和为$0$，就是找到任意一个在$B$中的前缀和。
如果整个序列的异或和不为$0$，假如找到一个合法前缀和等于整个序列的前缀和，那么一定有答案。因为这个位置到最后这段区间的异或和为$0$，于是这段区间就可以直接和上一段合并。
（后面值考虑异或和不为$0$的情况。因为后面可以发现能够直接用它维护的东西，求出异或和为$0$的答案。）
于是我们要找到这些合法的前缀和，只要有一个等于整个序列的异或和就完事了。
对于每种取值，我们求出它第一次出现的位置。这样意味着往后转移能有更多的选择（强调：我们现在的目的是尽可能找到所有合法的前缀和）。
“转移”：枚举$B$集合中的一个数，异或一下求出一个新的前缀和，试着在后面找。

做一个优先队列BFS，先处理位置小的。试着去转移到另一个前缀和，如果转移到了，加入队列，并且标记；如果转移不到，还是要标记。因为现在位置是最小的，所以现在找不到，以后也找不到。
于是现在的问题是如何在某个位置后面快速地找到一个数（前缀和）且使其尽量靠前。并且支持区间异或一个数。
这个东西显然就是个大分块。每个块维护bitset用以查询，整块打标记散块暴力，这样就可以在$O(Q{2}^k\sqrt{n})$的时间复杂度内解决了。

还有什么用哈希和vector上二分，或双关键字基数排序之后二分的方法，看起来感觉可能快一些，但实际并没有什么卵用，反而常数大会拖慢程序速度。因为总体复杂度还是$O(Q{2}^k\sqrt{n})$，只是局部找数字的时候快一点。然而预处理这些东西的时间远远超过它们能够优化的时间……
看似优美的想法使我浪费两个小时卡常数

---

代码

using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <cassert>
#include <ctime>
#define N 100010
#define B 1000
#define M (N/B+5)
int input(){
	char ch=getchar();
	while (ch<'0' || ch>'9')
		ch=getchar();
	int x=0;
	do{
		x=x*10+ch-'0';
		ch=getchar();
	}
	while ('0'<=ch && ch<='9');
	return x;
}
int n,m;
int a[N],ps[N],all;
int k,b[5];
int bel[N],L[N],R[N];
bitset<1<<20> s[M];
int tag[M];
void build(int x){
	bitset<1<<20> *sx=&s[x];
	sx->reset();
	for (int i=L[x];i<=R[x];++i)
		sx->set(ps[i]);
}
void change(int l,int c){
	int x=bel[l];
	for (int i=l;i<=R[x];++i)
		ps[i]^=c;
	build(x);
	for (int i=bel[l]+1;i<=m;++i)
		tag[i]^=c;
}
int find(int x,int y,int v){
	v^=tag[x];
	for (int i=y;i<=R[x];++i)
		if (ps[i]==v)
			return i;
	return 0;
}
int query(int l,int v){
	if (s[bel[l]][v^tag[bel[l]]]){
		int tmp=find(bel[l],l,v);
		if (tmp)
			return tmp;
	}
	for (int i=bel[l]+1;i<=m;++i)
		if (s[i][v^tag[i]])
			return find(i,L[i],v);
	return 0;
}
bool vis[1<<20];
int h[100],nh;
bool cmph(int son,int fa){return son>fa;}
int q[100];
int BFS(){
	int cnt=0;
	nh=1;
	h[0]=0;
	vis[0]=1;
	q[++cnt]=0;
	while (nh){
		int x=h[0];
		pop_heap(h,h+nh--);
		for (int i=0;i<k;++i){
			int v=(ps[x]^tag[bel[x]])^b[i];
			if (vis[v]==0){
				vis[v]=1;
				q[++cnt]=v;
				int tmp=query(x+1,v);
				if (tmp){
					if (v==all){
						for (int j=1;j<=cnt;++j)
							vis[q[j]]=0;
						return 1;
					}
					h[nh++]=tmp;
					push_heap(h,h+nh,cmph);
				}
			}
		}
	}
	for (int i=1;i<=cnt;++i)
		vis[q[i]]=0;
	return 0;
}
int main(){
	freopen("contest.in","r",stdin);
	freopen("contest.out","w",stdout);
	n=input();
	int Q=input();
	for (int i=1;i<=n;++i){
		a[i]=input();
		ps[i]=a[i]^ps[i-1];
	}
	all=ps[n];
	for (int i=1;i<=n;i+=B){
		++m;
		L[m]=i;
		for (int j=i;j<i+B && j<=n;++j)
			bel[j]=m,R[m]=j;
	}
	for (int i=1;i<=m;++i)
		build(i);
	while (Q--){
		int op=input();
		if (op==1){
			int x=input(),c=input();
			c=a[x]^c;
			a[x]^=c;
			all^=c;
			change(x,c);
		}
		else{
			k=input();
			for (int i=0;i<k;++i)
				b[i]=input();
			sort(b,b+k);
			k=unique(b,b+k)-b;
			if (all){
				if (BFS())
					printf("yes\n");
				else
					printf("no\n");
			}
			else{
				int ans=0;
				for (int i=0;i<k && ans==0;++i)
					ans|=query(1,b[i]);
				if (ans)
					printf("yes\n");
				else
					printf("no\n");
			}
		}
	}
//	printf("time=%d %d\n",tim,clock());
	return 0;
}