求最长公共子序列（DP算法）

思路：

s1 = <x1,x2,x3...xn>

s2 = <y1,y2,y3...ym>

若它们的最长公共子序列为 s3 = <g1,g2,g3...gk>

则：

若xn == ym:

gk为xn和ym的最长公共子序列LCS的最后一个元素&&g(k-1)为xn-1和ym-1的LCS的元素

此时：lcs[xn][ym] = lcs[xn-1][ym-1] + 1.(即在len(zk-1)的长度的基础上长度 + 1)

若xn != ym:

则gk为（xn,ym-1）或 (xn-1,ym)的最长子序列的一个元素

此时：lcs[xn][ym] = max(len(lcs[xn][ym-1]) , len(lcs[xn-1][ym]) )

综上利用动态规划方程：

#include<iostream>
#include<string>
#include<cstring>
using namespace std;
string s1,s2;
int dp[1001][1001];
void Initial(){
	cin>>s1>>s2;
	memset(dp,0,sizeof(dp));
}
void DP(){
	for(int i = 0;i < s1.size();i++)
	for(int j = 0;j < s2.size();j++){
		if(s1[i] == s2[j]){
			if(i == 0 || j == 0)dp[i][j] = 1;
			else dp[i][j] = dp[i - 1][j - 1] + 1;
		}
		else dp[i][j] = max(dp[i - 1][j],dp[i][j - 1]);
	}
}
void Show(){
	cout<<"the longest num of substr is: "<<dp[s1.size() - 1][s2.size() - 1];
}
int main(){
	Initial();
	DP();
	Show();
}