最大子矩阵和 （DP + 降维处理）

这个题目是这样的：
一个M*N的矩阵，找到此矩阵的一个子矩阵，并且这个子矩阵的元素的和是最大的，输出这个最大的值。

例如：3*3的矩阵：

-1 3 -1
2 -1 3
-3 1 2

和最大的子矩阵是：

3 -1
-1 3
1 2

输入

第1行：M和N，中间用空格隔开（2 <= M,N <= 500)。
第2 - N + 1行：矩阵中的元素，每行M个数，中间用空格隔开。(-10^9 <= M[i] <= 10^9)

输出

输出和的最大值。如果所有数都是负数，就输出0。

输入示例

3 3
-1 3 -1
2 -1 3
-3 1 2

输出示例

7

我想起科幻小说《三体》中描绘了恢弘壮丽的“降维攻击”的场景：
“歌者”随手抛下了一张“二向箔”，整个银河系的三维空间奔腾汹涌地流入二向箔，塌缩成一个二维平面，三维结构被碾压在二维平面之上。同时，这一降维过程是全息的，所有的三维信息被保留在碾压后的二维空间里。这种致命的攻击令攻击者和被攻击者同归于尽，玉石俱焚，其结局黑暗得令人窒息。

我们知道，在线性空间的最大子段和，已有线性时间的DP算法。有没有可能，把二维的最大矩阵和转化为一维的最大子段和呢？

答案是肯定的，算法有着似三体描述的未来武器般的威力

对于一个子矩阵，有四个属性：起始行位置 i , 结束行位置 j, 起始列 x, 结束列 y

我们可以枚举i， j， 然后对后列进行降维压缩

再对降维后的一维数组进行最大子序列和动态规划就行了

#include <cstdio>#include <cstring>#include <algorithm>using namespace std;typedef long long ll;const ll INF = 0x3f3f3f3f3f3f3f3f;const int SIZE = 500 + 10;int gra[SIZE][SIZE];ll c[SIZE], dp[SIZE];int main() {  //freopen("in.txt", "r", stdin);  int m, n;  scanf("%d%d", &m, &n);  for (int i = 0; i < n; ++i) {    for (int j = 0; j < m; ++j) {      scanf("%d", &gra[i][j]);    }  }  ll ans = -INF;  for (int i = 0; i < n; ++i) {    memset(c, 0, sizeof(c));    for (int j = i; j < n; ++j) {      //对i到j行矩阵进行降维操作      for (int k = 0; k < m; ++k) {        c[k] += gra[j][k];      }      //对降维后的c[k]进行最大子序列和的动态规划      dp[0] = c[0];      if (ans < dp[0]) { ans = dp[0]; }      for (int k = 1; k < m; ++k) {        dp[k] = max(dp[k - 1] + c[k], c[k]);        if (ans < dp[k]) {ans = dp[k];}      }    }  }  printf("%lld\n", ans);  return 0;}