bzoj3744 Gty的妹子序列

我是萌萌的传送门

感觉这题还是不错的……虽然其实算是比较水的题= =

首先分块，令f[i][j]表示第i块到第j块的逆序对数，询问的时候直接计算不完整块与完整块以及不完整块之间的逆序对。

不完整块之间的逆序对直接树状数组暴力求，至于不完整块和完整块的逆序对，我是令s[i]表示前i块的权值前缀和，这样单次查询O(1)，可以减小一点常数，代价是空间稍微费了点……

预处理O(nsqrt(n)logn)，单次查询O(sqrt(n)logn)，空间O(nsqrt(n))，好吧我懒得算如何调节块大小来降低复杂度了，于是就随便找了个233当的块大小= =

细节并不多，然而因为各种脑残错误调了半天才敢交，又因为没看见要离散化和没加强制在线RE两次，真是废了……

/**************************************************************
    Problem: 3744
    User: hzoier
    Language: C++
    Result: Accepted
    Time:9688 ms
    Memory:49700 kb
****************************************************************/
 
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=50010,B=233,maxb=245;
void add(int,int);
int query(int);
int a[maxn],b[maxn],id[maxn],L[maxb]={0},R[maxb]={0},cntb,f[maxb][maxb]={{0}},s[maxb][maxn]={{0}},c[maxn]={0};
int n,m,l,r,lb,rb,ans=0;
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        cntb=id[i]=(i-1)/B+1;
        if(!L[id[i]])L[id[i]]=i;
        R[id[i]]=i;
    }
    copy(a+1,a+n+1,b+1);
    sort(b+1,b+n+1);
    for(int i=1;i<=n;i++)s[id[i]][a[i]=lower_bound(b+1,b+n+1,a[i])-b]++;
    for(int i=1;i<=cntb;i++){
        for(int j=1;j<=n;j++)s[i][j]+=s[i-1][j];
        for(int j=i;j<=cntb;j++){
            for(int k=L[j];k<=R[j];k++){
                f[i][j]+=query(a[k]+1);
                add(a[k],1);
            }
            f[i][j]+=f[i][j-1];
        }
        memset(c,0,sizeof(c));
    }
    for(int i=1;i<=cntb;i++)for(int j=1;j<=n;j++)s[i][j]+=s[i][j-1];
    scanf("%d",&m);
    while(m--){
        scanf("%d%d",&l,&r);
        l^=ans;r^=ans;
        if(id[l]>=id[r]-1){
            ans=0;
            for(int i=l;i<=r;i++){
                ans+=query(a[i]+1);
                add(a[i],1);
            }
            for(int i=l;i<=r;i++)add(a[i],-1);
        }
        else{
            lb=id[l]+1;rb=id[r]-1;
            ans=f[lb][rb];
            for(int i=l;i<L[lb];i++){
                ans+=query(a[i]+1)+s[rb][a[i]-1]-s[lb-1][a[i]-1];
                add(a[i],1);
            }
            for(int i=R[rb]+1;i<=r;i++){
                ans+=query(a[i]+1)+(s[rb][n]-s[rb][a[i]])-(s[lb-1][n]-s[lb-1][a[i]]);
                add(a[i],1);
            }
            for(int i=l;i<L[lb];i++)add(a[i],-1);
            for(int i=R[rb]+1;i<=r;i++)add(a[i],-1);
        }
        printf("%d
",ans);
    }
    return 0;
}
void add(int x,int d){
    while(x){
        c[x]+=d;
        x&=x-1;
    }
}
int query(int x){
    int ans=0;
    while(x<=n){
        ans+=c[x];
        x+=x&-x;
    }
    return ans;
}