hdu 1080 Human Gene Functions dp

题意：

给定两个基因，对比两个基因的相似度。即给定一个表，行列项均为ACTG-，表示基因，表内数字表示行和列对应基因的相似度，-为空。先给定两个字符串表示基因a和b，长度分别为n和m，我们可以在字符串中间加入-，使得两个字符串长度相同；现在问通过加入-，最大的基因相似度是多少？

题解：

dp题目，dp[i][j]表示用了a串用了前i个字符，b串用了前j个字符的情况下，得到的最大基因值。答案就是dp[n][m]。

由于‘-’不能和‘-’配对，所以可以得出dp公式：dp[i][j]=max(dp[i-1][j-1]+get_num(a[i-1],b[j-1]),dp[i-1][j]+get_num(a[i-1,]'-'),dp[i][j-1]+get_num('-',b[j-1]));

代码：

#include <iostream>
#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <ctime>
#include <vector>
#include <cmath>
#include <map>
#include <cstdlib>
using namespace std;

const int maxn=110;
int dp[maxn][maxn];
int e[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,-100000}};
map<char,int>mm;
char a[maxn],b[maxn];
void init()
{
    mm['A']=0;mm['C']=1;mm['G']=2;mm['T']=3,mm['-']=4;
}
int get_num(char a,char b)
{
    return e[mm[a]][mm[b]];
}
int main()
{
    init();
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int i,j,k,n,m;
        scanf("%d%s%d%s",&n,a,&m,b);
        memset(dp,0,sizeof(dp));
        for(i=1;i<=n;i++)
            dp[i][0]=dp[i-1][0]+get_num(a[i-1],'-');
        for(i=1;i<=m;i++)
            dp[0][i]=dp[0][i-1]+get_num('-',b[i-1]);
        //printf("%s %c\n",a,a[3]);
        //printf("%d %d\n",get_num('-',b[1]),dp[0][2]);
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=m;j++)
            {
                //printf("%d\n",j);
                dp[i][j]=max(dp[i-1][j]+get_num(a[i-1],'-'),dp[i][j-1]+get_num('-',b[j-1]));
                dp[i][j]=max(dp[i][j],dp[i-1][j-1]+get_num(a[i-1],b[j-1]));
            }
        }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}