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;
}



posted @ 2012-03-26 09:14  书山有路,学海无涯  阅读(180)  评论(0编辑  收藏  举报