[BZOJ1084][SCOI2005]最大子矩阵

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
char buf[10000000], *ptr = buf - 1;
inline int readint(){
    int f = 1, n = 0;
    char ch = *++ptr;
    while (ch < ‘0‘ || ch > ‘9‘){
        if(ch == ‘-‘) f = -1;
        ch = *++ptr;
    }
    while (ch <= ‘9‘ && ch >= ‘0‘){
        n = (n << 1) + (n << 3) + ch - ‘0‘;
        ch = *++ptr;
    }
    return f * n;
}
const int maxn = 100 + 10, INF = 0x7f7f7f7f;
int N, M, K;
int s[maxn], sum[maxn][3];
int f[maxn][11], g[maxn][maxn][11];
void work1(){
    s[0] = 0;
    for(int i = 1; i <= N; i++) s[i] = s[i - 1] + readint();
    for(int i = 0; i <= N; i++)
        for(int k = 1; k <= K; k++)
            f[i][k] = -INF;
    for(int i = 1; i <= N; i++)
        for(int k = 1; k <= K; k++){
            f[i][k] = f[i - 1][k];
            for(int j = 0; j < i; j++)
                f[i][k] = max(f[i][k], f[j][k - 1] + s[i] - s[j]);
        }
    printf("%d\n", f[N][K]);
}
void work2(){
    sum[0][1] = sum[0][2] = 0;
    for(int i = 1; i <= N; i++){
        sum[i][1] = sum[i - 1][1] + readint();
        sum[i][2] = sum[i - 1][2] + readint();
    }
    for(int i = 0; i <= N; i++)
        for(int j = 0; j <= N; j++)
            for(int k = 1; k <= K; k++)
                g[i][j][k] = -INF;    
    for(int i = 1; i <= N; i++)
        for(int j = 1; j <= N; j++)
            for(int k = 1; k <= K; k++){
                g[i][j][k] = max(g[i - 1][j][k], g[i][j - 1][k]);
                for(int L = 0; L < i; L++)
                    g[i][j][k] = max(g[i][j][k], g[L][j][k - 1] + sum[i][1] - sum[L][1]);
                for(int L = 0; L < j; L++)
                    g[i][j][k] = max(g[i][j][k], g[i][L][k - 1] + sum[j][2] - sum[L][2]);
                if(i == j)
                    for(int L = 0; L < i; L++)
                        g[i][j][k] = max(g[i][j][k], g[L][L][k - 1] + sum[i][1] + sum[i][2] - sum[L][1] - sum[L][2]);
            }
    printf("%d\n", g[N][N][K]);
}
int main(){
    fread(buf, sizeof(char), sizeof(buf), stdin);
    N = readint();
    M = readint();
    K = readint();
    if(M == 1) work1();
    else work2();
    return 0;
}