exLCS

exLCS

时间限制:1s

空间限制:512MB

题目描述:

给出两个仅有小写字母组成的字符串str1 和str2,试求出两个串的最长公共子序列。

公共子序列定义如下:

若有a1 < a2 < … < ak 和b1 < b2 < … < bk,满足

str1[ai] = str2[bi];

则称找到了一个长度为k 的公共子序列。

输入格式:

第一行一个字符串str1。

第二行一个字符串str2。

输出格式:

一行,一个整数,表示str1 与str2 的最长公共子序列的长度。

样例读入:

abdcba

abafdsfa

样例输出:

4

样例解释:

如果字符串从0 开始标号,可以验证{an} = {0,1,2,5},{bn} = {0,1,4,7} 是满足要求的方案。

数据范围:

对于10% 的数据,保证jstr1j ⩽ 10; jstr2j ⩽ 10

对于30% 的数据,保证jstr1j ⩽ 20; jstr2j ⩽ 30

对于60% 的数据,保证jstr1j ⩽ 1000; jstr2j ⩽ 1000

对于100% 的数据,保证jstr1j ⩽ 1000; jstr2j ⩽ 106

题解

由于len(s1)很小,最长长度一定小于len(s1)。

考虑变换状态定义,f[i][j]表示匹配到s1[i],公共子序列长度为j时,s2匹配到的最小的位置。

再记录next[i][j]表示当前在s2的i位置,下一个j字母出现的位置。

转移方程:

f[i+1][j]=min(f[i+1][j],f[i][j])

f[i+1][j+1]=min(f[i+1][j+1],next[f[i][j]+1][a[i+1]-'a'])

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#define LL long long
using namespace std;
const int N=1005,M=1000005,inf=0x3f3f3f3f;
int n,m,ans;
int f[N][N],nxt[M][26];
char a[N],b[M];
void solve(){
    f[0][1]=nxt[0][a[0]-'a'];
    for(int i=0;i<=n;i++)f[i][0]=-1;
    for(int i=0;i<n-1;i++)
        for(int j=0;j<=n&&f[i][j]<inf;j++){
            f[i+1][j]=min(f[i+1][j],f[i][j]);
            if(j<n)f[i+1][j+1]=min(f[i+1][j+1],nxt[f[i][j]+1][a[i+1]-'a']);
        }
}
int main(){
    scanf("%s%s",a,b);
    n=strlen(a);
    m=strlen(b);
    for(int i=0;i<=n;i++)
        for(int j=0;j<=n;j++)
            f[i][j]=inf;
    for(int i=0;i<26;i++)
        nxt[m][i]=inf;
    for(int i=m-1;i>=0;i--){
        memcpy(nxt[i],nxt[i+1],sizeof(nxt[i]));
        nxt[i][b[i]-'a']=i;
    }
    solve();
    for(int i=n;i>0;i--)
        if (f[n-1][i]<inf){
            ans=i;
            break;
        }
    printf("%d
",ans);
    return 0;
}
原文地址:https://www.cnblogs.com/chezhongyang/p/7648806.html