bzoj2821 作诗(Poetize)

2821: 作诗(Poetize)

Time Limit: 50 Sec   Memory Limit: 128 MB
Submit: 2321   Solved: 658
[ Submit][ Status][ Discuss]

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

很神奇的分块题！！！

为什么说它神奇呢？因为我的程序在本机过了，26s，然而交到网站上就TLE了......并不知道为什么(大概RP不好吧)，希望有好心的小伙伴帮我看看，谢谢啦！！

这道题是强制在线的。我们将整个序列分块，然后预处理两个量：s[i][j]表示第i种汉字在前j块中的出现次数，ans[i][j]表示第i块到第j块出现偶数次的汉字的数目。

对于每个询问，先将整块的ans[i][j]加到答案中，再把多余部分统计一遍数目，和块内的进行比较，并对答案做出调整。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define maxn 100005
using namespace std;
int n,m,c,block,num,pre=0,l,r;
int a[maxn],ans[320][320],s[320][maxn],pos[maxn],f[maxn];
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int main()
{
	n=read();c=read();m=read();
	F(i,1,n) a[i]=read();
	block=int(sqrt(n));num=n/block;
	F(i,1,n) pos[i]=(i-1)/block+1;
	F(i,1,num)
	{
		int cnt=0;
		memset(f,0,sizeof(f));
		F(j,(i-1)*block+1,n)
		{
			f[a[j]]++;
			if (!(f[a[j]]&1)) cnt++;
			else if (f[a[j]]>2) cnt--;
			if (j%block==0) ans[i][j/block]=cnt;
		}
	}
	F(i,1,n) s[pos[i]][a[i]]++;
	F(i,2,num) F(j,1,c) s[i][j]+=s[i-1][j];
	F(i,1,m)
	{
		int x=read(),y=read();
		l=(x+pre)%n+1;r=(y+pre)%n+1;
		if (l>r) swap(l,r);
		int pl=pos[l],pr=pos[r];
		pre=0;
		memset(f,0,sizeof(f));
		if (pl==pr)
		{
			F(j,l,r)
			{
				f[a[j]]++;
				if ((f[a[j]]&1)&&f[a[j]]>2) pre--;
				if (!(f[a[j]]&1)) pre++;
			}
		}
		else
		{
			pre=ans[pl+1][pr-1];
			F(j,l,pl*block)
			{
				f[a[j]]++;
				int tmp=s[pr-1][a[j]]-s[pl][a[j]];
				if (tmp&1){if (f[a[j]]&1) pre++;else pre--;}
				else if (tmp){if (f[a[j]]&1) pre--;else pre++;}
				else{if ((f[a[j]]&1)&&f[a[j]]>2) pre--;if (!(f[a[j]]&1)) pre++;}
			}
			F(j,(pr-1)*block+1,r)
			{
				f[a[j]]++;
				int tmp=s[pr-1][a[j]]-s[pl][a[j]];
				if (tmp&1){if (f[a[j]]&1) pre++;else pre--;}
				else if (tmp){if (f[a[j]]&1) pre--;else pre++;}
				else{if ((f[a[j]]&1)&&f[a[j]]>2) pre--;if (!(f[a[j]]&1)) pre++;}
			}
		}
		printf("%d\n",pre);
	}
}