POJ1080解题报告
依旧是如何找到递归表达形式。这个我目前也没有更好的方法,只能遇到一道题目,就解决一道,想不到就抄别人的,变成自己的,积累经验。 这道题目中提到的对齐方法,可以任意的加空格,用‘-’表示。那这个情况就比较多了,要怎么处理呢。总结递归表达形式如下,对于任意第一个前i个字符和第二个前j个字符匹配的最大匹配值dp[i][j]可取一下三个之中最大:a[i]第一个字符串的前i个字符,b[i]第二个字符串的前j个字符。
- dp[i - 1][j - 1] + socre-matrix[a[i]][b[j]]
- dp[i - 1][j] + score-matrix[a[i]]['-']
- dp[i][j - 1] + score-matrix['-'][a[j]]
代码如下:
#include
#include
using namespace std;
#define A 0
#define C 1
#define G 2
#define T 3
int scoreMatrix[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 a[101], b[101], dp[101][101];
int len1, len2;
char c;
int char2index(char c) {
switch (c) {
case 'A': return A;
case 'C': return C;
case 'G': return G;
case 'T': return T;
}
}
int threeMax(int i, int j, int k) {
return max(max(i, j), k);
}
int main() {
int t;
cin >> t;
while (t--) {
cin >> len1;
for (int i = 1; i <= len1; ++i) {
cin >> c;
a[i] =char2index(c);
}
cin >> len2;
for (int i = 1; i <= len2; ++i) {
cin >> c;
b[i] = char2index(c);
}
dp[0][0] = 0;
for (int i = 1; i <= len1; ++i)
dp[0][i] = dp[0][i -1] + scoreMatrix[4][b[i]];
for (int i = 1; i <= len2; ++i)
dp[i][0] = dp[i - 1][0] + scoreMatrix[a[i]][4];
for (int i = 1; i <= len1; ++i)
for (int j = 1; j <= len2; ++j) {
dp[i][j] = threeMax(dp[i - 1][j - 1] + scoreMatrix[a[i]][b[j]], dp[i - 1][j] + scoreMatrix[a[i]][4], dp[i][j - 1] + scoreMatrix[4][b[j]]);
}
cout << dp[len1][len2] << endl;
}
return 0;
}