uva 1368 - DNA Consensus String

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

GejinZ

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

阅读量997

点赞数

CC 4.0 BY-SA版权

分类专栏：算法竞赛题解 acm-icpc竞赛时光记文章标签： uva

本文链接：https://blog.youkuaiyun.com/u013382399/article/details/37964271

算法竞赛题解同时被 2 个专栏收录

161 篇文章

订阅专栏

acm-icpc竞赛时光记

113 篇文章

订阅专栏

本文介绍了一个用于DNA序列分析的算法实现，通过统计多种碱基出现的频率来确定最优的碱基选择，并计算相应的代价。该算法能有效处理大规模的数据输入，并提供了一种优化的解决方案。

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

题目上给的样例input格式存在一定问题，行末多空格。

在这上面纠结了几次……

#include<stdio.h>
#include<string.h>
#include<iostream>
const int maxn = 1100;
int A[maxn];
int C[maxn];
int G[maxn];
int T[maxn];
char ans[maxn];
int ans_;
int m,n;

void judge(int pos)
{
    ans[pos]='A';
    int max_=A[pos];
    if(max_<C[pos])
    {
        ans[pos]='C';
        max_=C[pos];
    }
    if(max_<G[pos])
    {
//        printf("****\n");
        ans[pos]='G';
        max_=G[pos];
    }
    if(max_<T[pos])
    {
        ans[pos]='T';
        max_=T[pos];
    }

//    printf("ans[%d] = %c\n",pos,ans[pos]);
    if(ans[pos]=='A')
    {
        ans_+=m-A[pos];
    }
    if(ans[pos]=='C')
    {
        ans_+=m-C[pos];
    }
    if(ans[pos]=='G')
    {
        ans_+=m-G[pos];
    }
    if(ans[pos]=='T')
    {
        ans_+=m-T[pos];
    }
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        memset(A,0,sizeof(A));
        memset(C,0,sizeof(C));
        memset(G,0,sizeof(G));
        memset(T,0,sizeof(T));
        char ch;
        scanf("%d%d",&m,&n);
        getchar();
        for(int i=0;i<m;i++)
        {
            for(int j=0;j<n;j++)
            {
                scanf("%c",&ch);
                if(ch=='A')
                {
                    A[j]++;
                }
                if(ch=='C')
                {
                    C[j]++;
                }
                if(ch=='G')
                {
                    G[j]++;
                }
                if(ch=='T')
                {
                    T[j]++;
                }
            }
            getchar();
        }
        //输出查错
//        for(int i=0;i<m;i++)
//        {
//            for(int j=0;j<n;j++)
//            {
//                printf("%c",s[i][j]);
//            }
//            printf("\n");
//        }
//        for(int i=0;i<n;i++)
//            printf(" %d",A[i]);
//        printf("\n");
        //------------------------------------
        memset(ans,0,sizeof(ans));
        ans_=0;
        for(int i=0;i<n;i++)
        {
            judge(i);
        }
        ans[n]=0;
        puts(ans);
        printf("%d\n",ans_);
    }
    return 0;
}