Poj 1080 Human Gene Functions

题目大意：已知ATCG四种不同的碱基之间的相似度。现在给出两个碱基序列，可以在两个串间插入‘-’即空白，求两串的最大相似度。

思路：有点类似于LCS，状态转移方程为dp[i][j]=getMax(dp[(i-1)][j-1]+table[str1[i]][str2[j]],dp[(i-1)][j]+table[str1[i]][4],dp[i][j-1]+table[4][str2[j]])。

#include <stdio.h>
#include <string.h>
#include <memory.h>
int table[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 dp[105][105];
int str1[105];
int str2[105];
int m,n;
int change(char c) {
	if (c=='A')
		return 0;
	else if (c=='C')
		return 1;
	else if (c=='G')
		return 2;
	else if (c=='T')
		return 3;
	else 
		return 4;
}
int getMax(int a,int b,int c) {
	int temp;
	if (a>b)
		temp=a;
	else 
		temp=b;
	if (c>temp)
		temp=c;
	return temp;
}
int solve_dp() {
	int i,j,ctr;

	dp[0][0]=0;
	for (i=1;i<=m;i++)
		dp[i][0]=dp[i-1][0]+table[str1[i]][4];
	for (i=1;i<=n;i++)
		dp[0][i]=dp[0][i-1]+table[4][str2[i]];
	for (i=1;i<=m;i++) {
		for (j=1;j<=n;j++) {
			dp[i][j]=getMax(dp[(i-1)][j-1]+table[str1[i]][str2[j]],dp[(i-1)][j]+table[str1[i]][4],dp[i][j-1]+table[4][str2[j]]);
		}
	}
	ctr=dp[m][n];
	return ctr;
}
int main()
{
	int t,i,j;
	int result;
	char temp;

	scanf("%d",&t);
	while (t--) {
		scanf("%d",&m);
		getchar();
		for (i=1;i<=m;i++) {
			scanf("%c",&temp);
			str1[i]=change(temp);
		}
		scanf("%d",&n);
		getchar();
		for (i=1;i<=n;i++) {
			scanf("%c",&temp);
			str2[i]=change(temp);
		}
		result=solve_dp();
		printf("%d\n",result);
	}
	return 0;
}