UVA - 1368 DNA Consensus String :简单字符串处理

最新推荐文章于 2021-12-30 15:07:01 发布

_N_ull_

最新推荐文章于 2021-12-30 15:07:01 发布

阅读量346

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/qq_17249047/article/details/51615551

算法竞赛入门经典（刘汝佳）代码抄写&&泛读专栏收录该内容

23 篇文章

订阅专栏

本文介绍了一种用于DNA序列比对的算法实现，通过统计每条输入DNA序列中A、G、C、T四种碱基的数量，并找出每种碱基在对应位置上出现次数最多的一个，从而构建出一个最相似的参考序列。

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

思路：用4个数组存放每个字符串中对应元素下标的A G C T 的个数

比如，所给样例①中 A[0]=1，G[0]=0，C[0]=0, TT[0]=4；

这样，在 A G C T 对应下标寻找出现次数最多的，同时为最终字符串贡献一个字符；

由于按照ACGT的顺序比较保证了出现次数相同的情况下字典序最小的那个字符作为最终输出字符串的贡献字符；

                    #define LOCAL
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
using namespace std;
const int Size = 1000;
int A[Size+5];
int G[Size+5];
int C[Size+5];
int TT[Size+5];

int main()
{
#ifdef LOCAL
    freopen("in1.txt","r",stdin);
#endif
    int T, m,n;
    string S;
    scanf("%d",&T);
    while(T--){
        scanf("%d %d",&m,&n);
        memset(A,0,sizeof(A));
        memset(G,0,sizeof(G));
        memset(C,0,sizeof(C));
        memset(TT,0,sizeof(TT));
        int AnsN=0; string AnsS="";
        int M=m;
        while(m--){
            cin>>S;
            for(int i=0; i<n; i++){
                if(S[i]=='A') A[i]++;
                if(S[i]=='G') G[i]++;
                if(S[i]=='C') C[i]++;
                if(S[i]=='T') TT[i]++;
            }
        }
        int Max=-1;
        char MaxX;
        for(int i=0; i<n; i++){
            Max = A[i]; MaxX = 'A';
            if(C[i]>Max) { Max=C[i]; MaxX='C'; }
            if(G[i]>Max) { Max=G[i]; MaxX='G'; }
            if(TT[i]>Max) { Max=TT[i]; MaxX='T'; }
            AnsN+=(M-Max); AnsS+=MaxX;
        }
        cout<<AnsS<<endl;
        printf("%d\n",AnsN);
    }//while
    return 0;
}

额···表述能力果然很烂····