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

二维最长公共子序列

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std ;

char s1[105],s2[105];

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}};
//dp[i][j]表示长i和j的字符串的最长公共子序列 
int dp[105][105];
int mp[2][105];
int cal(char x){
    if(x=='A')return 0;
    if(x=='C')return 1;
    if(x=='G')return 2;
    if(x=='T')return 3; 
}
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        int n,m;
        scanf("%d%s%d%s",&n,s1,&m,s2);
        for(int i=0;i<2;i++){
            if(!i){
                for(int j=0;j<n;j++)
                    mp[0][j]=cal(s1[j]);
            }
            else{
                for(int j=0;j<m;j++)
                    mp[1][j]=cal(s2[j]);
            }
        }
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++)
            dp[i][0]=dp[i-1][0]+tab[mp[0][i-1]][4];
        for(int i=1;i<=m;i++)
            dp[0][i]=dp[0][i-1]+tab[mp[1][i-1]][4];
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                dp[i][j]=max(dp[i-1][j]+tab[mp[0][i-1]][4],max(dp[i][j-1]+tab[mp[1][j-1]][4],dp[i-1][j-1]+tab[mp[0][i-1]][mp[1][j-1]]));
            }
        }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}
View Code