[Ynoi2011]D2T1

题目大意:

给定一个数列$a$,有以下几种询问:

1. 给定$x$,在序列末尾插入$x$。
2. 给定$l,r$,输出$sumlimits_{i=l}^r a_i$。
3. 给定$x$,将数列中的所有数异或$x$。
4. 将当前数列从小到大排序。

解题思路:

考虑某一时刻数列的状态,一定是前面一段有序,后面一段无序。

首先考虑第三个操作,显然用一个变量记一下全局异或的值$X$就可以了。对于第一个操作,先把$x$异或上$X$再插入。

计算贡献的话,可以按位考虑,记录一下这一位有多少个即可。

如果没有排序操作,直接记一下每一位的前缀和即可。

考虑排序操作,前面有序的我们不去动它,后面无序的,直接把它插在正确的位置即可。显然一个数只会被插入一次。用01Trie或线段树维护一下区间各个位的出现次数即可做到区间查询(可以用前缀相减来搞)。

然而操作3会改变已经排完序的数的大小关系。我们考虑一下性质。

考虑最高位,如果它异或上了1,则这个序列的后一部分(最高位为1的那些数)会整体变到前面来。而每个部分内部的相对顺序不变。然后对每一部分,再对次高位考虑。依次类推。

对应到数据结构内部,我们发现,如果某一位异或上了1,相当于这一层的节点的左右儿子交换了一下。

我们不需要真的交换,只需要在询问的时候判断一下往哪边走即可。

这样就非常简单了,代码也很好写。

时间复杂度$O(nlog^2 A)$。

C++ Code:

#include<iostream>
#include<cstring>
const int N=2e5+5;
typedef long long LL;
int n,p,a[N],pre[N][31],A,m,cnt=1,ch[N<<4][2],wg[N<<4][31],sz[N<<4];
int STD;
inline void insert(const int&X){
	int u=1;
	for(int i=30;~i;--i){
		++sz[u];
		for(int j=0;j<31;++j)
		wg[u][j]+=X>>j&1;
		const int sn=X>>i&1;
		if(!ch[u][sn])ch[u][sn]=++cnt;
		u=ch[u][sn];
	}
	++sz[u];
	for(int j=0;j<31;++j)
	wg[u][j]+=X>>j&1;
}
LL query(int k){
	LL ret=0;
	int val=0;
	int u=1;
	for(int i=30;~i;--i){
		const int _0=STD>>i&1;
		if(sz[ch[u][_0]]>=k)u=ch[u][_0],val|=_0<<i;else{
			int x=sz[ch[u][_0]];
			k-=x;
			for(int j=0;j<31;++j){
				int cnt=wg[ch[u][_0]][j];
				if(A>>j&1)cnt=x-cnt;
				ret+=(LL)cnt<<j;
			}
			u=ch[u][!_0];
			val|=!_0<<i;
		}
	}
	val^=A;
	return ret+(LL)val*k;
}
inline LL calc(int l,int r){
	const int len=r-l+1;
	LL ret=0;
	for(int i=0;i<31;++i){
		int cnt=pre[r][i]-pre[l-1][i];
		if(A>>i&1)cnt=len-cnt;
		ret+=(LL)cnt<<i;
	}
	return ret;
}
inline LL solve(int l,int r){return query(r)-query(l-1);}
int main(){
	std::ios::sync_with_stdio(0),std::cin.tie(0),std::cout.tie(0);
	std::cin>>n;p=1;
	for(int i=1;i<=n;++i){
		std::cin>>a[i];
		for(int j=0;j<31;++j)
		pre[i][j]=pre[i-1][j]+(a[i]>>j&1);
	}
	std::cin>>m;
	while(m--){
		int op;
		std::cin>>op;
		switch(op){
			case 1:{
				int x;
				std::cin>>x;
				x^=A;
				a[++n]=x;
				for(int i=0;i<31;++i)
				pre[n][i]=pre[n-1][i]+(a[n]>>i&1);
				break;
			}
			case 2:{
				int l,r;
				std::cin>>l>>r;
				LL ans;
				if(r<p)ans=solve(l,r);else
				if(l>=p)ans=calc(l,r);else
				ans=solve(l,p-1)+calc(p,r);
				std::cout<<ans<<'
';
				break;
			}
			case 3:{
				int x;
				std::cin>>x;
				A^=x;
				break;
			}
			case 4:{
				for(int i=p;i<=n;++i)
				insert(a[i]);
				STD=A;
				p=n+1;
				memset(pre[n],0,sizeof*pre);
				break;
			}
		}
	}
	return 0;
}
原文地址:https://www.cnblogs.com/Mrsrz/p/10755999.html