2281. 巫师的总力量和(单调栈+前缀和的前缀和)

2281. 巫师的总力量和

题目描述

在这里插入图片描述

示例

在这里插入图片描述

在这里插入图片描述

数据范围

  • 1<= strength.length <= 105
  • 1 <= strength[i] <= 109

思路

参考灵茶山艾府

  • 以每个点作为最小值,计算包含该最小值的所有区间所造成的贡献。
  • 设左右边界 [ L , R ] {[L,R]} [L,R](即为包含该元素为最小值的左右区间端点),为了避免重复计算,考虑左侧严格小于当前元素的最近位置 L − 1 {L-1} L1,以及右侧 小于等于当前元素的最近元素位置 R + 1 {R + 1} R+1
  • 设当前枚举的巫师的能力值为 v {v} v,那么他对答案产生的贡献是 v {v} v 乘上在左右边界 [ L , R ] {[L,R]} [L,R]内的所有包含 v {v} v 的子数组的元素和的和。

如何计算子数组的元素和的和

  • 不妨固定右端点,移动左端点来考虑。

设子数组右端点为 r {r} r,左端点为 l {l} l,当前枚举的元素下标为 i {i} i,那么有 l < = i < = r {l<=i<=r} l<=i<=r

s t r e n g t h {strength} strength 数组的前缀和为 s {s} s,其中 s [ i ] = ∑ j = 0 i − 1 s t r e n g t h [ j ] {s[i] = \sum\limits_{j=0}^{i-1}strength[j]} s[i]=j=0i1strength[j],因此子数组 [ l , r ] {[l,r]} [l,r]的元素和可以表示为

s [ r + 1 ] − s [ l ] {s[r+1]-s[l]} s[r+1]s[l]

在范围 [ L , R ] {[L,R]} [L,R] 内的所有子数组的元素和的和可以表示为

∑ r = i + 1 R + 1 ∑ l = L i s [ r ] − s [ l ] = ∑ r = i + 1 R + 1 ( ( i − L + 1 ) ⋅ s [ r ] − ∑ l = L i s [ l ] ) = ( i − L + 1 ) ⋅ ∑ r = i + 1 R + 1 s [ r ] − ( R − i + 1 ) ⋅ ∑ l = L i s [ l ] \begin{aligned} & \sum\limits_{r=i+1}^{R+1}\sum\limits_{l=L}^{i}s[r]-s[l] \\ &= \sum\limits_{r=i+1}^{R+1}\big((i-L+1) · s[r] - \sum\limits_{l=L}^{i}s[l]\big) \\ &=(i-L+1) · \sum\limits_{r=i+1}^{R+1}s[r] - (R-i+1) ·\sum\limits_{l=L}^{i}s[l] \end{aligned} r=i+1R+1l=Lis[r]s[l]=r=i+1R+1((iL+1)s[r]l=Lis[l])=(iL+1)r=i+1R+1s[r](Ri+1)l=Lis[l]

所以还需计算出前缀和 s {s} s 的前缀和 s s {ss} ss, 其中 s s [ i ] = ∑ j = 0 i − 1 s [ j ] {ss[i] = \sum\limits_{j=0}^{i-1}s[j]} ss[i]=j=0i1s[j],上式即为

( i − L + 1 ) ⋅ ( s s [ R + 2 ] − s s [ i + 1 ] ) − ( R − i + 1 ) ⋅ ( s s [ i + 1 ] − s s [ L ] ) (i-L+1) ·(ss[R+2]-ss[i+1])-(R-i+1)·(ss[i+1]-ss[L]) (iL+1)(ss[R+2]ss[i+1])(Ri+1)(ss[i+1]ss[L])

在乘上 v {v} v 即为当前巫师贡献,累加所有贡献即为答案。


注意的点

在这里插入图片描述


代码

class Solution {
public:
    const int mod = 1e9 + 7;
    #define ll long long

    int totalStrength(vector<int>& a) {
        
        int n = a.size();
        vector<int> L(n, -1), R(n, n);
        stack<int> s; 
        // 直接求出左右
        for (int i = 0; i < n; i++) {
            while(!s.empty() && a[s.top()] >= a[i]) R[s.top()] = i, s.pop();
            if(!s.empty()) L[i] = s.top();
            s.push(i); 
        }

        long long pre = 0;
        vector<int> ppre(n + 2);
        for (int i = 1; i <= n; i++) {
            pre += a[i - 1];
            ppre[i + 1] = (ppre[i] + pre) % mod;
        }

        int res = 0;
        for (int i = 0; i < n; i++) {
            ll l = L[i] + 1, r = R[i] - 1;
            ll ans = ((i - l + 1) * (ppre[r + 2] - ppre[i + 1]) - \
                      (r - i + 1) * (ppre[i + 1] - ppre[l])) % mod;
            res = (res + ans * a[i]) % mod;
        }
        return (res + mod) % mod; // 防止负数
    }
};
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

ღCauchyོꦿ࿐

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值