LCS-最长公共子序列
有时候我们会遇到要求两个字符串的公共最长序列,这就是LCS问题,它通常是用动态规划来求解.
解析
假设现在有两个字符串A, B,字符串A的长度为
n
n
n, 字符串B的长度为
m
m
m, 既然是
d
p
dp
dp, 那么就设
d
p
[
i
]
[
j
]
dp[i][j]
dp[i][j]表示字符串从0的到
i
i
i位置和字符串从0到j位置的最长公共子序列, 那么就可以讨论
A
[
i
]
=
=
B
[
j
]
,
d
p
[
i
]
[
j
]
=
d
p
[
i
−
1
]
[
j
−
1
]
+
1
A[i]==B[j] , dp[i][j]=dp[i-1][j-1]+1
A[i]==B[j],dp[i][j]=dp[i−1][j−1]+1
A
[
i
]
!
=
B
[
j
]
,
d
p
[
i
]
[
j
]
=
m
a
x
(
d
p
[
i
−
1
]
[
j
]
,
d
p
[
i
]
[
j
−
1
]
)
A[i]!=B[j], dp[i][j]=max(dp[i-1][j],dp[i][j-1])
A[i]!=B[j],dp[i][j]=max(dp[i−1][j],dp[i][j−1])
这个式子还是比较好理解的,当两个位置的字符相等时,他们的最长公共子序列长度等于前一个位置+1,这个是显然的。当两个位置的字符不相等时,他们的公共子序列的长度就为两个前面一个的最大值。在计算过程中需要枚举字符串A的每一位去和B的每一位比较,要单独处理0位置的值。
说了这么多,看道题吧,传送门:
https://www.nowcoder.com/practice/b43fb39898f448e39adbcffde5ff0dfc?tpId=90&tqId=30810&tPage=2&rp=2&ru=/ta/2018test&qru=/ta/2018test/question-ranking
这道题的意思就是把字符串分成两部分,然后在分别求最长的公共子序列,那么就可以考虑把字符串分成两部分在分别 d p dp dp, 可以先预处理下,就是把只出现一次的字符去掉,因为这道题规模太小了,不这样做照样可以A。
AC代码
// 小学生一发的刷题之路
// 3月3号-平方串;
//
//
//
#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#include <queue>
#include <deque> //双向队列;
#include <cmath>
#include <set>
#include <stack>
#include <map>
#include <vector>
#include <cstdlib>
#include <iomanip>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
const double PI=acos(-1.0);
const double eps=1e-10;
const int maxn=1e5+5;
const int maxm=1e3+5;
const ll mod=1e9+7;
const int INF=1e8;
template<class T>
inline void read(T &ret){ //快速输入模版;
ret=0;
int f=1;
char c=getchar();
while(c<'0'||c>'9'){
if(c=='-') f=-1;
c=getchar();
}
while(c>='0'&&c<='9'){
ret=ret*10+c-'0';
c=getchar();
}
ret*=f;
}
template <class T>
inline void out(T ret){ //快速输出模版;
if(ret>9)
{
out(ret/10);
}
putchar(ret%10+'0');
}
string str;
int dp[100][100];
int sovle(int m,int n){
for(int i=0;i<=m;i++)
for(int j=m+1;j<n;j++){
if(str[i]==str[j]){
if(i==0||j==m+1){
dp[i][j]=1;
}else{
dp[i][j]=dp[i-1][j-1]+1;
}
}else{
if(i==0&j==m+1){
dp[i][j]=0;
}else if(i==0){
dp[i][j]=dp[i][j-1];
}else if(j==m+1){
dp[i][j]=dp[i-1][j];
}else{
dp[i][j]=max(dp[i][j-1],dp[i-1][j]);
}
}
}
return dp[m][n-1];
}
int main()
{
cin>>str;
int n=str.length();
int ans=0;
for(int i=0;i<n;i++){
int tmp=sovle(i,n);
ans=max(ans,tmp);
}
printf("%d\n",2*ans);
return 0;
}
新的开始,每天都要快乐哈!