P1140 相似基因

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

#include<iostream>
#include<algorithm>
using namespace std;
const int N=1e2+50;
int dp[N][N];
int la,lb,a[N],b[N];
int v[6][6]={
  {0,0,0,0,0,0},
  {0,5,-1,-2,-1,-3},
  {0,-1,5,-3,-2,-4},
  {0,-2,-3,5,-2,-2},
  {0,-1,-2,-2,5,-1},
  {0,-3,-4,-2,-1,0}
};
int main()
{
	for(int i=1;i<=N;i++)
	 for(int j=1;j<=N;j++)
	  dp[i][j]=-1e6;
	cin>>la;
	for(int i=1;i<=la;i++)
	{
		char x;
		cin>>x;
		if(x=='A') a[i]=1;
		if(x=='C') a[i]=2;
		if(x=='G') a[i]=3;
		if(x=='T') a[i]=4;
	}
	cin>>lb;
	for(int i=1;i<=lb;i++)
	{
		char x;
		cin>>x;
		if(x=='A') b[i]=1;
		if(x=='C') b[i]=2;
		if(x=='G') b[i]=3;
		if(x=='T') b[i]=4;
	}
	for(int i=1;i<=la;i++)
	 dp[i][0]=dp[i-1][0]+v[a[i]][5];
	for(int i=1;i<=lb;i++)
	 dp[0][i]=dp[0][i-1]+v[5][b[i]];
	for(int i=1;i<=la;i++)
	 for(int j=1;j<=lb;j++)
	 {
	 	dp[i][j]=max(dp[i][j],dp[i-1][j]+v[a[i]][5]);
	 	dp[i][j]=max(dp[i][j],dp[i][j-1]+v[5][b[j]]);
	 	dp[i][j]=max(dp[i][j],dp[i-1][j-1]+v[a[i]][b[j]]);
	 }
	cout<<dp[la][lb];
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值