最长公共子串问题

程序员代码面试指南：最长公共子串问题

动态规划：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;
}