BZOJ_2821_作诗(Poetize)_分块

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

---

用一个前缀桶记录一下前缀块中每个数出现的次数。

然后处理出ans[i][j]表示从i块到j块这部分的答案。

查询时一个块内直接暴力，

否则用和统计ans[i][j]差不多的方法统计零散部分的答案。

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define N 100050
#define M 318
int n,m,C,pos[N],L[M],R[M],ans[M][M],h[M][N],size,block,a[N],w[N],b[N];
int solve(int l,int r) {
    int p=pos[l],q=pos[r],re=0,i;
    if(p==q) {
        for(i=l;i<=r;i++) {
            w[a[i]]++;
            if(w[a[i]]!=1) {
                if(w[a[i]]%2==0) re++;
                else re--;
            }
        }
        for(i=l;i<=r;i++) w[a[i]]=0;
    }else {
        re=ans[p+1][q-1];
        b[0]=0;
        for(i=l;i<=R[p];i++) b[++b[0]]=a[i];
        for(i=L[q];i<=r;i++) b[++b[0]]=a[i];
        for(i=1;i<=b[0];i++) {
            if(w[b[i]]==0) {
                w[b[i]]=h[q-1][b[i]]-h[p][b[i]]+1;
                if(w[b[i]]!=1) {
                    if(w[b[i]]%2==0) re++;
                    else re--;
                }
            }else {
                w[b[i]]++;
                if(w[b[i]]%2==0) re++;
                else re--;
            }
        }
        for(i=1;i<=b[0];i++) w[b[i]]=0;
    }
    return re;
}
int main() {
    scanf("%d%d%d",&n,&C,&m);
    int i,j,k; size=sqrt(n); block=n/size;
    for(i=1;i<=block;i++) {
        L[i]=R[i-1]+1; R[i]=i*size;
        for(j=1;j<=C;j++) h[i][j]=h[i-1][j];
        for(j=L[i];j<=R[i];j++) {
            scanf("%d",&a[j]); pos[j]=i; h[i][a[j]]++;
        }
    }
    if(R[block]!=n) {
        block++; L[block]=R[block-1]+1; R[block]=n;
        for(i=1;i<=C;i++) h[block][i]=h[block-1][i];
        for(i=L[block];i<=n;i++) {
            scanf("%d",&a[i]); pos[i]=block; h[block][a[i]]++;
        }
    }
    for(i=1;i<=block;i++) {
        for(j=i;j<=block;j++) {
            ans[i][j]=ans[i][j-1];
            for(k=L[j];k<=R[j];k++) {
                if(w[a[k]]==0) {
                    w[a[k]]=h[j-1][a[k]]-h[i-1][a[k]]+1;
                    if(w[a[k]]!=1) {
                        if(w[a[k]]%2==0) ans[i][j]++;
                        else ans[i][j]--;
                    }
                }
                else {
                    w[a[k]]++;
                    if(w[a[k]]%2==0) ans[i][j]++;
                    else ans[i][j]--;
                }
            }
            for(k=L[j];k<=R[j];k++) w[a[k]]=0;
        }
    }
    int lstans=0,x,y;
    while(m--) {
        scanf("%d%d",&x,&y);
        x=(x+lstans)%n+1;
        y=(y+lstans)%n+1;
        if(x>y) swap(x,y);
        lstans=solve(x,y);
        printf("%d\n",lstans);
    }
}