一二计划（Day 10）求最大子矩阵和

思路

看到这思路时我被惊讶到了一下，着实想不到，而且有点dp的味道
首先我们知道子矩阵和有行与列组成而且是随机的，而我们的目的是找出子矩阵和的最大值，那么就有一个问题怎么快速计算子矩阵和
————————————————————————
假设我们是N*M的矩阵，而我们要将N行它压缩成一个个一维数组
用求最大连续子序列和的操作，来找出此时最大子矩阵的和
————————————————————————

最大连续子序列和

最大连续子序列和就是给你一个一维数组，找出下标连续且和最大的值
代码为

for(int i=1;i<=n;i++)
{
    dp[i]=max(a[i],dp[i-1]+a[i]);
    //实际意义为
    //if(dp[i-1]<0) dp[i]=a[i]
    //if(dp[i-1]>0) dp[i]=dp[i-1]+a[i]
}

————————————————————————
知道了怎么求最大连续子序列和，也就知道了怎么找子矩阵列的和
接下来我们要做的就是将原矩阵的N行变成 多个一维数组
例如原矩阵为
第一行 0 -2 -7 0
第二行9 2 -6 2
第三行-4 1 -4 1
第四行-1 8 0 -2
我们先将 第一行与第二行相加 也是是为 9 0 -6 2
再根据最大连续子序列和 发现是{9，0}或{9}
也就是如果在第一第二行面 最大的子矩阵为（用加粗表示）
第一行 0 -2 -7 0
第二行 9 2 -6 2
第三行 -4 1 -4 1
第四行 -1 8 0 -2
或
第一行 0 -2 -7 0
第二行 9 2 -6 2
第三行 -4 1 -4 1
第四行 -1 8 0 -2
接下来依次枚举第一行到第三行 （到四行），然后第二行到第三行第四行，后面同理。
每次得出来得最大值与上次得出来得最大值进行比较
代码

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<string.h>
typedef long long ll;
ll a[110][110],b[110][110],dp[110];
ll mx=0;
int n;
//我们是枚举每次子矩阵和的行，然后求出在这种行的要求情况下最大的子矩阵和（也就是求列）
void solve(int j)//求最大连续子序列和
{
    memset(dp,0,sizeof(dp));
    int i;
    for(i=1;i<=n;i++)
    {
        dp[i]=fmax(b[j][i],dp[i-1]+b[j][i]);
        mx=fmax(mx,dp[i]);//这次的子矩阵最大和与之前子矩阵的最大和进行比较。
    }
}
int main()
{
    while(~scanf("%d",&n))
    {
        int i,j,k;
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=n;j++)
            {
                scanf("%lld",&a[i][j]);
            }
        }
        for(i=1;i<=n;i++)//枚举从第i行开始加
        {
            memset(b,0,sizeof(b));
            for(j=i;j<=n;j++)//加到第j行
            {
                for(k=1;k<=n;k++)//计算i行~j行上每列上值的总和，也就是将i~j行变成一行，然后用最大连续子序列和求出子矩阵的最大和
    			//即这个是用来计算子矩阵的行数，大连续子序列和函数用来计算最大子矩阵的列数
                {
                    b[j][k]=a[j][k]+b[j-1][k];
                }
                solve(j);
            }
        }
        printf("%lld\n",mx);
    }
}

10000小时计划
101h