bzoj 3809 莫队

收获：

　　1、分块时顺便记录每个位置所属的块，然后一次排序就OK了。

　　2、要权衡在“区间移动”与“查询结果”之间的时间，莫队算法一般区间移动频率远大于查询结果，所以我们选择的辅助数据结构时就要注意了，我最开始写的是值域线段树，自己生成的极限数据要1m8s，改成树状数组后要24s，还是过不了，hzwer只要13s，细看hzwer的代码，发现Ta用的是分块，O(1)修改O(n^0.5)查询，改成分块后的确快多了。

　　3、块的大小理论最优值是O(n*m^-0.5)，最开始设成这个交上去35卡过，改成hzwer的n/2后29s，所以做题时要自己试一下。

 /**************************************************************
     Problem: 3809
     User: idy002
     Language: C++
     Result: Accepted
     Time:29336 ms
     Memory:26988 kb
 ****************************************************************/

 #include <cstdio>
 #include <cctype>
 #include <cstring>
 #include <cmath>
 #include <algorithm>
 #define maxn 100010
 #define maxm 1000010
 #define lowbit(i) ((i)&(-(i)))
 using namespace std;

 int n, m;
 int lx[maxn], rx[maxn], mccno[maxn], stot;
 int w[maxn];
 int ans[maxm];
 int cnt[maxn], siz[maxn], cmc[maxn];

 struct Qu {
     int l, r, a, b, id;
     bool operator<( const Qu & c ) const {
         return mccno[l]<mccno[c.l] || ( mccno[l]==mccno[c.l] && r<c.r );
     }
 };
 Qu qu[maxm];

 void add( int pos ) {
     siz[pos]++;
     if( siz[pos]== ) cmc[mccno[pos]]++;
 }
 void del( int pos ) {
     siz[pos]--;
     if( siz[pos]== ) cmc[mccno[pos]]--;
 }
 int query( int l, int r ) {
     int lm = mccno[l], rm = mccno[r];
     int rt = ;
     if( lm==rm ) {
         for( int j=l; j<=r; j++ )
             if( siz[j] ) rt++;
         return rt;
     }
     for( int b=lm+; b<rm; b++ )
         rt += cmc[b];
     for( int j=l; j<=rx[lm]; j++ )
         if( siz[j] ) rt++;
     for( int j=lx[rm]; j<=r; j++ )
         if( siz[j] ) rt++;
     return rt;
 }
 void partition() {
     int len = (int)ceil(sqrt(n/)+);
     stot = n/len;
     rx[] = ;
     for( int i=; i<=stot; i++ ) {
         lx[i] = rx[i-]+;
         rx[i] = rx[i-]+len;
     }
     if( rx[stot]!=n ) {
         stot++;
         lx[stot] = rx[stot-]+;
         rx[stot] = n;
     }
     for( int i=; i<=stot; i++ )
         for( int j=lx[i]; j<=rx[i]; j++ )
             mccno[j] = i;
 }
 void work() {
     sort( qu+, qu++m);

     int lf = qu[].l;
     int rg = qu[].r;
     for( int j=lf; j<=rg; j++ )
         add(w[j]);
     ans[qu[].id] = query( qu[].a, qu[].b );
     for( int q=; q<=m; q++ ) {
         while( qu[q].l<lf ) add(w[--lf]);
         while( qu[q].l>lf ) del(w[lf++]);
         while( qu[q].r>rg ) add(w[++rg]);
         while( qu[q].r<rg ) del(w[rg--]);
         ans[qu[q].id] = query( qu[q].a, qu[q].b );
     }
 }
 int main() {
     scanf( "%d%d", &n, &m );
     for( int i=; i<=n; i++ )
         scanf( "%d", w+i );
     for( int i=; i<=m; i++ ) {
         scanf( "%d%d%d%d", &qu[i].l, &qu[i].r, &qu[i].a, &qu[i].b );
         qu[i].id = i;
     }
     partition();
     work();
     for( int i=; i<=m; i++ )
         printf( "%d\n", ans[i] );
 }