poj 1080

状态转移方程如下：

res[i][j]=max{res[i-1][j-1]+f(str1[i],str2[j]),res[i-1][j]+f(str1[i],’-‘),res[i][j-1]+f(‘-‘,str2[j])};

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
char str1[101],str2[101];
int d[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 res[101][101];
int max(int a,int b,int c)
{
int x;
    x=a > b ? a : b;
return x > c ? x : c;
}
int hash(char c)
{
if(c=='A') return 0;
if(c=='C') return 1;
if(c=='G') return 2;
if(c=='T') return 3;
return 4;
}
int work(int n,int m)
{
int i,j,s,t;
for(i=1;i<=m;i++)
        res[0][i]=res[0][i-1]+d[hash(str2[i])][4];
for(i=1;i<=n;i++)
        res[i][0]=res[i-1][0]+d[hash(str1[i])][4];
for(i=1;i<=n;i++)
    {
        s=hash(str1[i]);
for(j=1;j<=m;j++)
        {
            t=hash(str2[j]);
            res[i][j]=max(res[i-1][j-1]+d[t][s],res[i][j-1]+d[t][4],res[i-1][j]+d[s][4]);
        }
    }
return res[n][m];
}
int main()
{
int m,n,Case;
    scanf("%d",&Case);
while(Case--)
    {
        scanf("%d %s",&n,str1+1);
        scanf("%d %s",&m,str2+1);
        memset(res,0,sizeof(res));
        printf("%d\n",work(n,m));
    }
return 0;
}

转载于:https://www.cnblogs.com/yu-chao/archive/2012/03/26/2417371.html