最长公共子序列问题——LCS算法

最长公共子序列问题——LCS算法

问题描述:

给出两个字符串A B，求A与B的最长公共子序列（子序列不要求是连续的）。比如两个串为：

abcicba
abdkscab

ab是两个串的子序列，abc也是，abca也是，其中abca是这两个字符串最长的子序列。

输入格式:

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

输出格式:

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

输入样例:

abcicba
abdkscab

输出样例:

abca

这题是典型的动态规划题型

●明确原始问题 C[n,m]：X[1,n]和Y[1,m]的最长公共子序列长度

●递推关系确立

推导过程：点击这里

伪码描述：

●Longest-Common-Subsequence(X,Y)

输入两个序列X,Y
输出X和Y的最长公共子序列
n<-length(X)
m<-length(Y)
//初始化
新建二维数组C[0..n,0..m]和追踪二维数组rec[0..n,0..m]
for i<-0 to n do
    C[i,0]<-0
end
for j<-0 to m do
    C[0,i]<-0
end
//动态规划
for i<-1 to n do
    for j<-1 to m do
        //当X[i]=Y[i]时
        if X[i]=Y[i] then
            C[i,j]<-C[i-1,j-1]+1
            rec[i,j]<-"LU"
        end
        else if C[i-1,j]≥C[i,j-1] then
            C[i,j]<-c[i-1,j]
            rec[i,j]<-"U"
        end
        else
            C[i,j]<-C[i,j-1]
            rec[i,j]<-"L"
        end
    end
end
return C,rec

●Print-LCS(rec,X,i,j)

输入追踪数组rec,序列X,当前位置i和j
输出X[1..i]和Y[1..j]的最长公共子序列
if i=0 or j=0 then
    //终止递归
    return NULL
end
    //三种递归式的情况
if rec[i,j]="LU" then
    Print-LCS(rec,X,i-1,j-1)
    print x[i]
end
else if rec[i,j]="U" then
    Print-LCS(rec,X,i-1,j)
end
else
    Print-LCS(rec,X,i,j-1)
end

AC代码:

#include<stdio.h>
#include<string.h>

using namespace std;
const int MAXSIZE=1001;
int C[MAXSIZE][MAXSIZE];
char rec[MAXSIZE][MAXSIZE];
char X[MAXSIZE];
char Y[MAXSIZE];

void Initialization(int n,int m){//初始化
    for(int i=0;i<=n;i++)
        C[i][0]=0;
    for(int i=0;i<=m;i++)
        C[0][i]=0;
}
void Lcs(int n,int m){//LCS算法
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if(X[i]==Y[j]){
                C[i][j]=C[i-1][j-1]+1;
                rec[i][j]='W';
            }else if(C[i-1][j]>=C[i][j-1]){
                C[i][j]=C[i-1][j];
                rec[i][j]='U';
            }else{
                C[i][j]=C[i][j-1];
                rec[i][j]='L';
            }
        }
    }
}
void Print_Lcs(int n,int m){//根据rec数组输出
    int i=n;
    int j=m;
    if(i==0||j==0)return;
    if(rec[i][j]=='W'){
        Print_Lcs(i-1,j-1);
        printf("%c",X[i]);
    }else if(rec[i][j]=='U'){
        Print_Lcs(i-1,j);
    }else{
        Print_Lcs(i,j-1);
    }
}
int main(){ 
    scanf("%s%s",X+1,Y+1);
    size_t n=strlen(X+1);
    size_t m=strlen(Y+1);
    Initialization(n,m);
    Lcs(n,m);
    Print_Lcs(n,m);
}