[DP_BASIC][LCS]Human Gene Functions

/*This Code is Submitted by Iris for Problem 1316 at 2013-07-26 12:19:13*/
#include <iostream>
#include <algorithm>
using namespace std;
//这个问题用DP的话,需要考虑
//dp[i][j]表示长为i,j时的最大相似度
//dp[i][j] = max{dp[i-1][j-1]+cal[i][j],
//               dp[i][j-1]+cal['-'][j],
//               dp[i-1][j]+cal[i]['-']}.
/**通过分析非常小规模的问题可以看出来,
这种转移方程可以表示出插入空格的情况**/
/**目前只能说,通过猜想状态转移方程(套模板)并且分析小规模问题
观察方程运转情况来证明其正确性,对其正确性的证明还没有更深的理解**/
int a[110],b[110],v[110][110];
int r[5][5]={{0,-3,-4,-2,-1},
             {-3,5,-1,-2,-1},
             {-4,-1,5,-3,-2},
             {-2,-2,-3,5,-2},
             {-1,-1,-2,-2,5}};
int change(char c){
    switch (c){
    case'A': return 1;
    case'C': return 2;
    case'G': return 3;
    case'T': return 4;
    default: return 0;
    }
}
int main(){
    int n,st1len,st2len,i,j;
    char ch;
    cin>>n;
    while(n--)
    {
        cin>>st1len;
        for(i=1;i<=st1len;i++){
            cin>>ch;//cin will ignore leading blanks
            a[i]=change(ch);
        }
        cin>>st2len;
        for(i=1;i<=st2len;i++){
            cin>>ch;
            b[i]=change(ch);
        }
        v[0][0]=0;
        for(i=1;i<=st1len;i++) v[i][0]=r[a[i]][0]+v[i-1][0];//str1 全部和blank匹配的结果
        for(i=1;i<=st2len;i++) v[0][i]=r[0][b[i]]+v[0][i-1];//str2 全部和blank匹配的结果
        for(i=1;i<=st1len;i++)
            for(j=1;j<=st2len;j++)
                v[i][j]=max(v[i-1][j-1]+r[a[i]][b[j]],
                max(v[i-1][j]+r[a[i]][0],v[i][j-1]+r[0][b[j]]));
            cout<<v[st1len][st2len]<<endl;
    }
    return 0;
}