最长公共子序列

      比较经典的算法问题。

      不再赘述，dp[i][j] = max(dp[i-1][j-1] + (A[i]==B[j]?1:0),  dp[i-1][j], dp[i][j-1])；

      下面的cpp程序是返回最长的公共子序列而不是序列长度。

      

vector<int> LCS(vector<int> &A,vector<int> &B)
{
	int m=A.size(), n=B.size();
	//dp[i][[j]为子序列A_i、B_j的最长公共子序列
	vector<vector<vector<int> > > dp(m,vector<vector<int> >());
	for(int i=0;i<m;i++) dp[i].insert(dp[i].begin(),n,vector<int>());
	for(int j=0;j<n;j++){
		if(A[0]==B[j]) dp[0][j].push_back(A[0]);
	}
	for(int i=1;i<m;i++){
		if(A[i]==B[0]) dp[i][0].push_back(B[0]);
	}
	for(int i=1;i<m;i++){
		for(int j=1;j<n;j++){
			int L1 = dp[i-1][j].size(),
				L2 = dp[i-1][j-1].size(),
				L3 = dp[i][j-1].size();
			if(A[i]==B[j]) L2++;
			int L = max(max(L1,L2),L3);
			if(L==L1) dp[i][j] = dp[i-1][j];
			else if(L==L2){ 
				dp[i][j] = dp[i-1][j-1];
				if(A[i]==B[j]) dp[i][j].push_back(A[i]);
			}
			else dp[i][j] = dp[i][j-1];
		}
	}
	return dp[m-1][n-1];
}