1265:【例9.9】最长公共子序列
【题目描述】
一个给定序列的子序列是在该序列中删去若干元素后得到的序列。确切地说,若给定序列X=<x1,x2,…,xm>X=<x1,x2,…,xm>,则另一序列Z=<z1,z2,…,zk>Z=<z1,z2,…,zk>是X的子序列是指存在一个严格递增的下标序列<i1,i2,…,ik><i1,i2,…,ik>,使得对于所有j=1,2,…,k有:
Xij=ZjXij=Zj
例如,序列Z=<B,C,D,B>是序列X=<A,B,C,B,D,A,B>的子序列,相应的递增下标序列为<2,3,5,7>。给定两个序列X和Y,当另一序列Z既是X的子序列又是Y的子序列时,称Z是序列X和Y的公共子序列。例如,若X=<A,B,C,B,D,A,B>和Y=<B,D,C,A,B,A>,则序列<B,C,A>是X和Y的一个公共子序列,序列 <B,C,B,A>也是X和Y的一个公共子序列。而且,后者是X和Y的一个最长公共子序列.因为X和Y没有长度大于4的公共子序列。
给定两个序列X=<x1,x2,…,xm>X=<x1,x2,…,xm>和Y=<y1,y2…yn>Y=<y1,y2….yn>.要求找出X和Y的一个最长公共子序列。
【输入】
共有两行。每行为一个由大写字母构成的长度不超过1000的字符串,表示序列X和Y。
【输出】
第一行为一个非负整数。表示所求得的最长公共子序列的长度。若不存在公共子序列.则输出文件仅有一行输出一个整数0。
【输入样例】
ABCBDAB
BDCABA
【输出样例】
4
【提示】
最长公共子串(Longest Common Substirng)和最长公共子序列(Longest Common Subsequence,LCS)的区别为:子串是串的一个连续的部分,子序列则是从不改变序列的顺序,而从序列中去掉任意的元素而获得新的序列;也就是说,子串中字符的位置必须是连续的,子序列则可以不必连续。字符串长度小于等于1000。
最长公共子序列是简单动态规划的一种题型。
其思路也是打表+缩小规模的递归表达。
思路:分情况讨论。
设有两个字符串,A(1,2,3,……,x)共有x个字符,B(1,2,3,……,y)共有y个字符。设函数L(x,y)是求字符串A和字符串B分别从头到x,y字符的最长公共子序列的长度。
1、从两个字符串的最后一个考虑。
a)如果Ax==By,两个字符串的最后一个字符相等,则公共子序列中一定有这个字符。L(x,y)=L(x-1,y-1)+1。
b)如果Ax!=By,两个字符串的最后一个字符不相等,则这两个字符至少会有一个不在最长公共子序列中,L(x,y)=max( L(x-1,y), L(x,y-1) )
2、边界条件,x=0或者y=0时,(下标小于最前面的字符时),返回0。
代码:
#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
//思路:最长公共子序列 如果两个子序列的结尾一样,那一定子序列包括最后那个
//如果不一样,一定至少一个不会有用,所以 max(c[x-1][y],c[x][y-1])。
int c[1000][1000];
string a,b;
int fag=<