poj1080

类似于最长公共子序列

dp[i][j]表示字符窜1 的第  i  个字符匹配  字符串 2 的第 j 个位置的最大值

如果 1 中第 i 个字符和 2 中第j个字符匹配 dp[i][j]=dp[i-1][j-1]+匹配值

不匹则用 * 来表示，一个在字符串1中的补*，一个在字符串2中补*

1中dp[i][j-1]+匹配值   2中dp[i][j-1]+匹配值（想不懂的可以用笔画星星补一下就知道 了）

#include<iostream>
#include<cstring>
#include<cstdio>
#include<map>
using namespace std;
int val[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[111][111];
int T,n1,n2;
map<char,int>m;
int main()
{
    cin>>T;
    m['A']=0;m['C']=1;m['G']=2;m['T']=3;m['*']=4;
    while(T--)
    {
        string s1,s2;
        cin>>n1>>s1>>n2>>s2;
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n1;i++)
        {
            int x=m[s1[i-1]];
            int y=m['*'];
            dp[i][0]=dp[i-1][0]+val[x][y];
        }
        for(int i=1;i<=n2;i++)
        {
            int x=m['*'];int y=m[s2[i-1]];
            dp[0][i]=dp[0][i-1]+val[x][y];
        }
        for(int i=1;i<=n1;i++)
            for(int j=1;j<=n2;j++)
            {
                int a=dp[i-1][j-1]+val[m[s1[i-1]]][m[s2[j-1]]];
                int b=dp[i][j-1]+val[m['*']][m[s2[j-1]]];
                int c=dp[i-1][j]+val[m[s1[i-1]]][m['*']];
                dp[i][j]=max(a,max(b,c));
            }
        cout<<dp[n1][n2]<<endl;

    }
    return 0;
}