hdu 1080

题目：http://acm.hdu.edu.cn/showproblem.php?pid=1080、

模型很像最长公共子序列类的DP。状态转移方程是  dp [i][j]=  max(dp[i-1][j]+val[i]['-'],dp[i][j-1]+val['-'][j],dp[i-1][j-1]+val[i][j]). （val[i][j]为题中权值）

dp[i][j]为长度为i的s1串与长度为j的s2串匹配得到的最大权值。  那么对于当前s1串上的s1[i]字符可以对应的有两种方案： ⑴‘-’符号和⑵ s[j].   同样。 s[j]也可以对应‘-’符号。

下面是AC代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn = 110;
char s1[maxn],s2[maxn];
int len1,len2;
int dp[maxn][maxn];
int table[6][6]={
    0,0,0,0,0,0,
    0,5,-1,-2,-1,-3,
    0,-1,5,-3,-2,-4,
    0,-2,-3,5,-2,-2,
    0,-1,-2,-2,5,-1,
    0,-3,-4,-2,-1,-1000000,
};
int getindex(char a){
    if(a=='A') return 1;
    else if(a=='C') return 2;
    else if(a=='G') return 3;
    else if(a=='T') return 4;
    else return 5;
}
int val(char a,char b){
    int c=getindex(a);  int d=getindex(b);
    return table[c][d];
}
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        memset(dp,0,sizeof(dp));
        scanf("%d %s",&len1,s1+1);
        scanf("%d %s",&len2,s2+1);
        int max_len = max(len1,len2);
        for(int i=1;i<=len1;i++){
            dp[i][0]=dp[i-1][0]+val(s1[i],'-');
        }
        for(int i=1;i<=len2;i++){
            dp[0][i]=dp[0][i-1]+val('-',s2[i]);
        }
     //   printf("%d\n",val('A','A'));
        for(int i=1;i<=len1;i++){
           for(int j=1;j<=len2;j++){
        dp[i][j]=max(dp[i-1][j]+val(s1[i],'-'),dp[i][j-1]+val('-',s2[j]));
        dp[i][j]=max(dp[i-1][j-1]+val(s1[i],s2[j]),dp[i][j]);
           }
        }
//        for(int i=1;i<=len1;i++){
//           for(int j=1;j<=len2;j++){
//            printf("%d ",dp[i][j]);
//            }
//            printf("\n");
//        }
        printf("%d\n",dp[len1][len2]);
    }
    return 0;
}
/*
2
7 AGTGATG
5 GTTAG
7 AGCTATT
9 AGCTTTAAA

*/