storm

努力学习编译知识,恶补计算机体系结构!

基因匹配问题

就是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;
}


阅读更多
文章标签: system c
个人分类: 算法
想对作者说点什么? 我来说一句

基因算法设计与应用-基因算法

2009年07月08日 389KB 下载

语法分析-if-else的匹配问题

2010年06月02日 15KB 下载

栈的经典应用-括号匹配问题

2009年10月22日 986B 下载

没有更多推荐了,返回首页

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭