[BZOJ3295][Cqoi2011]动态逆序对

题意

---

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

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

---

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

#include <cstdio>#include <cstdlib>#define N 100010typedef 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 100010typedef 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;}