动态规划模型——最长公共子序列（LCS）

给定两个序列，求一个最长的子序列（子序列可以不连续）

//最长公共子序列做法
#include<cstdio>
#include<vector>

using namespace std;

int LCS(vector<int>& A,vector<int>& B){//A为模板串,A,B下标从1开始
	vector<vector<int> > dp(A.size(),vector<int>(B.size(),0));//dp[i][j] 表示A[i]和B[j]之前的LCS长
	for(int i=1;i<A.size();i++){
		for(int j=1;j<B.size();j++){
			if(A[i]==B[j]){
				dp[i][j]=dp[i-1][j-1]+1;
			}else{
				if(dp[i-1][j]>dp[i][j-1])
					dp[i][j]=dp[i-1][j];
				else
					dp[i][j]=dp[i][j-1];
			}
		}
	}
	return dp[A.size()-1][B.size()-1];
}

原模型是要求子序列和模板序列一一对应的。如果允许子序列中出现重复元素，需要对原状态转移方程做修改。原本要求一一对应的话A[i]与B[j]相同时的子问题是继承与dp[i-1][j-1]。但是B[j]或者A[i]可能是dp[i-1][j]或dp[i][j-1]中子序列的重复元素。故也需要像不相等时一样继承dp中的较大值再加1。

#include<cstdio>
#include<vector>
#include<algorithm>
using namespace std;
int LCS(vector<int>& A,vector<int>& B){//A为模板串,A,B下标从1开始
	vector<vector<int> > dp(A.size(),vector<int>(B.size(),0));//dp[i][j] 表示A[i]和B[j]之前的LCS长
	for(int i=1;i<A.size();i++){
		for(int j=1;j<B.size();j++){
			if(A[i]==B[j]){
				dp[i][j]=max(dp[i-1][j],dp[i][j-1])+1;
			}else{
				dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
			}
		}
	}
	return dp[A.size()-1][B.size()-1];
}