最长公共字串
目前,在各大网站或者是博客园,有着不少的关于这个问题的讲解,但是目前觉得都不是很清晰明白!现在普遍的解法都是用动态规划的方式来做这道笔试题。或者是经典的面试题。
这里我做一个自认为比较清楚的讲解。说到动态规划很多人都不会感到很陌生。但是有多少知道动态规划什么?精髓在哪里? 这里仅仅提供一些个人见解,首先,动态规划不是一种解题的方式,它只是一种以空间换时间的优化方式,同时可能会对代码有一个比较好的优化方式。也就是说它并不是一种算法。如果有人告诉你这个题用动态规划去做,请抽他,因为他这样说等于啥都没有说!
动态规划的精髓在于状态转移方程。我们首先应该找出最简单的一种情况,然后再去找到 从一个状态转移到另外一个状态的对应关系。可能还是有很多人不知道这是什么意思。
举个简单的例子,从数字5到数字6需要加一。很好理解吧!
现在我们来看这个题。首先,这个题和最长公共子序列是不一样的。就是最长公共子序列是可以不连续的,但是最长公共字串必须是连续的。
假设需要求的字符串为 str1 , str2 .函数 f(m,n) 求分别以str1[m] , str2[n] 结尾的公共字符串长度。
首先,我们应该考虑 f(m,n),在二维数组的(m,n)出代表的意义,这样才可以更好的找出状态转移方程!这里该位置代表着一 str1[m],str2[n]结尾的公共字串的长度是多少。然后我们考虑当前m,n位置处,如果str1[m]==str2[n],那么公共字串,长度肯定就等于str1[m-1],str2[n-1] + 1,因为前面如果是字串就加1,如果不是就正常m,n处是开头也是结尾,然后如果str1[m]!=str2[n],就说明如果加上最后一个m,n处就 不是公共字串了,所以这里是0,
最后就有了状态转移方程! if(str1[m] == str2[n]) { f(m,n) = f(m-1,n-1) + 1;}else{ f(m,n) = 0;}
剩下的就是代码了!
#include<iostream>
#include<string.h>
using namespace std;
int c[10000][10000];
char str1[10000];
char str2[10000];
void func(int m,int n){
if((m<0)||(n<0)) return ;
for(int i=0;i<m;i++)
for(int j=0;j<n;j++) c[i][j]=0;
int besti=0,bestj=0;
int count=0;
for(int i=0;i<m;i++)
for(int j=0;j<n;j++)
if(str1[i]==str2[j]) {
if((i==0)||(j==0)) c[i][j]=1; //增加判断是否为第一个,第一个不能再向下
else c[i][j]=c[i-1][j-1] + 1 ;
}
else c[i][j]=0;
/*
for(int i=0;i<m;i++){
for(int j=0;j<n;j++)
cout<<c[i][j]<<' ';
cout<<endl;
}
*/
for(int i=0;i<m;i++)
for(int j=0;j<n;j++)
if(c[i][j]>count) { count=c[i][j]; besti=i; bestj=j; } //查找最长子字符串
cout<<count<<endl;
if(count==0) cout<<endl; //公共子字符串长度为1,输出空行
else
for(int i=besti-count+1;i<=besti;i++) cout<<str1[i]; //否则输出靠X左边(如果多个)子字符串
}
int main() {
int m=0;
int n=0;
cin>>str1;
cin>>str2;
m=strlen(str1);
n=strlen(str2);
func(m,n);
return 0;
}