51nod 1006 最长公共子序列Lcs（经典动态规划）

传送门

Description

给出两个字符串A B，求A与B的最长公共子序列（子序列不要求是连续的）。

 

比如两个串为：

 

abcicba

abdkscab

 

ab是两个串的子序列，abc也是，abca也是，其中abca是这两个字符串最长的子序列。

Input

第1行：字符串A
第2行：字符串B
(A,B的长度 <= 1000)

Output

输出最长的子序列，如果有多个，随意输出1个。

Sample Input

abcicba
abdkscab

Sample Output

abca

思路　

记：
Xi=﹤x1，⋯，xi﹥即X序列的前i个字符 (1≤i≤m)（前缀）

Yj=﹤y1，⋯，yj﹥即Y序列的前j个字符 (1≤j≤n)（前缀）

假定Z=﹤z1，⋯，zk﹥∈LCS(X , Y) 。

• 若xm=yn（最后一个字符相同），则不难用反证法证明：该字符必是X与Y的任一最长公共子序列Z（设长度为k）的最后一个字符，即有zk = xm = yn 且显然有Zk-1∈LCS(Xm-1 , Yn-1)即Z的前缀Zk-1是Xm-1与Yn-1的最长公共子序列。此时，问题化归成求Xm-1与Yn-1的LCS（LCS(X , Y)的长度等于LCS(Xm-1 , Yn-1)的长度加1）。

• 若xm≠yn，则亦不难用反证法证明：要么Z∈LCS(Xm-1, Y)，要么Z∈LCS(X , Yn-1)。由于zk≠xm与zk≠yn其中至少有一个必成立，若zk≠xm则有Z∈LCS(Xm-1 , Y)，类似的，若zk≠yn 则有Z∈LCS(X , Yn-1)。此时，问题化归成求Xm-1与Y的LCS及X与Yn-1的LCS。LCS(X , Y)的长度为：max{LCS(Xm-1 , Y)的长度, LCS(X , Yn-1)的长度}。

由于上述当xm≠yn的情况中，求LCS(Xm-1 , Y)的长度与LCS(X , Yn-1)的长度，这两个问题不是相互独立的：两者都需要求LCS(Xm-1，Yn-1)的长度。另外两个序列的LCS中包含了两个序列的前缀的LCS，故问题具有最优子结构性质考虑用动态规划法。

也就是说，解决这个LCS问题，你要求三个方面的东西：

• • 1、LCS(Xm-1，Yn-1)+1；
  • 2、LCS(Xm-1，Y)，LCS(X，Yn-1)；
  • 3、max{ LCS(Xm-1, Y)，LCS(X, Yn-1) }。

最长公共子序列的结构

最长公共子序列的结构有如下表示：

设序列X=< x1, x2, …, xm >和Y=< y1, y2, …, yn >的一个最长公共子序列Z=< z1, z2, …, zk >，则：

• 若xm=yn，则zk=xm=yn且Zk-1是Xm-1和Yn-1的最长公共子序列；
• 若xm≠yn且zk≠xm ，则Z是Xm-1和Y的最长公共子序列；
• 若xm≠yn且zk≠yn ，则Z是X和Yn-1的最长公共子序列。
  其中Xm-1 = < x1, x2, …, xm-1 >，Yn-1 = < y1, y2, …, yn-1 >，Zk-1 = < z1, z2, …, zk-1 >。

 我们定义c[i, j]表示X_i和Y_i的LCS的长度。如果i = 0或j = 0，即一个序列长度为0，那么LCS的长度为0.根据LCS问题的最优子结构性质，可得到如下公式：

		0	若i = 0 或 j = 0
c[i, j]	=	c[i - 1,j - 1] + 1	若i, j > 0 且 xi = yi
		max(c[i, j - 1],c[i - 1, j])	若i, j > 0 且 xi≠yi

 

+ View Code?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30

#include<stdio.h> #include<string.h> const int maxn = 1005; int dp[maxn][maxn] = {0}; int main() {     char a[maxn],b[maxn],lcs[maxn];     int i,j;     scanf("%s",a);scanf("%s",b);     int lena = strlen(a),lenb = strlen(b);     for (i = 1;i <= lena;i++)     {         for (j = 1;j <= lenb;j++)         {             if (a[i-1] == b[j-1])   dp[i][j] = dp[i-1][j-1]+1;             else    dp[i][j] = dp[i][j-1]>dp[i-1][j]?dp[i][j-1]:dp[i-1][j];         }     }     i = lena,j = lenb;     int len = dp[lena][lenb];     lcs[len] = '\0';     while (dp[i][j])     {         if (dp[i][j] == dp[i-1][j]) i--;         else if (dp[i][j] == dp[i][j-1])    j--;         else lcs[--len] = a[i-1],i--,j--;     }     printf("%s\n",lcs);     return 0; }