lg P1140 相似基因

最新推荐文章于 2025-01-21 23:42:28 发布

J__lxj

最新推荐文章于 2025-01-21 23:42:28 发布

阅读量111

点赞数

CC 4.0 BY-SA版权

分类专栏：动态规划文章标签：动态规划

本文链接：https://blog.youkuaiyun.com/weixin_44652585/article/details/89364231

动态规划专栏收录该内容

13 篇文章

订阅专栏

本文介绍了一种用于DNA序列匹配的动态规划算法，通过将字母转换为编号，建立价值表，并初始化f数组，实现了序列匹配的价值最大化。算法考虑了三种匹配情况，包括A和空匹配、B和空匹配以及A和B之间的匹配。

问题描述：戳这里
就是找匹配，使得价值最大。
注意对f数组的初值，和预处理。
1，把字母换做是编号。
2，把对应的价值表做成一个x数组
3，剩余初值设为负无穷，因为价值里有负值。
4，边界 0,0 即A 和全是空的匹配， B 和全是空的匹配。
5，三种选择情况，A 和空匹配，B 和空匹配，A和B 匹配。
考虑以上情况后，此题可以轻松解决。
code：

#include<bits/stdc++.h>
 using namespace std;
const int x[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;
int a[500],b[500],f[500][500];
int main()
{
  cin>>n; 
  for(int i=1;i<=n;++i) 
  {
    char c;
    cin>>c;
    if(c=='A') a[i]=0;
    if(c=='C') a[i]=1;
    if(c=='G') a[i]=2;
    if(c=='T') a[i]=3;
  }
  cin>>m;
  for(int i=1;i<=m;++i)
  {
    char c;
    cin>>c;
    if(c=='A') b[i]=0;
    if(c=='C') b[i]=1;
    if(c=='G') b[i]=2;
    if(c=='T') b[i]=3;	
  }
  for(int i=1;i<=n;++i)
    for(int j=1;j<=m;++j) f[i][j]=-1e7;
  for(int i=1;i<=n;++i) f[i][0]=f[i-1][0]+x[a[i]][4];
  for(int i=1;i<=m;++i) f[0][i]=f[0][i-1]+x[b[i]][4];
  for(int i=1;i<=n;i++)
  
   for(int j=1;j<=m;j++)
   {
   	f[i][j]=max(f[i][j],f[i-1][j]+x[a[i]][4]);
   	f[i][j]=max(f[i][j],f[i][j-1]+x[4][b[j]]);
   	f[i][j]=max(f[i][j],f[i-1][j-1]+x[a[i]][b[j]]);
   }
   cout<<f[n][m];
  return 0;
}