POJ 1050 最大子矩阵和

POJ 1050 最大子矩阵和

题目：http://poj.org/problem?id=1050
题意：给一个N*N的矩阵，求一个子矩阵，使得该子矩阵中的数字之和为所有子矩阵中最大的。
思路：这里给出了两种思路：
1.枚举所有可能的组合。
2.dp,其实这是根据一维的最大子段和改编的题目

枚举：
首先计算出所有的A[0][0]至A[i][j]这样的子矩阵的和，当然不是一个个的计算，有一点点小的技巧，具体参见代码。
然后利用四重循环枚举所有可能的矩阵的左上角顶点(i,j)和右下角顶点(ii,jj),则这个子矩阵的和即为A[ii][jj]-A[i-1][jj]-A[ii][j-1]+A[i-1][j-1].这样找出所有的和中最大的值即为所求。
代码：

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int maxn=107;
int a[maxn][maxn];

int main()
{
    int n;
    scanf("%d",&n);
    memset(a,0,sizeof(a));
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    {
        scanf("%d",&a[i][j]);
        if(i==0) a[i][j]+=a[i][j-1];
        else if(j==0) a[i][j]+=a[i-1][j];
        else a[i][j]+=(a[i-1][j]+a[i][j-1]-a[i-1][j-1]);
    }
    //枚举所有可能
    int ans=a[1][1];
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    for(int ii=i;ii<=n;ii++)
    for(int jj=j;jj<=n;jj++)
        if(ii==i&&jj==j) continue;
        else ans=max(ans,a[ii][jj]-a[i-1][jj]-a[ii][j-1]+a[i-1][j-1]);
    printf("%d\n",ans);
}

dp思路：
其实就是每次选择两行i,j，并将第i至第j行间的元素全部压缩至一行中，这样问题就变成了一维数组的最大子段和问题，其状态转移方程为
p[i]=p[i-1]>0? p[i-1]+p[i]:p[i];
注意要枚举所有合理的i,j组合，故每次选择的i和j的顺序十分重要，此代码中的i,j枚举顺序可以避免重复运算，减少时间开销。
代码：

#include <cstring>
#include <iostream>
using namespace std;
const int maxn=107;
int a[maxn][maxn];

int main()
{
    int n;
    scanf("%d",&n);
    int ans=-128;
    int tmp;
    for(int i=1;i<=n;i++)
    {
        tmp=0;
        for(int j=1;j<=n;j++)
        {
            scanf("%d",a[i]+j);
            if(tmp<0) tmp=a[i][j];
            else tmp+=a[i][j];
            if(tmp>ans)
                ans=tmp;
        }
    }
    for(int i=1;i<n;i++)            //每次选择i,j两行之间的合并为一行,注意i,j选择的顺序很重要，这
        for(int j=i+1;j<=n;j++)     //样选择可以尽量减少运算
        {
            tmp=0;
            for(int k=1;k<=n;k++)
            {
                a[i][k]+=a[j][k];
                if(tmp<0) tmp=a[i][k];
                else tmp+=a[i][k];
                if(tmp>ans) ans=tmp;
            }
        }
    printf("%d\n",ans);
}