UVA 1368 - DNA Consensus String

本文介绍了一个简单的DNA序列匹配算法,该算法旨在从多个给定的DNA序列中找到一个与所有序列差异最小的新序列。通过统计每个位置上出现次数最多的碱基来构建目标序列,并计算总差异数。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

/*

字符串类的水题, 就是给m个DNA序列, 让你找出一个DNA序列 这个序列要求 跟给的m个序列不同之处最小,也就是最相似

*/

#include<cstdio>

#include<cstring>


char DNA[51][1010],s[1100];
int pos[4]= {0,2,6,19},m,n,d[1100];
int main()
{
    int t,pd[26];
    scanf("%d",&t);
    while(t--)
    {
        memset(s,0,sizeof(s));
        memset(d,0,sizeof(d));
        scanf("%d %d",&m,&n);
        for(int i = 0; i < m; i++)
            scanf("%s",DNA[i]);
        for(int j = 0; j < n; j++)
        {
            memset(pd,0,sizeof(pd));
            for(int i = 0; i < m; i++)
            {
                pd[DNA[i][j]-'A']++;
            }
            int max = 0;
            for(int i = 0; i < 4; i++)
            {
                if(max < pd[pos[i]])
                {
                    max = pd[pos[i]];
                    s[j] = pos[i]+'A';
                }
            }
            d[j] = max;
        }
        puts(s);
        int temp=0;
        for(int i = 0; i < n; i++)
            temp += d[i];
        printf("%d\n",m*n-temp);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值