【动态规划入门】最长公共子序列

一、题目描述

给定两个字符串s1,s2，返回这两个字符串的最长公共子序列的长度。如果不存在公共子序列，则返回0.

示例如下：

输入：
abcde
ace
输出：3

 

二、思路

  按照动态规划的解题步骤：

1. 确定dp数组的含义
   dp[i][j]:：长度为[0, i - 1]的字符串s1与长度为[0, j - 1]的字符串s2的最长公共子序列为dp[i][j]
2. 确定递推公式
   主要就是两大情况：s1[i - 1] 与 s2[j - 1]相同，s1[i - 1] 与 s2[j - 1]不相同
   如果s1[i - 1] 与 s2[j - 1]相同，那么找到了一个公共元素，所以dp[i][j] = dp[i - 1][j - 1] + 1;
   如果s1[i - 1] 与 s2[j - 1]不相同，那就看看s1[0, i - 2]与s2[0, j - 1]的最长公共子序列 和s1[0, i - 1]与s2[0, j - 2]的最长公共子序列，取最大的。
   即：dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
3. 数组初始化
   由于s1[0,i-1]和空串的最长公共子序列长度为0，因此dp[i][0]为0，同理，dp[0][j]也为0
4. 确定遍历顺序
   因为dp[i][j]可以从dp[i-1][j]，dp[i][j-1]以及dp[i-1][j-1]推出，因此有三个方向可以推出dp[i][j]，因此遍历顺序是从前到后，从上到下。

 

三、C++代码

#include<bits/stdc++.h>
using namespace std;

//最长公共子序列 
#define maxn 10010
int dp[maxn][maxn];   //dp[i][j]表示以下标i-1结尾的字符串s1和以下标j-1结尾的字符串s2的最长重复子数组的长度 


int main(){
	
	 
	 string s1,s2;
	 cin>>s1>>s2;
	 
	 int N=s1.size();
	 int M=s2.size();
	 
	 //dp数组初始化
	 for(int i=0;i<N;i++){
			dp[i][0]=0;
	 } 
	 for(int j=0;j<M;j++){
			dp[0][j]=0;
	 }
	 
	 //确定递推关系
	 for(int i=1;i<=N;i++){
			for(int j=1;j<=M;j++){
				if(s1[i-1]==s2[j-1]){
					dp[i][j]=dp[i-1][j-1]+1;
				}else{
					dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
				}
			}
	 } 
	 
	 cout<<dp[N][M];
	 
	 
	  
}