hdu 1080 Human Genne Function

这道题我不会。dp感觉老是自己写不出来啊，状态转移方程不好建啊，问题抽象。
我一定可以超越其他人，其他人混得比我好，只是因为他们比我早几年工作吧。沉下心来学习，不要被世俗所打扰，守护你的高傲。
用的是dp的思想。
和求最长LCS的思路大体相似。
对于一个序列s1,s2.
用数组dp[i][j] 表示 s1 的第i个字符和s2的第j个字符 比较时最大值 ，，这时分3种情况：

s1[i]==s2[j]   ：   那么dp[i][j]=dp[i-1][j-1]+val(s1[i],s2[j])s1[i]!=s2[j] ,那么，i 可以和‘-’ 配对，也可以选择和 j-1配对。

dp[i][j]=max(dp[i-1][j]+val(i,’-‘,dp[i][j-1]+val(‘-‘,j)

状态转移方程分析

 dp [i][j]=  max(dp[i-1][j]+val[i]['-'],dp[i][j-1]+val['-'][j],dp[i-1][j-1]+val[i][j]). （val[i][j]为题中权值）

参考博客： http://blog.csdn.net/w00w12l/article/details/7956797

下面是代码：

#include <iostream>#include<cstring>#include<cstdio>using namespace std;#define mm 105char s1[mm];char s2[mm];int dp[mm][mm];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,-10000};int getIndex(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 if(c=='-') return 4;}int val(char c,char d){    int i1=getIndex(c);    int i2=getIndex(d);    return score[i1][i2];}int MAX(int a,int b){    return b>a?b:a;}int main(){    int t;    cin>>t;    int len1,len2;    while(t--){        memset(dp,0,sizeof(dp));        scanf("%d%s",&len1,s1+1); //if index start from 0,then will need index -1                                   // so index start from 1        scanf("%d%s",&len2,s2+1);        for(int i=1;i<=len1;i++){            dp[i][0]=dp[i-1][0]+val(s1[i],'-'); // s1[i] math -        }        for(int j=1;j<=len2;j++){            dp[0][j]=dp[0][j-1]+val('-',s2[j]);        }        for(int i=1;i<=len1;i++){            for(int j=1;j<=len2;j++){                dp[i][j]=MAX(dp[i-1][j]+val(s1[i],'-'),dp[i][j-1]+val('-',s2[j]));                dp[i][j]=MAX(dp[i-1][j-1]+val(s1[i],s2[j]),dp[i][j]);            }        }        cout<<dp[len1][len2]<<endl;    }    return 0;}