poj 1080 Human Gene Functions（动态规划）

原文：http://blog.acmj1991.com/?p=704

题意：给你两字符串，让你按照题目给出表格的对应关系来进行积分计算，并求最大积分

思路：动态规划 dp[i][j] 表示第一个字符串的前i与第二个字符串前j 对应的最大积分

dp[i][j]=max(dp[i-1][j-1]+mapij,dp[i-1][j]+map0i,dp[i][j-1]+mapj0)

mapij 表示第一个字符串的第i个字符于第二个字符串的第j个字符缩对应的积分

#include<stdio.h>
#include<string.h>
 
#define maxN 110
int map[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 dp[maxN][maxN];
char str1[maxN],str2[maxN];
 
int get_id(char s)
{
    if(s=='A')return 0;
    else if(s=='C')return 1;
    else if(s=='G')return 2;
    else if(s=='T')return 3;
    else return 4;
}
int max(int x,int y){return x>y?x:y;}
 
int main()
{;
    int t,len1,len2;
    scanf("%d",&t);
    while(t--){
        scanf("%d%s",&len1,str1);
        scanf("%d%s",&len2,str2);
        for(int i=1;i<=len1;i++)
            dp[0][i]=dp[0][i-1]+map[4][get_id(str1[i-1])];
        for(int i=1;i<=len2;i++)
            dp[i][0]=dp[i-1][0]+map[4][get_id(str2[i-1])];
        for(int i=1;i<=len2;i++){
            for(int j=1;j<=len1;j++){
                int u=get_id(str1[j-1]),v=get_id(str2[i-1]);
                dp[i][j]=dp[i-1][j-1]+map[u][v];
                dp[i][j]=max(dp[i][j],dp[i][j-1]+map[u][4]);
                dp[i][j]=max(dp[i][j],dp[i-1][j]+map[v][4]);
            }
        }
        int maxx=-0x0fffffff;
        for(int i=1;i<=len1;i++){
            dp[len2+1][i]=dp[len2][i-1]+dp[0][len1]-dp[0][i-1];
            maxx=max(maxx,dp[len2+1][i]);
        }
        for(int i=1;i<=len2;i++){
            dp[i][len1+1]=dp[i-1][len1]+dp[len2][0]-dp[i-1][0];
            maxx=max(maxx,dp[i][len1+1]);
        }
        maxx=max(maxx,dp[len2][len1]);
        printf("%d\n",maxx);
    }
}