[SCOI2005][BZOJ 1084]最大子矩阵

#include<bits/stdc++.h>
using namespace std;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){x=(x<<3)+(x<<1)+ch-‘0‘;ch=getchar();}
    return x*f;
}
#define MN 105
#define inf 0x7f
int n,m,k,a[MN][3],ans;
int f[MN][15][6];
void rw(int &x,int y){if(y>x)x=y;}
int main(){
    n=read(),m=read(),k=read();
    register int i,j;
    for(i=1;i<=n;i++)for(j=1;j<=m;j++) a[i][j]=read();
    if(m==1){
        for(i=1;i<=n;i++) for(j=1;j<=k;j++){
            rw(f[i][j][0],f[i-1][j][0]);
            rw(f[i][j][0],f[i-1][j][1]);
            rw(f[i][j][1],f[i-1][j][1]+a[i][1]);
            rw(f[i][j][1],f[i-1][j-1][0]+a[i][1]);
        }
        printf("%d\n",max(f[n][k][1],f[n][k][0]));
    }
    else{
        memset(f,-inf,sizeof f);
        for(i=0;i<=n;i++)for(j=0;j<=k;j++) f[i][j][1]=0;
        for(i=1;i<=n;i++)for(j=1;j<=k;j++){
            for(int h=1;h<=5;h++) rw(f[i][j][1],f[i-1][j][h]);
            
            f[i][j][2]=max(f[i-1][j][2],f[i-1][j][5])+a[i][1];
            f[i][j][3]=max(f[i-1][j][3],f[i-1][j][5])+a[i][2];
            rw(f[i][j][2],f[i-1][j-1][4]+a[i][1]);
            rw(f[i][j][3],f[i-1][j-1][4]+a[i][2]);
            rw(f[i][j][2],max(f[i-1][j-1][1],f[i-1][j-1][3])+a[i][1]);
            rw(f[i][j][3],max(f[i-1][j-1][1],f[i-1][j-1][2])+a[i][2]);
            
            f[i][j][4]=max(f[i-1][j-1][1],f[i-1][j][4])+a[i][1]+a[i][2];
            rw(f[i][j][4],max(f[i-1][j-1][2],f[i-1][j-1][3])+a[i][1]+a[i][2]);
            rw(f[i][j][4],f[i-1][j-1][5]+a[i][1]+a[i][2]);
            
            f[i][j][5]=f[i-1][j][5]+a[i][1]+a[i][2];
            rw(f[i][j][5],max(f[i-1][j-1][2],f[i-1][j-1][3])+a[i][1]+a[i][2]);
            if(j>=2) rw(f[i][j][5],max(f[i-1][j-2][1],f[i-1][j-2][4])+a[i][1]+a[i][2]);
        }
        ans=-1e12;
        for(i=1;i<=5;i++) rw(ans,f[n][k][i]);
        printf("%d\n",ans);
    }
    return 0;
}