LA - 3602 - DNA Consensus String

本文探讨了一个关于ACGT序列的问题,目标是在给定的序列集合中找到一个新序列,使得其与所有已知序列在指定位置上的距离和达到最小。详细介绍了算法步骤,包括统计每个字符的出现频率并选择出现次数最多的字符作为新序列对应位置的字符。

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

题意:有m个长度为n的ACGT序列,求一个长度为n的序列,使得每一位与这m个序列相应位的距离和最小(相同为0,不同为1)(4 <= m <= 50, 4 <= n <= 1000)。

题目链接:https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&category=19&problem=1603

——>>统计每一位A、C、G、T出现的次数,取出现最多的。

#include <cstdio>
#include <cstring>

using namespace std;

const int maxn = 1000 + 10;

int main()
{
    int A[maxn], C[maxn], G[maxn], T[maxn], t, m, n, i, j;
    char s[maxn];
    scanf("%d", &t);
    while(t--){
        scanf("%d%d", &m, &n);
        memset(A, 0, sizeof(A));
        memset(C, 0, sizeof(C));
        memset(G, 0, sizeof(G));
        memset(T, 0, sizeof(T));
        for(j = 0; j < m; j++){
            scanf("%s", s);
            for(i = 0; i < n; i++){
                switch(s[i])
                {
                    case 'A':{
                        A[i]++;
                        break;
                    }
                    case 'C':{
                        C[i]++;
                        break;
                    }
                    case 'G':{
                        G[i]++;
                        break;
                    }
                    case 'T':{
                        T[i]++;
                        break;
                    }
                }
            }
        }
        char ret[maxn];
        int cnt = 0;
        for(i = 0; i < n; i++){
            int Max = -1, num;
            if(A[i] > Max){
                Max = A[i];
                ret[i] = 'A';
                num = A[i];
            }
            if(C[i] > Max){
                Max = C[i];
                ret[i] = 'C';
                num = C[i];
            }
            if(G[i] > Max){
                Max = G[i];
                ret[i] = 'G';
                num = G[i];
            }
            if(T[i] > Max){
                Max = T[i];
                ret[i] = 'T';
                num = T[i];
            }
            cnt += m - num;
        }
        for(i = 0; i < n; i++) putchar(ret[i]);
        printf("\n");
        printf("%d\n", cnt);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值