#include<stdio.h>
#define N 105
int n,map[N][N];
void input()
{
scanf("%d",&n);
for(int i = 1;i<=n;i++)
for(int j = 1;j<=n;j++)
{
scanf("%d",&map[i][j]);
}
}
int cal_array_max(int a[])
{
int max = a[1],temp = a[1];
for(int i = 2;i<=n;i++)
{
if(temp>0)
temp += a[i];
else
temp = a[i];
if(temp > max)
max = temp;
}
return max;
}
int result(int i,int j)
{
int temp[N] = {0};
for(int k = 1;k<=n;k++)
for(int l = i;l<=j;l++)
{
temp[k] += map[l][k];
}
return cal_array_max(temp);
}
void search()
{
int max = -1000,temp;
for(int i = 1;i<=n;i++)
for(int j = i;j<=n;j++)
{
temp = result(i,j);
if(temp>max)
max = temp;
}
printf("%d\n",max);
}
int main()
{
input();
search();
return 0;
}
还是动态规划的思想。
这道题比较巧妙。他不是直接使用动态规划的思想,而是先将问题转化,转化后再使用动态问题解决转化问题。
为了一个矩阵的最大子阵。解决这道题有点暴力。它以行为单位,result(i,j)表示第i行到第j行间(矩阵的高度必须是j-i+1)的最大子矩阵。即result(i,j)表示的最大子矩阵必须是严格的从第i行开始,在第j行结束,不能在第i和第j行的中间开始或结束。这样,我们最后只需要利用冒泡排序的思想,对起始行和结尾行进行一个遍历,来找到最大值。
而计算第i行到第j行间的最大子矩阵result(i,j),我们将二维的矩阵转化为一维的数列。这只需要我们将矩阵的每一列之和求出来,然后分别放在对应的数组里。
接下来只需要找到数组中数列的子列使的子列的和最大即可。这里,我们就用到了动态规划的方法来求。
我们建立决策:设temp[j]表示以a[j]结尾的子序列的最大和(必须包含a[j])。这个子问题并不是指前j个元素中的最大子列和。而是包含a[j]的最大连续子序列的和。我们只需要求出这样的序列,然后再求出temp数组中的最大值,既可求出最大连续子序列的和。
状态转移方程:temp[j] = max { a[j] , temp[j-1] + a[j] } 。(当temp[j-1]<0时,取a[j]单独。)
根本思想还是对第j个元素a[j]取还是不取得问题,如果取(即可以续接到前面的队伍中),则继续连接。如果不取,则断开。用这样的决策来保证连续。