程序员代码面试指南:最长公共子串问题
动态规划:dp[M][N]矩阵,d[i][j]表示必须以str1[i]和str2[j]为公共子串的最后一个字符时,公共子串的长度。
1. 第一列,如果str1[i]==str2[0],则dp[i][0]=1,否则dp[i][0]=0;
2. 第一行,同理。
3. 其他位置,如果str1[i] != str2[j],则dp[i][j]=0;如果str1==str2[0],dp[i][j]=dp[i-1][j-1]+1;
得到dp后,找到最大值,即找到了公共子串的最后一个字符,从该位置向前截取字符串即为最长公共子串。
C++实现如下:
#include <string>
using namespace std;
string getLCSTR(const string& str1, const string&str2)
{
string result;
if (str1.empty() || str2.empty())
{
return result;
}
int **dp = new int*[str1.length()];
for (int i = 0; i < str1.length(); i++)
{
dp[i] = new int[str2.length()];
}
for (int i = 0; i < str1.length(); i++)
{
dp[i][0] = str1.at(i) == str2.at(0) ? 1 : 0;
}
for (int j = 0; j < str2.length(); j++)
{
dp[0][j] = str1.at(0) == str2.at(j) ? 1 : 0;
}
for (int i = 1; i < str1.length(); i++)
{
for (int j = 1; j < str2.length(); j++)
{
dp[i][j] = str1.at(i) == str2.at(j) ? dp[i-1][j-1] + 1 : 0;
}
}
int end = 0;
int length = 0;
for (int i = 0; i < str1.length(); i++)
{
for (int j = 0; j < str2.length(); j++)
{
if (dp[i][j] > length)
{
length = dp[i][j];
end = i;
}
}
}
result = str1.substr(end - length + 1, length);
for (int i = 0; i < str1.length(); i++)
{
delete[] dp[i];
}
delete[] dp;
return result;
}
int main(int argc, char* argv[])
{
string str1 = "1AB2345D";
string str2 = "12345EF";
string common = getLCSTR(str1,str2);
printf("%s\n", common.c_str());
return 0;
}