动态规划(DP)，Human Gene Functions

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=1027

http://poj.org/problem?id=1080

解题报告：

1、类似于LCS

2、gene[i][j]表示str1[i-1]和str2[j-1]的分值串没有，则应该扣分

3、递推公式

temp1=gene[i-1][j-1]+score[_map[str1[i-1]]][_map[str2[j-1]]];

temp2=gene[i-1][j]+score[_map[str1[i-1]]][4];

temp3=gene[i][j-1]+score[4][_map[str2[j-1]]];

gene[i][j]=max(temp1,max(temp2,temp3));

4、在初始化边界条件时，认为一个字符串为空，则要扣分

#include <cstdio>
#include <algorithm>
#include <map>
#define NUM 105

using namespace std;

int score[][]= {{,-,-,-,-},{-,,-,-,-},{-,-,,-,-},{-,-,-,,-},{-,-,-,-,}};

map<char,int> _map;

char str1[NUM],str2[NUM];
int len1,len2;

int gene[NUM][NUM];///gene[i][j]表示基因子串str1[i-1]和str2[j-1]的分值.

int main()
{
    _map['A']=,_map['C']=,_map['G']=,_map['T']=,_map['-']=;
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%s",&len1,str1);
        scanf("%d%s",&len2,str2);
        ///初始化边界条件
        gene[][]=;
        for(int i=; i<=len1; i++)
            gene[i][]=gene[i-][]+score[_map[str1[i-]]][];
        for(int i=; i<=len2; i++)
            gene[][i]=gene[][i-]+score[][_map[str2[i-]]];
        int m1,m2,m3;
        for(int i=; i<=len1; i++)
        {
            for(int j=; j<=len2; j++)
            {
                m1=gene[i-][j]+score[_map[str1[i-]]][];///str1取i-1个字符，str2取'-';
                m2=gene[i][j-]+score[][_map[str2[j-]]];///str1取'-',str2取j-1个字符;
                m3=gene[i-][j-]+score[_map[str1[i-]]][_map[str2[j-]]];///str1取i-1个字符，str2取j-1个字符
                gene[i][j]=max(m1,max(m2,m3));
            }
        }
        printf("%d\n",gene[len1][len2]);
    }
    return ;
}