华为机试题63-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。

---

---

解题思路：

这题直接暴力求解就行。

输入字符串和子串长度后，根据子串长度，遍历每一个子串，并统计其含字符'C'和'G'的个数，同时用整型数组存储，整型数组的下标和子串首元素在原串中的下标一致，方便后续对ratio最高的子串进行打印。

代码如下：

#include <stdio.h>
#include <string.h>
#define    N    1000
int main()
{
    char str[N];
    int m,i,j,len,ra[N]={0};
    scanf("%s%d",str,&m);
    len=strlen(str);
    for(i=0;i<len-m;i++)            //i为每一个子串的首元素在原字符串中的下标
    {
        for(j=i;j<i+m;j++)            //遍历子串的每一个字符
        {
            if(str[j]=='C'||str[j]=='G')
                ra[i]++;
        }
    }
    j=0;
    for(i=0;i<len-m;i++)
    {
        if(ra[i]>ra[j])            //统计ratio最大的子串
            j=i;
    }
    for(i=j;i<j+m;i++)
        printf("%c",str[i]);
    return 0;
}