【动态规划】【SCOI2009】粉刷匠

【问题描述】
windy有 N 条木板需要被粉刷。
每条木板被分为 M 个格子。
每个格子要被刷成红色或蓝色。
windy每次粉刷，只能选择一条木板上一段连续的格子，然后涂上一种颜色。
每个格子最多只能被粉刷一次。
如果windy只能粉刷 T 次，他最多能正确粉刷多少格子？
一个格子如果未被粉刷或者被粉刷错颜色，就算错误粉刷。

【输入格式】
输入文件paint.in第一行包含三个整数，N M T。
接下来有N行，每行一个长度为M的字符串，'0'表示红色，'1'表示蓝色。

【输出格式】
输出文件paint.out包含一个整数，最多能正确粉刷的格子数。

【输入样例】
3 6 3
111111
000000
001100

【输出样例】
16

【数据规模和约定】
30%的数据，满足 1 <= N,M <= 10 ； 0 <= T <= 100 。
100%的数据，满足 1 <= N,M <= 50 ； 0 <= T <= 2500 。

此题考察动态规划。
状态：用f[t][i][j][k]表示t次涂色过后，涂到了第i个木板的第j格，且被涂为颜色k（由于不涂色等价于颜色涂错所以不用考虑不涂色的情况）的最大正确格数。

转移方程：f[t][i][j][k] = max{f[t-1][i-1][m][0..1]} + cost (j=1的情况)；
= max(f[t][i][j-1][k],f[t-1][i][j-1][1-k]) + cost (j>1的情况)。

提交情况：
提交5次。
第一次10分；
第二次10分；
第三次0分（刷的次数应当放在最外层，否则有后效性）；
第四次80分（更新ans出错）；
第五次100分。
Accode：

#include <cstdio>
#include <algorithm>

using std::max;
const int maxN = 60;
const int maxT = 2510;

int f[maxT][maxN][maxN][2];
bool mp[maxN][maxN];
int n, m, T, ans;

int main()
{
	freopen("paint.in", "r", stdin);
	freopen("paint.out", "w", stdout);
    scanf("%d%d%d", &n, &m, &T);
    for (int i = 1; i < n + 1; ++i)
    {
        getchar();
        for (int j = 1; j < m + 1; ++j)
            mp[i][j] = getchar() - '0';
    }
#define Rep(i, L, R) \
for (int i = (L); i < (R); ++i)
    Rep(t, 1, T + 1)
    Rep(i, 1, n + 1)
    Rep(j, 1, m + 1)
    Rep(k, 0, 2)
#undef Rep(i, L, R)
    {
        if (j == 1) f[t][i][j][k] = (k != mp[i][j])
                    + max(f[t - 1][i - 1][m][0],
                          f[t - 1][i - 1][m][1]);
        else f[t][i][j][k] = (k != mp[i][j])
             + max(f[t][i][j - 1][k],
                   f[t - 1][i][j - 1][!k]);
        if (f[t][i][j][k] > ans)
            ans = f[t][i][j][k];
	//ans必须取所有的ans中最大值。
    }
    printf("%d\n", ans);
	return 0;
}