J - DNA Consensus String

#include<stdio.h>
#include<string.h>
int main(void)
{
    int t;
    int n,m;
    char a;
    int b[1005][5];   //表示每列中A C G T的个数
    scanf("%d",&t);
    while(t--){
        memset(b,0,sizeof(b));
        scanf("%d%d",&n,&m);
        getchar();
        for(int i=0;i<n;i++){
            for(int j=0;j<m;j++){
                scanf("%c",&a);
                //记录每列A C G T出现的次数
                if(a=='A')b[j][0]++;       
                else if(a=='C')b[j][1]++;  
                else if(a=='G')b[j][2]++;
                else if(a=='T')b[j][3]++;
            }
            getchar();
        }
        int sum=0;
        for(int i=0;i<m;i++){
            int max=-1000000;
            int p=0;
            for(int j=0;j<4;j++){
                if(b[i][j]>max){  //max表示第i列中出现的字符次数最多的值是多少
                    max=b[i][j];
                    p=j;     //记录该字符
                }
            }
            sum+=n-max;    //加上每一列的距离的最小值
            if(p==0)printf("A");
            else if(p==1)printf("C");
            else if(p==2)printf("G");
            else if(p==3)printf("T");
        }
        printf("\n");
        printf("%d\n",sum);
    }
    return 0;
}

猜你喜欢

转载自blog.csdn.net/qq_43328040/article/details/88367872
DNA