BZOJ2821(作诗)

2821: 作诗(Poetize)

Time Limit: 50 Sec   Memory Limit: 128 MB
Submit: 1004   Solved: 321
[ Submit][ Status]

Description

神犇SJY虐完HEOI之后给傻×LYD出了一题：
SHY是T国的公主，平时的一大爱好是作诗。
由于时间紧迫，SHY作完诗之后还要虐OI，于是SHY找来一篇长度为N的文章，阅读M次，每次只阅读其中连续的一段[l,r]，从这一段中选出一些汉字构成诗。因为SHY喜欢对偶，所以SHY规定最后选出的每个汉字都必须在[l,r]里出现了正偶数次。而且SHY认为选出的汉字的种类数（两个一样的汉字称为同一种）越多越好（为了拿到更多的素材！）。于是SHY请LYD安排选法。
LYD这种傻×当然不会了，于是向你请教……
问题简述：N个数，M组询问，每次问[l,r]中有多少个数出现正偶数次。

Input

输入第一行三个整数n、c以及m。表示文章字数、汉字的种类数、要选择M次。
第二行有n个整数，每个数Ai在[1, c]间，代表一个编码为Ai的汉字。
接下来m行每行两个整数l和r，设上一个询问的答案为ans(第一个询问时ans=0)，令L=(l+ans)mod n+1, R=(r+ans)mod n+1，若L>R，交换L和R，则本次询问为[L,R]。

Output

输出共m行，每行一个整数，第i个数表示SHY第i次能选出的汉字的最多种类数。

Sample Input

5 3 5
1 2 2 3 1
0 4
1 2
2 2
2 3
3 5

Sample Output

2
0
0
0
1

【题解】分块解决，也看了其他的题解，我的无需二分，且理论复杂度在O(n*sqrt(n)),但有一定的常数。

步骤：

   一、预处理出第 i 块到第 j 块的正偶数答案。

        二、预处理出前 i 块数字 p 出现的个数。

三、处理是先通过预处理，O(1)计算l,r之间整块的答案数量。然后在判断两边零碎的数，具体判断详见代码。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>

#define maxn 101000

using namespace std;

int f[maxn][318],w[318][318],twice[maxn],S,a[maxn],cc[maxn];

inline void read(int &x)
{
	x=0;char ch;
	while(ch=getchar(),ch==' ' || ch=='\r' || ch=='\n');
	while(ch!=' ' && ch!='\r' && ch!='\n')x=x*10+(int)ch-48,ch=getchar();	
}

int main()
{
	memset(f,0,sizeof(f));
	memset(w,0,sizeof(w));
	memset(twice,0,sizeof(twice));
	
	int n,c,m,ans=0,l;
	read(n);read(c);read(m);
	S=sqrt(n)+1;l=1;

	for(int i=1;i<=n;i++)
	{		
		read(a[i]);
		if(((int)(i/S)+1)==l)f[a[i]][l]++;
		else{
		     	l++;for(int j=1;j<=c;j++)f[j][l]+=f[j][l-1];
				f[a[i]][l]++;
		 }
	}
	
	for(int i=1;i<=l;i++)
	{
		memset(cc,0,sizeof(cc));int now=0;
	 	for(int j=i;j<=l;j++)
	 	{
	 	 	for(int p=(j-1)*S;p<=j*S-1;p++)
	   			if(p){
					cc[a[p]]++;
					if(cc[a[p]]>1 && (cc[a[p]]&1))now--;
					 else if(!(cc[a[p]]&1))now++;
				}
			w[i][j]=now;
		}	
	}

	for(int l,r;m;m--){
		
		read(l);read(r);
		l=(l+ans)%n+1; r=(r+ans)%n+1;
		if (l>r) swap(l,r);
	    int head=l/S+2,tail=r/S;
	    
	    //把中间整块的求出累加； 
		if(head<=tail)ans=w[head][tail];
		 else ans=0;
		 
		//小范围边界； 
		int headt=min((int)(l/S+1)*S-1,n);if(head==tail+2)headt=r;
		int tailh=tail*S; if(head==tail+2)tailh=l;
		if(tailh==0)tailh=1;		

		//将数统计； 
		for(int j=l;j<=headt;j++)twice[a[j]]++;
		if(tailh>headt)//考虑同一块的情况； 
		  for(int j=tailh;j<=r;j++)twice[a[j]]++;
		
		for(int j=l;j<=headt;j++)
		 if(twice[a[j]]){
				int t=0;
				if((head-1)<=tail)t=f[a[j]][tail]-f[a[j]][head-1];
				if(t){
					if((t&1) && (!((t+twice[a[j]])&1)))ans++;
					 else if((!(t&1)) && ((t+twice[a[j]])&1))ans--;
				}else if(!(twice[a[j]]&1))ans++;
			twice[a[j]]=0;
		 }
		 
		for(int j=tailh;j<=r;j++)
		 if(twice[a[j]]){
				int t=0;
				if((head-1)<=tail)t=f[a[j]][tail]-f[a[j]][head-1];
				if(t){
					if((t&1) && (!((t+twice[a[j]])&1)))ans++;
					 else if((!(t&1)) && ((t+twice[a[j]])&1))ans--;
				}else if(!(twice[a[j]]&1))ans++;
			twice[a[j]]=0;
			}
		printf("%d\n",ans);
	}
	return 0;
}