to the max

这个严格意义上不能算是我写的吧；

因为刚开始我完全没有想出来这道题该怎么解，是看了网上的解题才知道该怎么解的，现在做题完全就是在被虐呀；

回到正文，其实这道题的解题思路就如下面这张图一样；

下面这张图是别人博客的，大家可以去看看他的，他那里有详细的解释：

http://blog.csdn.net/piaoyi0208/article/details/7728972

假设最大子矩阵的结果为从第r行到k行、从第i列到j列的子矩阵，如下所示(ari表示a[r][i],假设数组下标从1开始)：
  | a11 …… a1i ……a1j ……a1n |
  | a21 …… a2i ……a2j ……a2n |
  |  .     .     .    .    .     .    .   |
  |  .     .     .    .    .     .    .   |
  | ar1 …… ari ……arj ……arn |
  |  .     .     .    .    .     .    .   |
  |  .     .     .    .    .     .    .   |
  | ak1 …… aki ……akj ……akn |
  |  .     .     .    .    .     .    .   |
  | an1 …… ani ……anj ……ann |

#include<iostream>
#include<string.h>
using namespace std;
int dp[101][101][101],a[101][101];
int main()
{
	int n;
	while(cin>>n)
	{
		memset(dp,0,sizeof(dp));
		memset(a,0,sizeof(a));
		int i,j,k,max1=-128;
		for(i=1;i<=n;i++)//把二维数组压缩成一维数组，
		{
			for(j=1;j<=n;j++)
			{
				cin>>a[i][j];
				//dp[i][i][j]=a[i][j];
			}
		}
		for(i=1;i<=n;i++)//类似于最大字段和
		{
			for(j=1;j<=n;j++)
			{
				for(k=j+1;k<=n;k++)
				{
					dp[j][k][i]=dp[j][k-1][i]+a[k][i];
				}
			}
		}
		for(i=1;i<=n;i++)
		{
			for(j=i;j<=n;j++)
			{
				for(k=2;k<=n;k++)
				{
					if(dp[i][j][k-1]>0)
						dp[i][j][k]+=dp[i][j][k-1];
				}
				for(k=1;k<=n;k++)
				{
					if(dp[i][j][k]>max1)
						max1=dp[i][j][k];
				}
			}
		}
		cout<<max1<<endl;
	}
}

在写完了这道题后我又去南阳理工学院写了一道差不多的题，发现上面的代码还是有问题，但是他居然还是让我AC。。。。

看来是抽风了。。。。。。。。。

附上南阳理工学院OJ上的代码：（这道题的输出输出不太一样，还有的他的行列是不一样的）

 
#include<iostream>
#include<string.h>
using namespace std;
int dp[101][101][101],a[101][101];
int main()
{
	int n;
	cin>>n;
	while(n--)
	{
		int hang,lie;
		cin>>hang>>lie;
		memset(dp,0,sizeof(dp));
		memset(a,0,sizeof(a));
		int i,j,k,max1=-99999999;
		for(i=1;i<=hang;i++)
		{
			for(j=1;j<=lie;j++)
			{
				cin>>a[i][j];
				dp[i][i][j]=a[i][j];//这个地方是关键。。
			}
		}
		for(i=1;i<=lie;i++)
		{
			for(j=1;j<=hang;j++)
			{
				for(k=j+1;k<=hang;k++)
				{
					dp[j][k][i]=dp[j][k-1][i]+a[k][i];
					//cout<<dp[j][k][i]<<" ";
				}
			}
		}
		for(i=1;i<=hang;i++)
		{
			for(j=i;j<=hang;j++)
			{
				for(k=2;k<=lie;k++)
				{
					if(dp[i][j][k-1]>0)
						dp[i][j][k]+=dp[i][j][k-1];
				}
				for(k=1;k<=lie;k++)
				{
					if(dp[i][j][k]>max1)
						max1=dp[i][j][k];
				}
			}
		}
		cout<<max1<<endl;
	}
}
        

附上测试数据一组：

1

2     2

-1     2

如果是第一个代码输出的是2；（这个答案明显是错的呀O.O）；

第二个代码输出的是4：

为什么不一样呢，你们好好看看这两个代码不一样的地方吧