前缀与差分模版
前缀和与差分
前缀和
题目:
输入一个长度为 n 的整数序列。
输入一对 l,r。输出原序列中从第l 个数到第 r 个数的和。
例如:
- 整数序列为:
2 1 3 6 4
- 输入
1 2
输出3
- 输入
1 3
输出6
- 输入
2 4
输出10
思路:
- 原数组1,a[2],a[3]...
- 前缀和数组:s[1],s[2],s[3]...
- s[1]=a[1]
- s[2]=a[1]+a[2]
- s[3]=a[1]+a[2]+a[3]
- 前缀和数组
S[i]=a[1]+a[2]+a[3]...+a[i]
- 数组a[l,r]的和为
S[r]-S[l-1]=a[l]+...+a[r]
- 规定
S[0]=0,S[i]=S[i-1]+a[i]
代码:
#include<bits/stdc++.h>
using namespace std;
const int N=100010;
int n,l,r;
int a[N],s[N];
int main(){
scanf("%d",&n);
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];
scanf("%d%d",&l,&r);
printf("%d\n",s[r]-s[l-1]);
return 0;
}
如何只使用一个数组存储:
for(int i = 1; i <= n; ++i)
{
scanf("%d", &arr[i]);
arr[i]+=arr[i-1];
}
二维前缀和
思路:
二维前缀和预处理公式:s[i,j]=s[i,j-1]+s[i-1,j]-s[i-1,j-1]+arr[i,j]代表的是(1,1)到(i,j)的元素之和
矩形(x1,y1)到(x2,y2)的所有和为:s[x2,y2]-s[x2,y1-1]-s[x1-1,y2]+s[x1-1,y1-1]
代码:
#include<bits/stdc++.h>
using namespace std;
const int N=1010;
int n,m,q;
int a[N][N],s[N][N];
int main(){
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<=n;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;
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
printf("%d\n",s[x2][y2]-s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1]);//输出和
}
return 0;
}
差分
题目
输入一个长度为 n 的整数序列。接下来输入 m 个操作,每个操作包含三个整数 l,r,c,表示将序列中 [l,r]之间的每个数加上 c。请你输出进行完所有操作后的序列。
- 输入格式:第一行包含两个整数 n 和 m。第二行包含 n 个整数,表示整数序列。接下来 m 行,每行包含三个整数 l,r,c,表示一个操作。
- 输出格式:共一行,包含 n 个整数,表示最终序列。
输入样例:
6 3
1 2 2 1 2 1
1 3 1
3 5 1
1 6 1
输出样例:
3 4 5 3 4 2
思路
原数组
1
,s[2],s[3]
差分数组
1
,a[2],a[3]
a[1]=s[1]-s[0]
a[2]=s[2]-s[1]
a[3]=s[3]-s[2]
差分数组公式:a[i]=s[i]-s[i-1]
事实上:原数组为就是差分数组的前缀和数组
s[1]=a[1]
s[2]=a[1]+a[2]
s[i]=a[1]+...+a[i]
将[l,r]区间内的每一个数都加上一个常数C,那么对于差分数组的处理是a[1]+=c,a[r+1]-=c
代码
#include<bits/stdc++.h>
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(){
scanf("%d%d",&n,&n);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
for(int i=1;i<=n;i++) insert(i,i,a[i]);
while(m--){
int l,r,c;
scanf("%d%d%d",&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++) printf("%d",b[i]);
return 0;
}
二维差分
思路
将(x1,y1)到(x2,y2)内的每个元素都加上一个常数c,对于矩阵b来说有:
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[x2 + 1][y1] -= c;
b[x1][y2 + 1] -= c;
b[x2 + 1][y2 + 1] += c;
}
int main()
{
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 ++ )
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];
for (int i = 1; i <= n; i ++ )
{
for (int j = 1; j <= m; j ++ ) printf("%d ", b[i][j]);
puts("");
}
return 0;
}