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+S2为ABAAXGFAABXFGA,可以看到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;
}