原代码:
#include <stdio.h> #include <stdlib.h> #include <string.h> int main() { int T,m,n,i,j,max,t,count,k; char ch[100][1500],sh[1500]; scanf("%d",&T); while(T--) { count=0; scanf("%d %d",&m,&n); for(i=0;i<m;i++) scanf("%s",ch[i]); for(j=0;j<n;j++) { int c[4]={0}; for(i=0;i<m;i++) { if(ch[i][j]=='A') c[0]++; if(ch[i][j]=='C') c[1]++; if(ch[i][j]=='G') c[2]++; if(ch[i][j]=='T') c[3]++; } max=c[0]; for(k=1;k<4;k++) if(c[k]>max) { t=max;max=c[k];c[k]=max; } if(max==c[0]) sh[j]='A'; else if(max==c[1]) sh[j]='C'; else if(max==c[2]) sh[j]='G'; else if(max==c[3]) sh[j]='T'; } for(i=0;i<m;i++) for(j=0;j<n;j++) {if(ch[i][j]!=sh[j]) count++;} printf("%s\n%d\n",sh,count); } return 0; }
其运行结果:
假如上组示例字符串的长度长,那么由于sh挨个输出时,没有在该组应该停止的那个位置停止输出,因为它没有在sh[j]处遇到"\0",所以对于每组的sh[j]都要赋值其为0.
DNA Consensus String
最新推荐文章于 2022-07-27 22:31:25 发布