poj1080(dp)

题目链接:http://poj.org/problem?id=1080

题干要求:
给定两个字符串(ATCG)用-补齐，不同的字符对（--无定义）会有不同权值，要求权值和最大
（字符串长度不超过100）
考虑:

一开始的思路是找最大子列，但是思考以后有以下几个问题：

当有两个相同长度的最大子列时如何处理
当最大子列一个对应尾端一个对应非尾端，是否有可能添加-导致和更小，如果有怎么处理

那么简单的通过寻找最大子列并不能较好的解决问题，然后寻找其他最优子结构:

由于后项字符不影响前子字符串的结果，即满足最优子结构

记字符串X的长度为n，字符串Y的长度为m

考虑以下状态转移方程

如果X[i]==Y[j]

DP[i][j]=DP[i-1][j-1]+权值

如果不相等

DP[i][j]=max(DP[i][j-1]+权值,DP[i-1][j]+权值,DP[i-1][j-1]+权值)

解释一下相等情况下很好想,去掉末尾总和加5

不相等的情况分别对应(-,Y[j]) （X[i],-）(X[i],Y[j])   //其实还有(-,-)的可能但是无定义所以不考虑

观察到存在重叠子结构特征，故dp有效

附上ac代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<map>
#define MEMSET(a) memset(a,0,sizeof(a))
#define GI(a) Getindex[a]
using namespace std;
const int MAXN=105;
int V[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 n,m;
char X[MAXN];
char Y[MAXN];
int DP[MAXN][MAXN];
map<char,int> Getindex;
int dp(int i,int j){
    if(DP[i][j])
        return DP[i][j];
    if(i==0&&j==0)
        DP[i][j]=0;
    else if(i==0&&j)
        DP[i][j]=dp(i,j-1)+V[GI('-')][GI(Y[j])];
    else if(j==0&&i)
        DP[i][j]=dp(i-1,j)+V[GI(X[i])][GI('-')];
    else{
        if(X[i]==Y[j])
            DP[i][j]=dp(i-1,j-1)+5;
        else{
            int t1,t2,t3;
            t1=dp(i-1,j-1)+V[GI(X[i])][GI(Y[j])];
            t2=dp(i-1,j)+V[GI(X[i])][GI('-')];
            t3=dp(i,j-1)+V[GI('-')][GI(Y[j])];
            DP[i][j]=max(t1,max(t2,t3));
        }
    }
    return DP[i][j];
}
int main(){
    Getindex.clear();
    Getindex['A']=0;
    Getindex['C']=1;
    Getindex['G']=2;
    Getindex['T']=3;
    Getindex['-']=4;

    int T;
    scanf("%d",&T);
    while(T--){
        MEMSET(X);
        MEMSET(Y);
        MEMSET(DP);
        scanf("%d",&n);
        scanf("%s",&X[1]);
        scanf("%d",&m);
        scanf("%s",&Y[1]);
        dp(n,m);
        printf("%d\n",DP[n][m]);
    }
    return 0;
}