最长公共子串（DP）

用动态规划做了一个，自己测试了几个测试用例也通过了。

#include <iostream>
using namespace std;
int main()
{
	string a="jahrgiuaerhgiaeg";
	string b="hsiresahrgiuahgise";
	int l1=a.size();
	int l2=b.size();
    if(l1==0||l2==0)
    {
        cout<<"";
        return 0;
    }
	vector<vector<int>>dp(l2,vector<int>(l1,0));
	for(int i=0;i<l1;++i)
	{
		if(a[i]==b[0])
			dp[0][i]=1;
	}
	for(int i=0;i<l2;++i)
	{
		if(a[0]==b[i])
			dp[i][0]=1;
	}
	for(int i=1;i<l2;++i)
	{
		for(int j=1;j<l1;++j)
		{
			if(a[j]==b[i])
				dp[i][j]=dp[i-1][j-1]+1;
            cout<<dp[i][j]<<" ";
		
		}
        cout<<endl;
	}
	int max_l=0;
	int max_i=0;
	int max_j=0;

  for(int i=0;i<l1;++i)
  {
		for(int j=0;j<l2;++j)
		{
			if(dp[j][i]>max_l)
			{
                
				max_i=i;
				max_l=dp[j][i];
			}
				
		
		}
  }
    
    string res;
  while(max_l--)
  {
	res=a[max_i--]+res;
  }
    cout<<res;
  system("pasue");
  
}

输出为：

          思路还是比较简单的，首先得出两个字符串a、b的长度l1和l2，建立一个二维数组dp用于存储以a[j]和b[i]作为结尾的最长公共子串的长度，意思就是说dp[i][j]的意思是以a[j]和b[i]为结尾的最长公共子串的长度，如果说a[j]!=b[i]，也就是说不存在以两者为结尾的子串，所以dp[i][j]=0，如果a[j]==b[i]，那说明可以将a[j]和b[i]作为公共子串的最后一个字符，其长度就是dp[i-1][j-1] + 1.

    然后就是预处理的部分，dp矩阵的第一行，对于其中的某个位置(0,j)来说，如果a[j]==b[0]，那么dp[0][j]=1，否则为0.同理对dp第一列进行处理，如果a[0]==b[i]，那么dp[i][0]=1,否则为0。然后对除去第一行和第一列的二维dp数组进行遍历，利用上面的规则重置每个值，最后输出最大值就是最长公共子串。

    其实认真观察上面的结论和程序，会知道其实根本就不需要一个dp二维矩阵也可以完成最长公共子串的寻找，因为最后只需要一个最大公共子串的长度和它的坐标，然后就能得出结果，并且每一个dp中的元素的值只与它的左上方值有关，所以可以简化函数，让它的空间复杂度为O（1）。代码如下：

#include <iostream>
using namespace std;
int main()
{
    string a="jahrgiuaerhgiaeg";
    string b="hsiresahrgiuahgise";
    int l1=a.size();
    int l2=b.size();
    if(l1==0||l2==0)
    {
        cout<<"";
        return 0;
    }
    int row = 0, col = l2-1; 
    int max_l= 0,max_i = 0; //记录最大长度和a字符串的某个截至位置
    while(row < l1){
        int i = row, j = col;
        int num = 0;
        while(i < l1 && j < l2){ //因为是像右下方遍历，所以只会增加，只用小于最大值就ok
            if(a[i] == b[j])
                ++num;
            else
                num = 0;
            if(num > max_l){ //更新最大值
                max_l =num;
                max_i=i;
            }
            i++;
            j++;
        }
        if(col > 0) // 没到最左往左移动
            col--;
        else
            row++;  //到达最左往下移动
    }
    string res;
    while(max_l--)
    {
    res=a[max_i--]+res;
    }
    cout<<res;
    system("pasue");
  
}

输出结果为：

ahrgiua ，可以看出和上面一种方法得出的结果是一样的。