最长公共子序列（LCS）

最长公共子序列(Longest Common Sebsequence,LCS)的问题描述为：

给定两个字符串（或数字序列）A和B， 求一个字符串，使得这个字符串是A和B的最长公共部分（子序列可以不连续）

动态规划的做法

令dp[i][j]表示字符串A的i号位和字符串B的j号位之前的LCS长度，可以分为两种情况（下标从1开始）

（1）若A[i]=B[j]，则dp[i][j]=dp[i-1][j-1]+1;

(2)若不等于，则dp[i][j]=max(dp[i-1][j],dp[i][j-1])

边界上的状态

dp[i][0]=dp[0][j]=0(0<=i<=n,0<=j<=m)其中，n为a的长度，m为b的长度

这样，dp[i][j]只与之前的状态有关，由边界出发就可以得到整个dp数组，最终dp[n][m]就是问题需要的答案，时间复杂度为O(mn)

代码如下：

#include<bits/stdc++.h>
using namespace std;
const int N=100;
char A[N],B[N];
int dp[N][N];
int main(){
	int n;
	gets(A+1);//从下标为1开始读入
	gets(B+1);
	int lenA=strlen(A+1);//由于读入时下标从1开始，因此读取长度也从+1开始
	int lenB=strlen(B+1);
	//边界
	for(int i=0;i<=lenA;i++)
		dp[i][0]=0;
	for(int j=0;j<=lenB;j++)
		dp[0][j]=0; 
	//状态转移方程
	for(int i=1;i<=lenA;i++){
		for(int j=1;j<=lenB;j++){
			if(A[i]==B[j]){
				dp[i][j]=dp[i-1][j-1]+1;
			}else{
				dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
			}
		}
	} 
	//dp[lenA][lenB]是答案
	printf("%d\n",dp[lenA][lenB]);
	return 0; 
}