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