java题目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 \le n \le 1000 \1n1000  ,输入的字符串只包含 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。

 

 

 1 import java.io.*;
 2 import java.util.*;
 3 
 4 
 5 /*
 6 思路:从索引0开始遍历整个字符串,找出GC值最大的子字符串
 7 */
 8 public class Main {
 9     public static void main(String[] args) throws IOException {
10         Scanner sc = new Scanner(System.in);
11         while(sc.hasNext()) {
12             String s = sc.next();
13             int n =sc.nextInt();
14             int countGC = 0;
15             String result = null;
16             for(int i =0; i<=s.length()-n;i++){
17                 int tmpcountGc =0;
18                 String s1 = s.substring(i,i+n);  //取子字符串
19                 for(int j =0; j<s1.length();j++) {
20                     if(s1.charAt(j) == 'G' || s1.charAt(j) == 'C') {
21                         tmpcountGc++;  //子字符串的GC次数
22                     }
23                 }
24                 if( tmpcountGc > countGC) {  //如果当前GC大于上一次
25                     countGC = tmpcountGc;  //更新GC次数
26                     result = s1;  //更新要输出的字符串结果
27                 } else if(tmpcountGc == countGC) {   //如果当前GC 和上一次一样
28                     result = result;  //保持输出结果不变
29                 }
30             }
31             System.out.println(result);
32         }
33     }
34 }

 

posted @ 2022-03-12 18:10  海漠  阅读(202)  评论(0编辑  收藏  举报