BZOJ 2141 排队 (三维偏序CDQ+树状数组)

题目大意：略

洛谷传送门

和 [CQOI2015]动态逆序对 这道题一样的思路

一开始的序列视为$n$次插入操作

把每次交换操作看成四次操作，删除$x$，删除$y$，加入$x$，加入$y$

把每次操作都看成一个三元组 $<x,w,t>$ 操作位置，权值，以及操作时间

变成了一道三维偏序裸题

外层按操作时间$t$升序，回溯时按操作位置$x$排序

处理左区间对右区间的影响时，正反两次树状数组求逆序对即可

#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N1 30100
#define ll long long
#define dd double
#define inf 0x3f3f3f3f3f3f3f3fll
using namespace std;
 
int gint()
{
    int ret=0,fh=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')fh=-1;c=getchar();}
    while(c>='0'&&c<='9'){ret=ret*10+c-'0';c=getchar();}
    return ret*fh;
}
int n,m,ma,nn,K;
 
struct BIT{
int s[N1];
void update(int x,int w){for(int i=x;i<=ma;i+=(i&(-i))) s[i]+=w;}
int query(int x){int ans=0; for(int i=x;i;i-=(i&(-i))) ans+=s[i]; return ans;}
}b;
struct node{int p,x,w,t,ans;}a[N1],tmp[N1];
int c[N1],q[N1],que[N1],tl;
 
void CDQ(int L,int R)
{
    if(R-L<=1) return;
    int M=(L+R)>>1;
    CDQ(L,M); CDQ(M,R);
    int i=M-1,j=R-1,cnt=0,x;
    while(i>=L&&j>=M)
    {
        if(a[i].x>a[j].x){
            b.update(a[i].w,a[i].p);
            que[++tl]=i; i--;
        }else{
            a[j].ans+=b.query(a[j].w-1);
            j--;
        }
    }
    while(i>=L) {i--;}
    while(j>=M) {a[j].ans+=b.query(a[j].w-1); j--;}
    while(tl) {x=que[tl--]; b.update(a[x].w,-a[x].p);}
     
    i=L,j=M,cnt=0;
    while(i<M&&j<R)
    {
        if(a[i].x<=a[j].x){
            b.update(a[i].w,a[i].p);
            que[++tl]=i; tmp[++cnt]=a[i]; i++;
        }else{
            a[j].ans+=b.query(ma)-b.query(a[j].w);
            tmp[++cnt]=a[j]; j++;
        }
    }
    while(i<M) {tmp[++cnt]=a[i]; i++;}
    while(j<R) {a[j].ans+=b.query(ma)-b.query(a[j].w); tmp[++cnt]=a[j]; j++;}
    while(tl) {x=que[tl--]; b.update(a[x].w,-a[x].p);}
 
    for(i=L;i<R;i++) a[i]=tmp[i-L+1];
}
int cmp(node s1,node s2){return s1.t<s2.t;}
 
int main()
{
    scanf("%d",&n);
    int i,j,x,y; ll ans=0;
    for(i=1;i<=n;i++) c[i]=a[i].w=gint();
    sort(c+1,c+n+1);
    ma=unique(c+1,c+n+1)-(c+1);
    for(i=1;i<=n;i++) a[i].w=lower_bound(c+1,c+ma+1,a[i].w)-c,q[i]=a[i].w;
    for(i=1;i<=n;i++) a[i].p=1,a[i].t=i,a[i].x=i;
    for(i=n;i>=1;i--) ans+=b.query(a[i].w-1),b.update(a[i].w,1);
    memset(b.s,0,sizeof(b.s));
    scanf("%d",&m); nn=n;
    for(i=n+1;i<=n+m;i++)
    {
        x=gint(); y=gint();
        a[++nn].p=-1,a[nn].t=nn,a[nn].x=x,a[nn].w=q[x];//,a[nn].id=i-n;
        a[++nn].p=-1,a[nn].t=nn,a[nn].x=y,a[nn].w=q[y];//a[nn].id=i-n;
        swap(q[x],q[y]);
        a[++nn].p=1,a[nn].t=nn,a[nn].x=x,a[nn].w=q[x];//a[nn].id=i-n;
        a[++nn].p=1,a[nn].t=nn,a[nn].x=y,a[nn].w=q[y];//a[nn].id=i-n;
    }
    CDQ(1,nn+1);
    sort(a+1,a+nn+1,cmp);
    printf("%lld
",ans);
    for(i=n+1;i<=nn;i+=4)
    {
        ans+=-a[i].ans-a[i+1].ans+a[i+2].ans+a[i+3].ans;
        printf("%lld
",ans);
    }
    return 0;
}