「Luogu4158」[SCOI2009]粉刷匠

「Luogu4158」[SCOI2009]粉刷匠

problem

题目描述

\(windy\)有 \(N\) 条木板需要被粉刷。 每条木板被分为 \(M\) 个格子。 每个格子要被刷成红色或蓝色。

\(windy\)每次粉刷，只能选择一条木板上一段连续的格子，然后涂上一种颜色。 每个格子最多只能被粉刷一次。

如果\(windy\)只能粉刷 \(T\) 次，他最多能正确粉刷多少格子？

一个格子如果未被粉刷或者被粉刷错颜色，就算错误粉刷。

输入输出格式

输入格式：

第一行包含三个整数，\(N\) \(M\) \(T\)。

接下来有\(N\)行，每行一个长度为\(M\)的字符串，'0'表示红色，'1'表示蓝色。

输出格式：

包含一个整数，最多能正确粉刷的格子数。

输入输出样例

输入样例#1：

3 6 3
111111
000000
001100

输出样例#1：

16

说明

\(30\%\)的数据，满足 \(1 \le N,M \le 10\) ； \(0 \le T \le 100\) 。

\(100\%\)的数据，满足 \(1 \le N,M <= 50\) ； \(0 \le T \le 2500\) 。

Solution

蒟蒻看到这题想了\(n\)多种完全不正确的处理方法，如果是在考场上估计已经光速凉凉了\(QAQ\)

对于每一行，我们可以把它分成若干个颜色不同的连续段（对应若干次颜色不同的粉刷），从左到右考虑

注意到题目：

一个格子如果未被粉刷或者被粉刷错颜色，就算错误粉刷。

意思就是不刷白不刷

如果是颜色不对，多刷这一格不会比不刷差，就不需要考虑不刷的情况

所以我们有：

状态：设\(dp_{i,j,k,0/1}\)表示当前枚举到第\(i\)行，第\(j\)个格子，已经涂了\(k\)次（分成了\(k\)段），当前格子涂的颜色是\(0/1\)

方程：

\[dp_{i,j,k,x}=max(dp_{i,j-1,k,x},dp_{i,j-1,k-1,x\space xor\space 1})+[x=color_{i,j}] \]

(方括号是艾弗森括号，当其中的条件为真时值为\(1\)，否则为\(0\))

这个转移应该很好理解吧

---

处理完每一行之后，\(max(dp_{i,m,k,0},dp_{i,m,k,1})\)即为第\(i\)行分配\(k\)次粉刷次数所能正确粉刷的最大格子数

这样，每一行作为一组，原问题就转化为一个分组背包问题，在此不再赘述

Code

代码中所用的字母与上述转移方程略有不同，请注意辨别

#include <cstdio>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#define maxn 55
#define maxm 55
#define maxt 2505
using namespace std;
typedef long long ll;

template <typename T> void read(T &t)
{
	t=0;int f=0;char c=getchar();
	while(!isdigit(c)){f|=c=='-';c=getchar();}
	while(isdigit(c)){t=t*10+c-'0';c=getchar();}
	if(f)t=-t;
}

int n,m,t;
int col[maxn][maxm];
int dp[maxn][maxm][maxm][2],tdp[maxt];

int main()
{
	read(n),read(m),read(t);
	for(register int i=1;i<=n;++i)
	{
		char c[maxm];
		scanf("%s",c+1);
		for(register int j=1;j<=m;++j)
			col[i][j]=c[j]-'0';
	}
	for(register int li=1;li<=n;++li)
		for(register int i=1;i<=m;++i)
			for(register int k=1;k<=min(t,m);++k)
				for(register int x=0;x<=1;++x)
					dp[li][i][k][x]=max(dp[li][i-1][k][x],dp[li][i-1][k-1][x^1])+(x==col[li][i]);
	for(register int i=1;i<=n;++i)
		for(register int j=t;j>=0;--j)
			for(register int k=1;k<=min(j,m);++k)
				tdp[j]=max(tdp[j],tdp[j-k]+max(dp[i][m][k][0],dp[i][m][k][1]));
	printf("%d",tdp[t]);
	return 0;
}