POJ1080

跟最长公共子序列很像

#include <iostream>
#include <cstdio>
#include <string.h>
using namespace std;
char a[105],b[105];
int dp[105][105];
char s[8]="ACGT-";
int score[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 getn(char c){
    for(int i=0;i<5;i++){
        if(s[i]==c)return i;
    }
}
int main(){
    int t,x,y;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&x);
        scanf("%s",a+1);
        scanf("%d",&y);
        scanf("%s",b+1);
        dp[0][0]=0;
        for(int i=1;i<=x;i++)dp[i][0]=score[getn(a[i])][4]+dp[i-1][0];
        for(int i=1;i<=y;i++)dp[0][i]=score[4][getn(b[i])]+dp[0][i-1];
        for(int i=1;i<=x;i++){
            for(int j=1;j<=y;j++){
                int temp1=dp[i-1][j-1]+score[getn(a[i])][getn(b[j])];                //相当于LCS的a[i]=b[j];

                int temp2=dp[i][j-1]+score[4][getn(b[j])];                            //相当于    a[i]    '-'        匹配
                                                                                    //            b[j-1]    b[j]

                int temp3=dp[i-1][j]+score[getn(a[i])][4];                            //相当于    a[i-1]    a[i]    匹配
                                                                                    //            b[jj]    '-'
                dp[i][j]=max(temp1,max(temp2,temp3));
            }
        }
        printf("%d
",dp[x][y]);
    }
    return 0;
}