poj1080 Human Gene Functions

#include<iostream>#include<cstring>#include<cstdio>#define max(a,b) (a>b?a:b)using namespace std;char str1[101],str2[101];char map[101];int gene[101][101];int main(){map['A']=0;map['C']=1;map['G']=2;map['T']=3;map['-']=4;int score[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 icase;int i,j,k;int first,second;scanf("%d",&icase);for(k=0;k<icase;k++){scanf("%d%s",&first,str1);scanf("%d%s",&second,str2);gene[0][0]=0;for(i=1;i<=first;i++)gene[i][0]=gene[i-1][0]+score[map[str1[i-1]]][4];for(i=1;i<=second;i++)gene[0][i]=gene[0][i-1]+score[4][map[str2[i-1]]];int m1,m2,m3;for(i=1;i<=first;i++)for(j=1;j<=second;j++){m1=gene[i-1][j]+score[map[str1[i-1]]][4];m2=gene[i][j-1]+score[4][map[str2[j-1]]];m3=gene[i-1][j-1]+score[map[str1[i-1]]][map[str2[j-1]]];gene[i][j]=max(m1,max(m2,m3));}printf("%d\n",gene[first][second]);}return 0;}