一维差分、二维差分算法

797. 差分
输入一个长度为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
输入样例：
6 3
1 2 2 1 2 1
1 3 1
3 5 1
1 6 1
输出样例：
3 4 5 3 4 2

前缀和：a₁、a₂…a_n

差分：b₁、b₂…b_n

a是b的前缀和，b是a的差分，举个例子：
a₃=b₁+b₂+b₃
a_n=b₁+…+b_n

我们想要对a数组的[ L,r ]进行操作，比如加上c，意思就是，a_L+c，a_L+1+c…a_r+c

如果我们直接循环对a进行操作，时间复杂度为O(n)，而如果我们借用差分数组(即b数组)来操作的话，时间复杂度只有O(1)

怎样借助b数组呢？首先，我们让b_L加上c，这时，a_L也会加上c，同理，a_L+1…a_n都会加上c

现在我们发现a[ L,n ]都已经加上c了，可我们想要的是a[ L,r ]区间加上c，那怎么办呢？再来操作一步，就是让b_r+1减c，那么a[ r+1，n ]先加c后减c，相当于不变，还是原来的a，而a[L,r]就是成功的加c了

总结一下就是只需要对b数组进行两部修改即可：
b[L]+=c;
b[r+1]-=c;

还有一点，就是，我们刚开始初始化a数组为0，那么b数组自然也为0，然后再通过insert操作构造差分数组b

#include<iostream>
using namespace std;
const int N=100010;

int n,m;
int a[N],b[N];

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

int main()
{
   cin>>n>>m;
   for(int i=1;i<=n;i++) cin>>a[i];
   
   for(int i=1;i<=n;i++) insert(i,i,a[i]);  //构造差分数组b[i] 
   
   while(m--)
   {
   	  int l,r,c;
   	  cin>>l>>r>>c;
   	  insert(l,r,c);
   }	
   
   for(int i=1;i<=n;i++) b[i]+=b[i-1]; //此时b[i]就是前缀和 
   
   for(int i=1;i<=n;i++) cout<<b[i]<<" ";
   return 0;
}

下面这样写更好理解…

#include<iostream>
using namespace std;
const int N = 1e5+10;
int n, m;
int a[N], b[N];
int main()
{
    cin >> n >> m;
    for(int i=1; i<=n; i++)
    {
        cin >> a[i];
        b[i] = a[i] - a[i-1];  //构造差分数组
    }
    
    while(m--)
    {
        int l, r, c;
        cin >> l >> r >> c;
        b[l] += c;
        b[r+1] -= c;
    }
    
    for(int i=1; i<=n; i++) 
    {
        b[i] += b[i-1]; 
        cout << b[i] <<" ";
    }        
}

---

798. 差分矩阵
输入一个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≤100000,
1≤x1≤x2≤n,
1≤y1≤y2≤m,
−1000≤c≤1000,
−1000≤矩阵内元素的值≤1000
输入样例：
3 4 3
1 2 2 1
3 2 2 1
1 1 1 1
1 1 2 2 1
1 3 2 3 2
3 1 3 4 1
输出样例：
2 3 4 1
4 3 4 1
2 2 2 2

二维差分与一位差分类似，需要对差分数组的四个数进行修改：
b[x1][y1]+=c;
b[x2+1][y1]-=c;
b[x1][y2+1]-=c;
b[x2+1][y2+1]+=c;

#include<iostream>
using namespace std;
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[x1][y2+1] -=c;
    b[x2+1][y1] -=c;
    b[x2+1][y2+1] +=c;
}
int main()
{
    cin.tie(0);
    cin >> n >> m >>q;
    
    for(int i=1; i<=n; i++)
     for(int j=1; j<=m; j++)
      {
          cin >> a[i][j];
          insert(i,j,i,j,a[i][j]);
      }
      
    while(q--)
    {
        int x1, y1, x2, y2, c;
        cin >> 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];
          cout << b[i][j] <<" ";
      }
      cout << endl;
    }     
}

这一步要是没看懂，就看下面这个博客… b[i][j]+=b[i-1][j]+b[i][j-1]-b[i-1][j-1]😉
求前缀和、求子矩阵的和详解

---

啊要哭了，辛辛苦苦写的博客刷新一下没了，又重新写了一遍😭😭