首页 > 其他 > 详细

LCS poj1080

时间:2018-11-14 17:58:18      阅读:128      评论:0      收藏:0      [点我收藏+]

题目链接:https://vjudge.net/problem/POJ-1080

参考博客:https://yq.aliyun.com/ziliao/372259

题意:给两个字符串,只含有‘A‘,‘G‘,‘C‘,‘T‘四个字符,现在根据题目给出的表里面的值,在两个字符串里面增加‘-‘,使最后的两个字符匹配的值最大,输出最大的值。

思路:开一个二维数组,dp[i][j]表示str1[i]和str2[j]之前的字符相匹配可以得到的最大价值,要确定两个要素,一个是递推公式,一个是初始化的值,首先先推递推公式:

如果str1[i]==str2[j],那么dp[i][j]=max(dp[i][j],dp[i-1][j-1]+5);

如果不相等,那么接下来可以分成三种情况:(假设用mp[a]表示字符a对应的数组下标)

1.str1[i]和一个‘-‘匹配,那么dp[i][j]=max(dp[i][j],dp[i-1][j]+value[ mp[ str1[i] ] ] [ mp[‘-‘] ]);   

2.str2[j]和一个‘-‘匹配,那么dp[i][j]=max(dp[i][j],dp[i][j-1]+value[ mp[ ‘-‘] ] [ mp[ str2[j] ] ]);

3.str1[i]直接和str2[j]匹配,那么dp[i][j]=max(dp[i][j],dp[i-1][j-1]+value[ mp[ str1[i] ] ] [ mp[ str2[j] ] ]);

然后就要确定初始值了,对于dp[0][0],肯定是等于0,而dp[i][0],就是dp[i][0]=dp[i-1][0]+value[ mp[ str1[i] ]][ mp[‘-‘] ];dp[0][j]也是一样,其他的dp[i][j]就赋无穷小就可以了。声明一下,这不是我想出来的,我脑瓜子还是不行。

我的代码:

#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<map>
#include<stack>
#include<cmath>
#include<vector>
#include<set>
#include<cstdio>
#include<string>
#include<deque> 
using namespace std;
typedef long long LL;
#define eps 1e-8
#define INF 0x3f3f3f3f
#define maxn 1005
map<char,int>mp;
int n,m,k,t;
int dp[205][205]; 
char str1[105],str2[105];
int value[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}};
void init()
{
    mp[A]=0;
    mp[C]=1;
    mp[G]=2;
    mp[T]=3;
    mp[-]=4;
}
int main()
{
    init();
    scanf("%d",&t);
    while(t--){
        scanf("%d %s",&n,str1+1);
        scanf("%d %s",&m,str2+1);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                dp[i][j]=-INF;
            }
        }
        dp[0][0]=0;//设为0 
        for(int i=1;i<=n;i++){//初始化 
            dp[i][0]=value[mp[str1[i]]][4]+dp[i-1][0];
        }
        for(int i=1;i<=m;i++){
            dp[0][i]=value[4][mp[str2[i]]]+dp[0][i-1];
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                if(str1[i]==str2[j])
                dp[i][j]=max(dp[i-1][j-1]+5,dp[i][j]);
                else
                {
                    dp[i][j]=max(dp[i-1][j]+value[mp[str1[i]]][4],dp[i][j]);
                    dp[i][j]=max(dp[i][j-1]+value[4][mp[str2[j]]],dp[i][j]);
                    dp[i][j]=max(dp[i-1][j-1]+value[mp[str1[i]]][mp[str2[j]]],dp[i][j]);
                }
            }
        }
        printf("%d\n",dp[n][m]);
    } 
    return 0;
}

 

LCS poj1080

原文:https://www.cnblogs.com/6262369sss/p/9959237.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!