关闭

bzoj[SCOI2005]最大子矩阵

标签: bzoj动态规划
217人阅读 评论(0) 收藏 举报
分类:

这题m只有2,所以一开始就想到了分情况。m=1的情况很容易想,就是普通动规枚举所有情况输出。m=2的情况一开始本来想先求出第一列,再得出第二列,这样只能过6个点,因为无法判断到两列应该分开当作两个矩阵,还是合成为一个矩阵(你可能会说,只要当一个矩阵就好,但在枚举中,有如i=1(行),k=2(矩阵数)的情况,这时必须分成两个)所以只能两个同步求,如果两行枚举到一样的就判断一下合成一行。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int a[110][5],f[110][110][20],sum[110][5];
int main()
{
    int n,m,K;
    scanf("%d%d%d",&n,&m,&K);
    memset(sum,0,sizeof(sum));
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            sum[i][j]=sum[i-1][j]+a[i][j];
        }
    memset(f,0,sizeof(f));
    if(m==1)
    {
        f[1][1][1]=a[1][1];
        for(int i=2;i<=n;i++)
            for(int k=1;k<=K;k++)
            {
                f[i][1][k]=f[i-1][1][k];
                for(int j=i-1;j>=k-1;j--)
                    f[i][1][k]=max(f[i][1][k],(sum[i][1]-sum[j][1])+f[j][1][k-1]);
            }
        printf("%d\n",f[n][1][K]);
    }
    else
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                for(int k=1;k<=K;k++)
                {
                    f[i][j][k]=max(f[i-1][j][k],f[i][j-1][k]);
                    for(int ii=i-1;ii>=k-1;ii--)f[i][j][k]=max(f[i][j][k],(sum[i][1]-sum[ii][1])+f[ii][j][k-1]);
                    for(int jj=j-1;jj>=k-1;jj--)f[i][j][k]=max(f[i][j][k],(sum[j][2]-sum[jj][2])+f[i][jj][k-1]);
                    if(i==j)
                    {
                        for(int t=i-1;t>=k-1;t--)
                            f[i][j][k]=max(f[i][j][k],(sum[i][1]-sum[t][1])+(sum[j][2]-sum[t][2])+f[t][t][k-1]);
                    }
                }
        printf("%d\n",f[n][n][K]);
    }
    return 0;
}
1
3

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:6147次
    • 积分:438
    • 等级:
    • 排名:千里之外
    • 原创:35篇
    • 转载:1篇
    • 译文:0篇
    • 评论:13条
    膜拜大(蒟)牛(蒻)
    game
    最新评论