2021HNCPC-E-差分,思维

2023大厂真题提交网址(含题解):

www.CodeFun2000.com(http://101.43.147.120/)

最近我们一直在将收集到的机试真题制作数据并搬运到自己的OJ上,供大家免费练习,体会真题难度。现在OJ已录入50+道2023年最新大厂真题,同时在不断的更新。同时,可以关注"塔子哥学算法"公众号获得每道题的题解。
在这里插入图片描述

题目大意:

求解每一个 f ( i ) f(i) f(i)代表 m a x ( a j , . . . , a i ) − ( j − i + 1 ) ≥ m max(a_j,...,a_i)-(j-i+1) \geq m max(aj,...,ai)(ji+1)m j j j的个数. j ≤ i j \leq i ji
n ≤ 1 e 6 n \leq 1e6 n1e6

思路:

关键:看到 m a x ( a j , . . . , a i ) max(a_j,...,a_i) max(aj,...,ai),想到枚举 a i a_i ai,考虑它所管辖的范围(单调栈求解).

PS:当有相同的值时,可以规定将区间贡献记在最右边的值上。

那么就转化成 ( j − i + 1 ) ≤ m a x ( a j , . . . , a i ) − m (j-i+1) \leq max(a_j,...,a_i)-m (ji+1)max(aj,...,ai)m.

然后枚举右半边,分情况讨论来区间加等差序列即可.

代码仓库:
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
#define pb push_back
#define mp make_pair
#define vi vector<int>
#define vll vector<ll>
#define fi first
#define se second
const int maxn = 1e6 + 5;
const int mod = 1e9 + 7;
int a[maxn] , l[maxn] , r[maxn] , b[maxn];
int n , m;
template <typename T>
void read(T & x){ x = 0;T f = 1;char ch = getchar();while(!isdigit(ch)){if(ch == '-') f = -1;
ch = getchar();}while (isdigit(ch)){x = x * 10 + (ch ^ 48);ch = getchar();}x *= f;}
template <typename T>
void write(T x){if(x < 0){putchar('-');x = -x;}if (x > 9)write(x / 10);putchar(x % 10 + '0');}
void add (int l , int r , int s , int d){
    if (l > r) return ;
    int t = s + (r - l) * d;
    b[l] += s;
    b[l+1] += d - s;
    b[r+1] -= d + t;
    b[r+2] += t;
}
void solve (){
    for (int i = 1 ; i <= n ; i++) b[i] = 0;
    for (int i = 1 ; i <= n ; i++){
        int g = a[i] - m;
        if (g <= 0) continue;
        int v = i - l[i] + 1;
        if (v >= g){
            add(i , min(i + g - 1 , r[i]) , g , -1);
        }else{
            int gap = g - v;
            int e1 = min(i + gap , r[i]);
            add(i , e1 , v , 0);
            add(e1 + 1 , min(r[i] , e1 + v - 1) , v - 1 , -1);
        }
    }
    for (int i = 1 ; i <= n ; i++) b[i] += b[i - 1];
    for (int i = 1 ; i <= n ; i++) b[i] += b[i - 1];
}
int main()
{
    while (~scanf("%d%d" , &n, &m)){
        for (int i = 1 ; i <= n ; i++) read(a[i]);
        stack<int> s;
        for (int i = 1 ; i <= n ; i++){
            while (s.size() && a[s.top()] <= a[i]) s.pop();
            if (s.size()) l[i] = s.top() + 1;
            else l[i] = 1;
            s.push(i);
        }
        while (s.size()) s.pop();
        for (int i = n ; i >= 1 ; i--){
            while (s.size() && a[i] > a[s.top()]) s.pop();
            if (s.size()) r[i] = s.top() - 1;
            else r[i] = n;
            s.push(i);
        }
        solve();
        for (int i = 1 ; i <= n ; i++) printf("%d " , b[i]);
        printf("\n");
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值