HDU 1080 Human Gene Functions DP -

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=1080

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxn=100+5;
const int INF=0x3f3f3f3f;
#define max3(a,b,c) max(a,max(b,c))
char s1[maxn],s2[maxn];
int d[maxn][maxn],len1,len2;
//d[i][j]表示s1前i个和s2前j个匹配的最大值
char str[]="ACGT-";
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 inline Value(char a,char b){
	int p1=find(str,str+5,a)-str;
	int p2=find(str,str+5,b)-str;
	return score[p1][p2];
}
int solve()
{
	memset(d,-INF,sizeof(d));

	d[0][0]=0;
	for(int i=1;i<=len1;i++){
		d[i][0]=d[i-1][0]+Value(s1[i],'-');
		d[0][i]=d[0][i-1]+Value(s2[i],'-');
	}
	for(int i=1;i<=len1;i++)
	for(int j=1;j<=len2;j++)
		d[i][j]=max3(d[i-1][j-1]+Value(s1[i],s2[j]),d[i-1][j]+Value(s1[i],'-'),d[i][j-1]+Value(s2[j],'-'));
					//s1[i]和s2[j] 匹配           s1[i]和'-'匹配              s2[j]和'-'匹配
	
	return d[len1][len2];
}
int main(int argc, char const *argv[])
{
	int T; cin>>T;
	while(T--)
	{
		scanf("%d %s",&len1,s1+1);
		scanf("%d %s",&len2,s2+1);
		
		cout<<solve()<<endl;
	}
	return 0;
}