POJ 1080 HUMAN GENE FUNCTION

fi][j] 表示 s1到i s2到j时能达到的最大的相似度。有三张情况。

s1[i] 与s2[j] 匹配。

s1[i]与‘-’匹配。

s2[j]与‘-’匹配。

取三者中最大的，还有就是初始化的时候需要注意并不是单纯的把全部都置为0就行了。空的时候是负的，就是说

for(int i = 1; i <= s1.size(); ++i)
            f[i][0] = f[i - 1][0] + score[m[s1[i - 1]]][4];
for(int j = 1; j <= s2.size(); ++j)
            f[0][j] = f[0][j - 1] + score[4][m[s2[j - 1]]];

 

这个初始化错了好多次。不仅仅是 f[0][j] = score[4][m[s2[j - 1]]];后面还要加上前面匹配的f[0][j - 1]这个部分！

 

304k  0MS

#include<iostream>
#include<map>
#include<string>
#include<cstring>
using namespace std;

int f[110][110];
string s1, s2;
int T;
map<char,int> m;
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, 0}};

int max(int a, int b, int c)
{
    int k = a > b ? a : b;
    return k > c? k : c;
}

int main()
{
    void ini();
    ini();
    cin >> T;
    while(T--)
    {
        int n;
        cin >> n >> s1;
        cin >> n >> s2;
        memset(f, 0, sizeof(f));
        for(int i = 1; i <= s1.size(); ++i)
            f[i][0] = f[i - 1][0] + score[m[s1[i - 1]]][4];
        for(int j = 1; j <= s2.size(); ++j)
            f[0][j] = f[0][j - 1] + score[4][m[s2[j - 1]]];
        for(int i = 1; i <= s1.size(); ++i)
        {
            for(int j = 1; j <= s2.size(); ++j)
                 f[i][j] = max(f[i - 1][j - 1] + score[m[s1[i - 1]]][m[s2[j - 1]]], f[i - 1][j] + score[m[s1[i - 1]]][4], f[i][j - 1] + score[4][m[s2[j - 1]]]);
        }
        cout << f[s1.size()][s2.size()] << endl;
    }
    return 0;
}

void ini()
{
    m['A'] = 0;
    m['C'] = 1;
    m['G'] = 2;
    m['T'] = 3;
}