P4059 [Code+#1]找爸爸（线性dp）

题目链接：点击这里

题目大意：
给出两个字符串，你需要在其中插入任意个空格，保证两个串长度相同。如果两个串的某同一位置都是字母则获得某给定价值，对于每个串的每个长度为 $k$ 的连续空格段要付出 $a+b(k-1)$ 的代价（ $a,b>0$ ），求可获得的最大价值

题目分析：
考虑设计状态，因为是两个串的匹配问题，状态肯定要有其两个串分别匹配到了什么地方，还要加一维的限制，本题的限制就是末尾是否有空格，故状态为 $dp[i][j][k]$ 表示两个串分别匹配到了 $i,j$ ，$k$ 代表是否末尾的空格位置，因为两个串都是空格的话显然会有额外花费，不会是最优状态；
所以 $k$ 要表示三种情况：两个串都没有空格，第一个串末尾是空格，第二个串末尾是空格
考虑状态的转移：
如果 $k=0$， $dp[i][j][0]$ 显然就是从 $dp[i-1][j-1][]$ 的最大值转移过来
$dp[i][j][1/2]$ 的转移方程如下：
$$dp[i][j][1]=max(dp[i][j-1][0]-a,dp[i][j-1][1]-b,dp[i][j-1][2]-a)$$
$$dp[i][j][2]=max(dp[i-1][j][0]-a,dp[i-1][j][1]-a,dp[i-1][j][2]-b)$$
解释一下方程：
$dp[i][j][1]$ 表示第一个串后面有空格，$dp[i][j-1][0]$ 表示没有空格所以 $-a$ 即可， $dp[i][j-1][1]$ 表示第一个串后面有空格了所以 $-b$ 即可， $dp[i][j-1][2]$ 表示第二个串后面有空格转移过来就得 $-a$
$dp[i][j][2]$ 转移方程含义同理
最后答案就是 $dp[n][m][]$ 的最大值

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 3e3+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
int n,m,a,b,x[maxn],y[maxn],d[4][4],dp[maxn][maxn][3];
char s[maxn];
map<char,int>mp;
int main()
{
	mp['A'] = 0;
	mp['T'] = 1;
	mp['G'] = 2;
	mp['C'] = 3;
	scanf("%s",s+1);
	n = strlen(s+1);
	for(int i = 1;i <= n;i++)
		x[i] = mp[s[i]];
	scanf("%s",s+1);
	m = strlen(s+1);
	for(int i = 1;i <= m;i++)
		y[i] = mp[s[i]];
	for(int i = 0;i < 4;i++)
		for(int j = 0;j < 4;j++)
			d[i][j] = read();
	a = read(),b = read();
	for(int i = max(n,m);i;i--)
	{
		dp[0][i][0] = dp[i][0][0] = dp[0][i][2] = dp[i][0][1] =-inf;
		dp[0][i][1] = dp[i][0][2] = -a-b*(i-1);
	}
	dp[0][0][1] = dp[0][0][2] = -inf;
	for(int i = 1;i <= n;i++)
		for(int j = 1;j <= m;j++)
		{
			dp[i][j][0] = *max_element(dp[i-1][j-1],dp[i-1][j-1]+3)+d[x[i]][y[j]];
			dp[i][j][1] = max({dp[i][j-1][0]-a,dp[i][j-1][1]-b,dp[i][j-1][2]-a});
			dp[i][j][2] = max({dp[i-1][j][0]-a,dp[i-1][j][1]-a,dp[i-1][j][2]-b});
		}
	printf("%d\n",*max_element(dp[n][m],dp[n][m]+3));
	return 0;
}