问题描述:
求一个M*N的矩阵的最大子矩阵和。
比如在如下这个矩阵中:
0 -2 -7 0
9 2 -6 2
-4 1 -4 1
-1 8 0 -2
拥有最大和的子矩阵为:
9 2
-4 1
-1 8
其和为15。
#include <iostream>
using namespace std;
int A[55][55];
int main(){
int n,m,ans;
cin >> n >> m;
ans = -1005;
for(int i=0;i<n;i++){
for(int j=0;j<m;j++){
cin >> A[i][j];
}
}
for(int i=0;i<n;i++){
for(int j=i;j<n;j++){
for(int l=0;l<m;l++){
for(int k=l;k<m;k++){
int tmp = 0;
for(int p=i;p<=j;p++){
for(int q=l;q<=k;q++){
tmp += A[p][q];
}
}
if(tmp > ans){
ans = tmp;
}
}
}
}
}
cout << ans << endl;
return 0;
}
上边代码的话就是在n行到m列中求出每次选定的 i行到j行 的 l列到k列 进而压缩成了
一维的最大子序和问题。
其实关键的是边界的问题,边界问题懂了,自然迎刃而解。
更更更关键的是先把最大子序和这个问题搞清楚,然后再来看这道最大子矩阵就会更助于理解了。
下面的代码只是把边界换了换,从而变成了在n行到m列中求出每次选定的 l列到k列 的 i行到j行。
#include <iostream>
using namespace std;
int ans[55][55];
int main() {
int n, m;
cin >> n >> m;
for (int i = 0; i < n; i++) {
for (int j = 0; j < m; j++) {
cin >> ans[i][j];
}
}
int Maxsum = INT_MIN;
for(int l = 0; l < m; l++){
for (int k = l; k < m; k++) {
for (int i = 0; i < n; i++) {
for (int j = i; j < n; j++) {
int Thissum = 0;
for (int p = l; p <= k; p++) {
for (int q = i; q <= j; q++) {
Thissum += ans[q][p];
}
}
if (Thissum > Maxsum) {
Maxsum = Thissum;
}
}
}
}
}
cout << Maxsum << endl;
return 0;
}