最大公共子序列

最大公共子序列
2.子问题的递归结构
设X={x1,x2,x3,...,xn}和Y={y1,y2,y3,...,yn}的最长公共子序列，可按以下方式递归计算：当xm=yn时，找出xm-1和yn-1的最长子序列，然后在其尾部加上xm（=yn）即可得X和Y的最长公共子序列。当xm<>yn时，必须解两个子问题，即找出xm-1和Y的一个最长公共子序列及X和Yn-1的一个最长公共子序列。这两个公共子序列中较长者即为X和Y的最长公共子序列。
由此递归结构容易看出最长公共子序列问题具有子问题重叠性质。例如，在计算X和Y的最长公共子序列时，可能要计算X和Yn-1及Xm-1和Y的最长公共子序列。而这两个子问题都包含一个公共子问题，即计算Xm-1和Yn-1的最长公共子序列。
首先建立子问题最优值的递归关系。用c[i][j]记录序列Xi和Yj的最长公共子序列的长度。其中，Xi={x1,x2,x3,...,xi}；Yj={y1,y2,y3,...,yj}。当i=0,j=0时，空序列是Xi和Yi的最长公共子序列，故此时c【i】【j】 = 0 。在其他的情况下，由最优子结构性质（当为题的最优解包含了其子问题的最优解时，称该问题具有最优子结构性质)可建立递归关系如下：

 0                         i= 0，j =0

c【i】【j】=               c【i-1】【j-1】+1          i，j>0;xi=yj

                          Max { c[i][j-1],c[i-1][j] }    i,j>0;xi<>yj

3.计算最优值。
    c【i】【j】存储Xi和Yi的最长公共子序列的长度，b【i】【j】记录c【i】【j】的值是哪一个子问题的解得到的，这在构造最长公共子问题序列时要用到。问题的最优值，记X和Y的最长公共子序列的长度记录于c【m】【n】中。
4.构造最长公共子序列
    首先从b 【 m 】【 n 】 开始，以其值在数组b中搜索。当b【i】【j】 =1 时 ， 表示Xi和Yj的最长公共子序列由Xi-1和Yj-1的最长公共子序列在尾部上加xi所得到的子序列.......……

算法设计例题：最长公共子序列（DP）

memory limit: 65536KB    time limit: 500MS

accept: 24    submit: 77

Description

一个给定序列的子序列是在该序列中删去若干元素后得到的序列。确切地说，若给定序列 X = { x1，x2，…，xm }，则另一序列Z={ z1，z2，…，zk }，X的子序列是指存在一个严格递增下标序列{ i1，i2，…，ik }，使得对于所有 j = 1，2，…，k ，有 zj = xij

给出两个字符序列X和Y，求出它们的最长公共子序列。

Input

输入的第一行为测试样例的个数T（ T < 40 ），接下来有T个测试样例。每个测试样例的第一行是字符串X，第二行是字符串Y。X和Y只包含大写字母，且长度不大于1000。

Output

对应每个测试样例输出一行，只有一个整数，表示字符串X和字符串Y的最长公共子序列的长度。

Sample Input

2
ABCDE
ACE
AAABBBCCC
AABBCC

Sample Output

3
6

#include"iostream"
#include<stdio.h>
using namespace std ;

class Element {
    public :

    int LcsLength (char x[],int xn , char y[], int yn , int b [1000][1000] )
    {
        int m = xn  ;
        int n = yn ;
        int c[1000][1000]  ;
        int i  ;
        int j;
        for (  i =1 ; i <= m ; i++ )
            c[i][0] = 0 ;
        for(i =1 ; i <=n ; i++ )
            c[0][i] = 0 ;
        for (  i =1 ; i <= m ; i++ )
            for(  j=1 ; j <= n ; j++ )
            {
                if( x[i] == y [j] )
                {
                     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 ;
                }
            }

            return c [ m ][ n ] ;
    }

    void Lcs(int i , int j , char x[] , int b [1000][1000])
    {
        if( i == 0 || j ==0 )
            return;
        if( b[i][j] == 1 )
        {

            Lcs ( i-1 , j-1 , x , b );
             cout<<x[i];
        }else if ( b[i][j] ==2 ){
            Lcs( i -1 , j , x , b );
        }else {
            Lcs( i  , j -1  , x , b );
        }
    }
};

int main()
{
    Element e ;
    int T  ;
    char x [1000] ;
    char y [1000] ;
    int xn  ;
    int yn  ;
    int b [1000][ 1000 ] ;
    cin>>T;
    while( T > 0 )
    {
        xn = 1 ;
        yn = 1 ;
        cin>>(x+1);
        cin>>(y+1);
        while(x[xn]!='\0')
        {
            xn++;
        }
        xn--;
        while(y[yn]!='\0')
        {
            yn++;
        }
        yn--;
        cout<<e.LcsLength(x,xn,y,yn,b)<<endl;
        T--;
    }
    return 0 ;
}