最大子矩阵

最近在看动归问题，碰到了最大子矩阵问题，本人愚钝，想了半天，也翻了半天博客，始终没找到（其实是没理解）简单好用的方法。于是干脆不看别人的了，自己思考一下。结果就想出了一个时间复杂度略高的方法，先将就这用，等慢慢熟练了，我想，简单方法应该一眼就能看懂了。我的思路是;定义一个二维数组dp[][],dp[i][j],表示i*j子矩阵的和。

例如：1 2 3

          1 2 3

          1 2 3

则dp[1][1]=1;dp[1][2]=3;dp[1][3]=6;dp[2][2]=6;代码如下：

<span style="font-size:18px;">for(int i=0;i<=n;i++)
{
	dp[0][i]=dp[i][0]=0;
}
for(int i=1;i<=n;i++)
{
	for(int j=1;j<=n;j++)
	{
		scanf("%d",&a[i][j]);
		dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]+a[i][j];
	}
}</span>

之后就是利用dp[][]寻找最大子矩阵了，我的想法比较暴力。因为是矩阵，所以i,j至少是2.拿上面的例子说事，当i==2时，没什么说的，最大肯定是第一和第二行的和。当i==3时。有可能是1,2,3行之和，也有可能是2,3行之和。不用考虑1,2行，因为在i==2时已经考虑过了。这样分析的话。还需要一个k。从0开始，因为矩阵行数大于等于2.所以k<=i-2;这样的话。利用dp[i][j]-dp[k][j];就能把所有符合条件的行的组合遍历一遍，从而找出最大。看代码：

<span style="font-size:18px;">int Max=0;
int x=0;
for(int i=2;i<=n;i++)
{
	for(int j=2;j<=n;j++)
	{
		for(int k=0;k<=i-2;k++)
		{
			x=dp[i][j]-dp[k][j];
			if(x>Max)
			<span style="white-space: pre;">	</span>Max=x;
		}
	}
}</span>

感觉这方法写起来和理解起来都不难，就是时间效率的问题。还是得再好好想想。