bzoj3295 [CQOI2011]动态逆序对

http://www.elijahqi.win/2018/01/15/bzoj3295-cqoi2011%e5%8a%a8%e6%80%81%e9%80%86%e5%ba%8f%e5%af%b9/
题目描述
对于序列A，它的逆序对数定义为满足i< j，且Ai>Aj的数对(i,j)的个数。给1到n的一个排列，按照某种顺序依次删除m个元素，你的任务是在每次删除一个元素之前统计整个序列的逆序对数。
输入输出格式
输入格式：

输入第一行包含两个整数n和m，即初始元素的个数和删除的元素个数。以下n行每行包含一个1到n之间的正整数，即初始排列。以下m行每行一个正整数，依次为每次删除的元素。

输出格式：

输出包含m行，依次为删除每个元素之前，逆序对的个数。

输入输出样例
输入样例#1： 复制
5 4
1
5
3
4
2
5
1
4
2
输出样例#1： 复制
5
2
2
1

样例解释
(1,5,3,4,2)(1,3,4,2)(3,4,2)(3,2)(3)。
说明
N<=100000 M<=50000
辣鸡蒟蒻我刚了一晚上的题以为凭借一己之力可以搞定最后还是看了会题解才搞定的 由于每次删数 那么不好办 我们不妨动态的从后往前看 变成一步一步往里面加数 针对每个这样的操作我构造一个三维偏序 这题一开始我一直在往二维数点方向搞 就有点偏了正确的应该往偏序方向走 (t1,x1,y1)(t2,x2,y2)针对这样两个偏序我们可以知道当t1< t2 x1< x2 y1>y2时构成一个逆序对 当t1< t2 x2>x1 y1< y2的时候同样也构成逆序对 那么我不妨构造的时候先以时间排好序 然后按照x的顺序搞一搞 在cdq里做两个归并排序 一个将x从小到大排序 然后每次树状数组统计比我y大的值即可 另一个将x从大到小排序 每次统计比我 y小的值的个数 最后累加起来是我的答案 然后针对答案数组倒序累加一下 因为现在分别是算出了加入这个数后的逆序对个数 所以我要的答案是一个后缀和 分别输出即可

#include<cstdio>
#include<algorithm>
#define N 110000 
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
struct node{
    int x,y,id;
}qr1[160000],qr2[160000],tmp[160000];
bool flag[N];
int pos[N],a[N],n,m,s[N];
inline void add(int x,int v){
    while(x<=n) s[x]+=v,x+=x&-x;
}
inline int query(int x){
    int tmp=0;while(x) tmp+=s[x],x-=x&-x;
    int tmp1=0;x=n;while(x) tmp1+=s[x],x-=x&-x;
    return tmp1-tmp;
}
inline int query1(int x){
    int tmp=0;while(x) tmp+=s[x],x-=x&-x;return tmp;
}
inline void clear(int x){
    for (;x<=n;x+=x&-x) if (!s[x]) return;else s[x]=0;
}
long long ans[N];
inline void cdq(int l,int r){
    if (l==r) return;int mid=l+r>>1;
    cdq(l,mid);cdq(mid+1,r);int h1=l,h2=mid+1,op=l;
    while(h1<=mid&&h2<=r){
        if (qr1[h1].x<qr1[h2].x){
            add(qr1[h1].y,1),tmp[op++]=qr1[h1++];
        }else{
            ans[qr1[h2].id]+=query(qr1[h2].y);tmp[op++]=qr1[h2++]; 
        }
    }
    while(h1<=mid) tmp[op++]=qr1[h1++];
    while(h2<=r) {ans[qr1[h2].id]+=query(qr1[h2].y);tmp[op++]=qr1[h2++];}
    for (int i=l;i<=r;++i){qr1[i]=tmp[i];clear(qr1[i].y);}
    h1=l,h2=mid+1,op=l;
    while(h1<=mid&&h2<=r){
        if (qr2[h1].x>qr2[h2].x){
            add(qr2[h1].y,1),tmp[op++]=qr2[h1++];
        }else{
            ans[qr2[h2].id]+=query1(qr2[h2].y-1);tmp[op++]=qr2[h2++];
        }
    }
    while(h1<=mid) tmp[op++]=qr2[h1++];
    while(h2<=r) ans[qr2[h2].id]+=query1(qr2[h2].y-1),tmp[op++]=qr2[h2++];
    for (int i=l;i<=r;++i) qr2[i]=tmp[i],clear(qr2[i].y);
}
int main(){
    freopen("bzoj3295.in","r",stdin);
    n=read();m=read();
    for (int i=1;i<=n;++i){
        a[i]=read();pos[a[i]]=i;
    }int t=n,mm=m;
    for (int i=1;i<=m;++i){
        qr1[t].id=i;int x=read();qr1[t].x=pos[x];qr1[t].y=x;flag[pos[x]]=1;
        qr2[t].x=pos[x];qr2[t].y=x;qr2[t].id=i;t--;
    }
    for (int i=1;i<=n;++i){
        if (!flag[i]) {
            qr1[t].x=i,qr1[t].y=a[i];qr1[t].id=++m;
            qr2[t].x=i,qr2[t].y=a[i];qr2[t].id=m;t--;
        }
    }
    cdq(1,n);
    for (int i=m-1;~i;--i) ans[i]+=ans[i+1];
    for (int i=1;i<=mm;++i) printf("%lld\n",ans[i]);
    return 0;
}