1768:最大子矩阵

1768:最大子矩阵

传送门

类似于求最大子序列和的思想，关键在于把一列(处于i、j行之间的列）看作一个数
用一个双重循环列举不同的行块，将行块的处理放在f函数中，等同于对普通一维数列的处理（即，将每列和视作一个数

最初原型的对一维数列a求最大子序列和，
dp[1]=a[1];
for(int i=2;i<=n;i++){
dp[i]=max(dp[i-1]+a[i],a[i]);
// dp[i]=max(dp[i-1]+a[i],dp[i-1]);突然迷糊了一下，最大子序列和 子序列必须要连续的呀，
//怎么能跳过a[i]不取还妄想带着a[i]之前的序列和继续呢
}

int a[MAX][MAX];
int col[MAX][MAX];//col[i][j]表示第j列上从第一行到第i行这i个数的和
int dp[MAX];

最后找dp数组中存储的最大子矩阵的和，一开始直接dp[n],离谱
再来找最后的dp数组中的最大值，也不可
在行块更新16次中要全程监控dp数组的最大值，
最大子矩阵可能出现在矩阵的任意部分

// cout<<dp[n];//求最长子序列的最大和 可不是 序列越长子序列和最大
//for(int i=0;i<=n;i++){
// cout<<dp[i]<<endl;
//}//dp数组是在不断更新了16次的，一位行块的上下界选取有16中取法，所以
//如果只是在这所有更新完之后的dp数组中找最大值就相当于默认最大子矩阵出现
//在最后一个列举的行块中，显然不对，在更新过程中要一直监控dp的值，
//最大子矩阵可能出现在矩阵的任意部分
//cout<<*max_element(dp,dp+MAX);

 #include<bits/stdc++.h>
using namespace std;
const int MAX=105;
	int res=1>>30; 
	int n;
	int a[MAX][MAX];
	int col[MAX][MAX];//col[i][j]表示第j列上从第一行到第i行这i个数的和 
	int dp[MAX];
void fun(int x,int y){
	for(int k=1;k<=n;k++){
		dp[k]=col[y][k]-col[x-1][k];//注意啦,减去的时x行以上的元素 
	//此时的dp数组就被初始化为了普通的一维数组，要做的就是求该一位数列的最大子序列和 	
	}
	for(int l=1;l<=n;l++){
		dp[l]=max(dp[l],dp[l-1]+dp[l]);
		res = max(res, dp[l]);
	} 
//这两个循环不能融为一个噢，dp数组一开始必须全部初始化，因为第二个循环
//不仅用到dp[l]，还要检索到dp中其他元素，比如dp[l-1]	
}
/*最初原型的对一维数列a求最大子序列和，
dp[1]=a[1];
for(int i=2;i<=n;i++){
	dp[i]=max(dp[i-1]+a[i],a[i]);
//	dp[i]=max(dp[i-1]+a[i],dp[i-1]);突然迷糊了一下，最大子序列和 子序列必须要连续的呀，
//怎么能跳过a[i]不取还妄想带着a[i]之前的序列和继续呢 
} 
*/ 
int main(){

	cin>>n;


	memset(col,0,sizeof(col));
//		memset(dp,0,sizeof(dp));
//for(int i=0;i<=n;i++){
//	dp[i]=1>>30;
//}

	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			cin>>a[i][j];
			col[i][j]=col[i-1][j]+a[i][j]; 
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			fun(i,j);
		}
	}
cout<<res;
    return 0;
}