[BZOJ3295][Cqoi2011]动态逆序对

题意

---

给定一个$N$个数序列，以及$M$个操作，每个操作删除一个数，求每次操作前数列中逆序对的个数

$(x,y)$满足$A[x]<A[y]$且$x>y$，即为逆序对。

---

第一想法是线段树套treap，交上去的时候发现treap没有打旋转……结果6.9s过了，打完旋转后9.6s……

#include <cstdio>
#include <cstdlib>
#define N 100010

typedef long long ll;

struct trp{
    trp *ls,*rs;
    int x,sz,rnd,ds;
}P[N<<5],*t=P,_nil,*nil;

struct xds{
    int l,r;
    trp *rt;
}T[N<<2];

int n,m,pst[N],Bit[N];
ll Ans;

inline void reaD(int &x){
    char Ch=getchar();x=0;
    for(;Ch>'9'||Ch<'0';Ch=getchar());
    for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=getchar());
}

inline trp *New(int x,int r){
    trp *w=t++;
    w->ls=w->rs=nil;w->ds=0;
    w->sz=1;w->x=x;w->rnd=r;
    return w;
}

inline void upd(trp *x){x->sz=x->ls->sz+x->rs->sz+1;}
inline void lrt(trp *&x){trp *y=x->ls;x->ls=y->rs;y->rs=x;upd(x);upd(y);x=y;}
inline void rrt(trp *&x){trp *y=x->rs;x->rs=y->ls;y->ls=x;upd(x);upd(y);x=y;}

void Inst(trp *&g,int x){
    if(g==nil){g=New(x,rand());return;}
    if(g->sz++,g->x>x){Inst(g->ls,x);/*if(g->rnd<g->ls->rnd)lrt(g);*/}//就是这个旋转......
    else{Inst(g->rs,x);/*if(g->rnd<g->rs->rnd)rrt(g);*/}
}

void build(int g,int l,int r){
    T[g].l=l;T[g].r=r;T[g].rt=nil;
    if(l==r) return;
    int mid=l+r>>1;
    build(g<<1,l,mid);build(g<<1|1,mid+1,r);
}

void InserT(int g,int p,int x){
    Inst(T[g].rt,x);
    if(T[g].l==T[g].r) return;
    int mid=T[g].l+T[g].r>>1;
    if(p<=mid) InserT(g<<1,p,x);
    else InserT(g<<1|1,p,x);
}

inline void insertb(int x){for(;x;x-=x&-x) Bit[x]++;}

inline int queryb(int x){
    int Res=0;
    for(;x<=n;x+=x&-x) Res+=Bit[x];
    return Res;
}

int w[30],wt;

inline int Qry(trp *g,int x,int k){
    int r=0;
    if(k==1)while(g!=nil)if(g->x>x) r+=g->sz-g->ls->sz,g=g->ls;else g=g->rs;
    else while(g!=nil)if(g->x<x) r+=g->sz-g->rs->sz,g=g->rs;else g=g->ls;
    return r;
}

int query(int g,int l,int r,int x,int k){
    if(T[g].l==l&&T[g].r==r) return Qry(T[g].rt,x,k);
    int mid=T[g].l+T[g].r>>1;
    if(r<=mid) return query(g<<1,l,r,x,k);
    if(l>mid) return query(g<<1|1,l,r,x,k);
    return query(g<<1,l,mid,x,k)+query(g<<1|1,mid+1,r,x,k);
}

void delt(trp *g,int x){
    g->sz--;if(g->x==x)return;
    if(g->x>x)delt(g->ls,x);else delt(g->rs,x);
}

void Dlt(int g,int p,int x){
    delt(T[g].rt,x);
    if(T[g].l==T[g].r) return;
    int mid=T[g].l+T[g].r>>1;
    if(p<=mid) Dlt(g<<1,p,x);
    else Dlt(g<<1|1,p,x);
}

int main(){
    _nil.ls=_nil.rs=nil;nil=&_nil;_nil.sz=_nil.x=_nil.ds=0;
    reaD(n);reaD(m);build(1,1,n);
    for(int i=1,x;i<=n;i++){
        reaD(x),pst[x]=i,InserT(1,i,x);
        Ans+=queryb(x);insertb(x);
    }
    for(int i=1,x;i<=m;i++){
        printf("%lld\n",Ans);reaD(x);
        Ans-=query(1,1,pst[x],x,1)+query(1,pst[x],n,x,2);
        Dlt(1,pst[x],x);
    }
    return 0;
}

看大神博客学了CDQ分治，跑得飞快……

#include <cstdio>
#include <algorithm>
#define N 100010

typedef long long ll;

struct Unit{
    int t,x,y;
    bool operator<(Unit a)const{return t<a.t;}
}A[N],B[N];

inline void reaD(int &x){
    char Ch=getchar();x=0;
    for(;Ch>'9'||Ch<'0';Ch=getchar());
    for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=getchar());
}

int n,m,pst[N],cnt,T[N],tt,Bit[N],t[N];
ll Ans[N];

inline void add(int x){
    for(;x;x-=x&-x)
        if(t[x]!=tt)t[x]=tt,Bit[x]=1;
        else Bit[x]++;
}

inline int query(int x){
    int res=0;
    for(;x<=n;x+=x&-x)
        if(t[x]==tt) res+=Bit[x];
    return res;
}

void S(int l,int r){
    if(l>=r) return;
    int mid=l+r>>1,l1=l,l2=mid+1;
    for(int i=l;i<=r;i++)
        if(A[i].x<=mid) B[l1++]=A[i];
        else B[l2++]=A[i];
    tt++;
    for(int i=mid+1,j=l;i<=r;i++){
        for(;B[j].t<B[i].t&&j<=mid;j++) add(B[j].y);
        Ans[B[i].t]+=query(B[i].y);
    }
    tt++;
    for(int i=l,j=mid+1;i<=mid;i++){
        for(;B[j].t<B[i].t&&j<=r;j++) add(B[j].y);
        Ans[B[i].t]+=(j-mid-1)-query(B[i].y);
    }
    for(int i=l;i<=r;i++)A[i]=B[i];
    S(l,mid);S(mid+1,r);
}

int w[30],wt;

int main(){
    reaD(n);reaD(m);cnt=n;
    for(int i=1,x;i<=n;i++) reaD(x),pst[x]=i;
    for(int i=1,x;i<=m;i++) reaD(x),T[x]=cnt--;
    for(int i=1;i<=n;i++)if(T[i]==0) T[i]=cnt--;
    for(int i=1;i<=n;i++) A[i]=(Unit){T[i],pst[i],i};
    std::sort(A+1,A+1+n);
    S(1,n);
    for(int i=1;i<=n;i++) Ans[i]+=Ans[i-1];
    for(int i=1;i<=m;i++) printf("%lld\n",Ans[n-i+1]);
    return 0;
}