最大子矩阵

最大子矩阵可以化解为最长子序列解决，给出例子：
4
0 -2 -7 0
9 2 -6 2
-4 1 -4 1
-1 8 0 -2
最大子矩阵的长宽没有限制，可以先假设最长子矩阵的行数为2，将原矩阵的两行两行相加，即：
一二行相加：9 0 -13 2
二三行相加：5 3 -10 3
三四行相加：-5 9 -4 3
相加后按照最长子序列来解，同理，如果最大子矩阵行数为3，即：
一二三行相加：5 1 -17 3
二三四行相加：4 11 -10 1
再假设最大子矩阵行数为四：
一二三四行相加：5 9 -17 1
有了这些行相加的序列，再加上最长子序列的算法：

for(int i=1;i<=n;i++){
    dp[i]=max(a[i],dp[i-1]+a[i]);
}

相结合，可以求出最大子矩阵的和，代码如下：

#include<bits/stdc++.h>
using namespace std;
#define mem(a,b) memset(a,b,sizeof(a))

int dp[105];
int a[105][105],b[105][105];
int n,maxn;

void solve(int j){//最长子序列的算法 
	mem(dp,0);
	for(int i = 1;i <= n;i++){
		dp[i] = max(b[j][i],dp[i-1] + b[j][i]);
		maxn = max(maxn,dp[i]);
	}
}

int main(){
	while(~scanf("%d",&n)){
		maxn = 0;
		for(int i = 1;i <= n;i++){
			for(int j = 1;j <= n;j++){
				scanf("%d",&a[i][j]);
			}
		}
		for(int i = 1;i <= n;i++){
	        mem(b,0);//b数组用来存行相加的结果 
	        for(int j=i;j<=n;j++){//从第i行开始加，最多加到第n行 
	            for(int k=1;k<=n;k++){
	                b[j][k]=a[j][k]+b[j-1][k];//最长子序列的dp 
	            }
	            solve(j);
	        }
	    }
	    printf("%d\n",maxn);
	}
} 

代码参考自https://blog.csdn.net/m0_37579232/article/details/89195533