bzoj千题计划241:bzoj3864: Hero meet devil
http://www.lydsy.com/JudgeOnline/problem.php?id=3864
题意:
给你一个DNA序列,求有多少个长度为m的DNA序列和给定序列的LCS为0,1,2....
求LCS方式:f[i][j]=max(f[i-1][j],f[i][j-1],f[i-1][j-1]*(s[i]==t[j]))
固定了i,相邻的j的f[i][j]值最多相差1
dp[i][j] 表示长度为i的DNA序列,将“f[ |S| ][j+1]是否比f[ |S| ][j] 大1” 这个状态压缩为j的方案数
若我们知道 状态j加上一个字母k可以到状态nxt[j][k]
那么dp[i+1][nxt[j][k]]+=dp[i][j]
关键是如何求得nxt[j][k]
再一次DP
枚举所有的状态i
令f[j] 表示加上字母k之前的LCS长度,g[j]表示加上字母k之后的LCS长度
g[j]=max(g[j-1],f[j])
如果加上的字母k和原序列第j个字母匹配 g[i]=max(g[j],f[j-1]+1)
g求完后,项邻的两个g要么相等,要么相差1
再把这个状态压缩起来即可
#include<cstdio>
#include<cstring>
#include<algorithm> using namespace std; const int mod=1e9+; int m; char ss[];
int L,S;
int s[]; int ch[]; int f[],g[];
int nxt[<<][]; int dp[][<<];
int ans[]; void pre()
{
int len; int c[];
for(int i=;i<S;++i)
{
memset(f,,sizeof(f));
for(int j=;j<=L;++j) f[j]=f[j-]+(i>>j-&);
for(int k=;k<;++k)
{
for(int j=;j<=L;++j)
{
g[j]=max(g[j-],f[j]);
if(s[j]==k) g[j]=max(g[j],f[j-]+);
}
nxt[i][k]=;
for(int j=;j<L;++j)
if(g[j+]-g[j]) nxt[i][k]+=<<j;
}
}
} int count(int x)
{
int sum=;
while(x)
{
sum+=x&;
x>>=;
}
return sum;
} void DP()
{
memset(dp,,sizeof(dp));
int now=,last=;
dp[][]=;
for(int i=;i<=m;++i)
{
memset(dp[now],,sizeof(dp[now]));
for(int j=;j<S;++j)
for(int k=;k<;++k)
{
dp[now][nxt[j][k]]+=dp[last][j];
dp[now][nxt[j][k]]-=dp[now][nxt[j][k]]>=mod ? mod : ;
}
swap(now,last);
}
memset(ans,,sizeof(ans));
int t;
for(int i=;i<S;++i)
{
t=count(i);
ans[t]+=dp[last][i];
ans[t]-=ans[t]>=mod ? mod : ;
}
for(int i=;i<=L;++i) printf("%d\n",ans[i]);
} int main()
{
ch['A'-'A']=;
ch['C'-'A']=;
ch['G'-'A']=;
ch['T'-'A']=;
int T;
scanf("%d",&T);
while(T--)
{
scanf("%s",ss+);
scanf("%d",&m);
L=strlen(ss+);
S=<<L;
for(int i=;i<=L;++i) s[i]=ch[ss[i]-'A'];
pre();
DP();
}
return ;
}
bzoj千题计划241:bzoj3864: Hero meet devil的更多相关文章
- bzoj千题计划300:bzoj4823: [Cqoi2017]老C的方块
http://www.lydsy.com/JudgeOnline/problem.php?id=4823 讨厌的形状就是四联通图 且左右各连一个方块 那么破坏所有满足条件的四联通就好了 按上图方式染色 ...
- bzoj千题计划196:bzoj4826: [Hnoi2017]影魔
http://www.lydsy.com/JudgeOnline/problem.php?id=4826 吐槽一下bzoj这道题的排版是真丑... 我还是粘洛谷的题面吧... 提供p1的攻击力:i,j ...
- bzoj千题计划280:bzoj4592: [Shoi2015]脑洞治疗仪
http://www.lydsy.com/JudgeOnline/problem.php?id=4592 注意操作1 先挖再补,就是补的范围可以包含挖的范围 SHOI2015 的题 略水啊(逃) #i ...
- bzoj千题计划177:bzoj1858: [Scoi2010]序列操作
http://www.lydsy.com/JudgeOnline/problem.php?id=1858 2018 自己写的第1题,一遍过 ^_^ 元旦快乐 #include<cstdio> ...
- bzoj千题计划317:bzoj4650: [Noi2016]优秀的拆分(后缀数组+差分)
https://www.lydsy.com/JudgeOnline/problem.php?id=4650 如果能够预处理出 suf[i] 以i结尾的形式为AA的子串个数 pre[i] 以i开头的形式 ...
- bzoj千题计划304:bzoj3676: [Apio2014]回文串(回文自动机)
https://www.lydsy.com/JudgeOnline/problem.php?id=3676 回文自动机模板题 4年前的APIO如今竟沦为模板,,,╮(╯▽╰)╭,唉 #include& ...
- bzoj千题计划292:bzoj2244: [SDOI2011]拦截导弹
http://www.lydsy.com/JudgeOnline/problem.php?id=2244 每枚导弹成功拦截的概率 = 包含它的最长上升子序列个数/最长上升子序列总个数 pre_len ...
- bzoj千题计划278:bzoj4590: [Shoi2015]自动刷题机
http://www.lydsy.com/JudgeOnline/problem.php?id=4590 二分 这么道水题 没long long WA了两发,没判-1WA了一发,二分写错WA了一发 最 ...
- bzoj千题计划250:bzoj3670: [Noi2014]动物园
http://www.lydsy.com/JudgeOnline/problem.php?id=3670 法一:KMP+st表 抽离nxt数组,构成一棵树 若nxt[i]=j,则i作为j的子节点 那么 ...
随机推荐
- 【绝对给力】Android开发免豆资料(教程+工具+源码)地址汇总
教程下载: [免费]android界面效果全汇总.pdf http://down.51cto.com/data/209179 Android终极开发教程[pdf高清版] http://down.51c ...
- 【原创】梵高油画用深度卷积神经网络迭代10万次是什么效果? A neural style of convolutional neural networks
作为一个脱离了低级趣味的码农,春节假期闲来无事,决定做一些有意思的事情打发时间,碰巧看到这篇论文: A neural style of convolutional neural networks,译作 ...
- Jmeter(二十二)_脚本上传Gitlab
Docker部署接口自动化持续集成环境第四步,代码上传到远程仓库! 接上文:Ubuntu部署jmeter与ant Gitlab在容器中部署好了之后,本地直接打开.我们可以在里面创建项目,上传脚本. 新 ...
- python 游戏(猜单词Hangman)
1.游戏思路和流程图 实现功能:随机一个单词让玩家猜测(后续难度实现修改为成语填空,成语必须要有提示,可修改猜的次数,增加连续猜成语,难度系数随着次数的增加而增加) 游戏流程图 2. 单词库和模块 i ...
- 高可用Kubernetes集群-14. 部署Kubernetes集群性能监控平台
参考文档: Github介绍:https://github.com/kubernetes/heapster Github yaml文件: https://github.com/kubernetes/h ...
- Magento 总结
ZEND EAV 速度 作者:李淼链接:https://www.zhihu.com/question/20656910/answer/25793452来源:知乎著作权归作者所有.商业转载请联系作者获得 ...
- 《Linux内核设计与实现》Chapter 1 读书笔记
<Linux内核设计与实现>Chapter 1 读书笔记 一.Unix的特点 Unix从Multics中产生,是一个强大.健壮和稳定的操作系统. 特点 1.很简洁 2.在Unix系统中,所 ...
- 《Linux内核分析与设计实现》读书笔记一
第一章 Linux内核简介 1.1 Unix的历史 Unix的特点: Unix很简洁,仅仅提供几百个系统调用并且有一个非常明确的设计目的: 在Unix中,所有的东西都被当做文件对待. Unix的内核和 ...
- JavaMail实现邮箱之间发送邮件功能
package com.minstone.message.util; import java.util.Date; import java.util.Properties; import javax. ...
- 80C51存储器与C51内存优化
80C51在物理结构上有四个存储空间:片内程序存储器.片外程序存储器.片内数据存储器和片外数据存储器.但在逻辑上,即从用户使用的角度上,80C51有三个存储空间:片内外统一编址的64KB的程序存储器地 ...