DAY4 前缀和与差分

本文中将要介绍前缀和与差分的简单知识，全部题目来自acwing，是对y总讲解的简单归纳。

题目：acwing795

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

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

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

输入格式

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

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

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

输出格式

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

数据范围

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

简单审题，会觉得题目不难，简单用两个for循环就能结束战斗，但是在提交的时候会显示超时。

两个for循环的实践复杂度为O（n^2），所以我们就要优化算法，降低时间复杂度，而前缀和的复杂度为O（n）。

#include<iostream>
using namespace std;
const int z = 1e6 + 10;
int sum[z];
int main() {
	int n, m; sum[0] = 0;
	cin >> n >> m;
	for (int i = 1; i <= n; i++) {
		int x;
		cin >> x;
		sum[i] = sum[i - 1] + x;
	}
	while (m--) {
		int l, r;
		cin >> l >> r;
		cout << sum[r] - sum[l - 1] << endl;
	}
	return 0;
}

本段代码只采用了一个for循环，用sum[i]来储存前i个数的和。在查询区间和的时候，答案就为sum[r]-sum[l-1](包括区间端点)。

前缀和的算法简言之就是降低时间复杂度，把存储前项改为存储前n项的和。

题目：acwing796

输入一个 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

#include<iostream>
using namespace std;
const int N = 1010;
int n, m, q;
int a[N][N], s[N][N];
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 = 0; j <= m; j++) {
			s[i][j] = s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1] + a[i][j];
		}}
	while (q--) {
			int x1, y1, x2, y2;
			cin >> x1 >> y1 >> x2 >> y2;
			cout << s[x2][y2] - s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1]<<endl;}
}

本题中s[i][j]代表以（0，0）为左上角坐标，（i，j）为右下角坐标为矩阵的和。

画一下图就知道s[i][j]的更新方式为是s[i][j] = s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1] + a[i][j]。

如果是以（x1,y1）为左上角坐标，（x2,y2）为右下角坐标，那么

s=s[x2][y2] - s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1]。简单画一下图就能分析出来

题目：acwing797

输入一个长度为 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

对于我一个之前没有学习过算法的萌新来说，直接用for循环把序列上的每一个数加上c。

在初步代码编写中我采用while循环中嵌套for循环，此时时间复杂度为O（n*m），对于大量数据输入会超时。此时就需要用差分的思想，降低时间复杂度为O（n+m）。

#include<iostream>
using namespace std;
const int n = 1e6 + 10;
int a[n]; int b[n];
int main() {
	int n, m;
	cin >> n >> m;
	for (int i = 1; i <= n; i++) cin >> a[i];
	b[0] = a[0] = 0;
	for (int i = 1; i <= n; 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++) {
		a[i] = a[i - 1] + b[i];
		cout << a[i ] << " ";
	}
}

用a[i]来存储数据，b[i]存储a[i]与a[i-1]的差。

值得注意的是while循环中对b[]的加减操作，由于b[i]为差，所以b[i]+c就相当于a[l]及其以后的数+c

同理可知b[r+1]-c就相当于a[r+1]及其以后的数-c。

此时时间复杂度由O（n*m）降低为O（n+m）。

题目：acwing798 

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

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

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

输入格式

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

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

接下来 q行，每行包含 55 个整数 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

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

本题中以（x1,y1）为左上角坐标，以（x2,y2）为右下角坐标。用a[i][j]来存储(i,j)的数值。

我们可以构建一个差分矩阵b，其定义为：b[i][j]=a[i][j]-a[i-1][j]-a[i][j-1]+a[i-1][j-1]。

insert函数实际上就是对b数组执行插入操作，等价于对a数组中的(x1,y1)到（x2,y2）中的元素都+c。在第一处insert函数就是把b数组构建成a的差分数组（相当于把（i，j）到（i，j）上加上a[i][j]）,第二处insert函数就是执行题目的意思。

值得注意的是b[i][j]的更新方式：b[i][j] += b[i - 1][j] + b[i][j - 1] - b[i - 1][j - 1](比较难理解)

也可以用其他的更新方式：a [i][j] = b[i][j]+ a[i-1][j] + a[i][j-1] - a[i-1][j-1]