hdu 1080

地址:http://acm.hdu.edu.cn/showproblem.php?pid=1080

题意:匹配两个人相似度。A,G,C,T,每两个都会有一个对应的值,给出两串基因,长度可以不一样,可以在基因中间加_使两串长度一样,然后有一个对应值,求最大对应值。

mark:LCS的升级版。dp[i][j] = max(dp[i-1][j]+tab[s[0][i-1]][4], dp[i][j-1]+tab[s[1][j-1]][4], dp[i-1][j-1]+tab[s[0][i-1]][s[1][j-1]]);

代码:

#include <stdio.h>
#include <string.h>
#include <stdlib.h>

char c[2][110];
int s[2][110];
int dp[110][110];
int len[2];
int tab[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};

int max(int a, int b, int c)
{
    a = a > b ? a : b;
    return a > c ? a : c;
}

int main()
{
    int t;
    int i,j;
    scanf("%d", &t);
    while(t--)
    {
        scanf("%d%s%d%s", len, c[0], len+1, c[1]);
        for(i = 0; i < 2; i++)
            for(j = 0; j < len[i]; j++)
            {
                if(c[i][j] == 'A') s[i][j] = 0;
                else if(c[i][j] == 'C') s[i][j] = 1;
                else if(c[i][j] == 'G') s[i][j] = 2;
                else s[i][j] = 3;
            }
        dp[0][0] = 0;
        for(i = 1; i <= len[0]; i++)
            dp[i][0] = dp[i-1][0]+tab[s[0][i-1]][4];
        for(i = 1; i <= len[1]; i++)
            dp[0][i] = dp[0][i-1]+tab[s[1][i-1]][4];
        for(i = 1; i <= len[0]; i++)
            for(j = 1; j <= len[1]; j++)
                dp[i][j] = max(dp[i-1][j]+tab[s[0][i-1]][4], dp[i][j-1]+tab[s[1][j-1]][4], dp[i-1][j-1]+tab[s[0][i-1]][s[1][j-1]]);
        printf("%d\n", dp[len[0]][len[1]]);
    }
    return 0;
}

转载于:https://www.cnblogs.com/andre0506/archive/2012/08/20/2648097.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值