洛谷P3759 - [TJOI2017]不勤劳的图书管理员

Description

给出一个$1..n(n\leq5\times10^4)$的排列$\{a_n\}$和数列$\{w_n\}(w_i\leq10^5)$，进行$m(m\leq5\times10^4)$次操作：

交换$a_{p_1},a_{p_2}$，并求$\sum_{i=1}^n \sum_{j=i+1}^n [a_i>a_j](w_{a_i}+w_{a_j})$。

Solution

树套树/CDQ分治，想锻炼一下代码能力所以写了树套树。

首先这是一个求逆序对的问题，那么我们考虑交换$a_{p_1},a_{p_2}$对答案有什么影响。易知只有$i\in[p_1+1,p_2-1]$对答案造成影响：

若$a_i<a_{p_1}$，那么答案减$a_i+a_{p_1}$；若$a_i>a_{p_1}$，那么答案加$a_i+a_{p_1}$。

若$a_i<a_{p_2}$，那么答案加$a_i+a_{p_2}$；若$a_i>a_{p_2}$，那么答案减$a_i+a_{p_2}$。

那么求出$\{c_1,s_1,c_2,s_2\}$分别表示区间$[p_1+1,p_2-1]$内小于/大于$x$的$i$的个数/$w_i$的和，然后就可以根据上两行计算答案的变化。而这可以用树套树来做。

时间复杂度$O(mlog^2n)$。

Code

//[TJOI2017]不勤劳的图书管理员

#include <algorithm>

#include <cstdio>

using std::swap;

typedef long long lint;

inline char gc()

{

    static char now[1<<16],*s,*t;

    if(s==t) {t=(s=now)+fread(now,1,1<<16,stdin); if(s==t) return EOF;}

    return *s++;

}

inline int read()

{

    int x=0; char ch=gc();

    while(ch<'0'||'9'<ch) ch=gc();

    while('0'<=ch&&ch<='9') x=x*10+ch-'0',ch=gc();

    return x;

}

const int N=5e4+10;

const int P=1e9+7;

int n,m,a[N],w[N];

namespace init

{

    int tr1[N],tr2[N];

    void add(int x,int v) {while(x<=n) tr1[x]++,tr2[x]=(tr2[x]+v)%P,x+=x&(-x);}

    int sum1(int x) {int r=0; while(x) r+=tr1[x],x-=x&(-x); return r;}

    int sum2(int x) {int r=0; while(x) r=(r+tr2[x])%P,x-=x&(-x); return r;}

    int calc()

    {

        int res=0;

        for(int i=1;i<=n;i++)

        {

            res=(res+(lint)sum1(n-a[i]+1)*w[a[i]])%P;

            res=(res+sum2(n-a[i]+1))%P;

            add(n-a[i]+1,w[a[i]]);

        }

        return res;

    }

}

const int N1=2e6;

struct info

{

    int c1,s1,c2,s2;

    info(int _c1=0,int _s1=0,int _c2=0,int _s2=0) {c1=_c1,s1=_s1,c2=_c2,s2=_s2;};

    friend info operator +(info x,info y)

    {

        int c1=x.c1+y.c1,s1=x.s1+y.s1,c2=x.c2+y.c2,s2=x.s2+y.s2;

        return info(c1,s1%P,c2,s2%P);

    }

};

namespace inTr

{

    int cnt,fa[N1],ch[N1][2],siz[N1]; int val[N1],sum[N1];

    inline void update(int p)

    {

        siz[p]=siz[ch[p][0]]+1+siz[ch[p][1]];

        sum[p]=((sum[ch[p][0]]+sum[ch[p][1]])%P+w[val[p]])%P;

    }

    inline int wh(int p) {return p==ch[fa[p]][1];}

    inline void rotate(int p)

    {

        int q=fa[p],r=fa[q],w=wh(p);

        fa[p]=r; if(r) ch[r][wh(q)]=p;

        fa[ch[q][w]=ch[p][w^1]]=q;

        fa[ch[p][w^1]=q]=p;

        update(q),update(p);

    }

    void splay(int &rt,int p)

    {

        for(int q=fa[p];q;rotate(p),q=fa[p]) if(fa[q]) rotate(wh(p)^wh(q)?p:q);

        update(rt=p);

    }

    int pre(int rt,int x)

    {

        int r=0;

        for(int p=rt;p;p=ch[p][val[p]<x]) if(val[p]<x) r=p;

        return r;

    }

    void ins(int &rt,int p)

    {

        if(!rt) {rt=p; return;}

        int q=rt;

        while(ch[q][val[q]<val[p]]) q=ch[q][val[q]<val[p]];

        fa[ch[q][val[q]<val[p]]=p]=q;

        update(q); splay(rt,p);

    }

    int del(int &rt,int x)

    {

        if(x==0||!rt) return 0;

        int p=pre(rt,x+1); if(val[p]!=x) return 0;

        splay(rt,p);

        int chCnt=(ch[p][0]>0)+(ch[p][1]>0);

        if(chCnt==0) {rt=0; return p;}

        if(chCnt==1) {rt=ch[p][ch[p][1]>0],fa[rt]=0; return p;}

        int q=ch[p][0]; while(ch[q][1]) q=ch[q][1];

        splay(rt,q);

        fa[ch[q][1]=ch[p][1]]=q;

        fa[p]=ch[p][0]=ch[p][1]=0;

        update(q);

        return p;

    }

    inline void change(int &rt,int x1,int x2)

    {

        int p=del(rt,x1); if(!p) p=++cnt;

        fa[p]=ch[p][0]=ch[p][1]=0,siz[p]=1;

        val[p]=x2,sum[p]=w[x2];

        ins(rt,p);

    }

    info query(int rt,int x)

    {

        int c1=0,s1=0;

        for(int p=rt;p;p=ch[p][val[p]<x])

            if(val[p]<x) c1+=siz[p]-siz[ch[p][1]],s1=(s1+sum[p]-sum[ch[p][1]]+P)%P;

        return info(c1,s1,siz[rt]-c1,sum[rt]-s1);

    }

}

namespace outTr

{

    #define Ls (p<<1)

    #define Rs (p<<1|1)

    int rt[N<<2];

    int in(int L,int R,int p1,int p2) {return ((L<=p2&&p2<=R)<<1)|(L<=p1&&p1<=R);}

    void change(int p,int L0,int R0,int p1,int p2)

    {

        int inP=in(L0,R0,p1,p2);

        if(inP==1) inTr::change(rt[p],a[p1],a[p2]);

        else if(inP==2) inTr::change(rt[p],a[p2],a[p1]);

        if(L0==R0) return;

        int mid=L0+R0>>1;

        if(in(L0,mid,p1,p2)) change(Ls,L0,mid,p1,p2);

        if(in(mid+1,R0,p1,p2)) change(Rs,mid+1,R0,p1,p2);

    }

    int optL,optR;

    info query(int p,int L0,int R0,int x)

    {

        if(optL<=L0&&R0<=optR) return inTr::query(rt[p],x);

        int mid=L0+R0>>1; info res=info(0,0,0,0);

        if(optL<=mid) res=res+query(Ls,L0,mid,x);

        if(mid<optR) res=res+query(Rs,mid+1,R0,x);

        return res;

    }

}

inline void change(int p1,int p2) {outTr::change(1,1,n,p1,p2);}

inline info query(int L,int R,int x)

{

    outTr::optL=L,outTr::optR=R;

    return outTr::query(1,1,n,x);

}

int main()

{

    n=read(),m=read();

    for(int i=1;i<=n;i++)

    {

        a[0]=read(),w[a[0]]=read();

        change(0,i); swap(a[0],a[i]);

    }

    int ans=init::calc();

    for(int i=1;i<=m;i++)

    {

        int p1=read(),p2=read();

        if(p1>p2) swap(p1,p2); else if(p1==p2) {printf("%d\n",ans); continue;}

        int x=a[p1],y=a[p2];

        info r1=query(p1+1,p2-1,x),r2=query(p1+1,p2-1,y);

        ans-=((lint)r1.c1*w[x]%P+r1.s1)%P; ans=(ans+P)%P;

        ans+=((lint)r1.c2*w[x]%P+r1.s2)%P; ans=ans%P;

        ans+=((lint)r2.c1*w[y]%P+r2.s1)%P; ans=ans%P;

        ans-=((lint)r2.c2*w[y]%P+r2.s2)%P; ans=(ans+P)%P;

        if(x<y) ans=(ans+w[x]+w[y])%P; else ans=(ans-w[x]-w[y]+P)%P;

        printf("%d\n",ans);

        change(p1,p2); swap(a[p1],a[p2]);

    }

    return 0;

}