最长公共子序列Lcs （51Nod - 1006）

 1 #include<cstdio>
 2 #include<cstdlib>
 3 #include<iostream>
 4 #include<cstring>
 5 #include<string>
 6 #include<cmath>
 7 #include<algorithm>
 8 using namespace std;
 9 int len1,len2;
10 char a[1001],b[1001];
11 int c[1001][1001];
12 int main()  
13 {
14     cin>>a>>b;
15     len1=strlen(a);
16     len2=strlen(b);
17     memset(c,0,sizeof(c));
18     int i,j;
19     for(i=1;i<=len1;i++)
20         for(j=1;j<=len2;j++)
21         {
22             if(a[i-1]==b[j-1]) c[i][j]=c[i-1][j-1]+1;
23             else c[i][j]=max(c[i-1][j],c[i][j-1]);
24         }
25     int x=len1;
26     int y=len2;
27     int k=c[x][y];//最长序列出现在右下角.（如果把算出来的c[i][j]画成表格的话）
28     char lcs[1001] = {‘\0‘};
29     while(i&&j)
30     {
31         if(a[i-1]==b[j-1]&&c[i][j]==c[i-1][j-1]+1)//此时i,j已经到了len1和len2。所以这里是倒着回去看的。
32         {//c[i][j]==c[i-1][j-1]+1证明是从左上角下来的也就是说,左上角位置的a[i-1]和b[j-1]末尾相同。
33             lcs[--k] = a[i-1];//先将k--再用是因为存数据的时候从0开始.
34             i--;
35             j--;
36         }
37         else if(a[i-1]!=b[j-1]&&c[i-1][j]>c[i][j-1])  i--;//看上方图就能明白。
38         else j--;
39     }
40     printf("%s\n",lcs);
41     return 0;  
42 }