这个题目是这样的:
一个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;
}