洛谷 P4059 [Code+#1]找爸爸

最新推荐文章于 2024-08-16 22:18:58 发布

北神不会咕

最新推荐文章于 2024-08-16 22:18:58 发布

阅读量146

点赞数

分类专栏：洛谷文章标签：算法 c++ 动态规划

本文链接：https://blog.csdn.net/KINGtxy520/article/details/126694219

版权

博客介绍了洛谷P4059题目的动态规划解决方案。面对n和m在可接受范围内的情况，博主通过开辟二维数组f[i][j]表示两个字符串分别考虑至第i位和第j位的状态。由于空格处理的复杂性，博主提出将空格的贡献转换为-B或-A，并采用三维数组f[i][j][0/1][0/1]来表示是否放置空格。虽然定义有些冗余，但对时间复杂度影响不大，转移状态时考虑放字符或空格进行更新。

摘要由CSDN通过智能技术生成

应该算是比较基础的一个DP，看到 n, m 的范围都在能开的下数组的范围内，果断开一个 f[i][j] 的数组表示一个串考虑到 i 位，另一个串考虑到 j 位，但是发现空格没有办法处理，如果加上末尾连续空格的长度则又要加上一维，时间复杂度有点。。。

所以应该将空格的长度转化为贡献，即如果前边有一个空格，那么空格的贡献是 -B 否则为 -A ，所以将定义改为 f[i][j][0/1][0/1] ，后边两维表示放没放空格，转移的时候考虑放字符还是放空格然后分别转移即可，写方程的时候才发现这个定义其实有些冗余，因为两个串的某个位置不可能都同时放空格，但是对时间复杂度影响不大，所以就没改。

#include <bits/stdc++.h>
#define N 3010
using namespace std;

int d[5][5], dp[N][N][2][2];
char s1[N], s2[N];
int a[N], b[N];

int main() {
	memset(dp, -0x3f, sizeof(dp));
	scanf("%s%s", s1 + 1, s2 + 1);
	int n = strlen(s1 + 1), m = strlen(s2 + 1);
	for(int i = 1; i <= n; i ++) {
		if(s1[i] == 'A') a[i] = 1;
		else if(s1[i] == 'T') a[i] = 2;
		else if(s1[i] == 'G') a[i] = 3;
		else if(s1[i] == 'C') a