动态规划-最长公共子序列

问题描述：

给定两个序列和，找出X和Y的一个最长公共子序列。

算法分析：

暴力法：

X有2^m个子序列，逐个匹配Y的子序列。

动态规划法：

(1)分析最优解的结构特征：

假设是X和Y的最长公共子序列。

证明：最优解包含局部最优解。

①，是和的最长公共子序列。

若不是，设M是，则，但，矛盾！

②，去掉，则是和的最长公共子序列。 

如果不是，设M是，则，在后+，则M也是和的最长公共子序列（因为），则，矛盾！

③，去掉，则是和的最长公共子序列。 

如果不是，设M是，则，在后+，则M也是和的最长公共子序列（因为），则，矛盾！

(2)建立最优值的递归式：c[i][j]是和的最长公共子序列长度。

若，则c[i][j]=c[i-1][j-1]+1 

若，比较和 与 和最长公共子序列，则c[i][j]=max{ c[i][j-1]，c[i-1][j] }

最长公共子序列长度递归式：

(3)自底向上计算最优值：

(4)构造最优解：

c[i][j]有3个来源，用辅助数组b[i][j]记录这3个来源。

c[i][j]=c[i-1][j-1]+1，则b[i][j]=1;

c[i][j]=c[i][j-1]，则b[i][j]=2; 

c[i][j]=c[i-1][j]，则b[i][j]=3; 

#include <iostream>
#include <cstring>
using namespace std;

#define N 1002
int c[N][N],b[N][N];
char s1[N],s2[N];
int len1,len2;
void LCSL(){
    int i,j;
    for(i=1;i<=len1;i++){
        for(j=1;j<=len2;j++){
            if(s1[i-1]==s2[j-1]){
                c[i][j]=c[i-1][j-1]+1;
                b[i][j]=1;
            }else{
                if(c[i][j-1]>=c[i-1][j]){
                    c[i][j]=c[i][j-1];
                    b[i][j]=2;
                }else{
                    c[i][j]=c[i-1][j];
                    b[i][j]=3;
                }
            }
        }
    }
}

void print(int i,int j){
    if(i==0||j==0) return;
    if (b[i][j]==1){
        print(i-1,j-1);
        cout<<s1[i-1];
    }else if(b[i][j]==2){
        print(i,j-1);
    }else{
        print(i-1,j);
    }
}

int main() {
    int i,j;
    cout<<"输入字符串s1:"<<endl;
    cin>>s1;
    cout<<"输入字符串s2:"<<endl;
    cin>>s2;
    len1=strlen(s1);
    len2=strlen(s2);
    for (i=0;i<=len1;i++){
        c[i][0]=0;
    }
    for (j=0;j<=len2;j++){
        c[0][j]=0;
    }
    LCSL();
    cout<<"s1和s2的最长公共子序列长度是： "<<c[len1][len2]<<endl;
    cout<<"s1和s2的最长公共子序列是： ";
    print(len1,len2);
    return 0;
}