DNA序列——来自华为OJ平台测试的基础题



/*
 * 一个DNA序列由A/C/G/T四个字母的排列组合组成。G和C的比例（定义为GC-Ratio）是序列中G和C两个字母的总的出现次数除以总的字母数目（也就是序列长度）。
 * 在基因工程中，这个比例非常重要。因为高的GC-Ratio可能是基因的起始点。
 * 给定一个很长的DNA序列，以及要求的最小子序列长度，研究人员经常会需要在其中找出GC-Ratio最高的子序列。
 * 知识点:  字符串
 * 题目来源:  内部整理
 * 练习阶段:  初级
 * 运行时间限制: 10Sec
 * 内存限制: 128MByte
 * 输入: 
 * 输入一个string型基因序列，和int型子串的长度
 * 输出: 
 * 找出GC比例最高的字串
 * 样例输入: 
 * AACTGTGCACGACCTGA
 * 5              
 * 样例输出: 
 * GCACG           
*/
#include <string>
#include <string.h>
#include <iostream>

using namespace std;

int main()
{
 std::string DNASequence;
 std::string NumberStr;
 
 
 int CharNum = 0;
 int MaxGCNum = 0;
 // 输入DNA序列
 getline(std::cin,DNASequence,'\n');

 int DNAlength = static_cast<signed>(DNASequence.size());

 if(DNAlength==0)
 {
  system("pause");
  return -1;
 }

 getline(std::cin,NumberStr,'\n');

 int Numlength = static_cast<signed>(NumberStr.size());

 if((Numlength > sizeof(int))||(Numlength == 0))
 {
  system("pause");
  return -1;
 }

 const int MinDNAlength = atoi(NumberStr.c_str());

 if((MinDNAlength<=0)|| (MinDNAlength>DNAlength))
 {
  system("pause");
  return -1;
 }

 if(DNAlength<MinDNAlength)
 {
  system("pause");
  return -1;
 }

 char* MaxGCDNA = new char[MinDNAlength+1];
 char* GCDNATemp = new char[MinDNAlength+1];
 int BeginNum = CharNum;
 
 int i=0;
 int Tempsize=0;
 while(CharNum<DNAlength)
 { 
  int MaxGCNumtemp = 0;
  Tempsize=0;
  
  for(i=BeginNum;i<BeginNum+MinDNAlength;i++)
  {
   if(i>=CharNum)
    CharNum++;

   char a = DNASequence.at(i);

   if((a=='G') || (a=='C'))
   {
    GCDNATemp[Tempsize] = a;
    MaxGCNumtemp++;
    Tempsize++;
   }
//   else if((a!='T') && (a!='A'))
//   {
//    system("pause");
//    return -1;
//   }
   else
   {
    GCDNATemp[Tempsize] = a;
    Tempsize++;
   }
  }
  GCDNATemp[Tempsize] = '\0';
  if(BeginNum==0)
  {
   strcpy_s(MaxGCDNA,MinDNAlength+1,GCDNATemp);
   MaxGCNum = MaxGCNumtemp;
  }
  if(MaxGCNumtemp>MaxGCNum)
  {
   strcpy_s(MaxGCDNA,MinDNAlength+1,GCDNATemp);
   MaxGCNum = MaxGCNumtemp;
  }
  BeginNum = (CharNum - MinDNAlength) + 1;

 } 

 for(int j=0;j<MinDNAlength;j++)
  std::cout<< MaxGCDNA[j];
 std::cout << endl;

 delete[] MaxGCDNA;
 delete[] GCDNATemp;

 system("pause");
 return 0;
}