4169: DNA序列---复制问题

DNA序列匹配算法

最新推荐文章于 2020-12-20 04:19:51 发布

cheerful_man

最新推荐文章于 2020-12-20 04:19:51 发布

阅读量830

点赞数

CC 4.0 BY-SA版权

分类专栏：咸鱼王臭垃圾还不刷题系列

本文链接：https://blog.youkuaiyun.com/foolishpichao/article/details/90271602

咸鱼王臭垃圾还不刷题系列专栏收录该内容

35 篇文章

订阅专栏

本文介绍了一种用于在DNA序列中查找特定模式的算法，通过正向和反向匹配，统计在一个较长的DNA序列中，特定模式出现的次数，特别关注半保留复制和逆复制过程。

题目描述

    核酸与蛋白质是生物体内两大重要的生物大分子，其中核酸是生物遗传信息从上一代传到下一代的的重要载体。核酸的基本组成单元是核苷酸，可将核苷酸分为脱氧核苷酸(DNA)和核糖核苷酸(RNA)。
    我们以DNA为例，所有的DNA核苷酸的磷酸基团和糖基团都是一样的，但其化学组成的碱基基团有四种类型，腺嘌呤（A）、鸟嘌呤（G）、胞嘧啶（C）、胸腺嘧啶（T）。因此，对DNA序列的单链研究可以转化为研究由A、G、C、T四个字母组成的字符串。
    DNA复制是指以原始的DNA分子为模板合成出相同分子的过程。在一些酶的作用下，DNA复制过程的有四种可能：完全复制、半保留复制、逆复制、逆半保留复制。
    给定一个带有遗传信息的模板mode，你能在DNA序列s中找到这样的复制过程最多有多少次？

输入

第一行：整数T ，表示以下有T组测试数据（ 1 ≤ T ≤ 4）
每组数据，占一行：
mode s （|mode |≤ 20， |s| ≤ 1000 ， mode与s之间只有一个空格）

输出

对每组测试数据，输出占一行。

样例输入

2
ACGA AGACGACGAGCAGCAACGA
AAA ACGATTAAAAAGTCTAAAA

样例输出

5
5

提示

1 半保留复制是指本次复制可以与前次复制部分重叠。例如：ACGACGA,算2次复制。类似逆半保留
复制，例如：AGCAGCA,算2次复制。
2 若模板是回文串，逆复制过程不再统计。

#include <iostream>
#include <string>
#include<algorithm>
using namespace std;
int ans;
bool chieck(string a,int len)
{
	int i,j;
	for(i=0,j=len-1;i<=(len-1)/2,j>(len-2)/2;i++,j--)
	{
		if(a[i]!=a[j])
		return false;
	}
	return true;
	
}

void jisuan(string mode,string s)
{
	string::size_type it;
	for(int i=0;i!=s.size();i++)
		{
			for(it=0;it<mode.size();it++)
			{
				if(s[i+it]!=mode[it])
				break;
			}
			if(it==mode.size())
			ans++;
		}
}
int main()
{
	int n;
	string mode,s;
	cin>>n;
	while(n--)
	{
		cin>>mode>>s;
		ans=0;
		if(chieck(mode,mode.size()))
		jisuan(mode,s);
	else
	{
		jisuan(mode,s);
		reverse(s.begin(),s.end());
		jisuan(mode,s);
    }
	
	cout<<ans<<endl;
    }
	return 0;
}
//ps：写的什么狗屎，垃圾思路，低等算法。