51nod (dp)

输入

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

输出

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

输入示例

abcicba
abdkscab

输出示例

abca

这一题是51nod教程上面的题。 在lcs里面更多的是求最长公共子序列的长度。这里是要求输出最长公共子序列。 我们知道长度是一定的，但是最长公共子序列可能有多个。  在dp的状态转移方程我们只能求出长度， 该怎么根据状态转移方程输出序列呢？ 

设ax 和 by分别是两个字符串末尾的字符。  那么字符串A和B的最长公共子序列怎么求？  如果 ax==by，那么

dp[x][y]=dp[x-1][y-1]+1。  如果ax != by，那么 dp[x][y] = max(dp[i-1][j],dp[i][j-1])。   在当 ax=by时，那么ax被记录在当前最长公共子序列中，不相等时，我们查找的是 dp[x-1][y]与 dp[x][y-1] 哪个长度1更长。  根据这个思路我们记录下每次对字符的选择，最后从ax ，ay开始往前回溯，打印出最长公共子序列。

代码如下：

[cpp]  view plain  copy

 

1. #include<cstdio>  
2. #include<cstring>  
3. #include<algorithm>  
4. using namespace std;  
5. char a[1010],b[1010];  
6. int dp[1010][1010],mark[1010][1010],lena,lenb;  
7.   
8. void lcs()  
9. {  
10.     int i,j;  
11.     memset(dp,0,sizeof(dp));  
12.     for(i=1;i<=lena;++i)  
13.         mark[i][0]=1;  
14.     for(i=1;i<=lenb;++i)  
15.         mark[0][i]=-1;   
16.     for(i=1;i<=lena;++i)  
17.     {  
18.         for(j=1;j<=lenb;++j)  
19.         {  
20.             if(a[i-1]==b[j-1])  
21.             {  
22.                 dp[i][j]=dp[i-1][j-1]+1;  
23.                 mark[i][j]=0;  
24.             }  
25.             else if(dp[i-1][j]>=dp[i][j-1])  
26.             {  
27.                 dp[i][j]=dp[i-1][j];  
28.                 mark[i][j]=1;  
29.             }  
30.             else  
31.             {  
32.                 dp[i][j]=dp[i][j-1];  
33.                 mark[i][j]=-1;  
34.             }  
35.         }  
36.     }  
37. }  
38.   
39. void output(int x,int y)  
40. {  
41.     if(!x&&!y)  
42.         return ;  
43.     if(mark[x][y]==0)  
44.     {  
45.         output(x-1,y-1);  
46.         printf("%c",a[x-1]);  
47.     }  
48.     else if(mark[x][y]==1)  
49.         output(x-1,y);  
50.     else  
51.         output(x,y-1);  
52. }  
53.   
54. int main()  
55. {  
56.     while(scanf("%s%s",a,b)!=EOF)  
57.     {  
58.         lena=strlen(a);  
59.         lenb=strlen(b);  
60.         lcs();  
61.         output(lena,lenb);  
62.         printf("\n");  
63.     }  
64.     return 0;  
65. }