[BJOI2017]开车

[BJOI2017]开车

直接做要用栈

修改？难以直接维护

统计边的贡献！

len*abs(pre)pre表示前缀car-stop

修改时候，整个区间的pre+1或者-1

分块，块内对pre排序并打标记

二分出0的位置，再根据pre的正负和本次是+1，-1，leni贡献+1还是-1的贡献

O(nsqrt(n)logn)

#include<bits/stdc++.h>
#define reg register int
#define il inline
#define fi first
#define se second
#define mk(a,b) make_pair(a,b)
#define numb (ch^'0')
using namespace std;
typedef long long ll;
template<class T>il void rd(T &x){
    char ch;x=;bool fl=false;
    while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
    for(x=numb;isdigit(ch=getchar());x=x*+numb);
    (fl==true)&&(x=-x);
}
template<class T>il void output(T x){if(x/)output(x/);putchar(x%+'');}
template<class T>il void ot(T x){if(x<) putchar('-'),x=-x;output(x);putchar(' ');}
template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');}

namespace Miracle{
const int N=*5e4+;
const int B=;
int n,m;
int L[N],R[N];
struct qs{
    int id,x;
}q[N];
int c[N],cnt;
int tot;
int pos[N],b[N];
int be[N],len[N];

vector<int>mem[B],sum[B];
int tag[B];
ll con[B];

int exi[N];//0: none 1:car -1:stop
int pre[N];
int buc[N];
bool cmp(int x,int y){
    return pre[x]<pre[y];
}
int li(int x){
    return lower_bound(c+,c+cnt+,x)-c;
}
int main(){
    rd(n);
    for(reg i=;i<=n;++i) rd(pos[i]),c[++cnt]=pos[i];
    for(reg i=;i<=n;++i) rd(b[i]),c[++cnt]=b[i];
    rd(m);
    for(reg i=;i<=m;++i){
        rd(q[i].id);rd(q[i].x);
        c[++cnt]=q[i].x;
    }
    sort(c+,c+cnt+);
    cnt=unique(c+,c+cnt+)-c-;

    int blo=sqrt(cnt)+;
    for(reg i=;i<=cnt;++i){
        be[i]=(i-)/blo+;
        if(!L[be[i]]) L[be[i]]=i;
        R[be[i]]=i;
        len[i]=c[i+]-c[i];
    }
    len[cnt]=;
    tot=be[cnt];

    for(reg i=;i<=n;++i){
        int po=lower_bound(c+,c+cnt+,b[i])-c;
        --exi[po];
        po=lower_bound(c+,c+cnt+,pos[i])-c;
        ++exi[po];
    }
    ll ans=;
    for(reg i=;i<=cnt;++i){
        pre[i]=pre[i-]+exi[i];
        ans+=(ll)len[i]*abs(pre[i]);
    }
    for(reg j=;j<=tot;++j){
        tag[j]=;
        for(reg i=L[j];i<=R[j];++i){
            mem[j].push_back(i);
            con[j]+=(ll)len[i]*abs(pre[i]);
        }
        sort(mem[j].begin(),mem[j].end(),cmp);
        sum[j].resize(mem[j].size());
        for(reg i=;i<(int)mem[j].size();++i){
            sum[j][i]=len[mem[j][i]];
            if(i) sum[j][i]+=sum[j][i-];
        }
    }
    printf("%lld\n",ans);
    for(reg o=;o<=m;++o){
        int l=li(pos[q[o].id]),r=li(q[o].x);
        exi[l]--;exi[r]++;

        int c=-;//l->r
        if(l>r) swap(l,r),c=;//r->l
        --r;//warinig!!!
        if(be[l]==be[r]){
            for(reg i=l;i<=r;++i) pre[i]+=c;
            ans-=con[be[l]];
            con[be[l]]=;
            mem[be[l]].clear();
            sum[be[l]].clear();
            int j=be[l];
            for(reg i=L[be[l]];i<=R[be[l]];++i){
                pre[i]+=tag[be[l]];
                mem[j].push_back(i);
            }
            sort(mem[j].begin(),mem[j].end(),cmp);
            sum[j].resize(mem[j].size());
            for(reg i=;i<(int)mem[j].size();++i){
                sum[j][i]=len[mem[j][i]];
                con[j]+=(ll)len[mem[j][i]]*abs(pre[mem[j][i]]);
                if(i) sum[j][i]+=sum[j][i-];
            }
            tag[be[l]]=;
            ans+=con[be[l]];
        }else{
            for(reg i=l;i<=R[be[l]];++i){
                pre[i]+=c;
            }
            int j=be[l];
            ans-=con[j];
            con[j]=;
            mem[j].clear();
            sum[j].clear();
            for(reg i=L[j];i<=R[j];++i){
                pre[i]+=tag[j];
                mem[j].push_back(i);
            }
            sort(mem[j].begin(),mem[j].end(),cmp);
            sum[j].resize(mem[j].size());
            for(reg i=;i<(int)mem[j].size();++i){
                sum[j][i]=len[mem[j][i]];
                con[j]+=(ll)len[mem[j][i]]*abs(pre[mem[j][i]]);
                if(i) sum[j][i]+=sum[j][i-];
            }
            tag[j]=;
            ans+=con[j];

            for(reg i=L[be[r]];i<=r;++i){
                pre[i]+=c;
            }
            j=be[r];
            ans-=con[j];
            con[j]=;
            mem[j].clear();
            sum[j].clear();
            for(reg i=L[j];i<=R[j];++i){
                pre[i]+=tag[j];
                mem[j].push_back(i);
            }
            sort(mem[j].begin(),mem[j].end(),cmp);
            sum[j].resize(mem[j].size());
            for(reg i=;i<(int)mem[j].size();++i){
                sum[j][i]=len[mem[j][i]];
                con[j]+=(ll)len[mem[j][i]]*abs(pre[mem[j][i]]);
                if(i) sum[j][i]+=sum[j][i-];
            }
            tag[j]=;
            ans+=con[j];

            for(reg j=be[l]+;j<=be[r]-;++j){
                pre[]=-tag[j];
                if(c==-){
                    int lp=upper_bound(mem[j].begin(),mem[j].end(),,cmp)-mem[j].begin();
                    --lp;//warning!! maybe -1
                    if(lp!=-){
                        ans+=(ll)sum[j][lp];
                        ans-=(ll)sum[j][mem[j].size()-]-sum[j][lp];
                        con[j]+=(ll)sum[j][lp];
                        con[j]-=(ll)sum[j][mem[j].size()-]-sum[j][lp];
                    }else{
                        ans-=(ll)sum[j][mem[j].size()-];
                        con[j]-=(ll)sum[j][mem[j].size()-];
                    }

                }else{
                    int lp=lower_bound(mem[j].begin(),mem[j].end(),,cmp)-mem[j].begin();
                    --lp;//warning!! maybe -1
                    if(lp!=-){
                        ans-=(ll)sum[j][lp];
                        ans+=(ll)sum[j][mem[j].size()-]-sum[j][lp];
                        con[j]-=(ll)sum[j][lp];
                        con[j]+=(ll)sum[j][mem[j].size()-]-sum[j][lp];
                    }else{
                        ans+=(ll)sum[j][mem[j].size()-];
                        con[j]+=(ll)sum[j][mem[j].size()-];
                    }
                }
                tag[j]+=c;
            }
            pre[]=;
        }
        pos[q[o].id]=q[o].x;
        printf("%lld\n",ans);
    }
    return ;
}

}
signed main(){
    Miracle::main();
    return ;
}

/*
   Author: *Miracle*
   Date: 2019/4/11 10:31:35
*/