BZOJ 1264[AHOI 2006 基因匹配]
来源:互联网 发布:c语言运算符优先级表 编辑:程序博客网 时间:2024/05/01 17:46
1264: [AHOI2006]基因匹配Match
基因匹配(match) 卡卡昨天晚上做梦梦见他和可可来到了另外一个星球,这个星球上生物的DNA序列由无数种碱基排列而成(地球上只有4种),而更奇怪的是,组成DNA序列的每一种碱基在该序列中正好出现5次!这样如果一个DNA序列有N种不同的碱基构成,那么它的长度一定是5N。 卡卡醒来后向可可叙述了这个奇怪的梦,而可可这些日子正在研究生物信息学中的基因匹配问题,于是他决定为这个奇怪星球上的生物写一个简单的DNA匹配程序。 为了描述基因匹配的原理,我们需要先定义子序列的概念:若从一个DNA序列(字符串)s中任意抽取一些碱基(字符),将它们仍按在s中的顺序排列成一个新串u,则称u是s的一个子序列。对于两个DNA序列s1和s2,如果存在一个序列u同时成为s1和s2的子序列,则称u是s1和s2的公共子序列。 卡卡已知两个DNA序列s1和s2,求s1和s2的最大匹配就是指s1和s2最长公共子序列的长度。 [任务] 编写一个程序: 从输入文件中读入两个等长的DNA序列; 计算它们的最大匹配; 向输出文件打印你得到的结果。
Input
输入文件中第一行有一个整数N,表示这个星球上某种生物使用了N种不同的碱基,以后将它们编号为1…N的整数。 以下还有两行,每行描述一个DNA序列:包含5N个1…N的整数,且每一个整数在对应的序列中正好出现5次。
Output
输出文件中只有一个整数,即两个DNA序列的最大匹配数目。
Sample Input
2
1 1 2 2 1 1 2 1 2 2
1 2 2 2 1 1 2 2 1 1
1 1 2 2 1 1 2 1 2 2
1 2 2 2 1 1 2 2 1 1
Sample Output
7
HINT
[数据约束和评分方法]
60%的测试数据中:1<=N <= 1 000
100%的测试数据中:1<=N <= 20 000
此题一眼看上去就是LCS,但是多了些判断
因为对于LCS,当A[i]==B[i]时,LCS的值+1,所以我们要做的就是找相同的值。
因为每一个数一定出现5次,所以可以记录下A数组中相同5个数的位置,然后扫一遍B,用B[i]来更新他们
即找到每个位置f[pos],然后找到f[1,....,pos-1]中的最大值+1即可,
根据别人的题解,要使用树状数组维护,不用的话我也不知道可不可以
下面是代码:
#include<cstdio>#include<algorithm>#define M 20020using namespace std;int n,ans,a[M*5],b[M*5],c[M*5],f[M*5],pos[M][6];void Update(int x,int y){ for(;x<=n*5;x+=x&-x) c[x]=max(c[x],y);}int Get_Ans(int x){ int re=0; for(;x;x-=x&-x) re=max(re,c[x]); return re;}int main(){ scanf("%d",&n); for(int i=1;i<=n*5;i++) { scanf("%d",&a[i]); pos[ a[i] ][ ++pos[a[i]][0] ]=i; } for(int i=1;i<=n*5;i++)scanf("%d",&b[i]); for(int i=1;i<=n*5;i++) { for(int j=5;j;j--) { int k=pos[b[i]][j]; f[k]=max( f[k] , Get_Ans(k-1)+1 ); Update(k,f[k]); ans=max(ans,f[k]); } } printf("%d\n",ans); return 0;}
1 0
- BZOJ 1264[AHOI 2006 基因匹配]
- AHOI 基因匹配 dp LCS LIS
- bzoj 1264: [AHOI2006]基因匹配Match
- bzoj 1264: [AHOI2006]基因匹配Match
- BZOJ 1264 AHOI2006 基因匹配Match 动态规划+树状数组
- BZOJ 1264: [AHOI2006]基因匹配Match|动态规划
- bzoj 1264 [AHOI2006]基因匹配Match(DP+树状数组)
- BZOJ 1264: [AHOI2006]基因匹配Match 树状数组,DP
- bzoj 1264: [AHOI2006]基因匹配Match (DP+树状数组)
- BZOJ 1264 [AHOI2006]基因匹配Match DP+BIT
- BZOJ 1264 Match 基因匹配 (dp 树状数组优化)
- bzoj 1264: [AHOI2006]基因匹配Match 树状数组
- bzoj 1264: [AHOI2006]基因匹配Match(树状数组)
- 1264: [AHOI2006]基因匹配Match
- BZOJ 1266 AHOI 2006 上学路线route 最小割
- BZOJ 1269 [AHOI 2006] STL(rope)\SPLAY 解题报告
- 【AHOI 2013】【BZOJ 3236】作业
- [BZOJ 3238][AHOI 2013]差异
- RabbitMQ消息队列(四):分发到多Consumer(Publish/Subscribe)
- Linux芯片级移植与底层驱动(基于3.7.4内核)
- Comet4J简介
- 利用libvirt给虚拟机添加virtio磁盘驱动
- Google Web Designer(HTML5网页可视化设计)
- BZOJ 1264[AHOI 2006 基因匹配]
- RabbitMQ消息队列(六):使用主题进行消息分发
- codeforces B
- lintcode:Minimum Adjustment Cost
- RabbitMQ消息队列(七):适用于云计算集群的远程调用(RPC)
- poj2367
- 2048游戏python版
- hihoCoder[Offer收割]编程练习赛3题目解析
- 2010年山东省第一届ACM大学生程序设计竞赛:Hello World!