andre_joy

导航

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

posted on 2012-08-20 21:11  andre_joy  阅读(216)  评论(0编辑  收藏  举报