问题 A: 求最长公共子串(串)
时间限制: 1 Sec 内存限制: 128 MB
提交: 34 解决: 19
题目描述
求采用顺序结构存储的串s和串t的一个最长公共子串,若没有则输出false,若最长的有多个则输出最先出现的那一串。
输入
输入两个字符串
输出
输出公共子串
样例输入
abcdef adbcef
样例输出
bc
经验总结
求最长公共子串(连续)的方法,可以用KMP,当然也可以用字符串hash,分别计算两个字符串的所有子串的hash值,然后一一对比,当两个字符串的hash值相等时,如果长度大于之前访问得到的公共子串长度的最大值,则更新最大值,并存储此子串的起始位置,最终得到最长的公共子串~
正确代码
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
#include <iostream>
#include <string>
using namespace std;
typedef long long LL;
const int maxn=1010;
const LL p=1e7+19;
const LL mod=1e9+7;
LL powP[maxn],H1[maxn]={0},H2[maxn]={0};
struct Node
{
int hashValue,length,start,end;
Node(int a,int b,int c,int d):hashValue(a),length(b),start(c),end(d){};
};
vector<Node> pr1,pr2;
void init(int len)
{
powP[0]=1;
for(int i=1;i<=len;++i)
{
powP[i]=(powP[i-1]*p)%mod;
}
}
void calH(LL H[],string &str)
{
H[0]=str[0];
for(int i=1;i<str.length();++i)
{
H[i]=(H[i-1]*p+str[i])%mod;
}
}
int calSingleSubH(LL H[],int i,int j)
{
if(i==0)
return H[j];
return ((H[j]-H[i-1]*powP[j-i+1])%mod+mod)%mod;
}
void calSubH(LL H[],int len,vector<Node> &p)
{
for(int i=0;i<len;++i)
for(int j=i;j<len;++j)
{
int value=calSingleSubH(H,i,j);
p.push_back(Node(value,j-i+1,i,j));
}
}
string getMax(string str1)
{
string str;
int ans=0;
for(int i=0;i<pr1.size();++i)
for(int j=0;j<pr2.size();++j)
{
if(pr1[i].hashValue==pr2[j].hashValue)
{
if(pr1[i].length>ans)
{
ans=pr1[i].length;
str=str1.substr(pr1[i].start,pr1[i].end);
}
}
}
return str;
}
int main()
{
string str1,str2;
getline(cin,str1);
getline(cin,str2);
init(max(str1.length(),str2.length()));
calH(H1,str1);
calH(H2,str2);
calSubH(H1,str1.length(),pr1);
calSubH(H2,str2.length(),pr2);
cout<<getMax(str1)<<endl;
return 0;
}