ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

DNA序列

2021-11-29 21:30:17  阅读:165  来源: 互联网

标签:子串 DNA int GC 序列 Ratio


DNA序列

描述
一个 DNA 序列由 A/C/G/T 四个字母的排列组合组成。 G 和 C 的比例(定义为 GC-Ratio )是序列中 G 和 C 两个字母的总的出现次数除以总的字母数目(也就是序列长度)。在基因工程中,这个比例非常重要。因为高的 GC-Ratio 可能是基因的起始点。

给定一个很长的 DNA 序列,以及限定的子串长度 N ,请帮助研究人员在给出的 DNA 序列中从左往右找出 GC-Ratio 最高且长度为 N 的第一个子串。
DNA序列为 ACGT 的子串有: ACG , CG , CGT 等等,但是没有 AGT , CT 等等

数据范围:字符串长度满足 1<=n<=1000 ,输入的字符串只包含 A/C/G/T 字母
输入描述:
输入一个string型基因序列,和int型子串的长度

输出描述:
找出GC比例最高的子串,如果有多个则输出第一个的子串

示例1
输入:
ACGT
2

输出:
CG

说明:
ACGT长度为2的子串有AC,CG,GT3个,其中AC和GT2个的GC-Ratio都为0.5,CG为1,故输出CG
示例2
输入:
AACTGTGCACGACCTGA
5

输出:
GCACG

说明:
虽然CGACC的GC-Ratio也是最高,但它是从左往右找到的GC-Ratio最高的第2个子串,所以只能输出GCACG。

滑动窗口

#include<iostream>
#include<string>
using namespace std;

int main(){
    string str;
    int N;
    while(cin>>str>>N){
        string ans="";
        int maxnum=0;
        for(int i=0;i<str.size()-N+1;++i){
            int cnt=0;
            for(int j=i;j<i+N;++j){
                if(str[j]=='C' || str[j]=='G') cnt++;
            }
            if(cnt>maxnum){
                maxnum=cnt;
                ans=str.substr(i,N);
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}

标签:子串,DNA,int,GC,序列,Ratio
来源: https://blog.csdn.net/Mwwwwwwww/article/details/121620015

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有