前缀和与差分

前缀和与差分

前缀和是什么？差分又是什么？

从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];
  }
}

以上全是一维的前缀和和差分

二维前缀和

根据容斥原理，
$su{m}_{i,j}=su{m}_{i-1,j}+su{m}_{i,j-1}-su{m}_{i-1,j-1}+a_{i,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\le 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;
}