【动态规划】POJ 1080 Human Gene Functions

这段时间要沉迷刷题一段时间了，就让CSDN陪我一起吧！

一、题目大意

这个题目的大致意思就是，给你两个字符串，每个字符串是由’A’、‘C’、‘G’、‘T’四个字符组成(长度不超过100)，可以在两个字符串中的任意位置增添任意的’-’，但增添后必须保证其长度相等，然后对应位置上的两个字符可以根据表去查找其对应的分数，把所有位置的分数相加，即得到最终的总分，题目要求求得总分的最大值。

二、题目思路以及AC代码

做完这道题，我才发现自己多么的菜，本来挺简单的一道题，想的极其复杂… 而且在看完题解后，又琢磨了半天才理解。

和其他的题目一样，最重要的还是要确定dp数组的定义，这里就把dp[i][j]定义为取第一个字符串的前i个字符和第二个字符串的前j个字符的情况下，分数的最大值。那这样如何进行状态转移呢？首先我给出状态转移的方程：

dp[i][j] = max(dp[i-1][j-1] + score[change(x[i])][change(y[j])],
　　　　　　dp[i-1][j] + score[change(x[i])][change(’-’))],
　　　　　　dp[i][j-1] + score[change(’-’)][change(y[j])])

可能有些人看到上面的递推公式有一点迷茫，别着急，马上就解释。首先，我们需要换一个思维来思考，不要往LCS上想，虽然他们感觉很像。首先考虑dp[i][j]肯定是由dp[i-1][j]、dp[i][j-1]、dp[i-1][j-1]这三个状态转换而来的，那上面的转换公式怎么理解呢？

dp[i-1][j]的意思就是在状态转换时，没有用到y[j]，因为下标j没有变化，所以是使用了’-’，而由于第一维由i-1变到了i，所以说明利用了x[i]，那么增加的score就是change(x[i])和change(’-’)的score.其他的两个递推公式也是如此。又由于是要最大分，所以要求他们三个的最大值。

这里主要的思想是你要将这两个字串想象成长度一样，并且都是由ACGT和-组成的字串，其实dp[i][j]就是要么让x[i]和y[j]匹配，要么让x[i]和’-‘匹配，这样j就不需要增加，要么让’-'和y[j]匹配，这样i就不需要增加，然后求三者的最大值，就可以进行状态转换了。

大致的思路讲解完毕，下面给出AC代码：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <iostream>
#include <algorithm>
using namespace std;
int dp[105][105];
int val[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,-5
};
char s1[105], s2[105];
int change(char x) {
	if (x == 'A')
		return 0;
	if (x == 'C')
		return 1;
	if (x == 'G')
		return 2;
	if (x == 'T')
		return 3;
	if (x == ' ')
		return 4;
}
int main() {                                     //dp[i][j]代表s1前i个与s2前j个匹配的最大值
	int t, l1, l2, i, j;
	scanf("%d", &t);
	while (t--) {
		scanf("%d", &l1);
		getchar();
		scanf("%s", s1 + 1);
		scanf("%d", &l2);
		getchar();
		scanf("%s", s2 + 1);
		memset(dp, 0, sizeof(dp));
		for (i = 1; i <= l1; i++)                      //理解意义是两个串的每一位进行对齐也就知道需要初始化
			dp[i][0] = dp[i - 1][0] + val[change(s1[i])][change(' ')];
		for (i = 1; i <= l2; i++)                      //s1与s2和空串匹配的值
			dp[0][i] = dp[0][i - 1] + val[change(s2[i])][change(' ')];
		for (i = 1; i <= l1; i++)
			for (j = 1; j <= l2; j++)                      //其实就是相当于s1串的一位与s2串的对齐剩下的与空格匹配
				dp[i][j] = max(dp[i - 1][j - 1] + val[change(s1[i])][change(s2[j])],
					max(dp[i - 1][j] + val[change(s1[i])][change(' ')],
						dp[i][j - 1] + val[change(s2[j])][change(' ')]));
		printf("%d\n", dp[l1][l2]);
	}
	return 0;
}

如果有什么问题，欢迎大家指正！！！