差分算法篇

引入 



AcWing借教室

        这道题题目较长,简单概括就是为了解决租借教室问题,我们需要处理接下来 n天的借教室信息,其中第 i天学校有 ri个教室可供租借。共有 m 份订单,每份订单用三个正整数描述,分别为 dj,sj,tj,表示某租借者需要从第 sj天到第 tj天租借教室(包括第 sj天和第 tj天),每天需要租借 dj个教室。我们按照订单顺序分配顺序,需要判断是否有订单无法满足,如果所有订单均可满足,则输出只有一行,包含一个整数 0。否则(订单无法完全满足)输出两行,第一行输出一个负整数 −1−1,第二行输出需要修改订单的申请人编号。

        现在我们来分析一下这道题,由于数据量庞大,并且要确定第一个不满足订单的位置,我们可以想到二分的方法(不了解的朋友可见http://t.csdnimg.cn/T6GS7),又因为每天的教室数目越到后面一定是递减的,所以,这里我们还需要用到另一个算法知识——差分

        下面我们来了解一下差分相关知识。

差分

        差分是一种常用的技巧,用于处理序列中元素之间的差值关系。差分的基本思想是在序列中相邻元素之间存储差值,通过这种方式可以在更快的时间复杂度下实现一些操作。

        在数学上,差分可以理解为序列中相邻元素的差值序列。差分序列的定义如下: 设原始序列为 a[1], a[2], ..., a[n],差分序列为 b[1], b[2], ..., b[n],其中 b[i] = a[i] - a[i-1]。通过差分操作,可以将原始序列中相邻元素之间的关系转化为差分序列中的元素。 差分是一种重要的处理序列的技巧,能够简化操作、提高效率,是算法设计中常用的方法之一。

举个栗子:

        假设有一个原始序列 a = {0, 0, 0, 0, 0, 0},初始化的差分数组 b = {0, 0, 0, 0, 0, 0}。 现在我们有一个操作序列如下: d = {1, -1, 2},表示依次对应的增量 l = {1, 3, 2},表示对应的左边界  r = {3, 5, 4},表示对应的右边界

        根据代码中的逻辑,对差分数组 b 进行操作:

差分数组中保存的值表示了原始序列中相邻元素之间的差值关系,用于记录对原始序列中区间的增量或减量。

具体来说,对于差分数组 b,如果 b[i] 表示原始序列中 a[i] 和 a[i-1] 之间的差值(即 a[i] - a[i-1]),那么通过差分数组可以方便地进行区间操作。差分数组中的值可以看作是对应位置处的增减量,通过对差分数组的累加可以快速得到对应位置处的值。

题目解析

那么通过二分与差分的两个概念,我们就可以写出下列代码:

1.若所有的订单都满足,则输出“0”,直接返回

 if(!check(m)) {
        puts("0");
        return 0;
    }

2. 通过二分算法,如果中间值订单不满足,则在左区间继续查找,right=mid; 反之在右区间, left=mid+1; 

 int l=1,r=m;
    while(l<r){
        int mid=(l+r) >>1;
        if(check(mid)) r=mid;
        else l=mid+1;
    }

3.之后,我们定义一个check函数,用于差分

bool check(int mid) // 定义布尔函数 check,参数为整型 mid
{
    memset(b, 0, sizeof b); // 将数组 b 的所有元素初始化为 0
    for(int i=1;i<=n;i++) b[i]=r[i]-r[i-1];
    for (int i = 1; i <= mid; i ++ )
    {
        b[s[i]] -= d[i];
        b[t[i] + 1] += d[i];
    }//差分

    for (int i = 1; i <= n; i ++ )
    {
        b[i] += b[i - 1];
        if (b[i] <0) return true;//我们定义的是不满足为真,所有这里为true
    }

    return false;
}

 下列是完整代码L;

#include<iostream>
#include <cstring>
#include <algorithm>
using namespace std;

typedef long long LL;
const int N=1000010;

int n,m;
int r[N];
int d[N],s[N],t[N];
LL b[N];//差分数组


bool check(int mid)
{
    memset(b, 0, sizeof b);
    for(int i=1;i<=n;i++) b[i]=r[i]-r[i-1];
    for (int i = 1; i <= mid; i ++ )
    {
        b[s[i]] -= d[i];
        b[t[i] + 1] += d[i];
    }

    for (int i = 1; i <= n; i ++ )
    {
        b[i] += b[i - 1];
        if (b[i] <0) return true;
    }

    return false;
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&r[i]);
    }
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&d[i],&s[i],&t[i]);
    }
    
    if(!check(m)) {
        puts("0");
        return 0;
    }
    
    int l=1,r=m;
    while(l<r){
        int mid=(l+r) >>1;
        if(check(mid)) r=mid;
        else l=mid+1;
    }
   
   
    printf("-1\n%d\n",r);
    return 0;
}

  • 25
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 2
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值