电子人的基因 Cyborg Genes

在这里插入图片描述
UVA10723
有两个字符串 S 1 , S 2 S_1,S_2 S1,S2,求能包含他们的最短串。显然 S 1 + S 2 S_1+S_2 S1+S2可以包含 S 1 S_1 S1,又包含 S 2 S_2 S2,但是却不一定是最短的。为什么呢?因为 S 1 , S 2 S_1,S_2 S1,S2中有公共的子序列,这些子序列实际上最后只需要出现一次,而非两次。在 S 1 + S 2 S_1+S_2 S1+S2中就把这些公共的子串被计算了两次。比如样例中的ABAAXGF和AABXFGA,最长公共子序列为ABXG,而 S 1 + S 2 S_1+S_2 S1+S2ABAAXGFAABXFGA,可以看到ABXG分别在 S 1 , S 2 S_1,S_2 S1,S2的部分中各出现了一次,而实际上我们只需要它出现一次即可。所以去掉一组ABXG得到AAFAABXFGA或ABAAXGFAFA都可以包含 S 1 , S 2 S_1,S_2 S1,S2,因此答案为 ∣ S 1 ∣ + ∣ S 2 ∣ − 最 长 公 共 子 序 列 长 度 |S_1|+|S_2|-最长公共子序列长度 S1+S2。因为得到的最短字符串为为 S 1 + S 2 − 最 长 公 共 子 序 列 S_1+S_2-最长公共子序列 S1+S2,而 S 1 , S 2 S_1,S_2 S1,S2是固定的,因此,解的数量为最长公共子序列的数量。
定义 d p 1 [ i ] [ j ] dp1[i][j] dp1[i][j]为考虑 S 1 S_1 S1 i i i个字符, S 2 S_2 S2 j j j个字符的最长公共子序列长度, d p 2 [ i ] [ j ] dp2[i][j] dp2[i][j]为最长公共子序列的数量。
初始化
d p 1 = 0 d p 2 [ i ] [ 0 ] = d p 2 [ 0 ] [ i ] = 1 , 即 此 时 最 长 公 共 子 序 列 唯 一 , 是 空 串 dp1=0\\dp2[i][0]=dp2[0][i]=1,即此时最长公共子序列唯一,是空串 dp1=0dp2[i][0]=dp2[0][i]=1
转移方程 d p 1 dp1 dp1常规 L C S LCS LCS转移方式, d p 2 dp2 dp2见代码
AC代码

#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
#include<map>
using namespace std;
string S1, S2;
int dp1[31][31], dp2[31][31];
int S1_Len, S2_Len;
void DP() {
	memset(dp1, 0x0, sizeof(dp1));
	memset(dp2, 0x0, sizeof(dp2));
	for (int i = 0; i <= S1_Len; ++i) {
		dp2[i][0] = 1;
	}
	for (int i = 0; i <= S2_Len; ++i) {
		dp2[0][i] = 1;
	}
	//开始转移
	for (int i = 1; i <= S1_Len; ++i) {
		for (int j = 1; j <= S2_Len; ++j) {
			//如果两个字符相等
			if (S1[i] == S2[j]) {
				dp1[i][j] = dp1[i - 1][j - 1] + 1;
				//新方案只有一种转移方式
				dp2[i][j] = dp2[i - 1][j - 1];
			}
			else if (dp1[i - 1][j] < dp1[i][j - 1]) {
				dp1[i][j] = dp1[i][j - 1];
				//新方案只有一种转移方式
				dp2[i][j] = dp2[i][j - 1];
			}
			else if (dp1[i - 1][j] > dp1[i][j - 1]) {
				dp1[i][j] = dp1[i - 1][j];
				//新方案只有一种转移方式
				dp2[i][j] = dp2[i - 1][j];
			}
			else {
				dp1[i][j] = dp1[i - 1][j];
				//新方案有两种种转移方式,方案数为两种转移啊方式之和
				dp2[i][j] = dp2[i - 1][j] + dp2[i][j - 1];
			}
		}
	}
}
int main() {
	int T;
	scanf("%d", &T);
	getchar();
	for (int Case = 1; Case <= T; ++Case) {
		getline(cin, S1);
		getline(cin, S2);
		S1_Len = S1.size();
		S2_Len = S2.size();
		S1 = " " + S1;
		S2 = " " + S2;
		DP();
		printf("Case #%d: %d %d\n", Case, S1_Len + S2_Len - dp1[S1_Len][S2_Len], dp2[S1_Len][S2_Len]);
	}
	return 0;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值