zoj1027

 

解析思路：

状态转移方程b[i][j]=max(b[i-1][j-1]+t[i][j],b[i][j-1]+t[4][j],b[i-1][j]+t[i][4]);i表示第一个字符串位子，j表示第二个字符串的位子；b[i][j]表示1到i，1到j位子的值

#include<stdio.h>#include<string.h>#include<math.h>#include<stdlib.h>int t[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 max(int a,int b,int c){  int temp;  temp=a>b?a:b;  return temp=c>temp?c:temp;}int main(){  int n;  scanf("%d",&n);  while(n--)  {    int i,j,k,m,temp1[200],temp2[200];    char s[2][200];    for(j=0;j<2;j++) scanf("%d %s",&m,s[j]);    for(k=1;k<=strlen(s[0]);k++)    {      if(s[0][k-1]=='A') temp1[k]=0;      else if(s[0][k-1]=='C') temp1[k]=1;      else if(s[0][k-1]=='G') temp1[k]=2;      else if(s[0][k-1]=='T') temp1[k]=3;    }    for(k=1;k<=strlen(s[1]);k++)    {      if(s[1][k-1]=='A') temp2[k]=0;      else if(s[1][k-1]=='C')  temp2[k]=1;      else if(s[1][k-1]=='G')  temp2[k]=2;      else if(s[1][k-1]=='T')  temp2[k]=3;    }    int p,q,best[200][200];    best[0][0]=0;    for(j=1;j<=strlen(s[0]);j++)      best[j][0]=best[j-1][0]+t[temp1[j]][4];    for(j=1;j<=strlen(s[1]);j++)      best[0][j]=best[0][j-1]+t[4][temp2[j]];    for(p=1;p<=strlen(s[0]);p++)    {      for(q=1;q<=strlen(s[1]);q++)      {        best[p][q]=max(t[temp1[p]][temp2[q]]+best[p-1][q-1],        best[p-1][q]+t[temp1[p]][4],best[p][q-1]+t[4][temp2[q]]);        }    }    printf("%d\n",best[p-1][q-1]);  }  return 0;}