【BZOJ2821】作诗(Poetize) 分块

时间:2024-11-11 15:36:26

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

  还没AC就写题解2333,原因是因为入BZ坑以来第二次遇到卡评测的事辣,虽然没法评测,遇到这种事莫名感兴趣XD。
  比较神的分块辣,比之前做的麻烦多了,一开始感觉弹飞绵羊思路挺有趣,一看这道题,强制在线,10^5,smg,看题解,脑冻好大啊,我自己肯定想不出,首先我们需要预处理出两个东西。
  1.f[i][j]:第i块到第j块中合法的个数。
  2.前缀和[i][j]:前i个汉字中j出现的个数,当然并不是用2维数组的来实现的,而是在结构体中以汉字为第一关键字,以所在序号为第二关键字排序,用的时候二分即可。
  这样,题目变得很明了了,查询时,对于同块或相邻块内的查询,直接暴力二分统计即可,若不在同一个块内,则开始直接ans=f[pos[l]+1][pos[r]-1](pos[]为所在块),也就是吧l,r之间的完整的块直接用f[][],剩下的暴力二分,注意细节。
  但是!!!跟着黄学长写,分成sqrt(n)块的我t了,要分成sqrt(n/logn)块,不知道哪里写残了,一定要写个分成sqrt(n)的版本出来。
  
 #include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#define N 100010
#define sN 1555
#define inf 0x7fffffff
using namespace std;
struct data{int v,p;}b[N];
int f[sN][sN],first[N],last[N],flag[N],temp[N],a[N],pos[N],L[sN],R[sN];
int n,m,block,cnt,c;
bool cmp(data a,data b) {if (a.v<b.v) return ; else if (a.v==b.v && a.p<b.p) return ; return ;}
inline int read()
{
char c;
int ans=;
while ((c=getchar())==' ' || c=='\n' || c=='\r');
ans=c-'';
while (isdigit(c=getchar())) ans=ans*+c-'';
return ans;
}
void part1_pre()
{
for (int i=;i<=cnt;i++)
{
int tot;
for (int j=L[i];j<=n;j++) temp[a[j]]=;
tot=;
for (int j=L[i];j<=n;j++)
{
if (!(temp[a[j]]&) && temp[a[j]]) tot--;
temp[a[j]]++;
if (!(temp[a[j]]&)) tot++;
f[i][pos[j]]=tot;
}
}
for (int i=;i<=n;i++)
{
b[i].v=a[i];
b[i].p=i;
}
sort(b+,b+n+,cmp);
for (int i=;i<=n;i++)
{
if (!first[b[i].v]) first[b[i].v]=i;
last[b[i].v]=i;
}
}
int findup(int x,int v)
{
int tem=,l=first[v],r=last[v];
while (l<=r)
{
int mid=(l+r)>>;
if (x<b[mid].p) r=mid-;
else {l=mid+;tem=mid;}
}
return tem;
}
int finddown(int x,int v)
{
int tem=inf,l=first[v],r=last[v];
while (l<=r)
{
int mid=(l+r)>>;
if (x>b[mid].p) l=mid+;
else {r=mid-;tem=mid;}
}
return tem;
}
int find(int x,int y,int v) {return max(findup(y,v)-finddown(x,v)+,);}
int query(int x,int y)
{
int ans=,t,t1;
if (pos[x]==pos[y] || pos[x]+==pos[y])
{
for (int i=x;i<=y;i++)
{
if (!flag[a[i]])
{
t=find(x,y,a[i]);
if (!(t&)) {ans++;flag[a[i]]=;}
}
}
for (int i=x;i<=y;i++) flag[a[i]]=;
}
else
{
int l=L[pos[x]+],r=R[pos[y]-];
ans=f[pos[x]+][pos[y]-];
for (int i=x;i<l;i++)
{
if (!flag[a[i]])
{
flag[a[i]]=;
t=find(x,y,a[i]);t1=find(l,r,a[i]);
if (!(t&))
{if ((t1&) || !t1) ans++;}
else
if (!(t1&) && t1) ans--;
}
}
for (int i=r+;i<=y;i++)
{
if (!flag[a[i]])
{
flag[a[i]]=;
t=find(x,y,a[i]);t1=find(l,r,a[i]);
if (!(t&))
{if ((t1&) || !t1) ans++;}
else
if (!(t1&) && t1) ans--;
}
}
for (int i=x;i<l;i++) flag[a[i]]=;
for (int i=r+;i<=y;i++) flag[a[i]]=;
}
return ans;
}
int main()
{
int ans=;
n=read();c=read();m=read();
block=sqrt((double)n/log((double)n)*log());
for (int i=;i<=n;i++) a[i]=read();
if (n%block) cnt=n/block+;
else cnt=n/block;
for (int i=;i<=n;i++) pos[i]=(i-)/block+;
for(int i=;i<=cnt;i++)
L[i]=(i-)*block+,R[i]=i*block;
R[cnt]=n;
part1_pre();
for (int i=;i<=m;i++)
{ int l=read(),r=read();
int x=(l+ans)%n+,y=(r+ans)%n+;
if (x>y) swap(x,y);
ans=query(x,y);
printf("%d\n",ans);
}
return ;
}

sqrt(n/logn)版本