前缀和计算

简介: 笔记

一维数组前缀和


14天阅读挑战赛


s[]数组用来储存


#include<iostream>
using namespace std;
#define N 100
int main()
{
  int arr1[N];
  int s[N];
  int n, m;
  scanf("%d %d", &n, &m);
  for (int i = 1; i <= n; ++i)
  scanf("%d", &arr1[i]);
  for (int i = 1; i <= n; ++i)
  s[i] = s[i - 1] + arr1[i];
  while (m--)
  {
  int l, r;
  scanf("%d %d", &l, &r);
  printf("%d", s[r] - s[l - 1]);
  }
}

l和r确定一个区间,打印这个区间内的数组和


二维数组前缀和


1.png2.png

#include<iostream>
using namespace std;
const int N = 1010;
int main()
{
  int n, m, q;
  scanf("%d%d%d", &n, &m, &q);//确定二维数组大小
  int arr1[N][N];
  int s[N][N];
  for (int i = 1; i <= n; i++)
  {
  for (int j = 1; j <= m; j++)
  {
    scanf("%d", &arr1[i][j]);
  }
  }
  for (int i = 1; i <= n; i++)
  {
  for (int j = 1; j <= m; j++)
  {
    s[i][j] = s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1] + arr1[i][j];//某点的区域面积,确定大区域
  }
  }
  while (q--)
  {
  int x1, x2, y1, y2;//某个特定的小区域(在大区域中的小区域)
  scanf("%d%d%d%d", &x1, &y1, &x2, &y2);//x2,y2是最右边的点
  printf("%d\n", s[x2][y2] - s[x1 - 1][y2] - s[x2][y1 - 1] + s[x1 - 1][y1 - 1]);
  }
  return 0;
}

s[x][y]表示的是红色区域的数字和,加了俩次的数字要减去

3.png

(x1,y1),(x2,y2)为某区域左上角和右下角坐标,求这俩个坐标围城的矩形里数字和


4.png



差分数组


5.png

思路:


首先给定一个原数组a:a[1], a[2], a[3],,,,,, a[n];


然后我们构造一个数组b : b[1] ,b[2] , b[3],,,,,, b[i];使得 a[i] = b[1] + b[2 ]+ b[3] +,,,,,, + b[i]


     也就是说,a数组是b数组的前缀和数组,反过来我们把b数组叫做a数组的差分数组。换句话说,每一个a[i]都是b数组中从头开始的一段区间和。


考虑如何构造差分b数组?


最为直接的方法


如下:


a[0 ]= 0;


b[1] = a[1] - a[0];


b[2] = a[2] - a[1];


b[3] =a [3] - a[2];


........


b[n] = a[n] - a[n-1];


问题:


给定区间[l ,r ],让我们把a数组中的[ l, r]区间中的每一个数都加上c,即 a[l] + c , a[l+1] + c , a[l+2] + c ,,,,,, a[r] + c;


a数组是b数组的前缀和数组,比如对b数组的b[i]的修改,会影响到a数组中从a[i]及往后的每一个数。


首先让差分b数组中的 b[l] + c ,a数组变成 a[l] + c ,a[l+1] + c,,,,,, a[n] + c;


然后我们打个补丁,b[r+1] - c, a数组变成 a[r+1] - c,a[r+2] - c,,,,,,,a[n] - c;


如果这里不打补丁,会影响这个区间后面的值

6.png

所以要打补丁

7.png



#include<iostream>
using namespace std;
const int N=100010;
    int n,m;
    int a[N];
    int b[N];
void insert(int l,int r,int c)
{
    b[l]+=c;
    b[r+1]-=c;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;++i)
    scanf("%d",&a[i]);
    for(int i=1;i<=n;++i)
    insert(i,i,a[i]);
    while(m--)
    {
        int l,r,c;
        scanf("%d%d%d",&l,&r,&c);
        insert(l,r,c);//插入后b直接就是差分数组
    }
    for(int i=1;i<=n;++i)
    {
        b[i]+=b[i-1];//求的是整个数组和
    }
      for(int i=1;i<=n;++i)
    {
        printf("%d ",b[i]);
    }
    return 0;
}

8.png


二维数组差分


9.png

数组b是数组a的差分,一维可以利用差分给其中一段数据加一个值,二维可以是某个矩阵加一个值,如给蓝框包围起来的矩阵加一个值

10.png

11.png

用bx1y1+=c,如果给蓝框部分加上C,则黄框-蓝框剩下的部分,也会加C,我们要减掉

12.png




黑色部分为bx2+1,y1-=C

13.png14.png

粉色部分为bx1,y2+1-=C


但右下角被减了俩次,我们加回来一次即可+x2+1,y2+1


差分数组bij满足性质:aij是bij的前缀和



const int N=1010;
int n, m, q;
int a[N][N], b[N][N];
void insert(int x1, int y1, int x2, int y2,int c)
{
  b[x1][y1] += c;
  b[x2 + 1][y1] -= c;
  b[x1][y2 + 1] -= c;
  b[x2 + 1][y2 + 1] += c;
}
int main()
{
  scanf("%d %d %d", &n, &m, &q);
  for (int i = 1; i <=n; ++i)
  {
  for (int j = 1; j <= m; ++j)
  {
    scanf("%d", &a[i][j]);
  }
  }
  for (int i = 1; i <=n; ++i)
  for (int j = 1; j <= m; ++j)
    insert(i, j, i, j, a[i][j]);
  while (q--)
  {
  int x1, y1, x2, y2, c;
  scanf("%d %d %d %d %d", &x1, &y1, &x2, &y2, &c);
  insert(x1, y1, x2, y2, c);
  }
  for (int i = 1; i <= n; ++i)
  for (int j = 1; j <= m; ++j)
    b[i][j] += b[i - 1][j] + b[i][j - 1] - b[i - 1][j - 1];
  for (int i = 1; i <= n; ++i)
  for (int j = 1; j <= m; ++j)
    printf("%d ", b[i][j]);
  return 0;
}


相关文章
|
编解码 算法 容器
AVFormatContext封装层:理论与实战(二)
AVFormatContext封装层:理论与实战(二)
118 0
Axure高保真原型设计:移动端多选图片上传
本文介绍了如何在Axure中利用中继器实现移动端应用的多选图片上传功能,适用于如微信、微博等社交平台。文章详细描述了主页面、相册页面和大图页面的制作步骤,并展示了如何通过中继器和交互设置实现图片的选择、上传及删除等功能。此教程有助于提升用户体验和应用功能性。
254 9
|
8月前
|
Ubuntu 机器人 开发工具
具身智能干货|ROS2理论与实践系列(一):ROS2概述与环境搭建
《ROS2理论与实践》系列课程主要由基础篇、核心篇、应用篇、进阶篇以及项目库五部分组成。本阶段为《ROS2理论与实践——核心篇》课程,核心篇课程设计以官方内容为标准,主要介绍ROS2中的通信机制与开发者工具,其中前者是整个ROS2框架中的核心构成,而后者则为开发者能够高效的构建应用程序提供有力支持。本阶段课程目的是帮助大家快速上手ROS2,为后续进阶奠定基础。
640 12
具身智能干货|ROS2理论与实践系列(一):ROS2概述与环境搭建
|
数据处理 调度 开发者
解密Python的异步编程:协程与事件循环的实战应用
在现代应用程序开发中,异步编程已经成为提高性能和响应速度的关键技术。Python的异步编程通过协程和事件循环提供了高效处理并发任务的能力。本文将深入探讨Python中异步编程的核心概念,包括协程的基本用法、事件循环的工作机制以及如何在实际项目中应用这些技术。通过对比同步和异步编程的性能差异,读者将能够理解异步编程的优势,并学会如何在Python中实现高效的异步任务处理。
水仙花数
水仙花数。
228 7
|
弹性计算 Oracle Ubuntu
阿里云免费试用云服务器上新了,最高配置4核16G,最长3个月免费试用
近日,阿里云免费试用云服务器又上新了,前几天还只有共享型s6实例和通用算力型u1实例可申请试用,现在又增加了通用型g7、计算型c7a、计算型c8y和计算型c7实例的云服务器参与免费试用,现在共计有11款不同配置和实例的云服务器可免费申请试用了。
1571 0
阿里云免费试用云服务器上新了,最高配置4核16G,最长3个月免费试用
|
传感器 自动驾驶 安全
无人驾驶出租车对网约车的影响将会是怎样的呢?
无人驾驶出租车对网约车的影响将会是怎样的呢?
|
存储 弹性计算 固态存储
阿里云服务器是如何收费的?阿里云服务器各收费项目收费标准参考
阿里云服务器收费标准包括实例价格、预留实例券价格、专有宿主机、块存储价格、存储容量单位包、带宽价格、快照服务价格等,云服务器价格主要由云服务器配置费用+磁盘价格+网络宽带价格,配置指的是云服务器的实例规格和cpu与内存配置,本文为大家分享一下2023年阿里云服务器所有收费项目的最新收费标准,以表格形式展示给大家,以供参考。
阿里云服务器是如何收费的?阿里云服务器各收费项目收费标准参考
|
JavaScript
html+css+js实现文本编辑器
html+css+js实现文本编辑器
174 0