uva1368 - DNA Consensus String

题目大意:给定m个长度为n的DNA序列,求一个DNA序列,使其到所有DNA序列的hamming距离尽量小,如有多组解,输出字典序最小的。(hamming距离定义为两个等长的字符串,字符不同的位置个数)。

题目思路:找hamming距离最小的字符串,即统计每个位置上A,C,T,G最多的字符组成的字符串。


#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define MAX 1010
#define INF 0xfffffff
char map[51][MAX];
char s[MAX];
char ans[]="ACGT";
int n,m,cnt;
int check(int x)
{
    int i,j,sum=-1;
    int a[5];
    memset(a,0,sizeof(a));
    for (i=0;i<n;i++)
    {
        if (map[i][x]=='A')
        {
           a[0]++;                   
        }
        else if (map[i][x]=='C')
        {
           a[1]++;        
        }
        else if (map[i][x]=='G')
        {
           a[2]++;       
        }
        else if (map[i][x]=='T')
        {
           a[3]++;     
        }
    }
    for (i=0;i<4;i++)
    {
        if (sum<a[i])
        {
           sum=a[i];
           j=i;             
        }
    }
    s[x]=ans[j];
    return n-sum;
}
int main()
{
    int T;
    scanf("%d",&T);
    while (T--)
    {
          int i,j;
          scanf("%d%d",&n,&m); 
          cnt=0; 
          for (i=0;i<n;i++)
          {
              scanf("%s",map[i]);    
          }
          for (i=0;i<m;i++)
          {
              cnt+=check(i);       
          }
          s[m]='\0';
          printf("%s\n%d\n",s,cnt);
    }
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值