LCS

LCS-最长公共子序列

有时候我们会遇到要求两个字符串的公共最长序列,这就是LCS问题,它通常是用动态规划来求解.

解析

假设现在有两个字符串A, B,字符串A的长度为 n n n, 字符串B的长度为 m m m, 既然是 d p dp dp, 那么就设 d p [ i ] [ j ] dp[i][j] dp[i][j]表示字符串从0的到 i i i位置和字符串从0到j位置的最长公共子序列, 那么就可以讨论
A [ i ] = = B [ j ] , d p [ i ] [ j ] = d p [ i − 1 ] [ j − 1 ] + 1 A[i]==B[j] , dp[i][j]=dp[i-1][j-1]+1 A[i]==B[j],dp[i][j]=dp[i1][j1]+1 A [ i ] ! = B [ j ] , d p [ i ] [ j ] = m a x ( d p [ i − 1 ] [ j ] , d p [ i ] [ j − 1 ] ) A[i]!=B[j], dp[i][j]=max(dp[i-1][j],dp[i][j-1]) A[i]!=B[j],dp[i][j]=max(dp[i1][j],dp[i][j1])
这个式子还是比较好理解的,当两个位置的字符相等时,他们的最长公共子序列长度等于前一个位置+1,这个是显然的。当两个位置的字符不相等时,他们的公共子序列的长度就为两个前面一个的最大值。在计算过程中需要枚举字符串A的每一位去和B的每一位比较,要单独处理0位置的值。
说了这么多,看道题吧,传送门:
https://www.nowcoder.com/practice/b43fb39898f448e39adbcffde5ff0dfc?tpId=90&tqId=30810&tPage=2&rp=2&ru=/ta/2018test&qru=/ta/2018test/question-ranking

这道题的意思就是把字符串分成两部分,然后在分别求最长的公共子序列,那么就可以考虑把字符串分成两部分在分别 d p dp dp, 可以先预处理下,就是把只出现一次的字符去掉,因为这道题规模太小了,不这样做照样可以A。

AC代码

//  小学生一发的刷题之路
//  3月3号-平方串;
//
//
//

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#include <queue>
#include <deque>                //双向队列;
#include <cmath>
#include <set>
#include <stack>
#include <map>
#include <vector>
#include <cstdlib>
#include <iomanip>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
const double PI=acos(-1.0);
const double eps=1e-10;
const int maxn=1e5+5;
const int maxm=1e3+5;
const ll mod=1e9+7;
const int INF=1e8;
template<class T>
inline void read(T &ret){       //快速输入模版;
    ret=0;
    int f=1;
    char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-') f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9'){
        ret=ret*10+c-'0';
        c=getchar();
    }
    ret*=f;
}
template <class T>
inline void out(T ret){     //快速输出模版;
    if(ret>9)
    {
        out(ret/10);
    }
    putchar(ret%10+'0');
}
string str;
int dp[100][100];

int sovle(int m,int n){
    
    for(int i=0;i<=m;i++)
        for(int j=m+1;j<n;j++){
            if(str[i]==str[j]){
                if(i==0||j==m+1){
                    dp[i][j]=1;
                }else{
                    dp[i][j]=dp[i-1][j-1]+1;
                }
            }else{
                if(i==0&j==m+1){
                    dp[i][j]=0;
                }else if(i==0){
                    dp[i][j]=dp[i][j-1];
                }else if(j==m+1){
                    dp[i][j]=dp[i-1][j];
                }else{
                    dp[i][j]=max(dp[i][j-1],dp[i-1][j]);
                }
            }
        }
    return dp[m][n-1];
}

int main()
{
    cin>>str;
    int n=str.length();
    int ans=0;
    for(int i=0;i<n;i++){
        int tmp=sovle(i,n);
        ans=max(ans,tmp);
    }
    printf("%d\n",2*ans);
    return 0;
}

新的开始,每天都要快乐哈!
在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值