BZOJ 3744: Gty的妹子序列 【分块求区间逆序对】
阅读原文时间:2021年04月20日阅读:1

题目分析:

n n   l o g n n\sqrt n~logn nn ​ logn的分块(+树状数组预处理)+ 主席树维护区间值求散块
n n l o g n n\sqrt {nlogn} nnlogn ​的分块+树状数组处理前i块小于j(及大于j)的数有多少个,相当于用分块代替了主席树。
n n n\sqrt n nn ​的分块,其实可以使上面的方法最后处理散块的时候用散块之间桶排归并,散块自身预处理前缀和/后缀和。

我打的第二种(求f[i][j]的时候循环顺序错了RE到自闭。。)。
Code:

#include<cstdio>
#include<cctype>
#include<algorithm>
#define LL long long
#define maxn 50005
#define S 250
using namespace std;
inline void read(int &a){
    char c;while(!isdigit(c=getchar()));
    for(a=c-'0';isdigit(c=getchar());a=a*10+c-'0');
}
int n,m,a[maxn],N,b[maxn],bel[maxn],arr[maxn];
int s1[maxn/S+3][maxn],s2[maxn/S+3][maxn],g[maxn/S+3],f[maxn/S+3][maxn/S+3];
inline void update(int i,int d){for(;i<=n;i+=i&-i) arr[i]+=d;}
inline int getsum(int i){int s=0;for(;i;i-=i&-i) s+=arr[i];return s;}
int solve(int x,int y){
    int ret=0;
    for(int i=y;i>=x;i--) ret+=getsum(a[i]-1),update(a[i],1);
    for(int i=y;i>=x;i--) update(a[i],-1);
    return ret;
}
int solve(int x,int y,int l,int r){
    int ret=0;
    for(int i=r;i>=l;i--) ret+=getsum(a[i]-1),update(a[i],1);
    for(int i=y;i>=x;i--) ret+=getsum(a[i]-1),update(a[i],1);
    for(int i=r;i>=l;i--) update(a[i],-1);
    for(int i=y;i>=x;i--) update(a[i],-1);
    return ret;
}
int main()
{
    read(n);
    for(int i=1;i<=n;i++) read(a[i]),b[i]=a[i],bel[i]=(i-1)/S;
    sort(b+1,b+1+n);
    N=unique(b+1,b+1+n)-b-1;
    for(int i=1;i<=n;i++) a[i]=lower_bound(b+1,b+1+N,a[i])-b;
    for(int i=1;i<=n;i++) s1[bel[i]][a[i]+1]++,s2[bel[i]][a[i]-1]++;
    for(int i=0;i*S<n;i++){
        for(int j=1;j<=N;j++) s1[i][j]+=s1[i][j-1];
        for(int j=1;j<=N;j++) s1[i][j]+=s1[i-1][j];
        for(int j=N;j>=1;j--) s2[i][j]+=s2[i][j+1];
        for(int j=N;j>=1;j--) s2[i][j]+=s2[i-1][j];
    }
    for(int i=0;i*S<n;i++){
        for(int j=min((i+1)*S,n);j>i*S;j--) g[i]+=getsum(a[j]-1),update(a[j],1);
        for(int j=min((i+1)*S,n);j>i*S;j--) update(a[j],-1);
        f[i][i]=g[i];
    }
    for(int i=bel[n];i>=0;i--)
        for(int j=i+1;j*S<n;j++){
            f[i][j]=f[i+1][j]+g[i];
            for(int k=i*S+1;k<=(i+1)*S;k++) f[i][j]+=s1[j][a[k]]-s1[i][a[k]];
        }
    read(m);int x,y,l,r,ans=0;
    while(m--){
        read(x),read(y),x^=ans,y^=ans,l=bel[x],r=bel[y];
        if(l+1>=r) printf("%d\n",ans=solve(x,y));
        else{
            ans=f[l+1][r-1];
            for(int i=x;i<=(l+1)*S;i++) ans+=s1[r-1][a[i]]-s1[l][a[i]];
            for(int i=r*S+1;i<=y;i++) ans+=s2[r-1][a[i]]-s2[l][a[i]];
            ans+=solve(x,(l+1)*S,r*S+1,y);
            printf("%d\n",ans);
        }
    }
}

手机扫一扫

移动阅读更方便

阿里云服务器
腾讯云服务器
七牛云服务器

你可能感兴趣的文章