前缀和 差分

一维前缀和

其实可以把它理解为数学上的一维数列的前n项和。
对于一个给定的数列 A, 它的前缀和数列S 是通过递推能求出来得

         S[i] = ∑A[j] ( j ∈ [1,i] )

部分和,即数列A中某个下标区间内的数的和,可以表示为前缀和相减的形式:

   sum(l,r) = ∑A[i] ( i ∈ [l,r] ) = S[r] - S[l-1]

模板：

S[i] = a[1] + a[2] + ... a[i]
a[l] + ... + a[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≤100000 ,
−1000≤数列中元素的值≤1000
输入样例：

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

输出样例

3
6
10

代码：

#include<iostream>
using namespace std;
const int N=100010;
int n,m;
int a[N],s[N]; 
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++)   s[i]=s[i-1]+a[i]; //前i项和
	while(m--)
	{
		int l,r;
		scanf("%d%d",&l,&r);
		printf("%d",s[r]-s[l-1]); 
	}
	return 0;
}

二维前缀和

与一维前缀和类似，设s[x][y]表示所有a[i][j]的。（1≤i≤x,1≤j≤y）
有一点像“矩形的面积”那样，把一整块区域的值都加起来。

红色区域之和可视为 总和 - 蓝色区域 - 绿色区域 + 黄色区域；
即 S[x2, y2] - S[x1 - 1, y2] - S[x2, y1 - 1] + S[x1 - 1, y1 - 1]

Sx2y2的求法：

Sij 为蓝色+绿色-重合部分+黄色；
S[i j] =S[i-1, j ]+S[i,j-1]-S[i-1,j-1]+a[i,j];
模板

S[i, j] = 第i行j列格子左上部分所有元素的和
以(x1, y1)为左上角，(x2, y2)为右下角的子矩阵的和为：
S[x2, y2] - S[x1 - 1, y2] - S[x2, y1 - 1] + S[x1 - 1, y1 - 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≤200000 ,
1≤x1≤x2≤n ,
1≤y1≤y2≤m ,
−1000≤矩阵内元素的值≤1000
输入样例：

3 4 3
1 7 2 4
3 6 2 8
2 1 2 3
1 1 2 2
2 1 3 4
1 3 3 4

输出样例：

17
27
21

代码：

#include<iostream>

using namespace std;
const int N=1010;
int n,m,q;
int  a[N][N],s[N][N];

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

一维差分

对于一个数列A，他的差分数列B定义为：
B[1] = A[1] ,B[i] = A[i] - A[i-1] ( 2<=i<=n)
容易发现，“前缀和” 和 “差分” 是一对互逆运算，差分序列B的前缀和序列就是原序列A，前缀和序列S的差分序列也是原序列A。
模板：

给区间[l, r]中的每个数加上c：B[l] += c, B[r + 1] -= c

例题：差分
输入一个长度为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的区间 [l,r] 加上 c（即把 Al，Al+1……Ar 都加上 c），其差分序列 B 的变化为 Bl 加 c，Br+1 减 c，其余位置不变。

代码：

#include<iostream>

using namespace std;
const int N=100010;
int a[N],b[N];

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

int main()
{
	int n,m;
	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]);
	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];
	
	for(int i=1;i<=n;i++)  cout<<b[i]<<" ";
	
	return 0;
 } 

二维差分

模板：

给以(x1, y1)为左上角，(x2, y2)为右下角的子矩阵中的所有元素加上c：
S[x1, y1] += c, S[x2 + 1, y1] -= c, S[x1, y2 + 1] -= c, S[x2 + 1, y2 + 1] += c

思路：

例题：差分矩阵
输入一个 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

代码

#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[x2+1][y1]-=c;
	b[x1][y2+1]-=c;
	b[x2+1][y2+1]+=c;
}

int main()
{
	cin>>n>>m>>q;
	
	for(int i=1;i<=n;i++)
	   for(int j=1;j<=m;j++)
	   cin>>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;
		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][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++)
	    cout<<b[i][j]<<" ";
	    
	
	return 0;
}