南阳理工：最长公共子序列

最长公共子序列

时间限制： 3000 ms  |  内存限制： 65535 KB

难度： 3

描述

咱们就不拐弯抹角了，如题，需要你做的就是写一个程序，得出最长公共子序列。
tip：最长公共子序列也称作最长公共子串(不要求连续)，英文缩写为LCS（Longest Common Subsequence）。其定义是，一个序列 S ，如果分别是两个或多个已知序列的子序列，且是所有符合此条件序列中最长的，则 S 称为已知序列的最长公共子序列。

输入

第一行给出一个整数N(0<N<100)表示待测数据组数
接下来每组数据两行，分别为待测的两组字符串。每个字符串长度不大于1000.

输出

每组测试数据输出一个整数，表示最长公共子序列长度。每组结果占一行。

样例输入

2  asdf  adfsd  123abc  abc123abc

样例输出

3  6

/*动态规划题目，分析如下

设序列X={x1,x2,x3......xm}和Y={y1,y2,y3,....yn}，要求最长公共子序列，可以按照下面方式递归计算，：当xm=yn时，找出Xm-1和Yn-1的最长公共子序列，然后在其尾部加上xm即可得到X和Y的最长公共子序列，当xm不等于yn的时候，必须解决两个子问题即找出Xm-1和Y的最长公共子序列和Yn-1和X的最长工公共子序列，然后这两个里面较长者为X和Y的最长公共子序列！

首先建立子问题的最优值递归关系。用C[i][j]表示Xi和Yj的最长公共子序列的长度。其中，Xi={x1,x2,x3......xi},Yj={y1,y2,y3,....yn}，当i=0或者j=0时空序列是Xi和Yj的最长公共子序列，故因此，C[i][j]=0，即有

代码如下：

 
C[i][j]=0，（i=0,j=0)
C[i][j]=c[i-1][j-1]+1,xi=yj
C[i][j]=max(C[i][j-1],C[i-1][j]),i,j>0,xi不等于yj

 
#include<stdio.h>
#include<string.h>
int c[1002][1002]={0},b[1002][1002]={0};
int main()
{
   int N,m,n,i,j;
   char x[1002],y[1002];
  scanf("%d",&N);
  while(N--)
  {
        scanf("%s",x);
        scanf("%s",y);
        m=strlen(x);
        n=strlen(y);
        for(i=1;i<=m;i++)
           for(j=1;j<=n;j++)
           {
               if(x[i-1]==y[j-1]) 
               {
                  c[i][j]=c[i-1][j-1]+1;
                  b[i][j]=1;
                 }
              else if(c[i-1][j]>=c[i][j-1])
              {
                   c[i][j]=c[i-1][j];
                   b[i][j]=2;
               }
               else
               {
                   c[i][j]=c[i][j-1];
                   b[i][j]=3;
              }
            }
          printf("%d\n",c[m][n]);
  }
  //system("pause");
  return 0;
}

利用上面的记录数组b[i][j]可以够找出最长公共子序列！

如下

#include<stdio.h>
#include<string.h>
int c[1002][1002]={0},b[1002][1002]={0},k;
char z[1002];
void lcs(int i,int j,char *x,char *y)
{
    if(i<0||j<0) return ;
    if(b[i][j]==1) 
    {
       //printf("%c",x[i-1]);
       z[k++]=x[i-1];
      lcs(i-1,j-1,x,y);
     }
   else if(b[i][j]==2) lcs(i-1,j,x,y);
   else lcs(i,j-1,x,y);
}     
int main()
{
   int N,m,n,i,j;
   char x[1002],y[1002];
  scanf("%d",&N);
  while(N--)
  {
        scanf("%s",x);
        scanf("%s",y);
        m=strlen(x);
        n=strlen(y);
        for(i=1;i<=m;i++)
           for(j=1;j<=n;j++)
           {
               if(x[i-1]==y[j-1]) 
               {
                  c[i][j]=c[i-1][j-1]+1;
                  b[i][j]=1;
                 }
              else if(c[i-1][j]>=c[i][j-1])
              {
                   c[i][j]=c[i-1][j];
                   b[i][j]=2;
               }
               else
               {
                   c[i][j]=c[i][j-1];
                   b[i][j]=3;
              }
            }
          printf("%d\n",c[m][n]);
          k=0;
          lcs(m,n,x,y);
          for(i=c[m][n]-1;i>=0;i--)
             printf("%c",z[i]);          
          printf("\n");
  }
  //system("pause");
  return 0;
}