[BZOJ2821]作诗(Poetize)

[BZOJ2821]作诗(Poetize)

题目大意:

一个长度为\(n(n\le10^5)\)的数列,\(m(m\le10^5)\)次询问,每次询问区间\([l,r]\)内出现次数为正偶数的数的个数。强制在线。

思路:

分块,\(f[i][j]\)表示第\(i\)块到第\(j\)块的答案,\(s[i][j]\)表示\(j\)在前\(i\)块出现次数。

时间复杂度\(\mathcal O(n\sqrt n)\)

源代码:

#include<cstdio>
#include<cctype>
#include<algorithm>
inline int getint() {
	register char ch;
	while(!isdigit(ch=getchar()));
	register int x=ch^'0';
	while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
	return x;
}
const int N=1e5+1,B=317;
int n,c,m,a[N],bel[N],f[B][B],s[B][N],cnt[N];
inline void init() {
	for(register int i=1;i<=n;i++) {
		bel[i]=i/B+1;
		s[bel[i]][a[i]]++;
	}
	for(register int i=1;i<=bel[n];i++) {
		for(register int j=1;j<=c;j++) {
			s[i][j]+=s[i-1][j];
		}
	}
	for(register int i=1;i<=bel[n];i++) {
		for(register int j=i;j<=bel[n];j++) {
			f[i][j]=f[i][j-1];
			for(register int k=std::max(B*(j-1),1);k<B*j&&k<=n;k++) {
				cnt[a[k]]++;
				f[i][j]+=cnt[a[k]]>0&&cnt[a[k]]%2==0;
				f[i][j]-=cnt[a[k]]>1&&cnt[a[k]]%2==1;
			}
		}
		for(register int k=std::max(B*(i-1),1);k<=n;k++) {
			cnt[a[k]]--;
		}
	}
}
inline int query(const int &l,const int &r) {
	if(bel[l]==bel[r]) {
		int ret=0;
		for(register int i=l;i<=r;i++) {
			cnt[a[i]]++;
			ret+=cnt[a[i]]>0&&cnt[a[i]]%2==0;
			ret-=cnt[a[i]]>1&&cnt[a[i]]%2==1;
		}
		for(register int i=l;i<=r;i++) {
			cnt[a[i]]--;
		}
		return ret;
	}
	int ret=f[bel[l]+1][bel[r]-1];
	for(register int i=l;bel[i]==bel[l];i++) {
		cnt[a[i]]=s[bel[r]-1][a[i]]-s[bel[l]][a[i]];
	}
	for(register int i=r;bel[i]==bel[r];i--) {
		cnt[a[i]]=s[bel[r]-1][a[i]]-s[bel[l]][a[i]];
	}
	for(register int i=l;bel[i]==bel[l];i++) {
		cnt[a[i]]++;
		ret+=cnt[a[i]]>0&&cnt[a[i]]%2==0;
		ret-=cnt[a[i]]>1&&cnt[a[i]]%2==1;
	}
	for(register int i=r;bel[i]==bel[r];i--) {
		cnt[a[i]]++;
		ret+=cnt[a[i]]>0&&cnt[a[i]]%2==0;
		ret-=cnt[a[i]]>1&&cnt[a[i]]%2==1;
	}
	for(register int i=l;bel[i]==bel[l];i++) {
		cnt[a[i]]=0;
	}
	for(register int i=r;bel[i]==bel[r];i--) {
		cnt[a[i]]=0;
	}
	return ret;
}
int main() {
	n=getint(),c=getint(),m=getint();
	for(register int i=1;i<=n;i++) {
		a[i]=getint();
	}
	init();
	for(register int i=0,ans=0;i<m;i++) {
		int l=(getint()+ans)%n+1,r=(getint()+ans)%n+1;
		if(l>r) std::swap(l,r);
		printf("%d\n",ans=query(l,r));
	}
	return 0;
}
posted @ 2019-01-07 16:31  skylee03  阅读(178)  评论(0编辑  收藏  举报