BZOJ_3744_Gty的妹子序列

BZOJ3744: Gty的妹子序列

https://lydsy.com/JudgeOnline/problem.php?id=3744

分析:

• 预处理出来每一块块首到所有位置的逆序对数。
• 查询时主席树上查即可。
• 卡常

代码：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <cmath>
using namespace std;
#define N 60050
typedef int ll;
int n,a[N],m;
int blo,size,L[N],R[N],pos[N];
int siz[N*23],ls[N*23],rs[N*23],cnt,V[N],ln,root[N];
ll ansblo[305][N];
void update(int l,int r,int x,int &p,int q) {
	p=++cnt; siz[p]=siz[q]+1; ls[p]=ls[q]; rs[p]=rs[q];
	if(l==r) return ;
	int mid=(l+r)>>1;
	if(x<=mid) update(l,mid,x,ls[p],ls[q]);
	else update(mid+1,r,x,rs[p],rs[q]);
}
int qmnx(int l,int r,int x,int p,int q) {
	if(l==r) return l<=x?siz[q]-siz[p]:0;
	int mid=(l+r)>>1;
	if(x<=mid) return qmnx(l,mid,x,ls[p],ls[q]);
	else return qmnx(mid+1,r,x,rs[p],rs[q])+siz[ls[q]]-siz[ls[p]];
}
int qmxx(int l,int r,int x,int p,int q) {
	if(l==r) return l>=x?siz[q]-siz[p]:0;
	int mid=(l+r)>>1;
	if(x<=mid) return qmxx(l,mid,x,ls[p],ls[q])+siz[rs[q]]-siz[rs[p]];
	else return qmxx(mid+1,r,x,rs[p],rs[q]);
}
ll work(int x,int y) {
	int p=pos[x],q=pos[y],i;
	ll re=0;
	if(p==q) {
		for(i=x+1;i<=y;i++) {
			re+=qmxx(1,ln,a[i]+1,root[x-1],root[i-1]);
		}
		return re;
	}
	re=ansblo[p+1][y];
	for(i=R[p];i>=x;i--) {
		re+=qmnx(1,ln,a[i]-1,root[i],root[y]);
	}
	return re;
}
int s[N];
void fix(int x){for(;x<=ln;x+=x&-x)s[x]++;}
int find(int x){int ret=0;for(;x;x-=x&-x)ret+=s[x];return ret;}
int main() {
	///////////////////////////////////////////////////////////////////////////////////////////////
	///*////////////////*/freopen("noon.in","r",stdin); freopen("noon.out","w",stdout);//////////////
	///////////////////////////////////////////////////////////////////////////////////////////////
	scanf("%d",&n);
	int i,j,x,y,k;
	size=sqrt(n);
	blo=(n+size-1)/size;
	for(i=1;i<=n;i++) scanf("%d",&a[i]),V[i]=a[i];
	sort(V+1,V+n+1);
	ln=unique(V+1,V+n+1)-V-1;
	for(i=1;i<=n;i++) a[i]=lower_bound(V+1,V+ln+1,a[i])-V;
	for(i=1;i<=n;i++) {
		update(1,ln,a[i],root[i],root[i-1]);
	}
	for(i=1;i<=blo;i++) {
		L[i]=R[i-1]+1; R[i]=min(i*size,n);
		for(j=L[i];j<=R[i];j++) {
			pos[j]=i;
		}
	}
	for(i=1;i<=blo;i++) {
		memset(s,0,sizeof(s));
		for(j=L[i];j<=n;j++) {
			ansblo[i][j]=ansblo[i][j-1]+j-L[i]-find(a[j]);
			fix(a[j]);
		}
	}
	scanf("%d",&m);
	ll ans=0;
	while(m--) {
		scanf("%d%d",&x,&y);
		x^=ans; y^=ans;
		ans=work(x,y);
		printf("%d\n",ans);
	}
}