poj——Human Gene Functions

最新推荐文章于 2020-03-02 17:08:49 发布

奈何陵琛

最新推荐文章于 2020-03-02 17:08:49 发布

阅读量450

点赞数

分类专栏：动态规划文章标签： poj dp ACM map

本文链接：https://blog.csdn.net/u012794562/article/details/38900269

版权

动态规划专栏收录该内容

9 篇文章 0 订阅

订阅专栏

题意：给你2个DNA序列，求它们的最大相似值，你可以在任一序列中加入空格“-”。例如：AGTGATG和GTTAG，根据对应的值可知GTTAG变为-GTTA-G时和第一个序列的相似值最大。

分析：dp[i][j]代表a序列的前i个和b序列的前j个配对的最大相似值;

对于a[i]来说，要么和b[j]配对，要么就和“-”配对；对于b[j]也是一样；

所以可以dp求解：

dp[i][j]=max(dp[i-1][j-1]+value[a[i]][b[j]],dp[i-1][j]+value[a[i]]['-'],dp[i][j-1]+value['-'][b[j]])

dp[len1][len2]即为所求。

#include<cstdio>
#include<map>
#include<cstring>
using namespace std;
int len1,len2,dp[101][101];
map<char,int>what;
char a[101],b[101];
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 deal()
{
	int i,j;
    dp[0][0]=0;
    for(i=1;i<=len1;i++)
        dp[i][0]=value[what[a[i]]][what['-']]+dp[i-1][0];
    for(j=1;j<=len2;j++)
        dp[0][j]=value[what['-']][what[b[j]]]+dp[0][j-1];
    for(i=1;i<=len1;i++)
        for(j=1;j<=len2;j++)
        {
            dp[i][j]=dp[i-1][j-1]+value[what[a[i]]][what[b[j]]];
            dp[i][j]=max(dp[i][j],dp[i-1][j]+value[what[a[i]]][what['-']]);
            dp[i][j]=max(dp[i][j],dp[i][j-1]+value[what['-']][what[b[j]]]);
        }
}
int main()
{
	int T,i,j;
	what['A']=0;what['C']=1;what['G']=2;what['T']=3;what['-']=4;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&len1);
		getchar();
		for(i=1;i<=len1;i++)
			scanf("%c",&a[i]);
		scanf("%d",&len2);
		getchar();
		for(i=1;i<=len2;i++)
			scanf("%c",&b[i]);
		deal();
		printf("%d\n",dp[len1][len2]);
	}
	return 0;
}