UVA - 1368 DNA Consensus String

最新推荐文章于 2024-04-30 19:34:36 发布

原创最新推荐文章于 2024-04-30 19:34:36 发布 · 187 阅读

0 ·

CC 4.0 BY-SA版权

UVA 专栏收录该内容

8 篇文章

订阅专栏

本文提供了一段C++代码，用于解决UVA-1368问题——寻找一组DNA序列的共识字符串，并计算变异总数。通过读取多组输入数据，程序能够找出每列中出现最多的碱基并输出最终的共识字符串及其变异数。

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

UVA - 1368 DNA Consensus String

#include"cstdio"
#include"string.h"
int m, n;
int sum;
char map[50][1000];
int cnt[4][1000];
char alp[] = { 'A','C','G','T' };
int main()
{
    int t; scanf("%d", &t);
    while (t--)
    {
        memset(cnt, 0, sizeof(cnt));
        sum = 0;
        scanf("%d%d", &m, &n);
        getchar();
        for (int i = 0; i < m; ++i)
        {
            for (int j = 0; j < n; ++j)
            {
                char c=fgetc(stdin);
                map[i][j] = c;
                if (c == 'A')++cnt[0][j];
                else if (c == 'C')++cnt[1][j];
                else if (c == 'G')++cnt[2][j];
                else ++cnt[3][j];
            }
            getchar();
        }
        for (int j = 0; j < n; ++j)
        {
            int maxi = 0;
            int max = cnt[0][j];
            for (int i = 1; i < 4; ++i)
            {
                if (max < cnt[i][j])
                {
                    maxi = i; max = cnt[i][j];
                }
            }
            sum += m - max;
            printf("%c", alp[maxi]);
        }
        printf("\n%d\n", sum);
    }
    return 0;
}