这题强呀……打了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 ;
}