差分差分矩阵

最新推荐文章于 2024-11-07 21:36:53 发布

你们姗姗来迟

最新推荐文章于 2024-11-07 21:36:53 发布

阅读量134

点赞数

文章标签：矩阵算法 c++

本文链接：https://blog.csdn.net/dcvpe/article/details/129699397

版权

文章介绍了如何利用差分技术处理整数序列和矩阵的更新操作。对于序列，每次将指定区间的每个元素加上一个常数，通过前缀和可以快速得到更新后的序列；对于矩阵，同样使用差分矩阵更新子矩阵元素，最后通过前缀和计算得出最终矩阵。

摘要由CSDN通过智能技术生成

1.差分

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

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

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

输入格式

第一行包含两个整数 nn 和 mm。

第二行包含 nn 个整数，表示整数序列。

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

输出格式

共一行，包含 nn 个整数，表示最终序列。

数据范围

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

输入样例：

6 3
1 2 2 1 2 1
1 3 1
3 5 1
1 6 1

输出样例：

3 4 5 3 4 2

#include<iostream>
using namespace std;

const int N=100010;
int a[N],b[N];
void insert(int x,int y,int c)
{
    b[x]+=c;
    b[y+1]-=c;
}
int main()
{
    int n,t;
    cin >> n >> t;
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);//输入
    for(int i=1;i<=n;i++) insert(i,i,a[i]);//初始化操作数组
    while(t--)
    {
        int x,y,num;
        scanf("%d %d %d",&x,&y,&num);
        insert(x,y,num);//差分
    }
    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;
}
/*
*主要利用前缀和求，当第i位的数+=C，i位以后的数也+=c，所以在j+1位-=c以后的数减去c与之前的c抵消;
由此在i-j区间内加得c。
*/

2.差分矩阵

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

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

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

输入格式

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

接下来 nn 行，每行包含 mm 个整数，表示整数矩阵。

接下来 qq 行，每行包含 55 个整数 x1,y1,x2,y2,cx1,y1,x2,y2,c，表示一个操作。

输出格式

共 nn 行，每行 mm 个整数，表示所有操作进行完毕后的最终矩阵。

数据范围

1≤n,m≤10001≤n,m≤1000,
1≤q≤1000001≤q≤100000,
1≤x1≤x2≤n1≤x1≤x2≤n,
1≤y1≤y2≤m1≤y1≤y2≤m,
−1000≤c≤1000−1000≤c≤1000,
−1000≤矩阵内元素的值≤1000−1000≤矩阵内元素的值≤1000

输入样例：

输出样例：

2 3 4 1
4 3 4 1
2 2 2 2

#include<iostream>
using namespace std;

const int N=1010;

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()
{
    int n,m,t;
    cin >> n >> m >> t;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            insert(i,j,i,j,a[i][j]);
        }
    }
    while(t--)
    {
        int x1,x2,y1,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]);
        }
        puts("");
    }
    return 0;
}