【莫队算法】bzoj3289 Mato的文件管理

莫队算法，离线回答询问，按一定大小(sqrt(n*log(n))左右)将答案分块，按 ①左端点所在块②右端点 双关键字排序。

然后暴力转移。

转移的时候用树状数组。

O(n*sqrt(n)*log(n))。

注意：①在一列数的后面添加一个数，逆序对数会增加 数列中比它大的数的个数。

②在一列数的后面删除一个数，逆序对数会减少 数列中比它大的数的个数。

③在一列数的前面添加一个数，逆序对数会增加 数列中比它小的数的个数。

④在一列数的前面删除一个数，逆序对数会减少 数列中比它小的数的个数。

 #include<cstdio>
 #include<algorithm>
 #include<cmath>
 using namespace std;
 struct Point{int x,y;}a[];
 struct ASK{int l,r,lb,p;}Q[];
 bool operator < (const ASK &a,const ASK &b){return a.lb!=b.lb ? a.lb<b.lb : a.r<b.r;}
 bool operator < (const Point &a,const Point &b){return a.x<b.x;}
 int n,m,sz,sum,num[],b[],cnt;
 int d[],ans[],res;
 inline int lowbit(const int &x){return x&(-x);}
 inline void add(int p,const int &x){while(p<=n){d[p]+=x;p+=lowbit(p);}}
 inline int getsum(int p){int res=;while(p){res+=d[p];p-=lowbit(p);}return res;}
 int Res,Num;char C,CH[];
 inline int G()
 {
     Res=;C='*';
     while(C<''||C>'')C=getchar();
     while(C>=''&&C<=''){Res=Res*+(C-'');C=getchar();}
     return Res;
 }
 inline void P(long long x)
 {
     Num=;if(!x){putchar('');puts("");return;}
     while(x>)CH[++Num]=x%,x/=;
     while(Num)putchar(CH[Num--]+);
     puts("");
 }
 void makeblock()
 {
     sz=sqrt((double)n*4.6); int L,R;
     for(sum=;sum*sz<n;sum++)
       {
         L=(sum-)*sz+; R=sum*sz;
         for(int i=L;i<=R;i++) num[i]=sum;
       }
     L=sz*(sum-)+; R=n;
     for(int i=L;i<=R;i++) num[i]=sum;
 }
 void LiSan()
 {
     int en=; sort(a+,a+n+); b[a[].y]=;
     for(int i=;i<=n;i++)
       {
         if(a[i].x!=a[i-].x) en++;
         b[a[i].y]=en;
       }
 }
 int main()
 {
     n=G();
     for(int i=;i<=n;i++)
       {
         a[i].x=G();
         a[i].y=i;
       }
     makeblock(); LiSan(); m=G();
     for(int i=;i<=m;i++)
       {
         Q[i].l=G(); Q[i].r=G();
         Q[i].lb=num[Q[i].l];
         Q[i].p=i;
       }
     sort(Q+,Q+m+);
     for(int i=Q[].l;i<=Q[].r;i++)
       {
         cnt++;
         add(b[i],);
         res+=(cnt-getsum(b[i]));
       } ans[Q[].p]=res;
     for(int i=;i<=m;i++)
       {
         if(Q[i].l<Q[i-].l)
           {
             for(int j=Q[i-].l-;j>=Q[i].l;j--)
               {
                 res+=getsum(b[j]-);
                 add(b[j],);
               }
             cnt+=(Q[i-].l-Q[i].l);
           }
         else
           {
             for(int j=Q[i-].l;j<Q[i].l;j++)
               {
                 res-=getsum(b[j]-);
                 add(b[j],-);
               }
             cnt-=(Q[i].l-Q[i-].l);
           }
         if(Q[i].r<Q[i-].r)
           {
             for(int j=Q[i-].r;j>Q[i].r;j--)
               {
                 res-=(cnt-getsum(b[j]));
                 cnt--;
                 add(b[j],-);
               }
           }
         else
           {
             for(int j=Q[i-].r+;j<=Q[i].r;j++)
               {
                 cnt++;
                 add(b[j],);
                 res+=(cnt-getsum(b[j]));
               }
           }
         ans[Q[i].p]=res;
       }
     for(int i=;i<=m;i++) P(ans[i]);
     return ;
 }