前缀和与差分
前缀和是什么?差分又是什么?
从1~k的和叫做k的前缀和
差分就是前缀和的逆运算
由此,可以写出基本的代码:
//这里是前缀和的初始化
inline void init() {
sum[1] = 1;
for (int i = 2; i <= n; ++i) {
sum[i] = sum[i - 1] + a[i];
}
}
//这里是前缀和查询操作
inline int find(int left,int right) {
return sum[right] - sum[left - 1];
}
//以下是差分操作
inline void init() {
for (int i = n; i >= 1; --i) {
b[i] = sum[i] - sum[i - 1];
}
}
以上全是一维的前缀和和差分
二维前缀和
根据容斥原理,
s
u
m
i
,
j
=
s
u
m
i
−
1
,
j
+
s
u
m
i
,
j
−
1
−
s
u
m
i
−
1
,
j
−
1
+
a
i
,
j
sum_{i,j} = sum_{i-1,j} + sum_{i,j-1} - sum_{i-1,j-1} + a_{i,j}
sumi,j=sumi−1,j+sumi,j−1−sumi−1,j−1+ai,j
problem&code
P1719 最大加权矩形
题目描述
为了更好的备战NOIP2013,电脑组的几个女孩子LYQ,ZSC,ZHQ认为,我们不光需要机房,我们还需要运动,于是就决定找校长申请一块电脑组的课余运动场地,听说她们都是电脑组的高手,校长没有马上答应他们,而是先给她们出了一道数学题,并且告诉她们:你们能获得的运动场地的面积就是你们能找到的这个最大的数字。
校长先给他们一个N*N矩阵。要求矩阵中最大加权矩形,即矩阵的每一个元素都有一权值,权值定义在整数集上。从中找一矩形,矩形大小无限制,是其中包含的所有元素的和最大 。矩阵的每个元素属于[-127,127],例如
0 –2 –7 0
9 2 –6 2
-4 1 –4 1
-1 8 0 –2
在左下角:
9 2
-4 1
-1 8
和为15。
几个女孩子有点犯难了,于是就找到了电脑组精打细算的HZH,TZY小朋友帮忙计算,但是遗憾的是他们的答案都不一样,涉及土地的事情我们可不能含糊,你能帮忙计算出校长所给的矩形中加权和最大的矩形吗?
输入格式
第一行:n,接下来是n行n列的矩阵。
输出格式
最大矩形(子矩阵)的和。
输入输出样例
输入
4
0 -2 -7 0
9 2 -6 2
-4 1 -4 1
-1 8 0 -2
输出
15
说明/提示
n
≤
120
n \leq 120
n≤120
很水,一看就知道是二维差分。
#include <bits/stdc++.h>
using namespace std;
const int inf = 2.1e9;
const int maxn = 150;
const int maxm = 20;
int n, m, t, matrix[maxn][maxn], ans = -inf, temp[maxn], dp[maxn];
template<typename Type>inline Type maxx(Type x, Type y) {return x < y ? y : x;}
inline void sum() {
memset(dp, 0, sizeof dp);
for (int i = 1; i <= n; ++i) {
dp[i] = maxx(dp[i], dp[i - 1] + temp[i]);
ans = maxx(ans, dp[i]);
}
}
inline void Do_Funcitional() {
for (int i = 1; i <= n; ++i) {
memset(temp, 0 ,sizeof temp);
for (int j = i; j <= n; ++j) {
for (int k = 1; k <= n; ++k) temp[k] += matrix[j][k];
sum();
}
}
}
signed main() {
scanf("%d", &n);
for (int i = 1; i <= n; ++i) {
for (int j = 1; j <= n; ++j) scanf("%d", &matrix[i][j]);
}
Do_Funcitional();
printf("%d\n", ans);
return 0;
}