给定一个n*m的矩阵A,求A中的一个非空子矩阵,使这个子矩阵中的元素和最大。
其中,A的子矩阵指在A中行和列均连续的一块。
输入
输入的第一行包含两个整数n, m,分别表示矩阵A的行数和列数。
接下来n行,每行m个整数,表示矩阵A。
1<=n, m<=500,A中每个元素的绝对值不超过5000。
输出
输出一行,包含一个整数,表示A中最大的子矩阵中的元素和。
样例输入 Copy
3 3 -1 -4 3 3 4 -1 -5 -2 8
样例输出 Copy
10
#include <iostream>
#include <bitset>
#include <algorithm>
#include <cmath>
#include <ctype.h>
#include <string>
#include <cstdio>
#include <set>
#include <cstdio>
#include <fstream>
#include <deque>
#include <vector>
#include <queue>
#include <map>
#include <stack>
#include <iomanip>
#define SIS std::ios::sync_with_stdio(false)
#define ll long long
#define INF 0x3f3f3f3f
const int mod = 1e9 + 7;
const double esp = 1e-5;
const double PI = 3.141592653589793238462643383279;
using namespace std;
const int N = 1e7 + 5;
const int maxn = 1 << 20;
ll powmod(ll a, ll b) { ll res = 1; a %= mod; while (b >= 0); for (; b; b >>= 1) { if (b & 1)res = res * a % mod; a = a * a % mod; }return res; }
ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }
/*void chafen(int l, int r, int k) {//差分函数
p[l] += k;
p[r + 1] -= k;
}*/
/*************************************************************/
int sum[505][505] = { 0 };
int a[505][505] = { 0 };
int main() {
int n, m;
int now = 0;
cin >> n >> m;
int ans = -INF;
for(int i=1;i<=n;i++)
for (int j = 1; j <= m; j++) {
cin >> a[i][j];
sum[i][j] = sum[i - 1][j] + a[i][j];//每一列的前缀和
ans = max(ans, a[i][j]);
}
if (ans <= 0)//全是负数
cout << ans << endl;
else {
for(int i=1;i<=n;i++)//i确定上边界
for (int j = 1; j <= i; j++) {//j确定下边界
now = 0;
for (int k = 1; k <= m; k++) {
now += (sum[i][k] - sum[j - 1][k]);
ans = max(ans, now);
if (now < 0)
now = 0;
}
}
cout << ans <<endl;
}
return 0;
}
/*
3 3
-1 -4 3
3 4 -1
-5 -2 8
*/