2021.12.10 P2516 [HAOI2010]最长公共子序列(动态规划+滚动数组)
2021.12.10 P2516 [HAOI2010]最长公共子序列(动态规划+滚动数组)
https://www.luogu.com.cn/problem/P2516
题意:
给定字符串 \(S\) 、 \(T\) ,都以 \(.\) 结尾,求 \(S\) 、 \(T\) 最长公共子序列的长度及个数。
分析:
一顿操作猛如虎,一看分数250……爆零了。原本就没准备拿几分,结果令人心塞。
第一问就是求最长公共子序列长度,数据范围比较小, \(O(n^2)\) 就行,上来就是一顿树状数组+LIS,忽视了里面有重复出现的字符。正解就是普普通通的DP,啥也没加,因为也加不了啥。
第二问实际上是第一问的深入思考(大雾)。
先来看求LCS的状态转移方程:
f[i-1][j-1]+1,&S_i==T_j\\
\max(f[i][j-1],f[i-1][j]),&S_i!=T_j
\end{cases}
\]
对于第一个转移方程,相当于是 \(S\) 的前 \(i\) 个字符全部出现在LCS中, \(T\) 的前 \(j\) 个字符同样全部出现在LCS中(情况一);
对于第二个转移方程,分两种情况:
情况a: \(S\) 的前 \(i\) 个字符全部出现在LCS中, \(T\) 的前 \(j-1\) 个字符全部出现在LCS中(情况二);
情况b: \(S\) 的前 \(i-1\) 个字符全部出现在LCS中, \(T\) 的前 \(j\) 个字符全部出现在LCS中(情况三);
但是还有特殊情况,即 \(S\) 的前 \(i-1\) 个字符全部出现在LCS中, \(T\) 的前 \(j-1\) 个字符同样全部出现在LCS中(情况四)。
\(g[i][j]\) 为求方案的数组。
g[i-1][j-1],&情况一\\
g[i][j-1],&情况二\\
g[i-1][j],&情况三\\
-g[i-1][j-1],&情况四
\end{cases}
\]
情况四是因为当在计算 \(f[i][j]\) 时,已经被 \(f[i][j-1]\) (按照情况二)和 \(f[i-1][j]\) (按照情况三)计算两遍了,所以减去。
至于其他情况为什么是加号而不是等号(我脑子抽抽,第一次写成等号,于是我挂了,shift!),那是因为这个位置会经过四重情况层层检验计算(大雾)。
代码如下:
100pts:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
#define IOS ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
using namespace std;
const int N=5e3+10;
const int mod=1e8;
int n,m,a[N],b[N],f[3][N],g[3][N];
int main(){
string s;cin>>s;
n=s.length()-1;
for(int i=1;i<=n;i++)a[i]=s[i-1]-'A'+1,g[0][i]=1;
cin>>s;
m=s.length()-1;
for(int i=1;i<=m;i++)b[i]=s[i-1]-'A'+1,g[0][i]=1;
g[1][0]=g[0][0]=1;
int k=1;
for(int i=1;i<=n;i++,k^=1){
for(int j=1;j<=m;j++){
g[k][j]=0;
if(a[i]==b[j])f[k][j]=f[k^1][j-1]+1;
else f[k][j]=max(f[k^1][j],f[k][j-1]);
if(a[i]==b[j])g[k][j]+=g[k^1][j-1];
if(f[k][j]==f[k][j-1])g[k][j]+=g[k][j-1];
if(f[k][j]==f[k^1][j])g[k][j]+=g[k^1][j];
if(f[k][j]==f[k^1][j-1])g[k][j]-=g[k^1][j-1];
g[k][j]=(g[k][j]%mod+mod)%mod;
}
}
cout<<f[k^1][m]<<endl<<g[k^1][m];
/*cout<<endl;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++)cout<<g[i][j]<<" ";
cout<<endl;
}*/
return 0;
}
0pts:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
#define IOS ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
using namespace std;
const int N=5e3+10;
int n,top,a[N],t[N],cnt[N];
int cap[30][N];
inline int lowbit(int x){
return x&-x;
}
inline void add(int x,int k){
for(int i=x;i<=n;i+=lowbit(i))t[i]=max(t[i],k);
}
inline int query(int x){
int fin=0;
for(int i=x;i>0;i-=lowbit(i))fin=max(fin,t[i]);
return fin;
}
inline void find(){
for(int i=1;i<=10;i++)cout<<query(i)<<" ";cout<<endl;
}
int main(){
string s;
cin>>s;
for(int i=1;i<s.length();i++){
int x=s[i-1]-'A'+1;
++cap[x][0];
cap[x][cap[x][0]]=++top;
}
for(int i=1;i<=26;i++)cap[i][0]=0;
cin>>s;
n=s.length()-1;
int maxn=0,num=0;
for(int i=1;i<s.length();i++){
int x=s[i-1]-'A'+1;
++cap[x][0];
if(!cap[x][cap[x][0]])cap[x][cap[x][0]]=++top;
a[i]=cap[x][cap[x][0]];
int len=query(a[i]-1)+1;
if(len==maxn)++num;
else if(len>maxn)maxn=len,num=1;
add(a[i],len);
//cout<<i<<" "<<a[i]<<" "<<len<<" "<<query(a[i]-1)<<endl;
//find();//
}
//for(int i=1;i<s.length();i++)cout<<a[i]<<" ";cout<<endl;
cout<<maxn<<endl<<num;
return 0;
}
2021.12.10 P2516 [HAOI2010]最长公共子序列(动态规划+滚动数组)的更多相关文章
- 洛谷P2516 [HAOI2010]最长公共子序列(LCS,最短路)
洛谷题目传送门 一进来就看到一个多月前秒了此题的ysn和YCB%%% 最长公共子序列的\(O(n^2)\)的求解,Dalao们想必都很熟悉了吧!不过蒟蒻突然发现,用网格图貌似可以很轻松地理解这个东东? ...
- 洛谷 P2516 [HAOI2010]最长公共子序列
题目传送门 解题思路: 第一问要求最长公共子序列,直接套模板就好了. 第二问要求数量,ans[i][j]表示第一个字符串前i个字符,第二个字符串前j个字符的最长公共子序列的数量 如果f[i][j]是由 ...
- P2516 [HAOI2010]最长公共子序列 题解(LCS)
题目链接 最长公共子序列 解题思路 第一思路: 1.用\(length[i][j]\)表示\(a\)串的前\(i\)个字符与\(b\)串的前\(j\)个字符重叠的最长子串长度 2.用\(num[i][ ...
- POJ 1159 Palindrome-最长公共子序列问题+滚动数组(dp数组的重复利用)(结合奇偶性)
Description A palindrome is a symmetrical string, that is, a string read identically from left to ri ...
- poj1159--Palindrome(dp:最长公共子序列变形 + 滚动数组)
Palindrome Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 53414 Accepted: 18449 Desc ...
- luogu P2516 [HAOI2010]最长公共子序列
传送门 首先那个\(O(n^2)\)的dp都会吧,不会自己找博客或者问别人,或是去做模板题(误) 对以下内容不理解的,强势推荐flash的博客 我们除了原来记录最长上升子序列的\(f_{i,j}\), ...
- 洛谷P2516 [HAOI2010]最长公共子序列
题目描述 字符序列的子序列是指从给定字符序列中随意地(不一定连续)去掉若干个字符(可能一个也不去掉)后所形成的字符序列.令给定的字符序列X="x0,x1,-,xm-1",序列Y=& ...
- Luogu P2516 [HAOI2010]最长公共子序列 DP
首先$LIS$显然:$f[i][j]=max(f[i][j-1],f[i-1][j],(a[i]==b[j])*f[i-1][j-1])$ 考虑如何转移数量: 首先,不管$a[i]$是否等于$b[j] ...
- P2516 [HAOI2010]最长公共子序列
传送门 看到数据范围,显然 $n^2$ 的 $dp$... 设 $f[i][j]$ 表示 $A$ 串考虑了前 $i$ 位,$B$ 串考虑了前 $j$ 位,最优情况下的方案数 但是好像没法判断转移来的是 ...
随机推荐
- python神器 Jupyter Notbook
python神器 Jupyter Notbook 简介 Jupyter Notebook是基于网页的用于交互计算的应用程序.其可被应用于全过程计算:开发.文档编写.运行代码和展示结果. Jupyter ...
- 6月7日 python 复习 collections
collections是Python内建的一个集合模块,提供了许多有用的集合类. namedtuple 我们知道tuple可以表示不变集合,例如,一个点的二维坐标就可以表示成: >>> ...
- Rsync未授权访问
1.漏洞名称 Rsync 未授权访问漏洞 2.漏洞原理 rsync是Linux下一款数据备份工具,支持通过rsync协议.ssh协议进行远程文件传输. 其中rsync协议默认监听873端口,如果目标开 ...
- Flutter入门教程(二)开发环境搭建
学习Flutter,首先需要搭建好Flutter的开发环境,下面我将一步步带领大家搭建开发环境并且成功运行flutter项目. Flutter环境配置主要有这几点: 系统配置要求 Java环境 Flu ...
- HashMap?ConcurrentHashMap?
Java团长在csdn 2019-02-06 19:58:00 104241 收藏 666 来源:crossoverJie(ID:crossoverJie) 前言 Map 这样的 Key ...
- Zookeeper 对于 Kafka 的作用是什么?
Zookeeper 是一个开放源码的.高性能的协调服务,它用于 Kafka 的分布式应用. Zookeeper 主要用于在集群中不同节点之间进行通信 在 Kafka 中,它被用于提交偏移量,因此如果节 ...
- spring集成mongodb简单使用和测试方式
@EnableMongoRepositories @ComponentScan(basePackages = "cn.example") @Configuration public ...
- 四种类型的数据节点 Znode?
1.PERSISTENT-持久节点 除非手动删除,否则节点一直存在于 Zookeeper 上 2.EPHEMERAL-临时节点 临时节点的生命周期与客户端会话绑定,一旦客户端会话失效(客户端与 zoo ...
- Elasticsearch 对于大数据量(上亿量级)的聚合如何实现?
Elasticsearch 提供的首个近似聚合是 cardinality 度量.它提供一个字段的基数, 即该字段的 distinct 或者 unique 值的数目.它是基于 HLL 算法的.HLL 会 ...
- 为什么HTTP/3要基于UDP?可靠吗?
目录 前言 为什么转用UDP? HTTP/3解决了那些问题? 队头阻塞问题 QPACK编码 Header 参考 推荐阅读: 计算机网络汇总 HTTP/3竟然是基于UDP的!开始我也很疑惑,UDP传输不 ...