2024.02.04二维前缀和差分

课堂内容

二维前缀和 s i , j = Σ 1 i Σ 1 j a x , y s_{i,j}=\Sigma_1^i\Sigma_1^ja_{x,y} si,j=Σ1iΣ1jax,y
构建: s i , j = a i , j + s i − 1 , j + s i , j − 1 − s i − 1 , j − 1 s_{i,j}=a_{i,j}+s_{i-1,j}+s_{i,j-1}-s_{i-1,j-1} si,j=ai,j+si1,j+si,j1si1,j1 (容斥原理)
查询: s x 2 , y 2 − s x 2 , y 1 − 1 − s x 1 − 1 , y 2 + s x 1 − 1 , y 1 − 1 s_{x_2,y_2} - s_{x_2,y_1 - 1} - s_{x_1 - 1,y_2} + s_{x_1 - 1,y_1 - 1} sx2,y2sx2,y11sx11,y2+sx11,y11
画图推导

#include <bits/stdc++.h>

using namespace std;
const int N = 2e3 + 10;
int a[N][N], s[N][N];
int n, m;

int main() {
    cin >> n >> m;
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            cin >> a[i][j];
            s[i][j] = s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1] + a[i][j];
        }
    }
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            cout << s[i][j] << " ";
        }
        cout << endl;
    }
    int q;
    cin >> q;
    while (q--) {
        int x1, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        cout << s[x2][y2] - s[x2][y1 - 1] - s[x1 - 1][y2] + s[x1 - 1][y1 - 1] << endl;
    }

    return 0;
}

前缀和的逆过程
二维差分 t i , j = a i , j − a i − 1 , j − a i , j − 1 + a i − 1 , j − 1 t_{i,j} = a_{i,j} - a_{i - 1,j} - a_{i,j - 1} + a_{i - 1,j - 1} ti,j=ai,jai1,jai,j1+ai1,j1
修改:
t x 1 , y 1 + = n u m t_{x_1,y_1} += num tx1,y1+=num (增加影响)
t x 1 , y 2 + 1 − = n u m t_{x_1,y_2 + 1} -= num tx1,y2+1=num (消除影响)
t x 2 + 1 , y 1 − = n u m t_{x_2 + 1,y_1} -= num tx2+1,y1=num (消除影响)
t x 2 + 1 , y 2 + 1 + = n u m t_{x_2 + 1,y_2 + 1} += num tx2+1,y2+1+=num (消除多了一块,要减掉)
得到修改后的数组: a i , j = Σ 1 i Σ 1 j t x , y = t i , j + t i − 1 , j + t i , j − 1 − t i − 1 , j − 1 a_{i,j}=\Sigma_1^i\Sigma_1^jt_{x,y}=t_{i,j}+t_{i-1,j}+t_{i,j-1}-t_{i-1,j-1} ai,j=Σ1iΣ1jtx,y=ti,j+ti1,j+ti,j1ti1,j1

#include <bits/stdc++.h>

using namespace std;
const int N = 2e3 + 10;
int a[N][N], s[N][N], t[N][N];
int n, m;

int main() {
    cin >> n >> m;
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            cin >> a[i][j];
            t[i][j] = a[i][j] - a[i - 1][j] - a[i][j - 1] + a[i - 1][j - 1];
            cout << t[i][j] << " ";
        }
        cout << endl;
    }

    int q;
    cin >> q;
    while (q--) {
        int x1, y1, x2, y2, num;
        cin >> x1 >> y1 >> x2 >> y2 >> num;
        t[x1][y1] += num;
        t[x1][y2 + 1] -= num;
        t[x2 + 1][y1] -= num;
        t[x2 + 1][y2 + 1] += num;
    }

    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            a[i][j] = t[i][j] + a[i - 1][j] + a[i][j - 1] - a[i - 1][j - 1];
            cout << a[i][j] << " ";
        }
        cout << endl;
    }

    return 0;
}

题目讲解

P3397 地毯
核心代码:

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

作业

https://www.luogu.com.cn/contest/154442#problems
外加P8218 P3406

  • 23
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值