HDU 1080 Human Gene Functions--DP--（最长公共子序列变形）

题意：两端基因序列匹配，四种不同的核苷酸TCGA匹配时有不同的分数，如T—G的分数是-2，还可以添加空格，空格与核苷酸匹配也有相应的分数，求最大的分数

分析：

      难点在空格的数量和位置不确定

      二维dp，dp[i][j]表示序列a的前i段和序列b的前j段匹配时的最大分数。接下来仔细分析当i和j匹配的情况：1.a[i]与b[j]匹配；2.a[i]与b[j-1]；3.a[i]与b[j+1]。

      所以方程：dp [i][j]=  max(dp[i-1][j]+mat[i]['#'],dp[i][j-1]+mat['#'][j],dp[i-1][j-1]+mat[i][j])

代码：

#include<iostream>
#include<string>
#include<cstring>
#define INF 0x3f3f3f3f
using namespace std;
int t,n,m,dp[200][200];
string a,b;
int mat[200][200];
int max(int a,int b)
{
	return a>b?a:b;
}
void f()
{
	mat['A']['A']=5;	mat['C']['C']=5;	mat['G']['G']=5;
	mat['T']['T']=5;	mat['A']['C']=-1;	mat['A']['G']=-2;
	mat['A']['T']=-1;	mat['A']['#']=-3;	mat['C']['A']=-1;
	mat['C']['G']=-3;	mat['C']['T']=-2;	mat['C']['#']=-4;
	mat['G']['A']=-2;	mat['G']['C']=-3;	mat['G']['T']=-2;
	mat['G']['#']=-2;	mat['T']['A']=-1;	mat['T']['C']=-2;
	mat['T']['G']=-2;	mat['T']['#']=-1;	mat['#']['A']=-3;
	mat['#']['C']=-4;	mat['#']['G']=-2;	mat['#']['T']=-1;
	mat['#']['#']=-INF;

}
int main()
{
    f();
	cin>>t;
	while(t--){
		cin>>n>>a;
		cin>>m>>b;
		dp[0][0]=0;
		for(int i=1;i<=n;i++) dp[i][0]=dp[i-1][0]+mat[a[i-1]]['#'];
		for(int i=1;i<=m;i++) dp[0][i]=dp[0][i-1]+mat['#'][b[i-1]];
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				dp[i][j]=max(dp[i-1][j-1]+mat[a[i-1]][b[j-1]],
				max(dp[i-1][j]+mat[a[i-1]]['#'],dp[i][j-1]+mat['#'][b[j-1]]));
			}
		}
		cout<<dp[n][m]<<endl;
	}
}