bzoj3295 [Cqoi2011]动态逆序对

Description

---

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

N<=100000 M<=50000

Solution

---

很裸的树套树呀，离线询问动态开点连离散都不用了
本来想打cdq分治但是时间不太够了

以后这种常数巨大的题目最好还是读入优化卡一卡

Code

---

#include <stdio.h>
#include <string.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)
#define lowbit(x) ((x)&(-(x)))

typedef long long LL;
const int N=200005;

struct treeNode {int l,r,sum;} t[N*101];

int root[N],tot;
int a[N],b[N],num[N],n,m;
int pos[N];
LL prt[N];

bool vis[N];

int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}

void modify(int &now,int tl,int tr,int x) {
    if (!now) t[now=++tot]=(treeNode) {0,0,0};
    t[now].sum++;
    if (tl==tr) return ;
    int mid=(tl+tr)>>1;
    if (x<=mid) modify(t[now].l,tl,mid,x);
    else modify(t[now].r,mid+1,tr,x);
}

int query(int now,int tl,int tr,int l,int r) {
    if (r<l) return 0;
    if (tl==l&&tr==r) return t[now].sum;
    int mid=(tl+tr)>>1;
    if (r<=mid) return query(t[now].l,tl,mid,l,r);
    if (l>mid) return query(t[now].r,mid+1,tr,l,r);
    int qx=query(t[now].l,tl,mid,l,mid);
    int qy=query(t[now].r,mid+1,tr,mid+1,r);
    return qx+qy;
}

int get(int x,int l,int r) {
    int ret=0;
    for (;x;x-=lowbit(x)) ret+=query(root[x],1,n,l,r);
    return ret;
}

int ins(int x,int v) {
    for (int i=x;i<=n;i+=lowbit(i)) modify(root[i],1,n,v);
    num[x]=v;
    int ret=get(x,v+1,n);
    ret+=get(n,1,v-1)-get(x,1,v-1);
    return ret;
}

int main(void) {
    n=read(); m=read();
    rep(i,1,n) {
        a[i]=read();
        pos[a[i]]=i;
    }
    rep(i,1,m) {
        b[i]=read();
        vis[b[i]]=1;
    }
    LL ans=0;
    rep(i,1,n) if (!vis[a[i]]) ans=ans+(LL)ins(i,a[i]);
    drp(i,m,1) {
        int delta=ins(pos[b[i]],b[i]);
        ans=ans+(LL)delta;
        prt[i]=ans;
    }
    rep(i,1,m) printf("%lld\n", prt[i]);
    return 0;
}