poj 1080 Human Gene Functions 算简单DP

题解：抄袭discuss，哈哈！

首先定义对X(xm) = {x1,x2,x3,...xm};Y(yn) = {y1,y2,y3,...yn};这两个按照题目的

要求得到最大值的排列M{...}所对应的最大值是f(m,n);
****************************************************
当xm != yn的时候
****************************************************
现在来思考一些结论:
1.
最大排列M{...}一定是
....xm
...._
由于以下不好排版,就程它为P(xm,_);
或者
...._
....yn
由于以下不好排版,就程它为P(_,yn);
或者
....xm
....yn
由于以下不好排版,就程它为P(xm,yn);
这样中的一种.因为xm,yn分别是X,Y的最后两个元素，在最大排列的上下两行中它们也必

然是最后的"非_"元素.

2.
如果最大值f(m,n)是以P(xm,_)这样的形式的,那么我们可以看到这个时候的yn是在那串

省略号里面的也就是说那串省略号是这样两个集合X(xm-1) = {x1,x2,x3,...xm-1};Y

(yn) = {y1,y2,y3,...yn}的一种按照题目要求的排列,现在我们称这两个集合的排列所

得到的值是w(m-1,n),同时定义t(xm,_)是"xm与_"相对应的值.那么就有f(m,n) = w(m-

1,n)+t(xm,_)我们现在证明w(m-1,n)是X(xm-1),Y(yn)集合排列的最大值:f(m-1,n).
由于如果存在一个w'(m-1,n) > w(m-1,n)那么w'(m-1,n)+t(xm,_) > w(m-1,n)+t(xm,_)

= f(m,n).也就是说存在一个比f(m,n)更大的值,矛盾.
所以w(m-1,n)是最大值f(m-1,n)

3.
同样可以证明如果最大值f(m,n)是以P(_,yn)这样的形式的,那么前面省略号所得到的值

w(m,n-1)也必然是最大值f(m,n-1)

4.
如果最大值是以P(xm,yn)结尾的,那么可以知道前面省略号所得到的值w(m-1,n-1)也必然

是最大值f(m-1,n-1).

从上面几个结论看过来,我们就可以知道当xm != yn时
f(m,n)一定是
f(m-1,n)+t(xm,_)
f(m,n-1)+t(xm,_)
f(m-1,n-1)+t(xm,yn)
中的最大的一个.

****************************************************************************
当xm == yn 时候
****************************************************************************
这个时候很简单，肯定是
f(m,n) = f(m-1,n-1) + t(xm,yn);

综合上面所有的东西,我们就完成了为动态规划提供子结构的过程(这样说是不是有语

病??),剩下就是LCS的编写工作了.

不过,有点要提醒一下,就是初始的时候,比如f(m,0)或者f(0,n)它们不是很简单的置成0

的,它们其中有一个不是空集
而是这样的f(m,0) = f(m-1,0)+t(m,_)     f(0,n) = f(0,n-1)+ t(_,n)
只有f(0,0) = 0.就是上下两个都是空集嘛.
小弟就是这里WA了.-_-|||

怡笑大方...

代码：

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <iostream>
#include <algorithm>
#include <cstdlib>
using namespace std;

const int MAXN = 108;
char str1[MAXN], str2[MAXN];
int t[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, -108}
              };
int dp[MAXN][MAXN];

int chge(char c) {
    if (c == 'A') return 0;
    if (c == 'C') return 1;
    if (c == 'G') return 2;
    if (c == 'T') return 3;
    if (c == '_') return 4;
    return -1;
}
int main() {
    int tt, lens1, lens2;
    scanf("%d", &tt);
    while(tt--) {
        scanf("%d%s", &lens1, str1);
        scanf("%d%s", &lens2, str2);
        for(int i = 0; i < lens1; i++) {
            int j = chge(str1[i]);
            dp[i+1][0] = dp[i][0]+t[j][4];
        }
        for(int i = 0; i < lens2; i++)
            dp[0][i+1] = dp[0][i]+t[chge(str2[i])][4];
        for(int i = 1; i <= lens1; i++)
        for(int j = 1; j <= lens2; j++) {
            dp[i][j] = dp[i-1][j-1]+t[chge(str1[i-1])][chge(str2[j-1])];
            if (str1[i-1] != str2[j-1]) {
                dp[i][j] = max(dp[i][j], dp[i-1][j]+t[chge(str1[i-1])][4]);
                dp[i][j] = max(dp[i][j], dp[i][j-1]+t[chge(str2[j-1])][4]);
            }
        }
        printf("%d\n", dp[lens1][lens2]);
    }
    return 0;
}