【莫队算法】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[50001];
struct ASK{int l,r,lb,p;}Q[50001];
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[50001],b[50001],cnt;
int d[50001],ans[50001],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=0;while(p){res+=d[p];p-=lowbit(p);}return res;}
int Res,Num;char C,CH[12];
inline int G()
{
    Res=0;C='*';
    while(C<'0'||C>'9')C=getchar();
    while(C>='0'&&C<='9'){Res=Res*10+(C-'0');C=getchar();}
    return Res;
}
inline void P(long long x)
{
    Num=0;if(!x){putchar('0');puts("");return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
    puts("");
}
void makeblock()
{
    sz=sqrt((double)n*4.6); int L,R;
    for(sum=1;sum*sz<n;sum++)
      {
        L=(sum-1)*sz+1; R=sum*sz;
        for(int i=L;i<=R;i++) num[i]=sum;
      }
    L=sz*(sum-1)+1; R=n;
    for(int i=L;i<=R;i++) num[i]=sum;
}
void LiSan()
{
    int en=1; sort(a+1,a+n+1); b[a[1].y]=1;
    for(int i=2;i<=n;i++)
      {
        if(a[i].x!=a[i-1].x) en++;
        b[a[i].y]=en;
      }
}
int main()
{
    n=G();
    for(int i=1;i<=n;i++)
      {
        a[i].x=G();
        a[i].y=i;
      }
    makeblock(); LiSan(); m=G();
    for(int i=1;i<=m;i++)
      {
        Q[i].l=G(); Q[i].r=G();
        Q[i].lb=num[Q[i].l];
        Q[i].p=i;
      }
    sort(Q+1,Q+m+1);
    for(int i=Q[1].l;i<=Q[1].r;i++)
      {
        cnt++;
        add(b[i],1);
        res+=(cnt-getsum(b[i]));
      } ans[Q[1].p]=res;
    for(int i=2;i<=m;i++)
      {
        if(Q[i].l<Q[i-1].l)
          {
            for(int j=Q[i-1].l-1;j>=Q[i].l;j--)
              {
                res+=getsum(b[j]-1);
                add(b[j],1);
              }
            cnt+=(Q[i-1].l-Q[i].l);
          }
        else
          {
            for(int j=Q[i-1].l;j<Q[i].l;j++)
              {
                res-=getsum(b[j]-1);
                add(b[j],-1);
              }
            cnt-=(Q[i].l-Q[i-1].l);
          }
        if(Q[i].r<Q[i-1].r)
          {
            for(int j=Q[i-1].r;j>Q[i].r;j--)
              {
                res-=(cnt-getsum(b[j]));
                cnt--;
                add(b[j],-1);
              }
          }
        else
          {
            for(int j=Q[i-1].r+1;j<=Q[i].r;j++)
              {
                cnt++;
                add(b[j],1);
                res+=(cnt-getsum(b[j]));
              }
          }
        ans[Q[i].p]=res;
      }
    for(int i=1;i<=m;i++) P(ans[i]);
    return 0;
}