BZOJ 2821 作诗(Poetize)

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

HINT

对于100%的数据，1<=n,c,m<=10^5

Source

By lydrainbowcat

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

分块+二分~

大体做法和蒲公英差不多，但是这道题lower_bound()换成了二分来查找，计算的时候要判断a和b之间是否有单独的块。

整块之外的部分要用暴力来更新，做法是求出整个a~b中的a[i]个数now1，再求出整块中的a[i]个数now2，然后综合判断是否要修改ans。注意题目里说的是正偶数，所以还要判断now是否为0。

但是我们注意到这里的暴力复杂度很高。所以我们开一个b[]数组来记录是否已经用a[i]更新过ans，最后的时候一定要把b[]清零。

分块大小的求法（来自黄学长）：

设分块大小为x，分块数n/x，预处理n/x*n

m与n同级，视为n个询问，每次询问二分x次n*x*logn（除非相同的数字很多，否则logn会很小）

n*（x*logn+n/x）

分块大小应该是sqrt(n/logn)

（代码看起来很长，实际上很好写~就是边界容易出错~）

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;

int n,c,m,a[100001],id[100001],cnt,jin,tot,f[1505][1505],num[100001],fir[100001],las[100001],ans;
bool b[100001];

struct node{
	int a,id;
}k[100001];

bool operator < (node u,node v)
{
	return u.a==v.a ? u.id<v.id:u.a<v.a;
}

void init()
{
	int ans;
	for(int i=1;i<=tot;i++)
	{
		ans=0;
		for(int j=(i-1)*jin+1;j<=n;j++) num[a[j]]=0;
		for(int j=(i-1)*jin+1;j<=n;j++)
		{
			num[a[j]]++;
			if(num[a[j]]&1 && (num[a[j]]-1)) ans--;
			else if(!(num[a[j]]&1) && num[a[j]]) ans++;
			f[i][id[j]]=ans;
		}
	}
	for(int i=1;i<=n;i++) k[i].id=i,k[i].a=a[i];
	sort(k+1,k+n+1);
	for(int i=1;i<=n;i++)
	{
		if(!fir[k[i].a]) fir[k[i].a]=i;
		las[k[i].a]=i;
	}
}

int finup(int u,int kk)
{
	int l=fir[kk],r=las[kk],ans=0;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(u<k[mid].id) r=mid-1;
		else ans=mid,l=mid+1;
	}
	return ans;
}

int findown(int u,int kk)
{
	int l=fir[kk],r=las[kk],ans=0x7fffffff;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(u>k[mid].id) l=mid+1;
		else r=mid-1,ans=mid;
	}
	return ans;
}

int findd(int u,int v,int k)
{
	return max(finup(v,k)-findown(u,k)+1,0);
}

int cal(int l,int r)
{
	int ans=0;
	if(id[l]==id[r] || id[l]+1==id[r])
	{
		for(int i=l;i<=r;i++)
		{
			if(b[a[i]]) continue;b[a[i]]=1;
			int now=findd(l,r,a[i]);
			if(!(now&1) && now) ans++;
		}
		for(int i=l;i<=r;i++) b[a[i]]=0;
	}
	else
	{
		ans=f[id[l]+1][id[r]-1];
		int ia=id[l]*jin+1,ib=(id[r]-1)*jin;
		for(int i=l;i<ia;i++)
		{
			if(b[a[i]]) continue;
			int now1=findd(l,r,a[i]),now2=findd(ia,ib,a[i]);
			if(now1&1 && now2 && !(now2&1)) ans--;
			else if(!(now1&1) && (now2&1 || !now2)) ans++;
			b[a[i]]=1;
		}
		for(int i=ib+1;i<=r;i++)
		{
			if(b[a[i]]) continue;
			int now1=findd(l,r,a[i]),now2=findd(ia,ib,a[i]);
			if(now1&1 && now2 && !(now2&1)) ans--;
			else if(!(now1&1) && (now2&1 || !now2)) ans++;
			b[a[i]]=1;
		}
		for(int i=l;i<ia;i++) b[a[i]]=0;
		for(int i=ib+1;i<=r;i++) b[a[i]]=0;
	}
	return ans;
}

int main()
{
	scanf("%d%d%d",&n,&c,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	jin=sqrt((double)n/log((double)n)*log(2));
	for(int i=1;i<=n;i++) id[i]=(i-1)/jin+1;
	tot=(n-1)/jin+1;
	init();
	while(m--)
	{
		int a,b;
		scanf("%d%d",&a,&b);
		a=(a+ans)%n+1;b=(b+ans)%n+1;
		if(a>b) swap(a,b);
		ans=cal(a,b);
		printf("%d\n",ans);
	}
	return 0;
}