【BZOJ3295】【CQOI2011】动态逆序对

cdq分治经典例题，然而智商掉线傻逼错误坑了两天

原题：

对于序列A，它的逆序对数定义为满足i<j，且Ai>Aj的数对(i,j)的个数。给1到n的一个排列，按照某种顺序依次删除m个元素，你的任务是在每次删除一个元素之前统计整个序列的逆序对数。

N<=100000 M<=50000

此题修改和询问绑定完全离线，可以直接倒序变成插入，然后就是三维数星星辣(ง •̀∀•́)ง

x为下标，y为值，z为时间轴，x排序，z_cdq分治，y树状数组

不用搞两次cdq分治，每次按x递增查找比y大的，然后按x递减查找比y小的

注意当x递增的时候查找的是比y大的，因为虽然x是递增，但是修改上去的都是比当前x小的

然后我先计算贡献再把小于mid的放左边大于mid的放右边然后进入下一层的做法是没有问题的

问题在于我傻逼题中删除的是值我删的是下标

注意longlong

代码：

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<queue>
 #include<vector>
 using namespace std;
 #define ll long long
 int read(){int z=,mark=;  char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')mark=-;  ch=getchar();}
     while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-'';  ch=getchar();}
     return z*mark;
 }
 struct cdd{int x,y,z;}a[],q[];
 int n,m,b[];  int N;
 bool usd[];  int tt=;
 ll e[];  int lbt[];
 ll ans[];
 int re[];
 void gtlbt(){for(int i=;i<=n;++i)  lbt[i]=i&-i;}
 void mdf(int x,int y){while(x<=n)  e[x]+=y,x+=lbt[x];}
 ll qr(int x){ll bwl=;  while(x)  bwl+=e[x],x-=lbt[x];  return bwl;}
 void cdq(int l,int r){
     if(l==r)  return ;
     int md=(l+r)>>;
     for(int i=l;i<=r;++i){
         if(a[i].z>md)  mdf(a[i].y,);
         else  ans[a[i].z]+=qr(n)-qr(a[i].y);
     }
     for(int i=l;i<=r;++i)if(a[i].z>md)  mdf(a[i].y,-);
     for(int i=r;i>=l;--i){
         if(a[i].z>md)  mdf(a[i].y,);
         else  ans[a[i].z]+=qr(a[i].y);
     }
     for(int i=l;i<=r;++i)if(a[i].z>md)  mdf(a[i].y,-);
     int t1=l,t2=md+;
     for(int i=l;i<=r;++i)  q[(a[i].z<=md?t1:t2)++]=a[i];
     for(int i=l;i<=r;++i)  a[i]=q[i];
     cdq(l,md),cdq(md+,r);
 }
 bool cmp(cdd x,cdd y){return x.x<y.x;}
 int main(){//freopen("ddd.in","r",stdin);
     memset(usd,,sizeof(usd));
     cin>>n>>m;  N=n-m;
     gtlbt();
     for(int i=;i<=n;++i)  b[i]=read(),re[b[i]]=i;
     for(int i=;i<=m;++i)  a[++tt].y=read(),a[i].x=re[a[i].y],a[i].z=i,usd[a[i].x]=true;
     for(int i=;i<=n;++i)if(!usd[i])  a[++tt].x=i,a[tt].y=b[i],a[tt].z=tt;
     sort(a+,a+n+,cmp);
     cdq(,n);
     for(int i=n;i>;--i)  ans[i-]+=ans[i];
     for(int i=;i<=m;++i)  printf("%I64d\n",ans[i]);
     return ;
 }