DNA序列

题目描述

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

给定一个很长的DNA序列，以及要求的最小子序列长度，研究人员经常会需要在其中找出GC-Ratio最高的子序列。

输入描述:

输入一个string型基因序列，和int型子串的长度

输出描述:

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

输入例子:

AACTGTGCACGACCTGA
5

输出例子:

GCACG

【代码】

1. #include<iostream>
2. #include<string>
3. using namespace std;
4.  
5. int  main()
6. {
7.     string str;
8.     int n;
9.     while(cin>>str>>n)
10.     {  
11.        int max=0;
12.        string s;
13.        string result;
14.        
15.         for(int i=0; i<=str.size()-n; i++)   //注意
16.         {
17.              s = str.substr(i,n);
18.   
19.               int ratio=0;
20.             for(int j=0; j<n; j++)
21.             {
22.                if(s[j]=='G' || s[j]=='C')
23.                     ratio++;
24.             }
25.            
26.             if(max < ratio)
27.               {
28.                    result = s;
29.                      max  = ratio;
30.               }

1.           }
2. 
   
3.       cout<<result<<endl;
4.    }
5.  return 0;
6. }