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; }