注意到m=1或者2,

当m=1时,是普通的最大连续字段和,只不过是k个:

设dp[i][j]表示前i个数中取出j个矩形的最大和

转移:

选:dp[i][j]=max{dp[l][j-1]+s[i]-s[l-1]}

不选:dp[i][j]=max(dp[i][j],dp[i-1][j])

复杂度O(n^2*K)

当m=2时,设f[i][j][k]表示第一列选到第i个数,第二列选到第j个数时,总共k个子矩形的答案

转移有4种情况

当这一位什么都不做的时候:f[i][j][k]=max(f[i-1][j][k],f[i][j-1][k])

当仅选取第一列的某段区间时:f[i][j][k]=max(f[l][j][k-1]+sum[i][1]-sum[l-1][1]) 1<=l<i

当仅选取第二列的某段区间时:f[i][j][k]=max(f[i][l][k-1]+sum[j][2]-sum[l-1][2]) 1<=l<j

当i==j时,可以选取两列一起的f[i][j][k]=max(f[l][l][k]+sum[i][1]+sum[i][2]-sum[l-1][1]-sum[l-1][2])

最后所有情况取max

复杂度O(n^3*K)

#include <cstdio>
using namespace std;
int n,m,k,x,y;
int dp[105][105][15];//dp[i][j][k] 表示 当前第一列选到第i个数,第二列选到第j个数,已经选出了k个子矩阵,的最大值
int s1[105],s2[105];//s1和s2数组分别表示第一列和第二列的前缀和数组
int max(int a,int b){return a>b?a:b;}

int main(){
    scanf("%d%d%d",&n,&m,&k);
    if(m==1){
        for(int i=1;i<=n;i++)scanf("%d",&x),s1[i]=s1[i-1]+x;
        for(int i=1;i<=k;i++)
            for(int j=1;j<=n;j++){
                dp[j][0][i]=dp[j-1][0][i];
                for(int skr=0;skr<j;skr++)dp[j][0][i]=max(dp[j][0][i],dp[skr][0][i-1]+s1[j]-s1[skr]);
            }
        printf("%d",dp[n][0][k]);
    }else{
        for(int i=1;i<=n;i++)scanf("%d%d",&x,&y),s1[i]=s1[i-1]+x,s2[i]=s2[i-1]+y;
        
        for(int i=1;i<=k;i++)
            for(int j=1;j<=n;j++)
                for(int sk=1;sk<=n;sk++){
                    dp[j][sk][i]=max(dp[j-1][sk][i],dp[j][sk-1][i]);
                    for(int st=0;st<j;st++)dp[j][sk][i]=max(dp[j][sk][i],dp[st][sk][i-1]+s1[j]-s1[st]);
                    for(int su=0;su<sk;su++)dp[j][sk][i]=max(dp[j][sk][i],dp[j][su][i-1]+s2[sk]-s2[su]);
                    if(j==sk)for(int sr=0;sr<j;sr++)dp[j][sk][i]=max(dp[j][sk][i],dp[sr][sr][i-1]+s1[j]+s2[sk]-s1[sr]-s2[sr]);
                }
    /*
    上面三个循环,
    第一个循环表示当前n*2矩阵中选k个矩阵,从1-k推,
    从dp数组的定义来看,第一个转移是,第一列的第j个数不选或者是第sk个数不选的最大值
    一次只操作一列!
    三重循环中的第一重循环 表示,操作第一列数,当前选到了第j个数,则可能第j到前面第st个数之间这个区域可以形成一个 构成答案中的一个矩阵
    第二重循环亦是如此
    最后如果j==sk,也就是第一列和第二列同时选取到同一行的数时,有可能第一列和第二列的第sr行到 第 j(sk)行 构成一个答案中的一个矩阵
    */
    printf("%d",dp[n][n][k]);
    }
    return 0;
}