POJ1080解题报告

依旧是如何找到递归表达形式。

这个我目前也没有更好的方法，只能遇到一道题目，就解决一道，想不到就抄别人的，变成自己的，积累经验。 这道题目中提到的对齐方法，可以任意的加空格，用‘-’表示。那这个情况就比较多了，要怎么处理呢。总结递归表达形式如下，对于任意第一个前i个字符和第二个前j个字符匹配的最大匹配值dp[i][j]可取一下三个之中最大：a[i]第一个字符串的前i个字符，b[i]第二个字符串的前j个字符。

1. dp[i - 1][j - 1] + socre-matrix[a[i]][b[j]]
2. dp[i - 1][j] + score-matrix[a[i]]['-']
3. 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;
}

转载于:https://www.cnblogs.com/sing1ee/archive/2012/02/01/2765002.html