BZOJ3745: [Coci2015]Norma【分治】

题目描述：

Mod 10⁹

题目分析：

分治真是强大。。

计算$[l,r]$的答案，只需要考虑怎么算跨过$mid$的答案，其余的分治即可。枚举左端点$i:mid\to l$，记$[i,mid]$的最小值为$mn$，最大值为$mx$，找到右边最后一个$\ge mn$的位置$j$，以及最后一个$\le mx$的位置$k$。

不妨设$j\le k$，这时答案分成了三个部分：

• 右端点$\le j$，$ans=mn\ast mx{\sum }_{p=mid+1}^{j}p-i+1$，可直接计算。
• $j<$右端点$\le k$，$ans=mn\ast {\sum }_{p=j+1}^{k}max[mid+1,p]\ast (p-i+1)=mn\ast {\sum }_{p=j+1}^{k}max[mid+1,p]\ast p-mn\ast {\sum }_{p=j+1}^k\ast max[mid+1,p]\ast (i-1)$
  预处理$\sum max[mid+1,p]\ast p$以及$\sum max[mid+1,p]$即可。
• $k<$右端点$\le r$，$ans={\sum }_{p=k+1}^{r}min[mid+1,p]\ast max[mid+1,p]\ast (p-i+1)$，同上面一样地拆分。
  预处理$\sum min[mid+1,p]\ast max[mid+1,p]\ast p$以及$min[mid+1,p]\ast max[mid+1,p]$即可。

随着$i:mid\to l$，$j$和$k$单调不减，所以复杂度是$O(r-l)$的，总复杂度就是$O(nlogn)$。

如果可以用$O(r-l)$或$O((r-l)\ast log)$的时间，来解决跨越$mid$的答案，这时就可以考虑分治。

Code：

#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
#define maxn 500005 
using namespace std;
const int mod = 1e9;
int n,a[maxn],ans;
int smn[maxn],SMN[maxn],smx[maxn],SMX[maxn],mnx[maxn],MNX[maxn];
void solve(int l,int r){
    if(l==r) {ans=(ans+1ll*a[l]*a[r])%mod;return;}
    int mid=(l+r)>>1;
    solve(l,mid),solve(mid+1,r);
    smn[mid]=SMN[mid]=smx[mid]=SMX[mid]=mnx[mid]=MNX[mid]=0;
    for(int i=mid+1,mx=0,mn=1e9;i<=r;i++){
        mx=max(mx,a[i]),mn=min(mn,a[i]);
        smn[i]=(smn[i-1]+mn)%mod,SMN[i]=(SMN[i-1]+1ll*mn*i)%mod;
        smx[i]=(smx[i-1]+mx)%mod,SMX[i]=(SMX[i-1]+1ll*mx*i)%mod;
        mnx[i]=(mnx[i-1]+1ll*mn*mx)%mod,MNX[i]=(MNX[i-1]+1ll*mn*mx%mod*i)%mod;
    }
    for(int i=mid,j=mid,k=mid,mx=0,mn=1e9;i>=l;i--){
        mx=max(mx,a[i]),mn=min(mn,a[i]);
        while(j<r&&a[j+1]>=mn) j++;
        while(k<r&&a[k+1]<=mx) k++;
        int x=min(j,k),y=max(j,k);
        ans=(ans+1ll*(mid+1-i+1+x-i+1)*(x-mid)/2%mod*mx%mod*mn)%mod;
        ans=(ans+MNX[r]-MNX[y]-1ll*(mnx[r]-mnx[y])*(i-1))%mod;
        if(j<k) ans=(ans+1ll*mx*(SMN[k]-SMN[j]-1ll*(smn[k]-smn[j])*(i-1)%mod))%mod;
        else ans=(ans+1ll*mn*(SMX[j]-SMX[k]-1ll*(smx[j]-smx[k])*(i-1)%mod))%mod;
    }
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    solve(1,n);
    printf("%d\n",(ans+mod)%mod);
}

还有大佬用线段树的做法。大致思路就是枚举右端点$r:1\to n$，计算左端点的答案，每次移动右端点都会对最小值或最大值进行区间修改，同样要维护$min[l,r]\ast max[l,r]\ast l$和$min[l,r]\ast max[l,r]$，然后线段树修改，大致思路就是这样。这里有一篇用线段树写的blog，不确定是不是一样的，姑且看一看。