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=,fh=;char c=getchar();
     while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
     while(c>=''&&c<=''){ret=ret*+c-'';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=; 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<=) return;
     int M=(L+R)>>;
     CDQ(L,M); CDQ(M,R);
     int i,j,x,cnt;

     i=M-,j=R-,cnt=;
     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==-) a[j].ans+=b.query(a[j].w-);
             j--;
         }
     }
     while(i>=L) {i--;}
     while(j>=M) {if(a[j].p==-) a[j].ans+=b.query(a[j].w-); j--;}
     while(tl) {x=que[tl--]; b.update(a[x].w,-a[x].p);}

     i=L,j=M,cnt=;
     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==-) 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==-) 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+];
 }
 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=; nn=n+m;
     for(i=;i<=n;i++) a[i].w=gint(),a[i].x=i,a[i].p=,a[i].t=i,pos[a[i].w]=i;
     for(i=n;i>=;i--) ans+=b.query(a[i].w),b.update(a[i].w,);
     memset(b.s,,sizeof(b.s));
     for(i=n+;i<=n+m;i++) a[i].w=gint(),a[i].x=pos[a[i].w],a[i].p=-,a[i].t=i;
     CDQ(,n+m+);
     sort(a+,a+n+m+,cmp);
     for(i=n+;i<=n+m;i++)
     {
         printf("%lld\n",ans);
         ans-=a[i].ans;
     }
     return ;
 }