习题3-7 DNA序列 UVa1368

最新推荐文章于 2021-03-16 22:17:30 发布

原创最新推荐文章于 2021-03-16 22:17:30 发布 · 556 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#NOIP #算法竞赛 #刘汝佳 #入门经典 #DNA序列

算法竞赛入门经典（第2版）第3章数组和字符串专栏收录该内容

18 篇文章

订阅专栏

本文介绍了一道算法竞赛题目UVa1368 DNA序列的解决方法，通过统计每种核苷酸出现次数找出最常出现的核苷酸，并计算次频繁核苷酸的总数。

算法竞赛入门经典（第2版）第3章数组和字符串
习题3-7 DNA序列 UVa1368
感悟。
1、上网站看英文原题。

2、此题比较简单，思路比较清晰，一气呵成，提交AC。

3、前面的训练帮助挺大的。

附上代码

环境Dev-cpp4.9.9.2

#include <stdio.h>
#include <string.h>

char s[50+10][1000+10];
char sout[1000+10];
int a[5];//a[0] A a[1] C a[2] G a[3] T
char sm[5]="ACGT";

void acount(char c){//统计相关字符个数
switch(c){
case 'A':
a[0]++;
break;
case 'C':
a[1]++;
break;
case 'G':
a[2]++;
break;
case 'T':
a[3]++;
break;
}
}

int getmax(){//获取最大字符对应的位置
int max=0;
int i;
for(i=0;i<4;i++){
if(a[max]<a[i]){
max=i;
}
}
return max;
}
int main(){
int T;
int i,j;
int sum;
int m,n;
int max;
scanf("%d",&T);
while(T--){
scanf("%d%d",&m,&n);
for(i=0;i<m;i++){
scanf("%s",s[i]);
}
sum=0;
for(j=0;j<n;j++){
memset(a,0,sizeof(a));
for(i=0;i<m;i++){
acount(s[i][j]);
}
max=getmax();
sout[j]=sm[max];
for(i=0;i<4;i++)
if(i!=max)
sum+=a[i];
}
sout[j]='\0';
printf("%s\n",sout);
printf("%d\n",sum);
}
return 0;
}