bzoj4184-shallot

题目

\(n\)次操作,每次向数集中插入或删除一个数,并询问此时的最大异或和。

分析

看到最大异或和就想到线性基,但是线性基是显然不支持删除的。
题目所需要的就是按时间顺序维护一个线性基。对于这一类插入,删除,整体查询的离线问题,我们可以用线段树(同理cdq分治)解决。
通过离线询问,我们可以记录每个数的插入和删除时间,用线段树维护时间线。对于一个数,我们在线段树上他插入时间和删除时间的区间中插入这个数,然后dfs线段树,在路上累计节点的信息,到达每个叶子节点的时候就计算并输出答案。这里的时间复杂度为\(O(f(n) \* n \* logn)\)
所以在这道题中,我们在线段树上在每个数出现和完全消失(所有相同的数都被删除)的时间段插入这个数,一次dfs,动态插入数到线性基中,到达叶子节点计算答案,这样我们以\(O(logn)\)的时间代价去除了删除操作。

代码

#include<cstdio>
#include<cctype>
#include<algorithm>
#include<map>
#include<list>
#include<cstring>
using namespace std;
int read() {
	int x=0,f=1;
	char c=getchar();
	for (;!isdigit(c);c=getchar()) if (c=='-') f=-1;
	for (;isdigit(c);c=getchar()) x=x*10+c-'0';
	return x*f;
}
const int maxn=5e5+10;
const int maxj=31;
int t[maxn*3];
struct lb {
	int p[maxj];
	void insert(int x) {
		for (register int j=maxj-1;~j;--j) if (x&(1<<j)) {
			if (!p[j]) {
				for (register int i=j-1;i>=0;--i) if (x&(1<<i)) x^=p[i];
				for (register int i=j+1;i<maxj;++i) if (p[i]&(1<<j)) p[i]^=x;
				p[j]=x;
				break;
			} else x^=p[j];
		}	
	}
	int get() {
		int ret=0;
		for (register int j=0;j<maxj;++j) ret^=p[j];
		return ret;
	}
};
void print(int x) {
	if (!x) {
		puts("0");
		return;
	}
	static char s[12];
	int tot=0;
	while (x) s[++tot]=x%10+'0',x/=10;
	while (tot--) putchar(s[tot+1]);
	puts("");
}
struct edge {
	int v,nxt;
} e[maxn<<4];
int tot=0;
void dfs(int now,int l,int r,lb b) {
	for (register int j=t[now];j;j=e[j].nxt) b.insert(e[j].v);
	if (l==r) {
		print(b.get());
		return;
	}
	int mid=(l+r)>>1;
	dfs(now<<1,l,mid,b);
	dfs(now<<1|1,mid+1,r,b);
}
void insert(int x,int L,int R,int l,int r,int d) {
	if (L==l && R==r) {
		e[++tot]=(edge){d,t[x]};
		t[x]=tot;
		return;
	}
	int mid=(L+R)>>1;
	if (r<=mid) insert(x<<1,L,mid,l,r,d); else 
	if (l>mid) insert(x<<1|1,mid+1,R,l,r,d); else {
		insert(x<<1,L,mid,l,mid,d);
		insert(x<<1|1,mid+1,R,mid+1,r,d);
	}
}
int a[maxn],bas[maxn],b[maxn],st[maxn];
int main() {
	#ifndef ONLINE_JUDGE
		freopen("test.in","r",stdin);
		freopen("my.out","w",stdout);
	#endif
	int n=read(),tk=0;
	for (register int i=1;i<=n;++i) {
		a[i]=read();
		if (a[i]>0) b[++tk]=a[i];
	}
	sort(b+1,b+tk+1);
	tk=unique(b+1,b+tk+1)-b-1;
	for (register int i=1;i<=n;++i) {
		if (a[i]>0) a[i]=lower_bound(b+1,b+tk+1,a[i])-b; else a[i]=-(lower_bound(b+1,b+tk+1,-a[i])-b);
		if (a[i]>0) {
			if (!bas[a[i]]) st[a[i]]=i;
			bas[a[i]]++;
		} else {
			a[i]=-a[i];
			bas[a[i]]--;
			if (!bas[a[i]]) insert(1,1,n,st[a[i]],i-1,b[a[i]]),st[a[i]]=0;
		}
	}
	for (register int i=1;i<=tk;++i) if (st[i]) insert(1,1,n,st[i],n,b[i]);
	dfs(1,1,n,lb());
}
posted @ 2017-04-17 20:40  permui  阅读(286)  评论(0编辑  收藏  举报