COCI2014/2015 Contest#2 F Solution

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
#include
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(aint &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<1,n)<return 0;
}

orz Leefir的数形结合套线段树，orz他的神奇的三角形，orz ShinFeb的5K的分类讨论