接着动态规划①的内容,动态规划基础的第二道题是最长公共子序列问题,通过洛谷上的模版题:最长公共子序列可以进一步试验自己的代码能力~
但是写的过程发现,洛谷的题目中有这样的一个数据说明:
这代表了什么呢?其实是它将题目复杂度分成了两个级别,例如这题有O(n²)和O(nlogn)
的解法,第一种只能过掉一半的数据,本人小白能力有限,就先介绍第一种写法,第二章nlogn的写法留个坑以后再补。。。
这个题目的O(n²)的做法也很好理解,就是两层循环,每层分别遍历一遍字符串(其实从这个地方就隐约感觉到有种重复计算的情况了,但具体怎么解决还没看懂= =),状态dp[i][j]代表了,第一个字符串第i个字符和第二个字符串第j个字符,看看他们之间的最长公共子序列长度是多少。假如遍历过程中,如果发现两个字符相等则很简单,长度直接+1(这个应该挺好理解的)。但如果两个字符不相等,那么这个状态怎么赋值呢?????
举个栗子:AEDBC和AEDC,当第一个i循环到“B”的位置,而第二个循环到“D”的位置时,我们要看它们两个各退一步的话,谁更大!如果第一个字符串退一步,即AED,第二个字符串保持不变即AED发现最长公共子序列为3。而如果第二个字符串退一步即AE,第一个字符串保持不变,即AEDB发现最长公共子序列为2.那么我们就取dp[i-1][j]即表示的是,虽然当前两个字符不能匹配上,但是多出来的那个两个字符可能可以跟之前没匹配的匹配上,取两者最大的值。所以有了状态转移方程:
最后贴出最终的O(n²)的代码:
#include<iostream>
using namespace std;
const int n=1001;
int dp[n][n];
int f[n],s[n];
int main(){
int r;
cin>>r;
for(int j=1;j<=r;j++){
cin>>f[j];
}
for(int j=1;j<=r;j++){
cin>>s[j];
}
for(int i=1;i<=r;i++){
for(int j=1;j<=r;j++){
if(f[i]==s[j]){
dp[i][j]=dp[i-1][j-1]+1;
}
else{dp[i][j]=max(dp[i-1][j],dp[i][j-1]);};
}
}
cout<<dp[r][r]<<endl;
return 0;
}