【高阶数据结构】线段树加乘(维护序列)详细解释乘与加懒标记

文章目录

  • 1.题目
    • [AHOI2009] 维护序列
  • 2.懒标记处理
    • 先加后乘的形式
      • 1. 先加后乘的操作
    • 先乘后加的形式
      • 2. 先乘后加的操作
        • **乘法操作**
        • **加法操作**
    • 懒标记的下传
  • 3.代码

1.题目

题目来源:https://www.luogu.com.cn/problem/P2023

[AHOI2009] 维护序列

题目背景

老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成。

题目描述

有一个长为 n n n 的数列 { a n } \{a_n\} {an},有如下三种操作形式:

  1. 格式 1 t g c,表示把所有满足 t ≤ i ≤ g t\le i\le g tig a i a_i ai 改为 a i × c a_i\times c ai×c ;
  2. 格式 2 t g c 表示把所有满足 t ≤ i ≤ g t\le i\le g tig a i a_i ai 改为 a i + c a_i+c ai+c ;
  3. 格式 3 t g 询问所有满足 t ≤ i ≤ g t\le i\le g tig a i a_i ai 的和,由于答案可能很大,你只需输出这个数模 p p p 的值。

输入格式

第一行两个整数 n n n p p p

第二行含有 n n n 个非负整数,表示数列 { a i } \{a_i\} {ai}

第三行有一个整数 m m m,表示操作总数。

从第四行开始每行描述一个操作,同一行相邻两数之间用一个空格隔开,每行开头和末尾没有多余空格。

输出格式

对每个操作 3,按照它在输入中出现的顺序,依次输出一行一个整数表示询问结果。


2.懒标记处理

先加后乘的形式

假设我们要在一个区间上做更新操作,区间内的某个数的值用 x x x 表示,addmul 分别代表加法因子和乘法因子。

1. 先加后乘的操作

先加后乘的更新过程是:
我们想在区间上的每个元素先加一个数 a a a,再乘以一个数 m m m,这个操作可以表示为:

( x + add ) ∗ mul (x + \text{add}) * \text{mul} (x+add)mul

  • 乘法更新
    假设当前要在区间上乘以 a a a,则操作变成:
    ( x + add ) ∗ mul ∗ a (x + \text{add}) * \text{mul} * a (x+add)mula
    新的乘法标记将变为 mul ∗ a \text{mul} * a mula,这是可以接受的。

  • 加法更新
    假设现在要在区间上加上 a a a,则变成:
    ( x + add ) ∗ mul + a (x + \text{add}) * \text{mul} + a (x+add)mul+a
    这个表达式不容易简化成一种标准形式。我们可以尝试将其转换为:
    ( x + add + a mul ) ∗ mul (x + \text{add} + \frac{a}{\text{mul}}) * \text{mul} (x+add+mula)mul

    然而,这样得到的 add 标记变成了 add + a mul \text{add} + \frac{a}{\text{mul}} add+mula,这个值可能是一个小数,很难表示或处理。因此,先加后乘的形式并不理想。

先乘后加的形式

2. 先乘后加的操作

另一种常见的更新方式是先乘后加,即首先进行乘法操作,然后再进行加法操作。我们可以表示为:

x ∗ mul + add x * \text{mul} + \text{add} xmul+add

乘法操作

如果我们在这个数上乘以 m m m,则更新如下:

( x ∗ mul + add ) ∗ m = x ∗ mul ∗ m + add ∗ m (x * \text{mul} + \text{add}) * m = x * \text{mul} * m + \text{add} * m (xmul+add)m=xmulm+addm

因此:

  • 新的乘法标记变成了 mul ∗ m \text{mul} * m mulm
  • 新的加法标记变成了 add ∗ m \text{add} * m addm
加法操作

如果我们在这个数上加上 a a a,则更新如下:

x ∗ mul + add + a x * \text{mul} + \text{add} + a xmul+add+a

这里:

  • 新的加法标记变为 add + a \text{add} + a add+a
  • 乘法标记保持不变。

懒标记的下传

考虑区间树的情况,假设父节点有乘法标记 m m m 和加法标记 a a a,其更新表达式为:

( x ∗ mul + add ) ∗ m + a = x ∗ mul ∗ m + add ∗ m + a (x * \text{mul} + \text{add}) * m + a = x * \text{mul} * m + \text{add} * m + a (xmul+add)m+a=xmulm+addm+a

  • 左右孩子节点的 sum 更新为:
    root.sum ∗ m + ( root.r − root.l + 1 ) ∗ a \text{root.sum} * m + (\text{root.r} - \text{root.l} + 1) * a root.summ+(root.rroot.l+1)a
    这是一个标准的加法和乘法更新,可以继续进行懒标记下传。

  • 乘法标记(mul)下传时,更新为:
    mul ∗ m \text{mul} * m mulm

  • 加法标记(add)下传时,更新为:
    add ∗ m + a \text{add} * m + a addm+a


3.代码

//为什么先加后乘的形式不可以
//我们要变成(x+add)*mul的形式
//假设现在要在这个区间上乘 a
//那么这个数就变成了 (x+add)*mul*a
//新的mul标记就变成了 mul*a 这个是可以的
//假设现在要在这个区间上加 a
//那么这个数就变成了 (x+add)*mul + a
//化成上面的形式 (x+add + a/mul)*mul
//显然新的add标记(add+ a/mul)可能是个小数,不好表示,故而这种方式不合适

//先乘后加形式
// x*mul +add的形式
// 在这个数上乘m
// (x*mul+add)*m
// x*mul*m + add*m
// 新的mul标记就变成了 mul*m
// 新的add标记就变成了 add*m
// 在这个数上加a
// x*mul + add + a
// mul标记不变
// 新的add标记就变成了 add + a
// pushdown的时候为什么l和r的懒标记怎么改
// 显然父亲结点的mul和add就是以先乘后加的形式下传
// 假设父亲结点为m和a
// (x*mul+add)*m+ a
// x*mul*m +add*m+a
// 左右孩子的 sum = (root.sum*m+(root.r-root.l+1)*add)
// mul : mul*m
// add : add*m+a

#include <cstdio>
#include <iostream>
using namespace std;
#define int long long
typedef long long ll;
using LL =long long;
const int N = 1e5 + 10;
int n, p, m;
int w[N];
struct Node{
    int l, r, sum, add, mul; 
} tr[4 * N];

void pushup(int u)
{
    tr[u].sum = (tr[u<<1].sum+tr[u<<1|1].sum)%p;
}

void cale(Node &root, int a, int m) 
{
    root.sum = (ll)((ll)(root.sum)*m +(ll)(root.r-root.l + 1)*a)%p;
    root.add = (ll)(root.add*m+a)%p;
    root.mul = (ll)root.mul*m%p;
}

void pushdown(int u)
{
    Node & root = tr[u],& left =tr[u<<1], &right =tr[u<<1|1];
    cale(left,root.add,root.mul);
    cale(right,root.add,root.mul);
    tr[u].add=0;
    tr[u].mul=1;
}

void build(int u, int l, int r)
{
    if(l==r){
        tr[u]={l,r,w[l],0,1};
    }
    else{
        tr[u]={l,r,0,0,1};
        int mid = l+r>>1;
        build(u<<1,l,mid);
        build(u<<1|1,mid+1,r);
        pushup(u);
    }
}

void modify(int u, int l, int r, int add, int mul)
{
    if(tr[u].l>=l&&tr[u].r<=r){
        cale(tr[u],add,mul);
    }
    else{
        pushdown(u);
        int mid =tr[u].l+tr[u].r>>1;
        if(l<=mid)modify(u<<1,l,r,add,mul);
        if(r >mid)modify(u<<1|1,l,r,add,mul);
        pushup(u);
    }
}

int query(int u, int l, int r)
{
    if(tr[u].l>=l &&tr[u].r<=r)return tr[u].sum;
    else{
        pushdown(u);
        int mid =tr[u].l+tr[u].r>>1;
        ll res =0;
        if(l<=mid)res += query(u<<1,l,r)%p;
        if(r >mid)res = (res+query(u<<1|1,l,r))%p;
        return res;
    }
}

signed main()
{
    cin>>n>>p;
    for(int i=1;i<=n;i++)cin>>w[i];
    build(1,1,n);
    cin>>m;
    while ( m -- )
    {
        int t, l, r, d;
        cin>>t>>l>>r;
        if ( t == 1 ) 
        {
            cin>>d;
            modify(1, l, r, 0, d);
        }
        else if ( t == 2 )
        {
            cin>>d;
            modify(1, l, r, d, 1);
        }
        else cout<<query(1, l, r)<<'\n';
    }
    return 0;
}


本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/954172.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

一文通透OpenVLA及其源码剖析——基于Prismatic VLM(SigLIP、DinoV2、Llama 2)及离散化动作预测

前言 当对机器人动作策略的预测越来越成熟稳定之后(比如ACT、比如扩散策略diffusion policy)&#xff0c;为了让机器人可以拥有更好的泛化能力&#xff0c;比较典型的途径之一便是基于预训练过的大语言模型中的广泛知识&#xff0c;然后加一个policy head(当然&#xff0c;一开…

《操作系统真象还原》第十三章——磁盘驱动程序

文件系统磁盘创建 创建磁盘 进入bochs安装目录&#xff0c;输入以下命令 ./bin/bximage 然后按照以下步骤创建硬盘 修改硬盘配置 vim boot.disk 添加以下代码行 ata0-slave: typedisk, path"hd80M.img", modeflat,cylinders162,heads16,spt63 完整配置如下 …

快速、可靠且高性价比的定制IP模式提升芯片设计公司竞争力

作者&#xff1a;Karthik Gopal&#xff0c;SmartDV Technologies亚洲区总经理 智权半导体科技&#xff08;厦门&#xff09;有限公司总经理 无论是在出货量巨大的消费电子市场&#xff0c;还是针对特定应用的细分芯片市场&#xff0c;差异化芯片设计带来的定制化需求也在芯片…

v-bind操作class

v-bind操作class 参考文献&#xff1a; Vue的快速上手 Vue指令上 Vue指令下 Vue指令的综合案例 指令的修饰符 文章目录 v-bind操作classv-bind对于样式控制的增强操作class案例(tab导航高亮)操作style操作style案例 结语 博客主页: He guolin-CSDN博客 关注我一起学习&#…

Kubernetes1.28 编译 kubeadm修改证书有效期到 100年.并更新k8s集群证书

文章目录 前言一、资源准备1. 下载对应源码2.安装编译工具3.安装并设置golang 二、修改证书有效期1.修改证书有效期2.修改 CA 证书有效期 三、编译kubeadm四、使用新kubeadm方式1.当部署新集群时,使用该kubeadm进行初始化2.替换现有集群kubeadm操作 前言 kubeadm 默认证书为一…

HarmonyOS NEXT应用开发边学边玩系列:从零实现一影视APP (三、影视搜索页功能实现)

在HarmonyOS NEXT开发环境中&#xff0c;我们可以使用nutpi/axios库来简化网络请求的操作。本文将展示如何使用HarmonyOS NEXT框架和nutpi/axios库&#xff0c;从零开始实现一个简单的影视APP&#xff0c;主要关注影视搜索页的功能实现。 为什么选择nutpi/axios&#xff1f; n…

高级运维:shell练习2

1、需求&#xff1a;判断192.168.1.0/24网络中&#xff0c;当前在线的ip有哪些&#xff0c;并编写脚本打印出来。 vim check.sh #!/bin/bash# 定义网络前缀 network_prefix"192.168.1"# 循环遍历1-254的IP for i in {1..254}; do# 构造完整的IP地址ip"$network_…

好用的php商城源码有哪些?

选择一个优秀的商城工具&#xff0c;能更好地帮助大家建立一个好用的商城系统。目前比较流行的都是开源PHP商城系统&#xff0c;那么现实中都有哪些好用的PHP商城源码值得推荐呢&#xff1f;下面就带大家一起来了解一下。 1.TigShop 【推荐指数】&#xff1a;★★★★★☆ 【推…

Docker Desktop 构建java8基础镜像jdk安装配置失效解决

Docker Desktop 构建java8基础镜像jdk安装配置失效解决 文章目录 1.问题2.解决方法3.总结 1.问题 之前的好几篇文章中分享了在Linux(centOs上)和windows10上使用docker和docker Desktop环境构建java8的最小jre基础镜像&#xff0c;前几天我使用Docker Desktop环境重新构建了一个…

Open FPV VTX开源之嵌入式OSD配置

Open FPV VTX开源之嵌入式OSD配置 1. 源由2. 安装3. 配置步骤一&#xff1a;备份/etc/telemetry.conf步骤二&#xff1a;修改/etc/telemetry.conf步骤三&#xff1a;配置时区步骤四&#xff1a;重启摄像头 4. 实测5. 参考资料 1. 源由 穿越机模拟图传延迟通常在10ms左右。 最…

数据平台浅理解

定义 数据平台架构是指用于收集、存储、处理和分析数据的一系列组件、技术和流程的整体架构设计。它就像是一个复杂的数据生态系统的蓝图&#xff0c;旨在高效地管理数据从产生源头到产生价值的整个生命周期。 主要层次 数据源层 这是数据的起点&#xff0c;包含各种类型的数据…

CSS3的aria-hidden学习

前言 aria-hidden 属性可用于隐藏非交互内容&#xff0c;使其在无障碍 API 中不可见。即当aria-hidden"true" 添加到一个元素会将该元素及其所有子元素从无障碍树中移除&#xff0c;这可以通过隐藏来改善辅助技术用户的体验&#xff1a; 纯装饰性内容&#xff0c;如…

【ArcGIS初学】产生随机点计算混淆矩阵

混淆矩阵&#xff1a;用于比较分类结果和地表真实信息 总体精度(overall accuracy) :指对角线上所有样本的像元数(正确分类的像元数)除以所有像元数。 生产者精度(producers accuracy) &#xff1a;某类中正确分类的像元数除以参考数据中该类的像元数(列方向)&#xff0c;又称…

认识机器学习中的结构风险最小化准则

上一篇文章我们学习了关于经验风险最小化准则&#xff0c;其核心思想是通过最小化训练数据上的损失函数来优化模型参数&#xff0c;从而提高模型在训练集上的表现。但是这也会导致一个问题&#xff0c;经验风险最小化原则很容易导致模型在训练集上错误率很低&#xff0c;但在未…

设计模式-工厂模式/抽象工厂模式

工厂模式 定义 定义一个创建对象的接口&#xff0c;让子类决定实列化哪一个类&#xff0c;工厂模式使一个类的实例化延迟到其子类&#xff1b; 工厂方法模式是简单工厂模式的延伸。在工厂方法模式中&#xff0c;核心工厂类不在负责产品的创建&#xff0c;而是将具体的创建工作…

Chatper 4: Implementing a GPT model from Scratch To Generate Text

文章目录 4 Implementing a GPT model from Scratch To Generate Text4.1 Coding an LLM architecture4.2 Normalizing activations with layer normalization4.3 Implementing a feed forward network with GELU activations4.4 Adding shortcut connections4.5 Connecting at…

Unity ShaderGraph中Lit转换成URP的LitShader

ShaderGraph中的LitShader如下&#xff1a; 在顶点和片元着色器暴露出了上图中的几个参数&#xff0c;要转换成URPLitShaderLab,首先要找到这几个参数&#xff0c;打开LitShader&#xff0c;找到第一个Pass&#xff0c;可以看到下图中的顶点和片元的定义函数&#xff0c;还有引…

uni-app的学习

uni-app 有着跨平台支持、丰富的插件和生态系统、高性能、集成开发工具HBuilderX的配合使用。允许使用者仅通过一套代码发布到多平台使用。 uni-app官网 uni-app 是一个适合开发跨平台移动应用和小程序的框架&#xff0c;能够大幅提高开发效率。 一、了解 1.1 工具准备 从Git…

USRP X310 Windows 烧录镜像

说明 USRP-X 系列设备包含两个用于两个以太网通道的 SFP 端口。由于 SFP 端口支持 1 千兆 (SFP) 和 10 千兆 (SFP) 收发器&#xff0c;因此 UHD 附带了多个 FPGA 图像&#xff0c;以确定上述接口的行为。 注意&#xff1a;Aurora 图像需要从 FPGA 源代码手动构建。 FPGA 图像…

Sprint Boot教程之五十八:动态启动/停止 Kafka 监听器

Spring Boot – 动态启动/停止 Kafka 监听器 当 Spring Boot 应用程序启动时&#xff0c;Kafka Listener 的默认行为是开始监听某个主题。但是&#xff0c;有些情况下我们不想在应用程序启动后立即启动它。 要动态启动或停止 Kafka Listener&#xff0c;我们需要三种主要方法…