UVA - 1368 DNA Consensus String

题目大意:计算这一组碱基序列的共同序列碱基以及出错数(不大会翻译)。一组碱基序列中的某一位出现最多次数的碱基为要输出的串在这一位的值,如果碱基出现次数相同则按字母顺序选择。要求输出的数字为最小修改的次数。

解题思路:先统计在这组串的某个位置四个碱基的出现次数,根据次数判断。

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
int main () {
  int t , m , n , coer , i , j , k , nn[4] , nt;
  char g[55][1005],costr[1005];
  cin >> t;
  while( t--) {
    scanf("%d%d",&m,&n);
    coer = 0;
    for( i = 0; i < m ; i++) 
       cin >> g[i];
       for( j = 0; j < n ; j++) {
          nn[0] = 0 ; nn[1] = 0 ; nn[2] = 0 ; nn[3] = 0;
         for( i = 0; i < m ; i++) {
           if( g[i][j] == 'A') nn[0]++;
           if( g[i][j] == 'C') nn[1]++;
           if( g[i][j] == 'G') nn[2]++;
           if( g[i][j] == 'T') nn[3]++;
           }
         int  cho = nn[0];
         nt = 0;
         for( k = 1; k < 4 ; k++) {
            if ( cho < nn[k] ) {   
               cho = nn[k];
               nt = k;
             }
         }
         coer +=( m-cho);
         switch(nt) {
           case 0: costr[j] = 'A' ;break;
           case 1: costr[j] = 'C' ;break;
           case 2: costr[j] = 'G' ;break;
           case 3: costr[j] = 'T' ;break;
         }
         costr[j+1] = '\0';
     } 
      cout << costr <<endl<<coer<<endl;

  } 
return 0;
}

看这道题的题目真被吓了一跳,其实并不难哈哈。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值