基因匹配问题

原创 2012年03月22日 19:45:11

就是dp, 数据处理的时候要注意。而且状态转换用了只用了一个sim[2][m],节省了空间。再者就是转换的时候有权值,加上权值比较即可。

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
using namespace std;
char weight[5][5] = 
{
  {5,-1,-2,-1,-3},
  {-1,5,-3,-2,-4},
  {-2,-3,5,-2,-2},
  {-1,-2,-2,5,-1},
  {-3,-4,-2,-1,0},
};
int LenA, LenB, i, j;
int sum;
bool Cur;
char F[256];
char Str[101];
char GeneA[101], GeneB[101];
short sim[2][101];
inline short max(short A, short B)
{
  return A > B ? A:B;
} 
int main()
{
    F['A'] = 0; F['C'] = 1; F['G'] = 2; F['T'] = 3;
    scanf("%d",&sum);
    while(sum--)
    {
       scanf("%d %s",&LenA, &Str); 
       for(i=0; i<LenA; i++)
       GeneA[i+1] = F[Str[i]];
       scanf("%d %s", &LenB,&Str);
       for(i=0; i<LenB; i++)
       GeneB[i+1] = F[Str[i]];
       sim[0][0] = 0;
       for(j=1; j<=LenB; j++)
       sim[0][j] = sim[0][j-1] + weight[4][GeneB[j]];
       for(i=1, Cur=1; i<=LenA; i++,Cur =!Cur)
       {
         sim[Cur][0] = sim[!Cur][0] + weight[GeneA[i]][4];
         for(j=1; j<=LenB; j++)
         {
           sim[Cur][j] = sim[!Cur][j-1] + weight[GeneA[i]][GeneB[j]];
           sim[Cur][j] = max(sim[Cur][j], sim[Cur][j-1] + weight[4][GeneB[j]]);
           sim[Cur][j] = max(sim[Cur][j], sim[!Cur][j] + weight[GeneA[i]][4]);
         }
       }
       printf("%d\n", sim[!Cur][LenB]);
    } 
   system("pause");
   return 0;
}


相关文章推荐

【AHOI2006】bzoj1264 基因匹配

动态规划+树状数组
  • sdfzyhx
  • sdfzyhx
  • 2017年03月19日 21:45
  • 162

BZOJ 1264 Match 基因匹配 (dp 树状数组优化)

1264: [AHOI2006]基因匹配MatchTime Limit: 10 Sec Memory Limit: 162 MB Description基因匹配(match) 卡卡昨天晚上做梦梦见...
  • w4149
  • w4149
  • 2017年07月02日 11:52
  • 140

基因数据处理12之samtool的tview来查看sam的匹配文件

基因数据处理12之samtool的tview来查看sam的匹配文件 具体的之前有文章讲过:http://blog.csdn.net/xubo245/article/details/50836185 ...

bzoj1264 [AHOI2006]基因匹配Match (树状数组优化DP)

bzoj1264 [AHOI2006]基因匹配Match 原题地址:http://www.lydsy.com/JudgeOnline/problem.php?id=1264题意: 给定n,给定两个长...
  • Bfk_zr
  • Bfk_zr
  • 2017年10月18日 21:10
  • 68

基因算法解决TSP问题

  • 2013年01月15日 22:29
  • 55KB
  • 下载

基因表达量计算与差异表达分析常见问题

问 1:在没有重复实验的情况下,用 RPKM 要怎么做检验呢? 答:如果要用泊松分布做差异分析模型的话,必须要用 reads count 的。只有 RPKM值的话,可以用 RPKM 的公式反推 re...

基因数据处理30之avocado运行avocado-cli中的avocado问题1和2

问题1:avocado中的run方法中: println("stats.coverage:" + stats.coverage)调用的是: lazy val coverage = Computi...

基因数据处理31之avocado运行avocado-cli中的avocado问题3-变异识别找不到RecordGroupSample(null)

读入的read为: val fqFile = "hs38DHSE1L100F1.sam"读取结果:cleanedReads.count:1 {"readNum": 0, "contig": {"con...

基因演算法解决函数极值问题-python版

基因演算法解决函数极值问题-python版算法 基因演算法,python,numpy

基因数据处理46之cloud-scale-bwamem安装(compile.pl安装没问题)

版本:https://github.com/ytchen0323/cloud-scale-bwamem/releases/tag/v0.2.11.需要设置spark路径: ...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:基因匹配问题
举报原因:
原因补充:

(最多只允许输入30个字)