动态逆序对 HYSBZ - 3295

https://www.lydsy.com/JudgeOnline/problem.php?id=3295

给一个原始序列 每次删一个数 问删这个数之前有多少逆序对

每个数从当前序列消失造成的影响就看 之前有多少数比当前数大 之后有多少数比当前数小 利用主席树做个前缀差即可 但是还有修改 所以要用动态主席树维护

也可以用CDQ分治来做 首先逆向考虑（正向也可以） 整个过程就变为不断向数列内插入数字了 每插入一个数 就看这个数和之前插入的数相比 位置在其之前的有多少比它大的 位置在其之后的有多少比它小的 这样就转换为三维偏序问题了

hysbz的评测机真的是恶心 用主席树一直迷之RE 把结构体里的long long改为int 然后把数组开大了一点就过了

 

主席树

#include <bits/stdc++.h>
using namespace std;
#define ll long long

struct node
{
    int l;
    int r;
    int val;
};

node tree[10000010];
int ary[100010],pos[100010],root[100010],groot[100010];
int n,q,num;

inline int lowbit(int x)
{
    return x&(-x);
}

int build(int l,int r)
{
    int cur,m;
    cur=num++;
    tree[cur].l=0,tree[cur].r=0,tree[cur].val=0;
    if(l==r) return cur;
    m=(l+r)/2;
    tree[cur].l=build(l,m);
    tree[cur].r=build(m+1,r);
    return cur;
}

int updateII(int rot,int tar,int val,int l,int r)
{
    int cur,m;
    cur=num++;
    tree[cur]=tree[rot];
    tree[cur].val+=val;
    if(l==r) return cur;
    m=(l+r)/2;
    if(tar<=m) tree[cur].l=updateII(tree[rot].l,tar,val,l,m);
    else tree[cur].r=updateII(tree[rot].r,tar,val,m+1,r);
    return cur;
}

void updateI(int tar,int val)
{
    int i;
    for(i=tar;i<=n;i+=lowbit(i)) groot[i]=updateII(groot[i],val,-1,1,n);
}

int queryII(int rot,int pl,int pr,int l,int r)
{
    int res;
    int m;
    if(pl<=l&&r<=pr) return tree[rot].val;
    res=0,m=(l+r)/2;
    if(pl<=m) res+=queryII(tree[rot].l,pl,pr,l,m);
    if(pr>m) res+=queryII(tree[rot].r,pl,pr,m+1,r);
    return res;
}

ll queryI(int tar,int val)
{
    ll res,tmp;
    int i;
    res=0;
    if(val<n)
    {
        tmp=queryII(root[tar-1],val+1,n,1,n);
        res+=tmp;
        for(i=tar-1;i>=1;i-=lowbit(i))
        {
            tmp=queryII(groot[i],val+1,n,1,n);
            res+=tmp;
        }
    }
    if(val>1)
    {
        tmp=queryII(root[n],1,val-1,1,n)-queryII(root[tar],1,val-1,1,n);
        res+=tmp;
        for(i=n;i>=1;i-=lowbit(i))
        {
            tmp=queryII(groot[i],1,val-1,1,n);
            res+=tmp;
        }
        for(i=tar;i>=1;i-=lowbit(i))
        {
            tmp=queryII(groot[i],1,val-1,1,n);
            res-=tmp;
        }
    }
    return res;
}

int main()
{
    ll ans;
    int i,val;
    scanf("%d%d",&n,&q);
    for(i=1;i<=n;i++)
    {
        scanf("%d",&ary[i]);
        pos[ary[i]]=i;
    }
    num=0;
    root[0]=build(1,n);
    for(i=1;i<=n;i++) root[i]=updateII(root[i-1],ary[i],1,1,n);
    ans=0;
    for(i=1;i<=n;i++)
    {
        if(ary[i]+1<=n) ans+=queryII(root[i],ary[i]+1,n,1,n);
    }
    for(i=1;i<=n;i++) groot[i]=root[0];
    while(q--)
    {
        printf("%lld\n",ans);
        scanf("%d",&val);
        updateI(pos[val],val);
        ans-=queryI(pos[val],val);
    }
    return 0;
}

 

CDQ分治

#include <bits/stdc++.h>
using namespace std;
#define ll long long

struct node1
{
    int tp;
    int x;
    int y;
    int val;
};

struct node2
{
    int l;
    int r;
    ll val;
};

node1 pre[500010],order[500010],tmp[500010];
node2 tree[400010];
ll ans[50010];
int ary[100010],val[100010],book[100010],pos[100010];
int n,q,tot;

void build(int l,int r,int cur)
{
    int m;
    tree[cur].l=l;
    tree[cur].r=r;
    tree[cur].val=0;
    if(l==r) return;
    m=(l+r)/2;
    build(l,m,2*cur);
    build(m+1,r,2*cur+1);
}

void update(int tar,int val,int cur)
{
    tree[cur].val+=val;
    if(tree[cur].l==tree[cur].r) return;
    if(tar<=tree[2*cur].r) update(tar,val,2*cur);
    else update(tar,val,2*cur+1);
}

ll query(int pl,int pr,int cur)
{
    ll res;
    if(pl<=tree[cur].l&&tree[cur].r<=pr) return tree[cur].val;
    res=0;
    if(pl<=tree[2*cur].r) res+=query(pl,pr,2*cur);
    if(pr>=tree[2*cur+1].l) res+=query(pl,pr,2*cur+1);
    return res;
}

void cdqI(int l,int r)
{
    int m,p,q,cnt,i;
    if(l==r) return;
    m=(l+r)/2;
    cdqI(l,m);
    cdqI(m+1,r);
    p=l,q=m+1,cnt=l;
    while(p<=m&&q<=r)
    {
        if(order[p].x<=order[q].x)
        {
            if(order[p].tp==1) update(order[p].y,order[p].val,1);
            tmp[cnt++]=order[p++];
        }
        else
        {
            if(order[q].tp==2&&order[q].y<n) ans[order[q].val]+=query(order[q].y+1,n,1);
            tmp[cnt++]=order[q++];
        }
    }
    while(q<=r)
    {
        if(order[q].tp==2&&order[q].y<n) ans[order[q].val]+=query(order[q].y+1,n,1);
        tmp[cnt++]=order[q++];
    }
    for(i=l;i<p;i++)
    {
        if(order[i].tp==1) update(order[i].y,-order[i].val,1);
    }
    while(p<=m) tmp[cnt++]=order[p++];
    for(i=l;i<=r;i++) order[i]=tmp[i];
}

void cdqII(int l,int r)
{
    int m,p,q,cnt,i;
    if(l==r) return;
    m=(l+r)/2;
    cdqII(l,m);
    cdqII(m+1,r);
    p=l,q=m+1,cnt=l;
    while(p<=m&&q<=r)
    {
        if(order[p].x>=order[q].x)
        {
            if(order[p].tp==1) update(order[p].y,order[p].val,1);
            tmp[cnt++]=order[p++];
        }
        else
        {
            if(order[q].tp==2&&order[q].y>1) ans[order[q].val]+=query(1,order[q].y-1,1);
            tmp[cnt++]=order[q++];
        }
    }
    while(q<=r)
    {
        if(order[q].tp==2&&order[q].y>1) ans[order[q].val]+=query(1,order[q].y-1,1);
        tmp[cnt++]=order[q++];
    }
    for(i=l;i<p;i++)
    {
        if(order[i].tp==1) update(order[i].y,-order[i].val,1);
    }
    while(p<=m) tmp[cnt++]=order[p++];
    for(i=l;i<=r;i++) order[i]=tmp[i];
}

int main()
{
    int i;
    scanf("%d%d",&n,&q);
    for(i=1;i<=n;i++)
    {
        scanf("%d",&ary[i]);
        pos[ary[i]]=i;
    }
    for(i=1;i<=q;i++)
    {
        scanf("%d",&val[i]);
        book[val[i]]=1;
    }
    tot=0;
    for(i=1;i<=n;i++)
    {
        if(book[i]==0)
        {
            tot++;
            order[tot].tp=1,order[tot].x=pos[i],order[tot].y=i,order[tot].val=1;
            tot++;
            order[tot].tp=2,order[tot].x=pos[i],order[tot].y=i,order[tot].val=q;
        }
    }
    for(i=q;i>=1;i--)
    {
        tot++;
        order[tot].tp=1,order[tot].x=pos[val[i]],order[tot].y=val[i],order[tot].val=1;
        tot++;
        order[tot].tp=2,order[tot].x=pos[val[i]],order[tot].y=val[i],order[tot].val=i;
    }
    for(i=1;i<=tot;i++) pre[i]=order[i];
    build(1,n,1);
    cdqI(1,tot);
    build(1,n,1);
    for(i=1;i<=tot;i++) order[i]=pre[i];
    cdqII(1,tot);
    for(i=q;i>=1;i--) ans[i]+=ans[i+1];
    for(i=1;i<=q;i++) printf("%lld\n",ans[i]);
    return 0;
}