
洛谷P2516 [HAOI2010]最长公共子序列 动态规划
发布日期:2021-05-08 02:34:17
浏览次数:20
分类:原创文章
本文共 1821 字,大约阅读时间需要 6 分钟。
分析
给出两个字符A和B,求A和B最长公共子序列长度以及最长的长度个数有几个。
###1. 首先求最长公共子序列
最长公共子序列问题:
###2. 求出了最长的长度,之后要同时把长度的数目计算清楚。
用cnt[i][j]表示A前i位,B前j位的最长公共子序列数目。
对于每个f[i][j]:
1. 若f[i][j]==f[i-1][j],说明状态是从f[i-1][j]转移来的,cnt[i][j]要加上cnt[i-1][j]
2. 若f[i][j]==f[i][j-1],说明状态是从f[i][j-1]转移来的,cnt[i][j]要加上cnt[i][j-1]
3. 若a[i]==b[j]且f[i][j]==f[i-1][j-1]+1,说明状态是从f[i-1][j-1]转移来的,cnt[i][j]要加上cnt[i-1][j-1]
4. 若a[i]!=b[j]且f[i][j]==f[i-1][j-1],那么状态cnt[i-1][j-1]会多加上一次,所以要让cnt[i][j]-=cnt[i-1][j-1]
C++ 代码
#include<bits/stdc++.h>using namespace std;const int N = 5010,mod = 1e8;int f[2][N],cnt[2][N]; //为防止内存溢出MLE,使用滚动数组,因为每个状态只与上一个状态有关char s[N],t[N];int main(){ scanf("%s%s",s+1,t+1); int n=strlen(s+1)-1,m=strlen(t+1)-1; for(int i=0;i<=m;i++) //初始化,长度为0最长子序列方案数为1(第一个序列长度为0) cnt[0][i]=1; cnt[1][0]=1; //初始化,长度为0最长子序列方案数为1(第二个序列长度为0) for(int i=1;i<=n;i++) { int now=i&1,pre=now^1; //当前状态为now,前一个状态为pre for(int j=1;j<=m;j++) { f[now][j]=max(f[pre][j],f[now][j-1]); //最长公共子序列前两种情况 cnt[now][j]=0; if(s[i]==t[j]) { f[now][j]=max(f[now][j],f[pre][j-1]+1); //最长公共子序列第三种种情况 if(f[now][j]==f[pre][j-1]+1) cnt[now][j]=cnt[pre][j-1]; //如果此时满足当前序列f[now][j]等于上一序列f[pre][j-1]+1,cnt当前状态则从上一个状态转移过来 } else { if(f[now][j]==f[pre][j-1]) cnt[now][j]-=cnt[pre][j-1]; //s[i]!=t[j]并且满足当前序列f[now][j]等于上一序列f[pre][j-1],cnt要减去上一状态的数目,防止重复 } if(f[now][j]==f[pre][j]) cnt[now][j]=(cnt[now][j]+cnt[pre][j])%mod; //状态是从f[i-1][j]转移来的,cnt[i][j]要加上cnt[i-1][j] if(f[now][j]==f[now][j-1]) cnt[now][j]=(cnt[now][j]+cnt[now][j-1])%mod; //状态是从f[i][j-1]转移来的,cnt[i][j]要加上cnt[i][j-1] } } printf("%d\n",f[n&1][m]); printf("%d\n",cnt[n&1][m]); return 0;}
发表评论
最新留言
关注你微信了!
[***.104.42.241]2025年03月22日 18时07分33秒
关于作者

喝酒易醉,品茶养心,人生如梦,品茶悟道,何以解忧?唯有杜康!
-- 愿君每日到此一游!
推荐文章
图像处理系列1.skimage
2019-03-05
好用的拼图小程序,图作妖
2019-03-05
读取二进制存储信息,将低位二进制存储转换为高位存储
2019-03-05
Hibernate操作Blob,将Blob转换为String
2019-03-05
Object Clone
2019-03-05
Javascript中String支持使用正则表达式的四种方法
2019-03-05
2021年判断浏览器最新写法,你都掌握了吗?
2019-03-05
简易版Http请求工具
2019-03-05
【法律】如何保障企业的合法权益:保密协议模板
2019-03-05
【IoT】 产品设计之结构设计:PMMA(亚克力)板、PC耐力板、PS有机板与MS板区别
2019-03-05
【IoT】蓝牙BLE基础:CC254x通信系列之模拟SPI协议
2019-03-05
【IoT】TI BLE CC2541 串口控制蓝牙详解
2019-03-05
【产品】项目管理的五个过程和九大知识领域之二
2019-03-05
【项目管理】项目管理流程浅析
2019-03-05
【企业】韬盛和夫六精进
2019-03-05
【Tool】如何使用 Uniflash 烧写 WIFI 芯片 CC3200
2019-03-05
html2canvas vue页面截图生成图片地址
2019-03-05
copy_{to, from}_user()的思考
2019-03-05
Web前端安全策略之CSRF的攻击与防御
2019-03-05