51Nod1577 异或凑数 线性基

原文链接https://www.cnblogs.com/zhouzhendong/p/51Nod1577.html

题意

  给定一个长度为 n 的序列。

  有 m 组询问,每一组询问给出 L,R,k ,询问 L,R 区间内是否能找出一些数,使它们 XOR 起来等于 k 。

  $n,m\leq 5\times 10^5,  0\leq a_i,k< 2^{30}$

题解

  由于 $n,m$ 同阶,所以以下时间复杂度描述时,对于 $n,m$ 不加区分。

  线性基合并是 $O(\log ^2 a_i)$ 的。

  直接线段树维护区间线性基或者 ST 表复杂度均为 $O(n\log ^3 a_i)$ 。

  CDQ分治时间复杂度为 $O(n\log ^2 a_i)$ 。

  以上算法均不能通过。

  考虑将询问离线,按照 R 从小到大排序。

  我们将 $a_i$ 从左到右依次加入。利用线性基维护尽量靠右的基向量即可(经典套路)。

  时间复杂度为 $O(n\log a_i)$ 。

代码

#include <bits/stdc++.h>
using namespace std;
int read(){
	int x=0;
	char ch=getchar();
	while (!isdigit(ch))
		ch=getchar();
	while (isdigit(ch))
		x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
	return x;
}
const int N=500005;
int n,m;
int a[N];
struct xxj{
	int v[30],p[30];
	void clear(){
		memset(v,0,sizeof v);
		memset(p,0,sizeof p);
	}
	void insert(int x,int y){
		for (int i=29;i>=0;i--)
			if (~x>>i&1)
				continue;
			else if (!v[i]){
				v[i]=x,p[i]=y;
				break;
			}
			else {
				if (y>p[i])
					swap(y,p[i]),swap(x,v[i]);
				x^=v[i];
			}
	}
	int query(int x,int y){
		for (int i=29;i>=0;i--)
			if (x>>i&1)
				if (!v[i]||p[i]<y)
					return 0;
				else
					x^=v[i];
		return 1;
	}
}xianxingji;
struct Query{
	int L,R,k,id,ans;
}q[N];
bool cmpR(Query a,Query b){
	return a.R<b.R;
}
bool cmpid(Query a,Query b){
	return a.id<b.id;
}
int main(){
	n=read();
	for (int i=1;i<=n;i++)
		a[i]=read();
	m=read();
	for (int i=1;i<=m;i++){
		q[i].L=read();
		q[i].R=read();
		q[i].k=read();
		q[i].id=i;
	}
	sort(q+1,q+m+1,cmpR);
	xianxingji.clear();
	for (int i=1,j=0;i<=m;i++){
		while (j<q[i].R)
			j++,xianxingji.insert(a[j],j);
		q[i].ans=xianxingji.query(q[i].k,q[i].L);
	}
	sort(q+1,q+m+1,cmpid);
	for (int i=1;i<=m;i++)
		puts(q[i].ans?"YES":"NO");
	return 0;
}

  

posted @ 2018-10-31 19:57  zzd233  阅读(361)  评论(0编辑  收藏  举报