C++ 算法学习——1.6 差分算法与二维差分算法

news2024/11/30 15:24:40
  1. 一维差分算法概述

    • 差分算法是一种用于计算序列中相邻元素之间差值的技术。
    • 在C++中,STL(标准模板库)提供了std::adjacent_difference函数来实现差分算法。
  2. std::adjacent_difference函数

    • std::adjacent_difference函数位于头文件<numeric>中。
    • 语法:std::adjacent_difference(first, last, result)
      • firstlast是定义了输入序列范围的迭代器。
      • result是存储结果的目标序列的起始位置

一维差分数组在处理连续序列中的区间更新和查询操作时非常有用。它可以帮助在常数时间复杂度内完成区间元素的增减操作,并支持常数时间复杂度的区间查询操作。以下是一维差分数组的一些主要用途和优点:

  1. 区间更新:对于原始数组中 [l, r] 区间内的元素增加 val,只需更新差分数组 diff 中 diff[l] += val 和 diff[r+1] -= val,然后再根据差分数组计算更新后的原始数组。

  2. 查询:原始数组中位置为L的元素等于差分数组中 [L]前所有元素的和 

  3. 应用领域:差分数组在处理数组序列的区间操作中广泛应用,比如处理频繁的区间修改、查询问题。

代码示例(构建): 

std::vector<int> differences(nums.size());

    std::adjacent_difference(nums.begin(), nums.end(), differences.begin());

  1. 二维差分算法概述

    • 二维差分算法是一种用于计算二维数组中差分的技术。
    • 通过预处理数组,可以在常数时间内计算子矩阵的增减值,而不必每次都一一计算。
  2. 实现方式

    • 我们可以使用差分数组diff表示相邻元素之间的差值关系。二维差分数组diff的构建规则如下:
  • diff[i][j] = matrix[i][j] - matrix[i-1][j] - matrix[i][j-1] + matrix[i-1][j-1]

代码示例(构建):

int **computeDifferenceArray(int **matrix, int rows, int cols) {
    int **diff = new int*[rows];
    for (int i = 0; i < rows; i++) {
        diff[i] = new int[cols]();
    }

    for (int i = 0; i < rows; i++) {
        for (int j = 0; j < cols; j++) {
            diff[i][j] = matrix[i][j]; // 复制原始数组的值

            // 处理边界情况
            if (i > 0) diff[i][j] -= matrix[i - 1][j]; // 上方元素
            if (j > 0) diff[i][j] -= matrix[i][j - 1]; // 左方元素
            if (i > 0 && j > 0) diff[i][j] += matrix[i - 1][j - 1]; // 左上方元素
        }
    }

    return diff;
}

 代码示例(计算)://如果对原数组(x1,y1)到(x2,y2)的矩形区域整体加上a.

void updateDiffArray(int** &diff, int x1, int y1, int x2, int y2, int a, int rows, int cols) {
    // Check if the indices are within bounds
    if (x1 >= 0 && x1 < rows && y1 >= 0 && y1 < cols &&
        x2 >= 0 && x2 < rows && y2 >= 0 && y2 < cols) {
        diff[x1][y1] += a;
        if (x2 + 1 < rows && y1 < cols) diff[x2 + 1][y1] -= a;
        if (x1 < rows && y2 + 1 < cols) diff[x1][y2 + 1] -= a;
        if (x2 + 1 < rows && y2 + 1 < cols) diff[x2 + 1][y2 + 1] += a;
    } 
}

最后通过二维前缀和还原到原数组即可。


 P1. 洛谷p2367语文成绩

#include<iostream>
#include<numeric>
#include<vector>
#include<cmath>
using namespace std;
long long ans=0;
long long result=10000;

int main()
{
    int n,q;cin>>n>>q;
    vector<int> nums,differs;
    for(int i=1;i<=n;i++)
    {
        int a;cin>>a;
        nums.push_back(a);
    }
    differs.resize(n);//很重要,不然调用函数会访问到未分配的内存,发生越界错误
    adjacent_difference(nums.begin(),nums.end(),differs.begin());
    //for(int i=0;i<n;i++) cout<<differs[i]<<" ";
    for(int i=1;i<=q;i++)
    {
        int left,right,value;
        cin>>left>>right>>value;
        differs[left-1]+=value;
        differs[right]-=value;
    }
    for(int i=0;i<n;i++)
    {
        ans+=differs[i];
        result=fmin(ans,result);
    }
    cout<<result;
    return 0;
}

P2. 洛谷p3397地毯

#include<iostream>
using namespace std;

int **computeDifferenceArray(int **matrix, int rows, int cols) {
    int **diff = new int*[rows];
    for (int i = 0; i < rows; i++) {
        diff[i] = new int[cols]();
    }

    for (int i = 0; i < rows; i++) {
        for (int j = 0; j < cols; j++) {
            diff[i][j] = matrix[i][j]; // 复制原始数组的值

            // 处理边界情况
            if (i > 0) diff[i][j] -= matrix[i - 1][j]; // 上方元素
            if (j > 0) diff[i][j] -= matrix[i][j - 1]; // 左方元素
            if (i > 0 && j > 0) diff[i][j] += matrix[i - 1][j - 1]; // 左上方元素
        }
    }

    return diff;
}

void updateDiffArray(int** &diff, int x1, int y1, int x2, int y2, int a, int rows, int cols) {
    // Check if the indices are within bounds
    if (x1 >= 0 && x1 < rows && y1 >= 0 && y1 < cols &&
        x2 >= 0 && x2 < rows && y2 >= 0 && y2 < cols) {
        diff[x1][y1] += a;
        if (x2 + 1 < rows && y1 < cols) diff[x2 + 1][y1] -= a;
        if (x1 < rows && y2 + 1 < cols) diff[x1][y2 + 1] -= a;
        if (x2 + 1 < rows && y2 + 1 < cols) diff[x2 + 1][y2 + 1] += a;
    } 
}

int** build2DPrefixSumArray(int** matrix, int rows, int cols) {
    // 分配内存来存储前缀和数组
    int** prefixSum = new int*[rows];
    for (int i = 0; i < rows; i++) {
        prefixSum[i] = new int[cols];
    }
 
    // 初始化前缀和数组
    for (int i = 0; i < rows; i++) {
        for (int j = 0; j < cols; j++) {
            prefixSum[i][j] = 0;
        }
    }
 
    // 计算第一行的前缀和
    for (int j = 1; j < cols; j++) {
        prefixSum[1][j] = matrix[1][j]+prefixSum[1][j-1];
    }
 
    // 计算第一列的前缀和
    for (int i = 2; i < rows; i++) {
        prefixSum[i][1] = prefixSum[i - 1][1] + matrix[i][1];
    }
 
    // 计算其他位置的前缀和
    for (int i = 2; i < rows; i++) {
        for (int j = 2; j < cols; j++) {
            prefixSum[i][j] = matrix[i][j] + prefixSum[i - 1][j] + prefixSum[i][j - 1] - prefixSum[i - 1][j - 1];
        }
    }
 
    return prefixSum;
}

int main()
{
    int n,q;cin>>n>>q;
    int** origin=new int*[n+1];
    for(int i=0;i<n+1;i++)
    {
        origin[i]=new int[n+1];
        fill(origin[i],origin[i]+n+1,0);
    }
    int** diff=computeDifferenceArray(origin,n+1,n+1);
    for(int i=1;i<=q;i++)
    {
        int x1,y1,x2,y2;
        cin>>x1>>y1>>x2>>y2;
        updateDiffArray(diff,x1,y1,x2,y2,1,n+1,n+1);
    }
    int** ans=build2DPrefixSumArray(diff,n+1,n+1);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        cout<<ans[i][j]<<" ";
        cout<<endl;
    }
    return 0;
}

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

基于MindSpore实现CycleGAN壁画修复

基于MindSpore实现CycleGAN壁画修复_哔哩哔哩_bilibili 本案例运行需要较大内存&#xff0c;建议在Ascend/GPU上运行。 模型介绍 模型简介 CycleGAN(Cycle Generative Adversarial Network) 即循环对抗生成网络&#xff0c;来自论文 Unpaired Image-to-Image Translation us…

指针赋值or常数赋值

int main (){int a 10;int b ;b a;int *c &a;int *d c; } 常数 a,b赋值&#xff1a; 都是将存储的值&#xff08;10&#xff09;赋值给别人。 指针赋值也是类似的&#xff1a; 指针存储的值&#xff08;&a&#xff09;为地址&#xff0c;就是把c指向的地址赋值给…

C语言 | Leetcode C语言题解之第458题可怜的小猪

题目&#xff1a; 题解&#xff1a; int poorPigs(int buckets, int minutesToDie, int minutesToTest){int base minutesToTest / minutesToDie 1;int pigs ceil(log(buckets)/log(base));return pigs; }

24-10-2-读书笔记(二十二)-《契诃夫文集》(一)上([俄] 契诃夫 [译] 汝龙)啊!真想生活。

文章目录 《契诃夫文集》&#xff08;一&#xff09;上&#xff08;[俄] 契诃夫 [译] 汝龙 &#xff09;早期生活——塔甘罗格&#xff08;人物家庭简介&#xff09;学生时期——莫斯科&#xff08;写作与学习&#xff09;流浪时期——哈萨林&#xff08;游历与流浪&#xff09…

VCSEL驱动电路

1.1 驱动电路 发射端可用MOS管控制VCSEL二极管负极方式发出脉冲光(正极对地)&#xff0c;具体作用过程如下&#xff1a; Step 1: MOS管断开, C2 电容充电(左侧HV)&#xff1b; Step 2: 信号控制MOS管打开&#xff1b; Step 3: MOS管打开后, C2电容左侧电压降为0V, 右侧变为…

当你系统有多个模块的时候,要设计统一入口页了。

本次给大家带来一批系统统一入口页的界面&#xff0c;这次都是科技感十足的界面。入口页将所有系统集中展示出来&#xff0c;并辅以其他设计元素进行氛围烘托&#xff0c;别看简单&#xff0c;但真的不好设计。

免费送源码:Java+B/S+MySQL 基于springboot网上书店管理系统 计算机毕业设计原创定制

基于springboot网上书店管理系统 摘 要 网上书店管理系统采用B/S结构、java开发语言、以及Mysql数据库等技术。系统主要分为管理员和用户两部分&#xff0c;管理员管理主要功能包括&#xff1a;首页、网站管理&#xff08;轮播图、网站公告&#xff09;人员管理&#xff08;管…

MATLAB - 机械臂手眼标定(眼在手内) - 估计安装在机器人上的移动相机的姿态

系列文章目录 前言 本示例展示了如何为装有手眼构型摄像头的机械臂或机械手执行和验证手眼校准。 一、概述 执行手眼校准有助于操作配备末端执行器&#xff08;简称 “手”&#xff09;的机械臂&#xff0c;该末端执行器依赖于摄像头提供的视觉数据。一旦完成了眼在手外的校准&…

数据结构阶段测试2的一点小补充

数据结构阶段测试2的一点小补充 1.已知⼩根堆为8,15,10,21,34,16,12&#xff0c;删除关键字8之后需重建堆&#xff0c;最后的叶⼦ 节点为() A. 34 B. 21 C. 16 D. 12 解题思路 向下调整算法删除堆顶元素 &#x1f4a1; 答案&#xff1a;C 删除堆顶元素的思路&#xff1a; …

环境对象刺激,recordnunber,记忆柱,记忆柱群,主注意对象,目的对象,状态中枢,奖惩预期,思维等等之间的联系

我们清醒状态下&#xff0c;随时都有目的&#xff0c;目的控制影响着我们思想行为的方向。目的是用对象来表征的&#xff08;目的对象&#xff09;&#xff0c;对象之所以能够表征目的&#xff0c;是因为对象能够被&#xff08;状态性&#xff09;赋值&#xff08;任何赋值都是…

UE4 材质学习笔记03(翻书(Flipbook)动画/环境混合)

一.FlipBook Animation 如果你想让游戏以每秒30帧的速度运行&#xff0c;所有内容都必须在33毫秒内渲染出来&#xff0c; 如果你想让游戏以每秒60帧的速度运行的话&#xff0c;必须在16毫秒内。 所以当一个效果需要很多细节的时候&#xff0c;往往会离线创建它&#xff0c;然…

【Qt】控件概述(4)—— 输出类控件

输出类控件 1. QLineEdit——单行输入框2. QTextEdit——多行输入框3. QComboBox——下拉框4. QSpinBox——微调框5. QDateEdit && QTimeEdit && QDateTimeEdit6 QDial——旋钮7. QSlider——滑动条 1. QLineEdit——单行输入框 QLineEdit是一个单行的输入框&…

BUU刷题-Pwn-shanghai2018_baby_arm(ARM_ROP_csu_init,ARM架构入门)

解题思路&#xff1a; 泄露或修改内存数据&#xff1a; 堆地址&#xff1a;无需栈地址&#xff1a;无需libc地址&#xff1a;无需BSS段地址&#xff1a;无需 劫持程序执行流程&#xff1a;ARM_ROP && mprotect函数(运行内存权限修改) && [[ARM_ROP_csu_init]…

【AI自然语言处理应用】通义晓蜜CCAI

通义晓蜜CCAI-对话分析AIO 对话分析AIO&#xff0c;即对话分析all-in-one API&#xff0c;是基于深度调优的对话大模型&#xff0c; 为营销服类产品提供智能化升级所需的生成式摘要总结、质检、分析等能力的官方应用。 面向对象&#xff1a;开发者、自研企业、传统呼叫中心采购…

02 nth_element 与第k小

题目&#xff1a; 方案一&#xff1a;sort排序 #include<bits/stdc.h> using namespace std;int main() {int n;int k;cin>>n>>k;int a[n]{0};for(int i0;i<n;i){cin>>a[i];}sort(a,an); cout<<a[k]<<endl;}方案二&#xff1a;…

【机器学习(十一)】糖尿病数据集分类预测案例分析—XGBoost分类算法—Sentosa_DSML社区版

文章目录 一、XGBoost算法二、Python代码和Sentosa_DSML社区版算法实现对比(一) 数据读入和统计分析(二)数据预处理(三)模型训练与评估(四)模型可视化 三、总结 一、XGBoost算法 关于集成学习中的XGBoost算法原理&#xff0c;已经进行了介绍与总结&#xff0c;相关内容可参考【…

leetcode面试题17.04:消失的数字(C语言版)

思路1 先排序&#xff0c;再依次查找&#xff0c;如果下一个值不等于前一个1&#xff0c;那么下一个值就是消失数字。 时间复杂度分析&#xff1a;冒泡排序的时间复杂度为O(N^2)&#xff0c;qsort排序时间复杂度为O(N*logN)。因此该思路不可行。 思路2 求和0到N&#xff0c;再减…

Python爬虫使用实例-mdrama

一个Python爬虫使用实例&#xff1a;主要用于下载指定的剧集音频。分别从网页和json文件中获取剧集的title和剧集中所存在音频的id&#xff0c;调用you-get&#xff0c;最后自动重命名下载文件夹为剧集名title。 目标网址&#xff1a; https://www.missevan.com/mdrama/其中为…

【C++】关键字+命名空间

大家好&#xff0c;我是苏貝&#xff0c;本篇博客带大家了解C的命名空间&#xff0c;如果你觉得我写的还不错的话&#xff0c;可以给我一个赞&#x1f44d;吗&#xff0c;感谢❤️ 目录 一. 关键字二. 命名空间2.1 命名空间的定义2.2 命名空间的使用a. 命名空间名称作用域限定…

R包的安装、加载以及如何查看帮助文档

0x01 如何安装R包 一、通过R 内置函数安装&#xff08;常用&#xff09; 1.安装CRAN的R包 install.packages()是一个用于安装 R 包的重要函数。 语法&#xff1a;install.packages(pkgs, repos getOption("repos"),...) 其中&#xff1a; pkgs&#xff1a;要安…