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

---

solution: 此题可以用分块来做。

/**************************************************************
    Problem: 2821
    User: Venishel
    Language: C++
    Result: Accepted
    Time:26868 ms
    Memory:14000 kb
****************************************************************/

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define Inf 0x7fffffff
#define N 100007
#define BK 1500
struct DATA{ int p, v ; }b[N];
inline bool operator<(DATA a,DATA b){
    return (a.v<b.v)||(a.v==b.v&&a.p<b.p);
}
int v[N], bl[N], tg[N], L[N], R[N], ft[N], lt[N], vis[N];
int f[BK][BK];
int n, m, blk, cnt;
void pre_wkg(){
    int tot;
    for (register int i = 1; i <= cnt; i++ ){
        for (register int j = L[i]; j <= n; j++ ) tg[ v[j] ] = 0;
        tot = 0;
        for (register int j = L[i]; j <= n; j++ ){
            if ( !( tg[v[j]]&1 ) && tg[v[j]] ) tot --;
            tg[v[j]] ++;
            if ( !( tg[v[j]]&1 ) ) tot++;
            f[i][ bl[j] ] = tot;
        }
    }
    for (register int i = 1; i <= n; i++ ) b[i].p = i, b[i].v = v[i];
    sort( b+1, b+n+1);
    for (register int i = 1; i <= n; i++ ) {
         if ( !ft[ b[i].v ] ) ft[ b[i].v ] = i;
         lt[ b[i].v ] = i;
    }
}
int ub(int x,int v){
    int l = ft[v], r = lt[v];
    int ret = 0;                        
    while(l<=r){
        int mid = (l+r)>>1;
        if(x<b[mid].p)r=mid-1;
        else {l=mid+1;ret=mid;}
    }
    return ret;
}
int lb(int x,int v){    
    int l=ft[v],r=lt[v];
    int ret=Inf;
    while(l<=r){
        int mid=(l+r)>>1;
        if(x>b[mid].p)l=mid+1;
        else {ret=mid;r=mid-1;}
    }
    return ret;
}
int find(int x,int y,int v){
    int t = max(ub(y,v) - lb(x,v)+1,0);
    return t;
}

int query(int x,int y){
    int ans=0,t1,t2,ya,bx=bl[x],by=bl[y];
    if( bx == by || bx + 1 == by ){
        for(int i = x ; i <= y; i++){
            ya=v[i];
            if(vis[ya])continue;
            int t = find( x, y, ya);
            if(!(t&1) ){ vis[ya]=1; ans++; }
        }
        for(int i = x; i <= y; i++ ) vis[v[i]] = 0;
    }
    else{
        int l = L[bx+1], r = R[by-1];
        ans = f[bx+1][by-1];
        for(int i = x; i < l; i++){
            ya = v[i]; 
            if( vis[ya] ) continue;
            t1 = find( x, y, ya), t2 = find( l, r, ya);
            if(!(t1&1)){ if((t2&1) || !t2 ) ans++; }
            else if( !(t2&1) && t2 )ans--;
            vis[ya] = 1;
        }
        for(int i = r + 1; i <= y; i++){
            ya = v[i];
            if( vis[ya] )continue;
            t1 = find( x, y, ya), t2 = find( l, r, ya);
            if( !( t1&1 ) ){ if( (t2&1) || !t2 ) ans++;}
            else if( !(t2&1) && t2 ) ans--;
            vis[ya] = 1;
        }
        for(int i = x; i < l; i++ ) vis[ v[i] ] = 0;
        for(int i = r+1; i <= y; i++ ) vis[v[i]]=0;
    }
    return ans;
}

int main(){
    int c;
    scanf( "%d%d%d", &n, &c, &m);
    for( int i = 1; i <= n; i++) scanf( "%d", &v[i]);
    blk = sqrt( (double) n / log((double)n) * log(2) );
    if ( n % blk ) cnt = n / blk + 1;
    else cnt = n / blk;
    for(int i = 1 ; i <= n; i++) bl[i] = ( i - 1 ) / blk + 1;
    for(int i = 1 ; i <= cnt; i++)
        L[i] = (i-1) * blk + 1, R[i] = i * blk;
    R[cnt] = n;
    pre_wkg();
    int ans = 0;
    for(int i = 1; i <= m; i++){
        int l, r;
        scanf( "%d%d", &l, &r);
        int x = ( l + ans ) % n + 1, y = ( r + ans ) % n + 1;
        if( x > y ) swap( x, y );
        ans=query( x, y);
         printf("%d\n", ans);
    }
}