一、前言
最近每一天都在恶补DP.......就连做梦都梦到别人跟我说:“你这个DP写得怕是有问题哦?” 我:“······”
但是尽管每天都在恶补DP,依旧深陷一个魔咒:看正解,海阔天空;自己写,波涛汹涌......
所以,一点点总结经典的DP模型和经典例题。/突然严肃
二、题目描述
给定字符串a,字符串b,求a和b的LCIS(最长公共上升子序列)——的长度。
(校内OJ-3183)
“求出这个序列”的问题见后_(:з」∠)_
三、分析
本题的关键字:最长+公共+上升+子序列(+长度)
可以看到,这个问题具有相当多的重叠子问题,LIS与LCS也是经典的dp问题。
——>于是,我们想到了DP!
DP核心一:定义状态——F[i][j]表示以a串的前i个字符b串的前j个字符且以b[j]为结尾构成的LCIS的长度
Ps.之前自己没有注意到“以b[j]结尾”,所以盯着方程半天没有反应.../共鸣失败,还好后来看到了,在这里重点强调一下
DP核心二:状态转移方程
我们来考察一下这个这个状态,思考这个状态能转移到哪些状态似乎有些棘手;/它转移给谁?
如果把思路逆转一下,考察这个状态的最优值依赖于哪些状态,就容易许多了。/谁转移到它?
首先,在a[i]!=b[j]的时候有F[i][j]=F[i-1][j]。
为什么呢?因为F[i][j]是以b[j]为结尾的LCIS,如果F[i][j]>0,那么就说明a[x](1<=x<=i)中,必然有一个字符a[k]等于b[j](如果F[i][j]等于0呢?那赋值与否都没有什么影响了)。
因为a[k]!=a[i],那么a[i]对F[i][j]没有贡献,于是我们不考虑它照样能得出F[i][j]的最优值。所以在a[i]!=b[j]的情况下必然有F[i][j]=F[i-1][j]。这一点参考LCS的处理方法。
那如果a[i]==b[j]呢?
首先,这个等于起码保证了长度为1的LCIS;然后我们还需要去找一个最长的且能让b[j]接在其末尾的LCIS。
之前最长的LCIS在哪呢?首先我们要去找的F数组的第一维必然是i-1。(为什么只能是i-1?①i已经被配对,无法再使用,②i-1必然比i-2要差)
第二维呢?那就需要枚举b[1]..b[j-1]了,因为你不知道这里面哪个最长且哪个小于b[j]。
这里还有一个问题,可不可能不配对呢?也就是在a[i]==b[j]的情况下,需不需要考虑第一种情况中F[i][j]=F[i-1][j]的决策呢?答案是不需要。
因为如果b[j]不和a[i]配对,那就是和之前的a[1]..a[i-1]配对(假设F[i-1][j]>0,等于0不考虑),这样必然没有和a[i]配对优越。
(为什么必然呢?因为b[j]和a[i]配对之后的转移是max(F[i-1][k])+1,而和之前的i`配对则是max(F[i`-1][k])+1,显然有F[i][j]>F[i`][j],i`>i)
于是我们得出了状态转移方程:
a[i]!=b[j]: F[i][j] = F[i-1][j]
a[i]==b[j]: F[i][j] = { max(F[i-1][k])+1 | 1<=k<=j-1&&b[j]>b[k] }
[二维DP+三层循环]
int mx=0;
for(int j=1;j<=n;j++)
for(int i=1;i<=m;i++)
{
if(b[j]!=a[i])
dp[i][j]=dp[i-1][j];
else//a[i]==b[j]
{
for(int k=1;k<j;k++)
if(dp[i][j]<dp[i-1][k]&&b[k]<b[j])
dp[i][j]=dp[i-1][k];
dp[i][j]++;
}
mx=max(mx,dp[i][j]);
}
不难看到,这是一个时间复杂度为O(n^3)的DP,离平方还有一段距离。
但是,这个算法最关键的是,如果按照一个合理的递推顺序,max(F[i-1][k])的值我们可以在之前访问F[i][k]的时候,通过维护更新一个max变量得到。
怎么得到呢?首先递推的顺序必须是状态的第一维在外层循环,第二维在内层循环。也就是算好了F[1][len(b)]再去算F[2][1]。
如果按照这个递推顺序,我们可以在每次外层循环的开始加上一个为0的max变量,然后开始内层循环。
当a[i]>b[j]的时候令max=F[i-1][j]。如果循环到了a[i]==b[j]的时候,则令F[i][j]=max+1。
Ps.自己之前有个疑问想了很久:为什么是“当a[i]>b[j]的时候令max=F[i-1][j]”,为什么不是a[j]<b[i]的时候咧?
后来问了某可爱的小伙伴才 大概 搞懂——
在刚刚讨论的 “a[i]==b[j]”一部分里面,
else
{
for(int k=1;k<j;k++)
if(dp[i][j]<dp[i-1][k]&&b[k]<b[j])
dp[i][j]=dp[i-1][k];
dp[i][j]++;
}
最后答案是F[len(a)][1]..F[len(a)][len(b)]的最大值。
[二维DP+两层循环]
for(int i=1;i<=m;i++){
int mx=0;//dp[i-1][1..n]中b[j]<a[i]的max,即a[i]可以接在后面的max
for(int j=1;j<=n;j++){
dp[i][j]=dp[i-1][j];
if(b[j]<a[i])mx=max(mx,dp[i-1][j]);
else if(b[j]==a[i])dp[i][j]=mx+1;
}
}
int mx=0;
for(int j=1;j<=n;j++)mx=max(mx,dp[m][j]);
算法的时间复杂度O(n^2)已经不能再优化,但是空间复杂度可以继续优化成一维。
可以看到,上式中,由于"dp[i][j]=dp[i-1][j]"的操作,条件语句中的dp值都是一样的,这个操作也只是将上一次的结果向下转移罢了。
下面是本题O(n^2)+一维的完整做法:
[一维DP+两层循环]
#include<cstdio>
#define M 5005
int a[M],b[M];
int dp[M];//dp[i][j],a{1~i},b{1~j}且以b[j]结尾的LCIS最大长度
int max(int a,int b){return a>b?a:b;}
int main(){
int m,n;
scanf("%d",&m);for(int i=1;i<=m;i++)scanf("%d",&a[i]);
scanf("%d",&n);for(int j=1;j<=n;j++)scanf("%d",&b[j]);
for(int i=1;i<=m;i++){
int mx=0;//dp[i-1][1..n]中b[j]<a[i]的max,即a[i]可以接在后面的max
for(int j=1;j<=n;j++)
if(b[j]<a[i])mx=max(mx,dp[j]);
else if(b[j]==a[i])dp[j]=mx+1;
}
int mx=0;
for(int j=1;j<=n;j++)mx=max(mx,dp[j]);
printf("%d\n",mx);
return 0;
}
最后附上一个 超级水 样例检验一下qwq~~
Sample Input
5
1 4 2 5 -12
4
-12 1 2 4
Sample Output
2
以上资料/博客参考:https://blog.csdn.net/Kanosword/article/details/51811084
四、拓展
以上的讨论只是针对如何求出LCIS的长度,但是总有那么一些 奇奇怪怪的 出题者微笑着问你——
你会求出 这个序列吗?不会啊?等WA吧 :)
所以接下来要 被迫 讨论一下如何求出 LCIS 这个序列!!
——To be continue...