bzoj2141 排队（树套树）

求动态逆序对，BIT套权值线段树即可。
空间和时间都是$O((n+m)log^2n)$的

最近怎么老是算错空间qaq

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 20010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,a[N],rt[N],owo=0,ans=0,aa[N];
struct node{
    int lc,rc,x;
}tr[6400000];
inline void ins(int &p,int l,int r,int x,int val){
    if(!p) p=++owo;tr[p].x+=val;
    if(l==r) return;int mid=l+r>>1;
    if(x<=mid) ins(tr[p].lc,l,mid,x,val);
    else ins(tr[p].rc,mid+1,r,x,val);
}
inline int query(int p,int l,int r,int x,int y){
    if(!p) return 0;if(x<=l&&r<=y) return tr[p].x;
    int mid=l+r>>1,res=0;
    if(x<=mid) res+=query(tr[p].lc,l,mid,x,y);
    if(y>mid) res+=query(tr[p].rc,mid+1,r,x,y);
    return res;
}
inline int ask1(int x){
    int res=0;for(;x;x-=x&-x) res+=rt[x];return res;
}
inline void add1(int x){
    for(;x<=m;x+=x&-x) rt[x]++;
}
inline void add(int x,int y,int val){
    for(;x<=n;x+=x&-x) ins(rt[x],1,m,y,val);
}
inline int ask(int x,int l,int r){
    if(l>r) return 0;int res=0;
    for(;x;x-=x&-x) res+=query(rt[x],1,m,l,r);return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    for(int i=1;i<=n;++i) aa[i]=a[i]=read();
    sort(aa+1,aa+n+1);m=unique(aa+1,aa+n+1)-aa-1;
    for(int i=n;i>=1;--i){
        a[i]=lower_bound(aa+1,aa+m+1,a[i])-aa;
        ans+=ask1(a[i]-1);add1(a[i]);
    }memset(rt,0,sizeof(rt));printf("%d\n",ans);
    for(int i=1;i<=n;++i) add(i,a[i],1);
    int ofo=read();
    while(ofo--){
        int x=read(),y=read();if(x>y) swap(x,y);
        if(a[x]<a[y]) ans++;
        if(a[x]>a[y]) ans--;
        ans+=ask(y-1,a[x]+1,m)-ask(x,a[x]+1,m);
        ans-=ask(y-1,1,a[x]-1)-ask(x,1,a[x]-1);
        ans+=ask(y-1,1,a[y]-1)-ask(x,1,a[y]-1);
        ans-=ask(y-1,a[y]+1,m)-ask(x,a[y]+1,m);
        add(x,a[x],-1);add(y,a[y],-1);swap(a[x],a[y]);
        add(x,a[x],1);add(y,a[y],1);printf("%d\n",ans);
    }return 0;
}