【BZOJ2821】作诗(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

　　还没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)版本