问题描述:
输入格式:
两行,每行一个字符串,分别表示一个DNA序列(每个字符串长度不超过1000)。
输出格式:
一个数,最长公共子序列元素的个数。
输入样例:
在这里给出一组输入。例如:
AGCT
ATT
输出样例:
在这里给出相应的输出。例如:
2
解析:首先要明白子序列的概念是什么,就是一个序列去掉0个或者多个元素后得到的结果,这里需要强调的就是,子序列不一定是连续的,有可能是不连续的,像题目中的最长公共子序列就是AT。想要用动态规划来解题,首先要找到递推关系式,然后根据地递推关系式自底向上的进行编程
这个就是该问题的递推关系式,其中C[i,j]表示长度为i和长度为j的两个序列的最长公共子序列的长度,它依赖于比它还要小的子问题。
有了递推关系式,就可以编程了
需要注意的地方已经注释了
#include<iostream>
using namespace std;
#define MAX 1002
string a,b;
int dp[MAX][MAX]={0};
int main(){
cin>>a>>b;
for(int i=1;i<=a.length();i++){//需要注意的是,这里的i和j都是从1开始的,因为最小的序列就是长度为1
for(int j=1;j<=b.length();j++){
if(a[i-1]==b[j-1])//这里i-1和j-1是比较两个字符串的第一个元素,之后是第二个...一直到a.length()-1个元素
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];
}
}
}
cout<<dp[a.length()][b.length()];
return 0;
}
}