首页 > 其他分享 >前缀和和差分归纳总结

前缀和和差分归纳总结

时间:2024-10-16 18:18:10浏览次数:9  
标签:y2 前缀 归纳 int cf 差分 x2 y1 x1

前缀和数组可以在O(1)的时间内求得某一区间中的所有数据的和

差分数组可以在O(1)的时间内对某一区间中的所有数据进行加减操作

原数组求差分及为差分数组,差分数组再求前缀和即为原数组

一维前缀和:

设原数组为a[N],前缀和数组为s[N],数组下标都从1开始存储

每个s[i]等于a[1] + …… + a[i]

对于前缀和数组s[i]

s[l - 1] = a[1] + a[2] + ... + a[l - 1]

s[l] = a[1] + a[2] + ... + a[l]

s[r] = a[1] + a[2] + ... + a[l] + ... + a[r]

那么

a[l ~ r] = s[r] - s[l - 1]

一维前缀和数组求法:

for (int i = 1; i <= n; i ++ )
	s[i] = s[i - 1] + a[i];

如果为了节省空间,输入的数据直接保存在了s数组中,那么有:

for (int i = 1; i <= n; i ++ ) s[i] += s[i - 1];

那么对于区间[l, r]中的数据的和,即为:s[r] - s[l - 1]


例题:前缀和

输入一个长度为 n的整数序列。

接下来再输入 m个询问,每个询问输入一对 l,r。

对于每个询问,输出原序列中从第 l个数到第 r个数的和。

输入格式

第一行包含两个整数 n和 m。

第二行包含 n个整数,表示整数数列。

接下来 m行,每行包含两个整数 l和 r,表示一个询问的区间范围。

输出格式

共 m行,每行输出一个询问的结果。

数据范围

1 ≤ l ≤ r ≤ n
1 ≤ n, m ≤ 1e5
−1000 ≤ 数列中元素的值 ≤ 1000

代码:

本题预处理前缀和数组的时间是O(n),求每个区间和的时间是O(1),一共的询问时间是O(m)

总时间复杂度为O(n + m)

#include <iostream>

using namespace std;

const int N = 1e5 + 10;

int n, m;
int s[N];

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

一维差分:

cf[i] = a[i] - a[i - 1];

对于每个cf[i]

cf[1] = a[1] - a[0];
cf[2] = a[2] - a[1];
cf[3] = a[3] - a[2];
...
cf[n] = a[n] - a[n - 1];

那么对cf数组求前缀和,就是原数组a

如果cf[i] += 1,那么a[i]以及之后的元素都会加上1

如果cf[j] -= 1,那么a[j]以及之后的元素都会减去1

因此如果我们想对区间[l, r]中的所有的元素都加上x,那么有:

cf[l] += x, cf[r + 1] -= x;

作用:快速将某一区间[l, r]内的所有元素都加上数x

时间复杂度由O(n)变为了O(1)

差分数组求法:

设原数组为a[N],差分数组为cf[N],数组下标都从1开始存储

先将差分数组初始化为0,memset(cf, 0, sizeof cf);

for (int i = 1; i <= n; i ++ )
	cf[i] = a[i] - a[i - 1];

也可以在输入元素时分别插入每个元素:

void insert(int l, int r, int x)
{
    cf[l] += x;
    cf[r + 1] -= x;
}

for (int i = 1; i <= n; i ++ )
{
	scanf("%d", &temp);
    insert(i, i, temp);
}

例题:差分

输入一个长度为 n的整数序列。

接下来输入 m个操作,每个操作包含三个整数 l,r,c,表示将序列中 [l,r]之间的每个数加上 c。

请你输出进行完所有操作后的序列。

输入格式

第一行包含两个整数 n和 m。

第二行包含 n个整数,表示整数序列。

接下来 m行,每行包含三个整数 l,r,c,表示一个操作。

输出格式

共一行,包含 n个整数,表示最终序列。

数据范围

1 ≤ n, m ≤ 100000
1 ≤ l ≤ r ≤ n
−1000 ≤ c ≤ 1000
−1000 ≤ 整数序列中元素的值 ≤ 1000

代码:

#include <iostream>

using namespace std;

const int N = 1e5 + 10;

int n, m;
int cf[N];

void insert(int l, int r, int x)
{
    cf[l] += x;
    cf[r + 1] -= x;
}

int main()
{
    scanf("%d%d", &n, &m);
    for (int i = 1; i <= n; i ++ )
    {
        int temp;
        scanf("%d", &temp);
        insert(i, i, temp);
    }
    
    while (m -- )
    {
        int l, r, x;
        scanf("%d%d%d", &l, &r, &x);
        insert(l, r, x);
    }
    
    // 求前缀和,将差分数组转换为原数组
    for (int i = 1; i <= n; i ++ ) cf[i] += cf[i - 1];
    
    for (int i = 1; i <= n; i ++ ) printf("%d ", cf[i]);
    
    return 0;
}

二维前缀和:

设原数组为a[N][N],前缀和数组为s[N][N],数组下标都从1开始存储

在这里插入图片描述

s[i][j]前缀和数组为从(0,0)到(i,j)区间中的元素和

s[i][j]数组中区间元素和比作区间面积:

边界为(x1, y1)(x2, y2)的区域中的元素和即为区域4的面积s4,

可以发现:s4 = (s1 + s2 + s3 + s4) - (s1 + s3) - (s1 + s2) + s1

s1 + s2 + s3 + s4 = s[x2][y2]

s1 + s3 = s[x2][y1 - 1]

s1 + s2 = s[x1 - 1][y2]

s1 = s[x1 - 1][y1 - 1]

那么可以得出:

二维区间和:

a[x1 ~ x2][y1 ~ y2] = s[x2][y2] - s[x2][y1 - 1] - s[x1 - 1][y2] + s[x1 - 1][y1 - 1];

二维前缀和数组求法:

间接求:

for (int i = 1; i <= n; i ++ )
    for (int j = 1; j <= m; j ++ )
        s[i][j] = s[i][j - 1] + s[i - 1][j] - s[i - 1][j - 1] + a[i][j];

直接求:

for (int i = 1; i <= n; i ++ )
    for (int j = 1; j <= m; j ++ )
    {
        scanf("%d", &s[i][j]);
        s[i][j] += s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1];
    }
例题:子矩阵的和

输入一个 n行 m列的整数矩阵,再输入 q个询问,每个询问包含四个整数 x1,y1,x2,y2,表示一个子矩阵的左上角坐标和右下角坐标。

对于每个询问输出子矩阵中所有数的和。

输入格式

第一行包含三个整数 n,m,q。

接下来 n行,每行包含 m个整数,表示整数矩阵。

接下来 q行,每行包含四个整数 x1,y1,x2,y2,表示一组询问。

输出格式

共 q行,每行输出一个询问的结果。

数据范围

1 ≤ n, m ≤ 1000
1 ≤ q ≤ 2e5
1 ≤ x1 ≤ x2 ≤ n
1 ≤ y1 ≤ y2 ≤ m
−1000 ≤ 矩阵内元素的值 ≤ 1000

代码:
#include <iostream>

using namespace std;

const int N = 1010;

int n, m, q;
int s[N][N];

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", &s[i][j]);
            s[i][j] += s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1];
        }
    
    while (q -- )
    {
        int x1, y1, x2, y2;
        scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
        
        int res = s[x2][y2] - s[x2][y1 - 1] - s[x1 - 1][y2] + s[x1 - 1][y1 - 1];
        printf("%d\n", res);
    }
    
    return 0;
}

二维差分:

设原数组为a[N][N],差分数组为cf[N][N],数组下标都从1开始存储

和分析一维差分的方法类似,结合二维前缀和数组理解

cf[x1][y1] += 1使得原数组中的(x1, n)(y1,m)区间内的数全部加1

为了只使得特定区间(x1, y1)(x2, y2)中的数据加1,还需要抵消其他区间的影响

使得cf[x2 + 1][y1] -= 1, cf[x1][y2 + 1] -= 1, cf[x2 + 1][y2 + 1] -= 1

二维差分数组求法:

void insert(int x1, int y1, int x2, int y2, int x)
{
    cf[x1][y1] += x;
    cf[x2 + 1][y1] -= x;
    cf[x1][y2 + 1] -= x;
    cf[x2 + 1][y2 + 1] += x;
}

for (int i = 1; i <= n; i ++ )
    for (int j = 1; j <= m; j ++ )
    {
        int t;
        scanf("%d", &t);
        insert(i, j, i, j, t);
    }

例题:差分矩阵

输入一个 n行 m列的整数矩阵,再输入 q个操作,每个操作包含五个整数 x1,y1,x2,y2,c,其中 (x1,y1)和 (x2,y2)表示一个子矩阵的左上角坐标和右下角坐标。

每个操作都要将选中的子矩阵中的每个元素的值加上 c。

请你将进行完所有操作后的矩阵输出。

输入格式

第一行包含整数 n,m,q。

接下来 n行,每行包含 m个整数,表示整数矩阵。

接下来 q行,每行包含 5个整数 x1,y1,x2,y2,c,表示一个操作。

输出格式

共 n行,每行 m个整数,表示所有操作进行完毕后的最终矩阵。

数据范围

1 ≤ n, m ≤ 1000
1 ≤ q ≤ 1e5
1 ≤ x1 ≤ x2 ≤ n
1 ≤ y1 ≤ y2 ≤ m
−1000 ≤ c ≤ 1000
−1000 ≤ 矩阵内元素的值 ≤ 1000

代码:

#include <iostream>

using namespace std;

const int N = 1010;

int n, m, q;
int cf[N][N];

void insert(int x1, int y1, int x2, int y2, int x)
{
    cf[x1][y1] += x;
    cf[x1][y2 + 1] -= x;
    cf[x2 + 1][y1] -= x;
    cf[x2 + 1][y2 + 1] += x;
}

int main()
{
    scanf("%d%d%d", &n, &m, &q);
    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
        {
            int t;
            scanf("%d", &t);
            insert(i, j, i, j, t);
        }
    
    while (q -- )
    {
        int x1, y1, x2, y2, x;
        scanf("%d%d%d%d%d", &x1, &y1, &x2, &y2, &x);
        insert(x1, y1, x2, y2, x);
    }
    
    for (int i = 1; i <= n; i ++ )
    {
        for (int j = 1; j <= m; j ++ )
        {
            // cf数组的前缀和即为原数组
            cf[i][j] += cf[i - 1][j] + cf[i][j - 1] - cf[i - 1][j - 1];
            printf("%d ", cf[i][j]);
        }
        puts("");
    }
    
    return 0;
}

标签:y2,前缀,归纳,int,cf,差分,x2,y1,x1
From: https://blog.csdn.net/a123_120/article/details/142988647

相关文章

  • 差分笔记
    差分笔记差分可看作前缀和的逆运算对于一个数组a[n]有:a[0]a[1]a[2]a[3]......a[n-2]a[n-1]构造一个差分数组b[n]使得其中每一项都为数组a每项的差:b[0]=a[0]b[1]=a[1]-a[0]......b[n-2]=a[n-2]-a[n-3]b[n-1]=a[n-1]-a[n-2]不难看出b的前缀和为a中的每一项a[n]......
  • 洛谷题单指南-字符串-P1470 [USACO2.3] 最长前缀 Longest Prefix
    原题链接:https://www.luogu.com.cn/problem/P1470题意解读:求s最长前缀长度,使得可以拆解成p集合中的字符串解题思路:动态规划:设s字符串下标从1开始,p集合用set<string>保存所有的元素状态表示:设f[i]表示前i个字符s[0~i-1]是否能拆解成p中的元素状态计算:对于j=i-1开始往后倒......
  • 三维前缀和
    二维前缀和模板:预处理:dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]+map[i][j];查询:dp[x2][y2]+dp[x1-1][y1-1]-dp[x1-1][y2]-dp[x2][y1-1]三维前缀和模板:预处理:pre[i][j][k]=mp[i][j][k]-pre[i-1][j-1][k]-pre[i-1][j][k-1]-pre[i][j-1][k-1]+pre[i][j][k-1]+pre[i][j-1][......
  • 高数导数积分知识点归纳
    高数公式函数1.基本公式:2.常见函数类型线性函数(a和b是常数):多项式函数(ai是常数):指数函数(a>0,且a≠1)对数函数(a>0, 且a≠1)三角函数:如正弦函数f(x)=sin(x),余弦函数f(x)=cos(x),正切函数f(x)=tan(x)如反正弦函数f(x)=arcsin(x),反余弦函数f(x)=arccos(x),反正切......
  • 每日OJ题_牛客_比那名居的桃子_滑动窗口/前缀和_C++_Java
    目录牛客_比那名居的桃子_滑动窗口/前缀和题目解析C++代码Java代码牛客_比那名居的桃子_滑动窗口/前缀和比那名居的桃子(nowcoder.com)描述:        小红有一天看到了一只桃子,由于桃子看上去就很好吃,小红很想把它吃掉。已知吃下桃子后,每天可以获得ai​的......
  • 【C++差分数组】P1672何时运输的饲料
    本文涉及知识点C++差分数组C++算法:前缀和、前缀乘积、前缀异或的原理、源码及测试用例包括课程视频P1672何时运输的饲料原文比较啰嗦,我简述一下:第x天运来F1(1<=F1<=1e6)千克的饲料,第D(1<=2e3)天还剩F2(1<=F2<=F1)千克饲料,某人养了C头牛,moves[i]={comi,leavei},表......
  • 前缀和笔记
    前缀和笔记对于一个一维数组a[m]其前i项和记作s[i]如果想要对a[m]中任意连续段的值进行求和,例如a[l]~a[r]则可使用前缀和数组进行o(n)计算inta[m],s[m];s[0]=0;//定义s[0]的值,防止边界问题for(inti=0;i<m;i++){ cin>>a[i]; s[i]=s[i-1]+a[i];}这样的话,s......
  • 二维前缀和 & 二维差分
    二维前缀和求二维前缀和后,能够实现\(O(1)\)求原数组二维区间和,但是不支持修改。lln,m,sum2[N][N],c[N][N];voidSum2_pre(){fr(i,1,n)fr(j,1,m)sum2[i][j]=sum2[i-1][j]+sum2[i][j-1]-sum2[i-1][j-1]+c[i][j];}llSum2(llx......
  • ARC136C Circular Addition [高维前缀和]
    Description给定一个长度为\(n\)的正整数序列\(A\),求有多少对\((i,j)\)使得\(A_i+A_j\)不发生进位操作。\(A_i<10^6\)。Solution显然对于每个\(A_i\),设\(B_i=999999-A_i\),那么\(A_i\)可以和所有位上的数都小于等于\(B_i\)对应位上的数的数匹配,考虑用桶加前缀......
  • 算法修炼之路之前缀和
    目录一:一维前缀和二:二维前缀和 三:LeetCodeOJ练习  1.第一题2.第二题 3.第三题 4.第四题5.第五题6.第六题一:一维前缀和这里通过例题来引出牛客_DP34【模板】前缀和画图分析:具体代码:#include<iostream>#include<vector>usingnamespacestd;int......