#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;
}
J - DNA Consensus String
猜你喜欢
转载自blog.csdn.net/qq_43328040/article/details/88367872
今日推荐
周排行