BZOJ 3295 [CQOI2011]动态逆序对 (三维偏序CDQ+树状数组)

题目大意：

题面传送门

还是一道三维偏序题

每次操作都可以看成这样一个三元组 $<x,w,t>$ ，操作的位置，权值，修改时间

一开始的序列看成n次插入操作

我们先求出不删除时的逆序对总数量，再统计每次删除元素时，减少的逆序对数量

然后就是三维偏序裸题了吧，第一维时间，第二维操作位置，第三维权值，用树状数组维护即可

由于逆序对可以在被删除元素的前面或者后面，所以在归并时需要正反遍历各统计一次

#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N1 150100
#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,nn,K;

struct node{
int p,t,x,w,ans;
}a[N1],tmp[N1];

struct BIT{
int s[N1];
void update(int x,int w){for(int i=x;i<=n;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;

int 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,j,x,cnt;
    
    i=M-1,j=R-1,cnt=0;
    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{
            if(a[j].p==-1) a[j].ans+=b.query(a[j].w-1);
            j--; 
        }
    }
    while(i>=L) {i--;}
    while(j>=M) {if(a[j].p==-1) 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{
            if(a[j].p==-1) a[j].ans+=b.query(n)-b.query(a[j].w);
            tmp[++cnt]=a[j]; j++;
        }
    }
    while(i<M) {tmp[++cnt]=a[i]; i++;}
    while(j<R) {if(a[j].p==-1) a[j].ans+=b.query(n)-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){
    if(s1.p!=s2.p) return s1.p>s2.p;
    return s1.t<s2.t;}
int pos[N1];

int main()
{
    scanf("%d%d",&n,&m);
    int i; ll ans=0; nn=n+m;
    for(i=1;i<=n;i++) a[i].w=gint(),a[i].x=i,a[i].p=1,a[i].t=i,pos[a[i].w]=i;
    for(i=n;i>=1;i--) ans+=b.query(a[i].w),b.update(a[i].w,1);
    memset(b.s,0,sizeof(b.s));
    for(i=n+1;i<=n+m;i++) a[i].w=gint(),a[i].x=pos[a[i].w],a[i].p=-1,a[i].t=i;
    CDQ(1,n+m+1);
    sort(a+1,a+n+m+1,cmp);
    for(i=n+1;i<=n+m;i++)
    {
        printf("%lld
",ans);
        ans-=a[i].ans;
    }
    return 0;
}