题目：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 ;
}