【题解】ZJOI2013蚂蚁寻路

　　这题强呀……打了10+30暴力之后苦想1h并不会做……于是去看题解。看题解的时候又莫名各种看错，结果看了好久才懂……记录一下血泪史吧。

　　这题不难发现走出来的图形就是一个高低高低的城堡型图案，命名为高峰跟低谷的话就是一共有k个低谷和k + 1个高峰，且交替出现。发现其实这个图形是由2 * k + 1个矩形所构成的，我们考虑将这样许多矩形看做转移的方向：f[i][j][p][h]代表以点（i, j）为左下角，当前枚举到的是第p个矩形，且高度为（i - h + 1）所能获得的最大权值。

　　转移方程：

　　　　\(f[i][j][p][h] = \left \{ f[i][j - 1][p][h], f[i][j - 1][p][h']] \right \} + s[j][i] - s[j][h - 1]; \)

　　其中第一部分代表和j - 1处于同一个矩形当中，第二部分则表示根据p的奇偶性选择<h或>h的转移。第二部分如果枚举，则复杂度太高了难以承受，所以用另一个数组g来保存，追加一维0/1来表示是>还是<中的最大值。然后发现其实 i 这一维是没有参与转移的，所以将这一维去掉，节约空间。期间有一个让我比较纠结的地方，就是

　　　　\( ans = max(f[i][j][k][i], g[i][j][k][i][0]); \)

　　起初并没有很理解为什么要将前一部分考虑进来。当高度等于i,不是说明这是一个低谷吗？为什么会符合条件呢？但其实有一种情况下是可以的：当没有低谷的时候，这样是一个符合条件的解。追加辅助数组的dp做的太少了，要多多加油呀ヾ(๑╹◡╹)ﾉ"

#include <bits/stdc++.h>
using namespace std;
#define INF 999999999
#define maxn 105
int n, m, k, ans = -INF, sum[maxn][maxn];
int f[maxn][maxn][maxn], g[maxn][maxn][maxn][];

int read()
{
    int x = , k = ;
    char c;
    c = getchar();
    while(c < '' || c > '') { if(c == '-') k = -; c = getchar(); }
    while(c >= '' && c <= '') x = x *  + c - '', c = getchar();
    return x * k;
}

void work()
{
    for(int p = ; p <= k; p ++)
        for(int h = ; h <= n; h ++)
        {
            f[][p][h] = -INF;
            g[][p][h][] = g[][p][h][] = -INF;
        }

    for(int i = ; i <= n; i ++)
        for(int j = ; j <= m; j ++)
        {
            for(int p = ; p <= k; p ++)
            {
                for(int h = ; h <= i; h ++)
                {
                    f[j][p][h] = max(f[j - ][p][h], g[j - ][p - ][h][p % ]);
                    f[j][p][h] += sum[j][i] - sum[j][h - ];
                }
                g[j][p][][] = -INF;
                for(int h = ; h <= i; h ++)
                    g[j][p][h][] = max(g[j][p][h - ][], f[j][p][h - ]);
                g[j][p][i][] = -INF;
                for(int h = i - ; h >= ; h --)
                    g[j][p][h][] = max(g[j][p][h + ][], f[j][p][h + ]);
            }
            ans = max(ans, max(f[j][k][i], g[j][k][i][]));
        }
}

int main()
{
    n = read(), m = read(), k = read();
    k = k *  + ;
    for(int i = ; i <= n; i ++)
        for(int j = ; j <= m; j ++)
        {
            int t = read();
            sum[j][i] += sum[j][i - ] + t;
        }
    work();
    printf("%d\n", ans);
    return ;
}