【51nod】---1006 最长公共子序列Lcs(动态规划&&字符串LCS)

题目链接这里呀
1006 最长公共子序列Lcs
基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注
给出两个字符串A B,求A与B的最长公共子序列(子序列不要求是连续的)。
比如两个串为:

abcicba
abdkscab

ab是两个串的子序列,abc也是,abca也是,其中abca是这两个字符串最长的子序列。
Input
第1行:字符串A
第2行:字符串B
(A,B的长度 <= 1000)
Output
输出最长的子序列,如果有多个,随意输出1个。
Input示例
abcicba
abdkscab
Output示例
abca

分析:注意哦,子序列不是子段,子序列可以是不连续的。
求最长公共子序列(LCS)用动态规划。
两字符串str1,str2,用dp[i][j]表示str1[1]~str1[i]和str2[1]~str2[j]两个子段最长公共子序列的长度。
若str1[i]==str2[j],那么此时的LCS长度就相当于在(i-1,j-1)的LCS长度上加1,即dp[i][j]=dp[i-1][j-1]+1;
若str1[i]!=str2[j],此时的LCS长度应取(i-1,j)和(i,j-1)的较大者,因为要求最长嘛!即dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
所以,dp[i][j]={ dp[i-1][j-1]+1; (str1[i]==str2[j])
max(dp[i-1][j],dp[i][j-1];(str1[i]!=str2[j])
}
dp[len1][len2]即为LCS长度;
若要求出LCS,还需要用到栈的数据结构哦。
用两个‘指针’ant1,ant2倒序指向两字符串中的字符,
若str1[ant1]==str1[ant2],则入栈(任一个就行哦),同时ant1–,ant2–;
若不相等,肯定首先字符不能入栈了哦,那么两个‘指针’怎么办呢?
两种情况:
如果dp[ant1-1][ant2]>dp[ant1][ant2-1],说明(ant1-1,ant2)这一部分的LCS比(ant1,ant2-1)的要长哦,那么当然就要研究看str1[ant1-1]和str1[ant2]是否相等咯,也就是ant1–,前移;
否则,ant2–。
LCS的全部字符入栈后,再依次输出。
代码如下:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX_N 1000
#include<stack>
int dp[MAX_N+11][MAX_N+11];
char str1[MAX_N+11],str2[MAX_N+11];
int main()
{
    scanf("%s %s",str1+1,str2+1);
    memset(dp,0,sizeof(dp));
    str1[0]='0';
    str2[0]='0';
    int len1=strlen(str1)-1;
    int len2=strlen(str2)-1;

    int ans=0;
    int i,j;
    for(i=1;i<=len1;i++)
    {
        for(j=1;j<=len2;j++)
        {
            if(str1[i]==str2[j])
                dp[i][j]=dp[i-1][j-1]+1;
            else
                dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
        }
        ans=max(ans,dp[i][len2]);//更新为当前最长,最后得出全部最长 
    }
    int ant1=len1;
    int ant2=len2;
    stack<char> stk;
    while(ant1>0&&ant2>0)
    {
        if(str1[ant1]==str2[ant2])//若对应两字符相等,说明是LCS组成部分,入栈 
        {
            stk.push(str1[ant1]);
            ant1--;
            ant2--;
        }
        else if(dp[ant1-1][ant2]>dp[ant1][ant2-1])
            ant1--;
        else
            ant2--;
    }
    while(!stk.empty() )
    {
        printf("%c",stk.top());
        stk.pop() ;
    }
    printf("\n");
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值