最长公共子序列问题(LCS) - 动态规划

参考自 <<

题目

题目描述

给出两个字符串，求两个字符串的最长公共字串。

输入

abcfbc abfcab
programming contest
abcd mnp

输出

4
2
0

题解

分析

可以看出：
F[i][j]=F[i-1][j-1]+1；（a[i]==b[j]）
F[i][j]=max(F[i-1][j],F[i][j-1])(a[i]!=b[j]);
n由于F(i,j)只和F(i-1,j-1), F(i-1,j)和F(i,j-1)有关, 而在计算F(i,j)时, 只要选择一个合适的顺序, 就可以保证这三项都已经计算出来了, 这样就可以计算出F(i,j). 这样一直推到f(len(a),len(b))就得到所要求的解了.

代码

#include<stdio.h>
#include<string.h>
#include<math.h>
#include <iostream>
#include <algorithm>
using namespace std;
int dp[1100][1100];
char s1[1000],s2[1000];
int main()
{
    while(~scanf("%s%s",s1,s2)){
        int len1=strlen(s1);
        int len2=strlen(s2);
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=len1;++i){
            for(int j=1;j<=len2;++j){
                if(s1[i-1]==s2[j-1])
                    dp[i][j]=dp[i-1][j-1]+1;
                else
                    dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
            }
        }
        printf("%d\n",dp[len1][len2]);
    }
    return 0;
}

拓展 - 打印序列 - 求序列路径

给出两个字符串，求两个字符串的最长公共字串及打印公共子串

输入

abcfbc abfcab
programming contest
abcd mnp

输出

LCS=4
a b c b
LCS=2
o n
LCS=0

代码

#include<stdio.h>
#include<string.h>
#include<math.h>
#include <iostream>
#include <algorithm>
using namespace std;

int dp[110][110],path[110][110];
char s1[110],s2[110];
void LCS(int len1,int len2)
{
    for(int i=1;i<=len1;i++){
        for(int j=1;j<=len2;j++){
            if(s1[i-1]==s2[j-1]){
                dp[i][j]=dp[i-1][j-1]+1;
                path[i][j]=1;
            }else if(dp[i][j-1]<=dp[i-1][j]){
                dp[i][j]=dp[i-1][j];
                path[i][j]=2;
            }else{
                dp[i][j]=dp[i][j-1];
                path[i][j]=3;
            }

        }
    }
    printf("LCS=%d\n",dp[len1][len2]);
}
void print(int x,int y)
{
    if(x==0||y==0) return;//注意边界
    if(path[x][y]==1){
        print(x-1,y-1);
        printf("%c ",s1[x-1]);
    }else if(path[x][y]==2){
        print(x-1,y);
    }else if(path[x][y]==3){
        print(x,y-1);
    }
}
int main()
{
    string s;
    while(~scanf("%s%s",s1,s2)){
        memset(dp,0,sizeof(dp));
        int len1=strlen(s1),len2=strlen(s2);
        LCS(len1,len2);
        print(len1,len2);
        cout<<endl;
    }
    return 0;
}

拓展 - 打印序列 - 求序列路径

给定两个字符串，输出两个字符串公共不重复部分。

输入

apple peach
ananas banana
pear peach

输出

appleach
bananas
pearch

代码

将两个字符串结合起来，他们的公共子串只输出一次。
根据LCS的原理，将每个字符都进行标记，看两个字符串中对应的字符究竟处于什么状态，然后输出，其标记为公共子串的字符只输出一次即可。
注意数组的初始化。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include <iostream>
#include <algorithm>
using namespace std;

int dp[110][110],path[110][110];
char s1[110],s2[110];
void LCS(int len1,int len2)
{
    for(int i=0;i<=len1;i++)
        path[i][0]=2;
    for(int i=0;i<=len2;i++)
        path[0][i]=3;

    for(int i=1;i<=len1;i++){
        for(int j=1;j<=len2;j++){
            if(s1[i-1]==s2[j-1]){
                dp[i][j]=dp[i-1][j-1]+1;
                path[i][j]=1;
            }else if(dp[i][j-1]<=dp[i-1][j]){
                dp[i][j]=dp[i-1][j];
                path[i][j]=2;
            }else{
                dp[i][j]=dp[i][j-1];
                path[i][j]=3;
            }

        }
    }
    /*for(int i=1;i<=len1;i++){
        for(int j=1;j<=len2;j++){
            printf("%d ",path[i][j]);
        }
        printf("\n");
    }
    printf("LCS=%d\n",dp[len1][len2]);
    */
}
void print(int x,int y)
{
    if(x==0&&y==0) return;//注意边界
    if(path[x][y]==1){
        print(x-1,y-1);
        printf("%c",s1[x-1]);
    }else if(path[x][y]==2){
        print(x-1,y);
        printf("%c",s1[x-1]);
    }else if(path[x][y]==3){
        print(x,y-1);
        printf("%c",s2[y-1]);
    }
}
int main()
{
    string s;
    while(~scanf("%s%s",s1,s2)){
        memset(dp,0,sizeof(dp));
        int len1=strlen(s1),len2=strlen(s2);
        LCS(len1,len2);
        print(len1,len2);
        cout<<endl;
    }
    return 0;
}