NC20242 [SCOI2005]最大子矩阵(多维DP)

题目链接

题意：
$给一个n\ast m的矩阵$
$找出k个不相交矩阵使和最大$
$求出这个最大和$
题解：
$n<=100,m<=2,k<=10$
$看到这个题先想到的就是一维的最大k个连续子序列和$
$如果是2个可以进行预处理求值$
$但是如果求k个子序列，那肯定需要状态转移$
$列一个dp[n][k]，前n个数分为k组$
$转移方程为$
$dp[i][k]=max(dp[i][k],dp[j][k-1]+sum[i]-sum[j])$
$这个dp下来需要O(n^{2}k)，但其实有更快的方法$
$但是在这里已经够了，更快的方法可以详见HDU1024$
$然后联想到这个题上，看怎么样能转移到这个题$

$首先看到这个数据范围，最先发现的肯定是这个m的范围$
$m<=2，最多只有两列$
$那说明，对于每一块矩阵，要不宽为1，要不宽为2$
$那我们完全可以求出两列的最大k个连续子序列和，然后进行合并$
$同时维护两列，那就可以列出dp[n][n][k]$
$在第一列选i个，第二列选j个，总共k组的最大值$
$对于每一列的求法，我们用上述的转移方程进行转移$
$现在需要考虑的是合并的做法$
$就是出现如果有两列可以并在一起比单独一列大的情况$
$这种情况肯定是在i，j相等的时候$
$因为是一个矩阵，所以找一下上边$
$转移方程为$
$dp[i][j][k]=max(dp[i][j][k],dp[p][p][k-1],s[i]-s[p]$
$此时的s表示的是两列共同的和，这时候就进行了对两列的合并$
$然后按照这个方法进行递推即可，详见代码$
AC代码

/*
    Author:zzugzx
    Lang:C++
    Blog:blog.csdn.net/qq_43756519
*/
#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define endl '\n'
#define SZ(x) (int)x.size()
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
const int mod=1e9+7;
//const int mod=998244353;
const double eps = 1e-10;
const double pi=acos(-1.0);
const int maxn=1e6+10;
const ll inf=0x3f3f3f3f;
const int dir[8][2]={{0,1},{1,0},{0,-1},{-1,0},{1,1},{1,-1},{-1,1},{-1,-1}};

int a[110][110];
int dp[110][110][20];
int s[110][3];
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int n,m,t;
    cin>>n>>m>>t;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            cin>>a[i][j];
            s[i][j]=s[i-1][j]+a[i][j];
        }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            for(int k=1;k<=t;k++){
                dp[i][j][k]=max(dp[i-1][j][k],dp[i][j-1][k]);
                for(int p=0;p<i;p++)
                    dp[i][j][k]=max(dp[i][j][k],dp[p][j][k-1]+s[i][1]-s[p][1]);
                for(int p=0;p<j;p++)
                    dp[i][j][k]=max(dp[i][j][k],dp[i][p][k-1]+s[j][2]-s[p][2]);
                if(i==j)
                    for(int p=0;p<i;p++)
                        dp[i][j][k]=max(dp[i][j][k],dp[p][p][k-1]+s[i][1]+s[j][2]-s[p][1]-s[p][2]);
            }
    cout<<dp[n][n][t];
    return 0;
}