bzoj 3745 [Coci2015]Norma——序列分治

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=3745

如果分治，就能在本层仅算过 mid 的区间了。

可以从中间到左边地遍历左边，给右边两个指针，表示第一个更新左边造成的最小值/最大值的位置。

两个位置共同的左边可以公式算长度，用左边的最值算；两个位置共同的右边可以预处理，处理出 算上长度(相对mid的)的最值乘积求和 与 不算长度的最值乘积求和（都是前缀），把前者加到答案里，后者乘上左边到mid的长度加到答案里即可；两个位置夹着的位置判断一下可以用左边的最大值还是最小值，所以要预处理右边最大值/最小值的算长度/不算长度前缀和，然后和公共右边一样的处理方法即可。

取模的地方要注意！可能加了一个东西，就不能再+mod、upd( )，而要直接upd( )。

注意把 l-1 位置的各种值赋成0。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=5e5+,mod=1e9;
int n,a[N],ans,s0[N],s1[N],s2[N];//s:可加入ans，相对长度
int ml0[N],ml1[N],ml2[N];//ml:单纯乘积相加
int rdn()
{
    int ret=;bool fx=;char ch=getchar();
    while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
    while(ch>=''&&ch<='') ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
    return fx?ret:-ret;
}
void upd(int &x){x-=(x>=mod?mod:);}
ll calc(int a){return (ll)a*(a+)>>1ll;}
int dis(int x,int y){return y-x+;}
void solve(int l,int r)
{
    if(l==r)
    {
        s2[l]=ml2[l]=(ll)a[l]*a[l]%mod;
        s1[l]=ml1[l]=s0[l]=ml0[l]=a[l];
        s2[l-]=s1[l-]=s0[l-]=ml2[l-]=ml1[l-]=ml0[l-]=;//
        ans+=s2[l]; upd(ans);
        return;
    }
    int mid=l+r>>;
    solve(l,mid); solve(mid+,r);

    int lo=a[mid],hi=a[mid],p0=mid+,p1=mid+;
    for(int i=mid,cd=;i>=l;i--,cd++)
    {
        hi=max(hi,a[i]); lo=min(lo,a[i]);
        while(a[p0]<=hi&&p0<=r) p0++;
        while(a[p1]>=lo&&p1<=r) p1++;
        int tl=min(p0,p1)-,tr=max(p0,p1);
//        printf("l=%d r=%d hi=%d lo=%d p0=%d p1=%d\n",l,r,hi,lo,p0,p1);

        int d=dis(mid+,tl);
        ans=(ans+( (ll)cd*d+calc(d) )%mod*lo%mod*hi)%mod;
//        printf("ans=%d ",ans);

        ans=(ans+s2[r]-s2[tr-])%mod+mod; upd(ans);
        ans=ans+(ll)(ml2[r]-ml2[tr-])*cd%mod;//cd not dis(i,tr-1)
        if(ans<) ans+=mod; else upd(ans);//if
//        printf("ans=%d ",ans);

        if(p1<p0)//最小值已更新
        {
            ans=ans+(ll)hi*(s1[tr-]-s1[tl])%mod;
            if(ans<)ans+=mod; else upd(ans);
            ans=(ans+(ll)(ml1[tr-]-ml1[tl])*cd%mod*hi%mod);//cd
            if(ans<) ans+=mod; else upd(ans);//if
        }
        if(p0<p1)//最大值已更新
        {
//            printf("tl=%d tr=%d dis(i,tl)=%d s0[%d]-s0[%d]=%d lo=%d mml=%d\n",
//            tl,tr,dis(i,tl),tr-1,tl,s0[tr-1]-s0[tl],lo,ml0[tr-1]-ml0[tl]);
            ans=ans+(ll)lo*(s0[tr-]-s0[tl])%mod;
            if(ans<)ans+=mod; else upd(ans);
            ans=(ans+(ll)(ml0[tr-]-ml0[tl])*cd%mod*lo%mod);
            if(ans<) ans+=mod; else upd(ans);//if
        }
    }

    hi=lo=a[l];
    s2[l]=ml2[l]=(ll)hi*lo%mod;  s1[l]=ml1[l]=lo; s0[l]=ml0[l]=hi;
    for(int i=l+,d=;i<=r;i++,d++)
    {
        hi=max(hi,a[i]); lo=min(lo,a[i]);
        s2[i]=s2[i-]+(ll)d*hi%mod*lo%mod; upd(s2[i]);
        s1[i]=s1[i-]+(ll)d*lo%mod; upd(s1[i]);
        s0[i]=s0[i-]+(ll)d*hi%mod; upd(s0[i]);

        ml2[i]=ml2[i-]+(ll)hi*lo%mod; upd(ml2[i]);
        ml1[i]=ml1[i-]+lo; upd(ml1[i]);
        ml0[i]=ml0[i-]+hi; upd(ml0[i]);
    }
    s2[l-]=s1[l-]=s0[l-]=ml2[l-]=ml1[l-]=ml0[l-]=;//
}
int main()
{
    n=rdn();
    for(int i=;i<=n;i++)a[i]=rdn();
    solve(,n);
    printf("%d\n",ans);
    return ;
}