最长公共子序列问题(LCS)

C++代码：

//LCS问题（动态规划）
string CommonSubsquence(string s1, string s2)
{
	int size1 = s1.size();//一维
	int size2 = s2.size();//二维

	size1++;
	size2++;

	//开辟数组
	int ** p = new int*[size2];
	p[0] = new int [size1*size2];
	memset(p[0],0,sizeof(int)*size1*size2);
	for(int i=1; i<=size2-1; i++)
		p[i] = p[i-1] + size1;

	//dp过程
	for(int j=1; j<size2; j++)
	{
		for(int i=1; i<size1; i++)
		{
			if(s1[i-1] == s2[j-1])
			{
				p[j][i] = p[j-1][i-1]+1;
			}
			else
			{
				p[j][i] = max(p[j-1][i],p[j][i-1]);
			}
		}
	}
	//输出结果
	for(int i=0; i<size1; i++)
	{
		for(int j=0; j<size2; j++)
		{
			cout<<p[j][i]<<' ';
		}
		cout<<endl;
	}

	//输出字符串
	string s;
	for(int i=size1-1,j=size2-1; i>0&&j>0;)
	{
		if(s1[i-1]==s2[j-1])
		{
			s.push_back(s1[i-1]);
			i--;
			j--;
		}
		else
		{
			//if(p[j][i-1]>=p[j-1][i]) "BDAB"
			if(p[j][i-1]>=p[j-1][i]) //"BCBA"
				i--;
			else
				j--;
		}
	}
	reverse(s.begin(),s.end());
	//释放数组
	delete[] p[0];
	delete p;

	return s;
}