[BZOJ3295] [Cqoi2011]动态逆序对

可用树套树可用CDQ分治（据说还有分块

基本上算是两种方法的版题。。。详解网上有很多，推荐   http://blog.csdn.net/u011542204/article/details/50571409

只写了树套树。。CDQ感觉和之前写的版差不多就懒得写了。。

上代码

#include<bits/stdc++.h>
using namespace std;
#define N 100005
#define X 5000005
#define lowbit(x) ((x)&-(x))
int num[N],pos[N],n,m,t[N],a1[N],a2[N],A[N],B[N],root[X],ls[X],rs[X],sz,sum[X];
long long ans;
int read(int& x){
    x=0; int f=1,a=getchar();
    while(a<'0' || a>'9') {if(a=='-') f=-1; a=getchar();}
    while(a>='0' && a<='9') x=x*10+a-'0',a=getchar(); x*=f;
}
int getans(int x){
    int ret=0;
    for(int i=x;i;i-=lowbit(i)) ret+=t[i]; return ret; 
}
void update(int& y,int l,int r,int k){
    if(!y) y=++sz;  sum[y]++; if(l==r) return;
    int mid=(l+r)>>1;
    if(k<=mid) update(ls[y],l,mid,k);
    else update(rs[y],mid+1,r,k);
} 
int askmore(int x,int y,int k){
    int ret=0; A[0]=B[0]=0; x--;
    for(int i=x;i;i-=lowbit(i)) A[++A[0]]=root[i];
    for(int i=y;i;i-=lowbit(i)) B[++B[0]]=root[i];
    int l=1,r=n;
    while(l<r){
        int mid=(l+r)>>1;
        if(k<=mid){
            for(int i=1;i<=A[0];i++) ret-=sum[rs[A[i]]],A[i]=ls[A[i]];
            for(int i=1;i<=B[0];i++) ret+=sum[rs[B[i]]],B[i]=ls[B[i]];
            r=mid;
        }
        else{
            for(int i=1;i<=A[0];i++) A[i]=rs[A[i]];
            for(int i=1;i<=B[0];i++) B[i]=rs[B[i]];
            l=mid+1;
        }
    }
    return ret;
}
int askless(int x,int y,int k){
    int ret=0; A[0]=B[0]=0; x--;
    for(int i=x;i;i-=lowbit(i)) A[++A[0]]=root[i];
    for(int i=y;i;i-=lowbit(i)) B[++B[0]]=root[i];
    int l=1,r=n;
    while(l<r){
        int mid=(l+r)>>1;
        if(k>mid){
            for(int i=1;i<=A[0];i++) ret-=sum[ls[A[i]]],A[i]=rs[A[i]];
            for(int i=1;i<=B[0];i++) ret+=sum[ls[B[i]]],B[i]=rs[B[i]];
            l=mid+1;
        }
        else{
            for(int i=1;i<=A[0];i++) A[i]=ls[A[i]];
            for(int i=1;i<=B[0];i++) B[i]=ls[B[i]];
            r=mid;
        }
    }
    return ret;
}
int main(){
    read(n); read(m);
    for(int i=1;i<=n;i++){
        read(num[i]); pos[num[i]]=i;
        a1[i]=getans(n)-getans(num[i]);                 //以该数为结尾的逆序对 
        for(int j=num[i];j<=n;j+=lowbit(j))  t[j]++;
        ans+=a1[i];
    }
    memset(t,0,sizeof(t));
    for(int i=n;i>=1;i--){
        a2[i]=getans(num[i]-1);
        for(int j=num[i];j<=n;j+=lowbit(j)) t[j]++;    //以该数为开头的逆序对 
    }
    for(int a,i=1;i<=m;i++){
        printf("%lld
",ans);
        read(a); a=pos[a];  
        ans-=(a1[a]+a2[a]-askmore(1,a-1,num[a])-askless(a+1,n,num[a]));    //删点转化为加点 
        for(int j=a;j<=n;j+=lowbit(j)) update(root[j],1,n,num[a]);
    }
    return 0;
}