标签：ch 洛谷 P3157 int tr pos ret tl CQOI2011

题目大意：

给定\(1\)到\(n\)的一个排列，按照给定顺序依次删除\(m\)个元素，计算每个元素删除之前整个序列的逆序对数量

基本套路：删边变加边

那么我们不就是求满足\(pos_i<pos_j,tim_i<tim_j,num_i>num_j\)的数量嘛

先按\(tim\)排序，然后归并\(pos_i\)，树状数组\(num_i\)

不过这道题我们需要正反跑两个\(cdq\)，因为我们需要分开统计\(pos_i<pos_j,num_i>num_j\)和\(pos_i>pos_j,num_i<num_j\)的贡献

但是我压缩到一个\(cdq\)里了\(emmmm\)

需要稍微注意的一点是，我们需要把答案累加的令一个\(ret_i\)数组中，其中\(ret_i\)表示在\(i\)时刻新产生了多少逆序对，最后还需要输出前缀和

不粘代码是不是太短了

#include<bits/stdc++.h>
using namespace std;
namespace red{
#define int long long
#define mid ((l+r)>>1)
#define lowbit(x) ((x)&(-x))
    inline int read()
    {
        int x=0;char ch,f=1;
        for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar());
        if(ch=='-') f=0,ch=getchar();
        while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
        return f?x:-x;
    }
    const int N=1e5+10;
    int n,m,idx,tot;
    int pos[N];
    int st[N];
    int ret[N];
    struct point
    {
        int x,id,tim;
        int val;
        inline bool operator < (const point &t) const
        {
            if(tim^t.tim) return tim<t.tim;
            return id<t.id;
        }
    }a[N<<2],t[N<<2];
    int tr[N<<1];
    inline void update(int x,int k)
    {
        for(int i=x;i<=n;i+=lowbit(i)) tr[i]+=k;
    }
    inline int query(int y)
    {
        int ret=0;
        for(int i=y;i;i-=lowbit(i))
            ret+=tr[i];
        return ret;
    }
    inline void cdq(int l,int r)
    {
        if(l==r) return;
        cdq(l,mid);
        cdq(mid+1,r);
        int tl=l,tr=mid+1,tot=l;
        while(tl<=mid&&tr<=r)
        {
            if(a[tl].id<=a[tr].id) update(a[tl].x,1),t[tot++]=a[tl++];
            else ret[a[tr].tim]+=query(n)-query(a[tr].x),t[tot++]=a[tr++];
        }
        while(tl<=mid) update(a[tl].x,1),t[tot++]=a[tl++];
        while(tr<=r) ret[a[tr].tim]+=query(n)-query(a[tr].x),t[tot++]=a[tr++];
        for(int i=l;i<=mid;++i) update(a[i].x,-1);
        tl=mid,tr=r;
        while(tl>=l&&tr>=mid+1)
        {
            if(a[tl].id>=a[tr].id) update(a[tl].x,1),--tl;
            else ret[a[tr].tim]+=query(a[tr].x-1),--tr;
        }
        while(tl>=l) update(a[tl].x,1),--tl;
        while(tr>=mid+1) ret[a[tr].tim]+=query(a[tr].x-1),--tr;
        for(int i=l;i<=mid;++i) update(a[i].x,-1);
        for(int i=l;i<=r;++i) a[i]=t[i];
    }
    inline void main()
    {
        n=read(),m=read();
        for(int x,i=1;i<=n;++i)
        {
            x=read();
            pos[x]=i;
            a[i].x=x;
            a[i].id=i;
            a[i].tim=1;
        }
        for(int x,tmp,i=1;i<=m;++i)
        {
            x=read();
            tmp=pos[x];
            a[tmp].tim=m-i+2;
        }
        sort(a+1,a+n+1);
        cdq(1,n);
        for(int i=1;i<=m+1;++i) ret[i]+=ret[i-1];
        for(int i=m+1;i>=2;--i) printf("%lld\n",ret[i]);
    }
}
signed main()
{
    red::main();
    return 0;
}

标签：ch,洛谷,P3157,int,tr,pos,ret,tl,CQOI2011
来源： https://www.cnblogs.com/knife-rose/p/12045586.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。