[SCOI2005] 最大子矩阵

本题地址

http://www.luogu.org/problem/show?pid=2331

题目描述

这里有一个n*m的矩阵，请你选出其中k个子矩阵，使得这个k个子矩阵分值之和最大。注意：选出的k个子矩阵不能相互重叠。
输入输出格式
输入格式：
第一行为n,m,k（1≤n≤100,1≤m≤2,1≤k≤10），接下来n行描述矩阵每行中的每个元素的分值(每个元素的分值的绝对值不超过32767)。

输出格式：
只有一行为k个子矩阵分值之和最大为多少。

输入输出样例
输入样例#1：
3 2 2
1 -3
2 3
-2 3
输出样例#1：
9

题解：

->这道省选题真的是= =语言简洁啊 一开始并没有读懂分值的意思（自动联想到2014普及T4关于矩阵分值的定义）
-> ->其实就是一个值嘛QAQ加个分字干嘛 然后并没有好好读题 m是1或者2 这个数据我只有%%% 然后就是面向数据乱搞嘛= =
-> ->m==1很好搞嘛 和上面那个T4一样的
-> ->dp(i,j)表示i个里选j个矩阵最大值 先dp(i,j)=dp(i-1,j)然后就dp(i,j)=max(dp(i,j),dp(k,j-1)+sum(i)-sum(k)) sum是前缀和
-> ->好像这样就可以到手30分了
-> ->m==2的时候嘛 发现对于每一行并没有什么影响 所以对于每行像上面一样DP
-> ->dp(i1,i2,j)表示上一列i1和下一列i2之前选j个矩阵
-> ->多出来的东西就是i1==i2的时候可以多出来一种同时选两列的方法（由于更新值的时候都是必定会选i1或者i2的，所以只存在i1==i2的时候能同时选两列，仔细想想上面的方程或者看代码）
-> ->对于多出来的情况，只要把两列当做一列看就好了（理解之后觉得作为DP专题最后一道有些水啊。。）
-> ->代码：（奇怪的宏只是为了代码看起来舒服一些，毕竟每次打dp要多打一个1或者2的感觉还是有点不爽）

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=100+10;
const int maxm=10+5;

int a[2][maxn];
int dp1[maxn][maxm],dp2[2][maxn][maxm];
int sum1[maxn],sum2[2][maxn];

int main()
{
    int n,m,x;cin>>n>>m>>x;
    if(m==1)
    {
        #define dp dp1
        #define sum sum1
        for(int i=1;i<=n;i++) 
        {
            scanf("%d",&a[1][i]);
            sum[i]=sum[i-1]+a[1][i];
        }
        dp[0][0]=0;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=i&&j<=x;j++)
            {
                dp[i][j]=dp[i-1][j];
                for(int k=j-1;k<=i;k++)
                {
                    dp[i][j]=max(dp[i][j],dp[k][j-1]+sum[i]-sum[k]);
                }
            }         
        printf("%d",dp[n][x]);
        #undef dp
        #undef sum
    }  
    else
    {
        #define dp dp2
        #define sum sum2
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d",&a[0][i],&a[1][i]);
            sum[0][i]=sum[0][i-1]+a[0][i];
            sum[1][i]=sum[1][i-1]+a[1][i];
        }  
        for(int i1=1;i1<=n;i1++)
            for(int i2=1;i2<=n;i2++)
            for(int j=1;j<=i1+i2&&j<=x;j++)
            {
                dp[i1][i2][j]=max(dp[i1-1][i2][j],dp[i1][i2-1][j]);
                for(int l=0;l<i1;l++)
                {
                    dp[i1][i2][j]=max(dp[i1][i2][j],dp[l][i2][j-1]+sum[0][i1]-sum[0][l]);
                }
                for(int l=0;l<i2;l++)
                {
                    dp[i1][i2][j]=max(dp[i1][i2][j],dp[i1][l][j-1]+sum[1][i2]-sum[1][l]);
                }
                if(i1==i2)
                {
                    for(int l=0;l< i1;l++)
                        dp[i1][i2][j]=max(dp[i1][i2][j],dp[l][l][j-1]+sum[0][i1]-sum[0][l]+sum[1][i2]-sum[1][l]);
                }
            } 
        printf("%d",dp[n][n][x]);
        #undef dp
        #undef sum
    }
    return 0;
}    

DP专题终于搞完啦= = 等什么时候变强了去把day2T2狠狠地干掉