Description
给出数列A,求 ∑ni=1∑nj=imaxjk=i∗minjk=i∗(j−i+1)
Solution
怎么看都像数据结构或者贪心的题目,然而正解却是分治。
这个分治就是分而治之,并没有归并。
首先对于区间
[l,r],mid=l+r/2
,我们在区间
[l,mid]
里枚举一个点,我们可以求出lMx=
maxmidi=l
以及IMn=
minmidi=l
(只要
O(n)
的从mid开始向l枚举即可),那么设当前左边枚举到
i
,然后对于区间
区间
[mid+1,L)
表示
maxLi=l
=lMx并且
minLi=1
=lMn
我们可以轻松的计算其值。
然后是区间
[L,R)
表示
maxLi=l
=lMx或者
minLi=1
=lMn,这两个不会同时满足因为同时满足的区间会并到第一种情况那里,然后无论是
maxLi=l
还是
minLi=1
都单调的。这个我们可以通过巧妙的前缀和搞定。
最后是区间
[R,r]
表示
maxLi=l
>lMx并且
minLi=1
#include<iostream>
#include<algorithm>
using namespace std;
const int INF=(unsigned)-1>>1,M=1e5+5,P=1e9;
int A[M],n,sMx[M],sMn[M],cMx[M],cMn[M],cM[M],sM[M];
inline void rd(int &a){
a=0;char c;
while(c=getchar(),!isdigit(c));
do a=a*10+(c^48);
while(c=getchar(),isdigit(c));
}
inline void Max(int &a,int b){if(a<b)a=b;}
inline void Min(int &a,int b){if(a>b)a=b;}
inline void Mod_add(int &a,int b){if((a+=b)>=P)a-=P;}
inline int dis(int l1,int l2){return 1ll*((l2-l1+1)*(l1+l2)>>1)%P;}
int solve(int l,int r){
if(l==r)return 1ll*A[l]*A[l]%P;
int mid=l+r>>1,ans;if((ans=solve(l,mid)+solve(mid+1,r))>=P)ans-=P;
for(int i=mid+1,rMx=0,rMn=INF,len;len=i-mid,i<=r;++i){
Max(rMx,A[i]),Min(rMn,A[i]),sMx[i]=rMx,sMn[i]=rMn;
cMx[i]=1ll*rMx*len%P,cMn[i]=1ll*rMn*len%P;
cM[i]=1ll*len*rMx*rMn%P,sM[i]=1ll*rMx*rMn%P;
}
for(int i=r-1;i>mid;--i){
Mod_add(cMx[i],cMx[i+1]),Mod_add(cMn[i],cMn[i+1]);
Mod_add(cM[i],cM[i+1]),Mod_add(sM[i],sM[i+1]);
Mod_add(sMx[i],sMx[i+1]),Mod_add(sMn[i],sMn[i+1]);
}
for(int i=mid,px=mid+1,pn=mid+1,lMx=0,rMx=0,lMn=INF,rMn=INF;i>=l;--i){
Min(lMn,A[i]),Max(lMx,A[i]);
for(;px<=r&&(rMx=max(A[px],rMx))<=lMx;++px);
for(;pn<=r&&(rMn=min(A[pn],rMn))>=lMn;++pn);
int L=min(pn,px),R=max(pn,px);
Mod_add(ans,1ll*lMx*lMn%P*dis(mid-i+2,L-i)%P);
Mod_add(ans,(cM[R]+1ll*sM[R]*(mid-i+1))%P);
if(L==pn)Mod_add(ans,((cMn[pn]-cMn[px]+1ll*(sMn[pn]-sMn[px])*(mid-i+1))%P+P)%P*lMx%P);
else Mod_add(ans,((cMx[px]-cMx[pn]+1ll*(sMx[px]-sMx[pn])*(mid-i+1))%P+P)%P*lMn%P);
}
return ans;
}
int main(){
cin>>n;for(int i=1;i<=n;++i)rd(A[i]);
cout<<solve(1,n)<<endl;
return 0;
}
orz Leefir的数形结合套线段树,orz他的神奇的三角形,orz ShinFeb的5K的分类讨论