[bzoj3744]Gty的妹子序列

题目大意

在线兹瓷区间求逆序对个数。

双倍经验

最简单的做法是预处理ans和sum这个不用说了具体可以看根号算法题库里的经典分块思路。
当然这道题有升级版Gty的文艺妹子序列，那道题多一个修改操作。
反正我直接改了一下那题的代码交了。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<ctime>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
const int maxn=50000+10,maxc=250+10;
int ans[maxc][maxc],num[maxc];
int sum[maxc][maxn],tree[maxn];
int belong[maxn],a[maxn],b[maxn],sta[maxn];
int i,j,k,l,r,s,t,n,m,now,top,c;
bool czy;
int lowbit(int x){
    return x&-x;
}
void change(int k,int t){
    while(k<=n){  
        tree[k]+=t;  
        k+=lowbit(k);  
    }  
}
int query(int k){
    int t=0;  
    while (k){  
        t+=tree[k];  
        k-=lowbit(k);  
    }  
    return t;  
}
void change2(int id,int k,int t){
    while(k<=belong[n]){  
        ans[id][k]+=t;  
        k+=lowbit(k);  
    }  
}
int query2(int id,int k){
    int t=0;  
    while (k){  
        t+=ans[id][k];  
        k-=lowbit(k);  
    }  
    return t;  
}
void change3(int id,int k,int t){
    while(k<=n){  
        sum[id][k]+=t;  
        k+=lowbit(k);  
    }  
}
int query3(int id,int k){
    int t=0;  
    while (k){  
        t+=sum[id][k];  
        k-=lowbit(k);  
    }  
    return t;  
}
int read(){
    int x=0;
    char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x;
}
void write(int x){
    while (x){
        sta[++top]=x%10;
        x/=10;
    }
    while (top){
        putchar(sta[top]+'0');
        top--;
    }
    putchar('\n');
}
int main(){
    //freopen("data6.in","r",stdin);freopen("3787.out","w",stdout);
    czy=1;
    n=read();
    c=floor(sqrt(n));
    fo(i,1,n) a[i]=read(),b[i]=a[i],belong[i]=(i-1)/c+1;
    sort(b+1,b+n+1);
    l=unique(b+1,b+n+1)-b-1;
    fo(i,1,n) a[i]=lower_bound(b+1,b+l+1,a[i])-b;
    fo(i,1,belong[n]-1){
        fo(j,(i-1)*c+1,i*c) change(a[j],1);
        fo(j,i+1,belong[n]){
            t=0;
            fo(k,(j-1)*c+1,min(j*c,n)) t+=query(n)-query(a[k]);
            change2(i,j,t);
        }
        fo(j,(i-1)*c+1,i*c) change(a[j],-1);
    }
    fo(i,1,belong[n])
        fo(j,1,min(i*c,n)) change3(i,a[j],1);
    fo(i,1,belong[n]){
        fo(j,(i-1)*c+1,min(n,i*c)){
            num[i]+=query(n)-query(a[j]);
            change(a[j],1);
        }
        fo(j,(i-1)*c+1,min(n,i*c)) change(a[j],-1);
    }
    m=read();
    while (m--){
        j=read();k=read();
        if (czy) j^=now,k^=now;
        now=0;
        l=belong[j];r=belong[k];
        if (r-l<=1){
            fo(i,j,k){
                now+=query(n)-query(a[i]);
                change(a[i],1);
            }
            fo(i,j,k) change(a[i],-1);
                //write(now);
            printf("%d\n",now);
            continue;
        }
        fo(i,l+1,r-1){
            now+=num[i];
            now+=query2(i,r-1)-query2(i,i);
        }
        fo(i,j,l*c){
            now+=query(n)-query(a[i]);
            now+=query3(r-1,a[i]-1)-query3(l,a[i]-1);
            change(a[i],1);
        }
        fo(i,(r-1)*c+1,k){
            now+=query(n)-query(a[i]);
            now+=query3(r-1,n)-query3(r-1,a[i])-query3(l,n)+query3(l,a[i]);
            change(a[i],1);
        }
        fo(i,j,l*c) change(a[i],-1);
        fo(i,(r-1)*c+1,k) change(a[i],-1);
        //write(now);
        printf("%d\n",now);
    }
    //printf("\n%d\n",clock());
}