BZOJ[1084]最大子矩阵

1084: [SCOI2005]最大子矩阵
Time Limit: 10 Sec Memory Limit: 162 MB
Description

这里有一个n*m的矩阵，请你选出其中k个子矩阵，使得这个k个子矩阵分值之和最大。注意：选出的k个子矩阵
不能相互重叠。
Input

第一行为n,m,k（1≤n≤100,1≤m≤2,1≤k≤10），接下来n行描述矩阵每行中的每个元素的分值(每个元素的
分值的绝对值不超过32767)。

Output

只有一行为k个子矩阵分值之和最大为多少。

Sample Input

3 2 2
1 -3
2 3
-2 3

Sample Output

9

m只有2，其实不是很大。。
这就比较简单了。。只用考虑两种情况。。
状态：f[i][j][k]表示第一列到了第i个格子，第二列到了第j个格子，已经选取了k个矩阵的最大得数。
先把现有的状态向后转移，转移成f[i·][j·][k]，然后取矩阵，向后转移，转移成f[i·][j·][k + 1]。这里要分三种情况讨论，转移第一列，转移第二列，和两列一起转移。

这样的话，就可以算出最大的值啦！！

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<set>
#include<map>
#include<sstream>
#include<iostream>
#include<algorithm>
using namespace std;
int a[1100][1100],s[1100][1100];
int f1[1100][1100],f2[110][110][110];
int mymax(int x,int y)
{
    return x>y?x:y;
}
int main()
{
    int n,m,t;
    scanf("%d%d%d",&n,&m,&t);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            s[j][i]=s[j][i-1]+a[i][j];
        }
    }
    int ans=0;
    if(n==1)
    {
        memset(f1,0,sizeof(f1));
        for(int i=0;i<=n;i++)
        {
            for(int k=0;k<=t;k++)
            {
                for(int j=0;j<=n;j++)
                {
                    f1[j][k]=mymax(f1[j][k],f1[i][k]);
                    f1[j][k+1]=mymax(f1[j][k+1],f1[i][k]+s[1][j]-s[1][i]); 
                }
            }
        }
        for(int i=0;i<=n;i++)
        {
            for(int j=0;j<=t;j++)
            {
                ans=mymax(ans,f1[i][j]);
            }
        }
        printf("%d\n",ans);
        return 0;
    }
    else
    {
        memset(f2,0,sizeof(f2));
        for(int i=0;i<=n;i++)  
        {
            for(int j=0;j<=n;j++)
            { 
                for(int k=0;k<=t;k++)
                {  
                    for(int kk=i+1;kk<=n;kk++)
                    {  
                        f2[kk][j][k]=mymax(f2[kk][j][k],f2[i][j][k]);  
                        f2[kk][j][k+1]=mymax(f2[kk][j][k+1],f2[i][j][k]+s[1][kk]-s[1][i]);  
                    }  
                    for(int kk=j+1;kk<=n;kk++)
                    {  
                        f2[i][kk][k]=mymax(f2[i][kk][k],f2[i][j][k]);  
                        f2[i][kk][k+1]=mymax(f2[i][kk][k+1],f2[i][j][k]+s[2][kk]-s[2][j]);  
                    }  
                    for(int kk=mymax(i,j)+1;kk<=n;kk++)
                    {  
                        f2[kk][kk][k]=mymax(f2[kk][kk][k],f2[i][j][k]);  
                        f2[kk][kk][k+1]=mymax(f2[kk][kk][k+1],f2[i][j][k]+s[1][kk]-s[1][mymax(i,j)]+s[2][kk]-s[2][mymax(i,j)]);  
                    }  
                }
            }
        }
        for(int i=0;i<=n;i++)
        {
            for(int j=0;j<=n;j++)
            {
                for(int k=0;k<=t;k++)
                {
                    ans=mymax(ans,f2[i][j][k]);
                }
            }
        }
        printf("%d\n",ans);
        return 0;
    }
    return 0;
}

开学了。。马上就要参加NOIP的考试了。。
希望大家都能获得好成绩！

版权声明：本文为博主原创文章，未经博主允许不得转载。