动态规划——最大连续子序列和&最大子矩阵

前言

最大连续子序列和是动态规划中最经典的问题之一。在一个给定的序列{A1，A2，...，An}中，找出一个连续的子序列{Ai，...，Aj}，使得这个连续的子序列的和最大，输出这个最大的序列和。

 

题目1

给出一个整数序列S，其中有N个数，定义其中一个非空连续子序列T中所有数的和为T的“序列和”。 对于S的所有非空连续子序列T，求最大的序列和。 变量条件：N为正整数，N≤1000000，结果序列和在范围（-2^63,2^63-1）以内。

输入描述:

第一行为一个正整数N，第二行为N个整数，表示序列中的数。

输出描述:

输入可能包括多组数据，对于每一组输入数据，
仅输出一个数，表示最大序列和。

示例1

输入

5
1 5 -3 2 4

6
1 -2 3 4 -10 6

4
-3 -1 -2 -5

输出

9
7
-1

思路

令dp[i]表示以A[i]作为末尾的连续序列的最大和。于是，最大连续子序列和便是数组dp中的最大值。
有两种情况：
①最大和的连续序列只有一个元素，即A[i]本身。
②最大和的连续序列有多个元素，即从前面的某个A[j]开始，一直到A[i]结束，也就是dp[i]=A[j]+...+A[i-1]+A[i]。而A[j]+...+A[i-1]=dp[i-1],即dp[i]=dp[i-1]+A[i]。
由于只有这两种情况，于是得到状态转移方程dp[i]=max{A[i],dp[i-1]+A[i]}。
只需将i从小到大枚举并依次遍历，即可得到整个dp数组。接着输出该数组中的最大值，即最大连续子序列的和。

#include<iostream>
#include<cstdio>
using namespace std;
const int MAXN=1000000;
long long A[MAXN];
long long dp[MAXN];
long long MaxSubsequence(int n)
{
    long long maximum=-999999;
    for(int i=0;i<n;i++)
    {
        if(i==0)
            dp[i]=A[i];
        else
            dp[i]=max(A[i],dp[i-1]+A[i]);
        maximum=max(maximum,dp[i]);
    }
    return maximum;
}
int main()
{
    int n;
    while(cin>>n)
    {
        for(int i=0;i<n;i++)
            cin>>A[i];
        long long answer=MaxSubsequence(n);
        cout<<answer<<endl;
    }
    return 0;
}

 

最大子矩阵

现在来了解最大连续子序列和的二维情况，即如何求一个二维矩阵的最大子矩阵和。

假设原二维矩阵中最大子矩阵所在的行是从i到j，那么只会出现下面这两种情况：

①当i=j时，求最大子矩阵和就转换成了求第i行元素的最大连续子序列和。

②当i！=j时，把从第i行到第j行的所有行的元素加起来，得到只有一行的一维数组，这个一维数组的最大连续子序列和便是最大子矩阵和。

（不必将原始矩阵中的第i行到第j行的所有行的元素次累加起来得到一维数组，而事先用一个辅助二维矩阵记录原始矩阵从上到下加起来的累加矩阵，于是求从第i行到第j行的一维数组只需将辅助矩阵进行一行减法便可得到，而不需要逐行进行多次累加。）

只需从小到大枚举并依次遍历i和j，接着输出所有子矩阵中的最大值，即求得最大子矩阵和。

题目2

已知矩阵的大小定义为矩阵中所有元素的和。给定一个矩阵，你的任务是找到最大的非空(大小至少是1 * 1)子矩阵。 比如，如下4 * 4的矩阵 0 -2 -7 0 9 2 -6 2 -4 1 -4 1 -1 8 0 -2 的最大子矩阵是 9 2 -4 1 -1 8 这个子矩阵的大小是15。

输入描述:

输入是一个N * N的矩阵。输入的第一行给出N (0 < N <= 100)。
再后面的若干行中，依次（首先从左到右给出第一行的N个整数，再从左到右给出第二行的N个整数……）给出矩阵中的N2个整数，整数之间由空白字符分隔（空格或者空行）。
已知矩阵中整数的范围都在[-127, 127]。

输出描述:

测试数据可能有多组，对于每组测试数据，输出最大子矩阵的大小。

示例1

输入

4
0 -2 -7 0
9 2 -6 2
-4 1 -4  1
-1 8  0 -2

输出

15

#include<iostream>
#include<cstdio>
using namespace std;
const int MAXN=100;
int matrix[MAXN][MAXN];//原始矩阵
int total[MAXN][MAXN];//辅助矩阵
int arr[MAXN];//一维数组
int dp[MAXN];
int MaxSubsequence(int n)
{
    int maximum=-128;
    for(int i=0;i<n;i++)
    {
        if(i==0)
            dp[i]=arr[i];
        else
            dp[i]=max(arr[i],dp[i-1]+arr[i]);
        maximum=max(maximum,dp[i]);
    }
    return maximum;
}
int MaxSubmatrix(int n)
{
    int maximal=0;
    for(int i=0;i<n;i++)
    {
        for(int j=i;j<n;j++)
        {
            for(int k=0;k<n;k++)//获得一维数组
            {
                if(i==0)
                    arr[k]=total[j][k];
                else
                    arr[k]=total[j][k]-total[i-1][k];
            }
            int current=MaxSubsequence(n);
            maximal=max(maximal,current);
        }
    }
    return maximal;
    
}
int main()
{
    int n;
    while(cin>>n)
    {
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                cin>>matrix[i][j];
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
            {
                if(i==0)
                    total[i][j]=matrix[i][j];
                else
                    total[i][j]=total[i-1][j]+matrix[i][j];
            }
        int answer=MaxSubmatrix(n);
        cout<<answer<<endl;
    }
    return 0;
}