COCI2014/2015 Contest#2 F Solution

来源：互联网发布：java 接口变量编辑：程序博客网时间：2024/04/30 09:03

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,r]，我们可以把它分段：
区间[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的分类讨论

0 0