最长公共子串

问题描述：已知strA和strB是两个字符串，例如 strA = "a12345a54312323b", strB ="b32321345a54321a"; 现在需要求得strA和strB的最长公共子串，这里子串是连续的，区别于子序列，例子中strA和strB的最长公共子串为"345a543"

  解体思路：

  假设已知以strA[i]和strB[j]为结尾字符的最长公共子串的长度为lcs[i,j],那么以strA[i+1]和strB[j+1]为结尾字符的最长公共子串的长度有两种情况：

  1.lcs[i+1,j+1] = 0， 此时strA[i+1] != strB[j+1]
  2.lcs[i+1,j+1] = lcs[i,j] + 1; 此时strA[i] = strB[j]
  由上面两个公式可知，可以用动态规划来解决这个问题
  
  由于只涉及到[i+1,j+1]和[i,j]，因此不用使用二维数组，只需一维数组就可以解决
  此时就有
  lcs[j+1] = lcs[j] + 1, strA[i] = strB[j]
  lcs[j+1] = 0,          strA[i] != strB[j]
  这里考虑边界条件lcs[0]会跟lcs[-1]相关，但这会出错，因此将lcs[j] = lcs[j-1] + 1改为lcs[j+1] = lcs[j] + 1
  
  还有一点需要注意：更新lcs时需要从后往前更新，这是因为如果从前往后更新的话，当lcs[j]需要使用lcs[j-1]时，由于lcs[j-1]已经被处理过了，那么lcs[j]

  所需要的lcs[j-1]就不是之前的lcs[j-1]了，这肯定会导致结果有问题

#include <iostream>
#include <string>
using namespace std;
int LongestCommonSubstring(const string &strA, const string &strB)
{
	int *lcs;                
	int lengthA, lengthB;
	int maxLen = 0, pos;
	lengthA = strA.length();
	lengthB = strB.length();

	lcs = new int[lengthB + 1];
	memset(lcs, 0, sizeof(int)* (lengthB + 1));


	for (int i = 0; i < lengthA; ++i)
	{
		for (int j = lengthB; j > 0; --j)
		{
			if (strA[i] == strB[j-1])		
			{
				lcs[j] = lcs[j - 1] + 1;
			}
			else
			{
				lcs[j] = 0;
			}
			
			if (lcs[j] > maxLen)		//求最长公共子序列
			{
				maxLen = lcs[j];
				pos = j - 1;			//记录最长公共子序列的尾元素在strB中的位置
			}
		}
		for (int j = 1; j <= lengthB; ++j)
			cout << lcs[j] << " ";
		cout << endl;
	}
	cout << strA << "和" << strB << "的最长公共子序列的长度为:" << maxLen << endl;
	cout << strA << "和" << strB << "的最长公共子序列为:" << endl;
	for (int k = pos - maxLen + 1; k <= pos; k++)
		cout << strB[k];
	cout << endl;

	return maxLen;
}

int main()
{
	string strA("a12345a54312323b");
	string strB("b32321345a54321a");
	int maxLen = LongestCommonSubstring(strA, strB);
	return 0;
}