bzoj 2821 分块

分块：

先预处理，将原序列分成长度为len的许多块，计算从第i块到第j块的答案，（可以做到O(n*n/len)）。

每次询问时，将询问的区间分成三部分，：左边，中间，右边，中间是尽量大的一个块区间，其答案已经计算得到，左右两边加起来最多有2*len个元素，暴力计算其对答案的影响。O(q*len*f(n))，f(n)是暴力加入一个元素的代价。

这道题f(n)是log(n)

总的复杂度：f(n) = O( n*n/len + q*len*log(n) )，

当len = n*(q*log(n))^-1/2时取最小值.

 /**************************************************************
     Problem: 2821
     User: idy002
     Language: C++
     Result: Accepted
     Time:24352 ms
     Memory:13728 kb
 ****************************************************************/

 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #include <algorithm>
 #define maxn 100010
 #define maxs 1510
 #define isok(a) ((a)>0&&!((a)&1))
 using namespace std;

 int n, c, m;
 int a[maxn], f[maxs][maxs];
 int lx[maxn], rx[maxn], mno[maxn], stot;
 int cnt[maxn], cur_ans;
 int vs[maxn], vt[maxn], vp[maxn];
 bool met[maxn];

 struct Pair {
     int w, pos;
     Pair(){}
     Pair( int w, int pos ) : w(w), pos(pos) {}
 };
 Pair vod[maxn];
 bool cmpwp( const Pair & a, const Pair & b ) {
     return a.w<b.w || (a.w==b.w && a.pos<b.pos);
 }
 bool cmpp( const Pair & a, const Pair & b ) {
     return a.pos<b.pos;
 }

 void part() {
     int len = (int)(double)sqrt((double)n/(log((double)n)/log()))+;
     stot = n/len + (n%len!=);
     for( int i=; i<=stot; i++ ) {
         lx[i] = rx[i-]+;
         rx[i] = rx[i-]+len;
         if( i==stot ) rx[i] = n;
         for( int j=lx[i]; j<=rx[i]; j++ )
             mno[j] = i;
     }
 }
 inline void update( int oldv, int newv ) {
     bool od, nw;
     od = isok(oldv);
     nw = isok(newv);
     if( od && !nw ) cur_ans--;
     else if( !od && nw ) cur_ans++;
 }

 void prep() {
     for( int i=; i<=n; i++ )
         vod[i] = Pair(a[i],i);
     sort( vod+, vod++n, cmpwp );
     for( int i=; i<=n; i++ ) {
         if( vod[i].w!=vod[i-].w )
             vs[vod[i].w] = vt[vod[i-].w] = i;
     }
     vt[vod[n].w] = n+;
     for( int i=; i<=n; i++ ) vp[i] = vod[i].pos;

     for( int i=; i<=stot; i++ ) {
         for( int j=i; j<=stot; j++ ) {
             for( int k=lx[j]; k<=rx[j]; k++ ) {
                 update( cnt[a[k]], cnt[a[k]]+ );
                 cnt[a[k]]++;
             }
             f[i][j] = cur_ans;
         }
         cur_ans = ;
         memset( cnt, , sizeof(cnt) );
     }
 }

 int qu_cnt( int w, int l, int r ) {
     return upper_bound( vp+vs[w], vp+vt[w], r )
         -  lower_bound( vp+vs[w], vp+vt[w], l );
 }
 int query( int l, int r ) {
     int ml=mno[l], mr=mno[r];
     int rt;
     if( ml==mr ) {
         for( int j=l; j<=r; j++ ) {
             update( cnt[a[j]], cnt[a[j]]+ );
             cnt[a[j]]+=;
         }
         rt = cur_ans;
         cur_ans = ;
         for( int j=l; j<=r; j++ )
             cnt[a[j]]--;
         return rt;
     }
     if( mno[l]==mno[l-] ) ml++;
     if( mno[r]==mno[r+] ) mr--;
     cur_ans = f[ml][mr];
     for( int j=l; j<lx[ml]; j++ ) {
         if( met[a[j]] ) continue;
         met[a[j]] = true;
         int t1 = qu_cnt(a[j],l,r);
         int t2 = qu_cnt(a[j],lx[ml],rx[mr]);
         update( t2, t1 );
     }
     for( int j=r; j>rx[mr]; j-- ) {
         if( met[a[j]] ) continue;
         met[a[j]] = true;
         int t1 = qu_cnt(a[j],l,r);
         int t2 = qu_cnt(a[j],lx[ml],rx[mr]);
         update( t2, t1 );
     }
     rt = cur_ans;

     cur_ans = ;
     for( int j=l; j<lx[ml]; j++ ) met[a[j]]=false;
     for( int j=r; j>rx[mr]; j-- ) met[a[j]]=false;
     return rt;
 }

 int main() {
     scanf( "%d%d%d", &n, &c, &m );
     for( int i=; i<=n; i++ ) scanf( "%d", a+i );
     part();
     prep();
     for( int i=,x=,l,r; i<=m; i++ ) {
         scanf( "%d%d", &l, &r );
         l = (l+x)%n+;
         r = (r+x)%n+;
         if( l>r ) swap(l,r);
         printf( "%d\n", x=query(l,r) );
     }
 }