【华为机试】HJ63 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比例最高的子串,如果有多个则输出第一个的子串

输入输出示例：

输入：

ACGT
2

输出：

CG

我的代码：

#include <iostream>
using namespace std;

int main() {
    string s;
    int n;
    cin >> s;
    cin >> n;
    int R = 0;
    int start = 0;
    //判别前n个
    for (int i = 0; i < n; i++) {
        if (s[i] == 'C' || s[i] == 'G') {
            R++;
        }
    }

    int Rt = R; //当前块的GC数
    for (int i = n - 1; i < s.length() - 1; i++) {
        //移动块并计算GC数
        if (s[i - n + 1] == 'C' || s[i - n + 1] == 'G') {
            Rt--;
        }

        if (s[i + 1] == 'C' || s[i + 1] == 'G') {
            Rt++;
        }
        // 更新最大块的数与位置
        if (Rt > R) {
            R = Rt;
            start = i - n + 2;
        }

    }
    for (int i = start; i < start + n; i++) {
        cout << s[i];
    }

}
// 64 位输出请用 printf("%lld")

我的代码分析：

监控一个大小为n的块在基因链上，用start来记录最大GC-radio的起始点，R记录最大GC数，遍历整个基因链，不断更新start和R，最后输出结果。

top1代码：

#include<iostream>
using namespace std;
int main(){
    string s;
    int n;
    string result;
    while(cin>>s>>n){
        float r=0;
        if(n==s.size()) cout<<s<<endl;
        for(int i=0;i<s.size()-n;i++){
            string temp=s.substr(i,n);
            int count=0;
            for(int k=0;k<n;k++){
              if(temp[k]=='G'||temp[k]=='C') count++;
            }
            if(r<(float)count/(float)n){
                r=(float)count/(float)n;
                //result="";
                result=temp;
            }
         }
        cout<<result;
    }
    return 0;
}

top1代码分析：

思路与我大致相同

总结：

无