您的位置:首页 > 其它

UVa1368 DNA Consensus string AC

2017-05-20 23:55 483 查看
#include "stdio.h"
#include "string.h"
#include "stdlib.h"
#define MAX(a, b) ((a)>(b)?(a):(b))
int main()
{
int T,m,n,Hamming=0,ACGTMAX=0;
char s[50][1000],ACGT[4];
scanf("%d",&T);
while (T--)
{
Hamming=0;
memset(ACGT,0,sizeof(ACGT));
scanf("%d%d",&m,&n);//DNA行数,DNA长度
getchar();//取回车
for (int i=0;i<m;i++)	gets(s[i]);//存入表中
for (int j=0;j<n;j++)
{
for (int i=0;i<m;i++)
{
switch (s[i][j])
{
case 'A':ACGT[0]++;
break;
case 'C':ACGT[1]++;
break;
case 'G':ACGT[2]++;
break;
case 'T':ACGT[3]++;
break;
}
}
ACGTMAX=MAX(MAX(MAX(ACGT[0], ACGT[1]), ACGT[2]),ACGT[3]);
if (ACGTMAX==ACGT[0])
{putchar('A');Hamming+=ACGT[1]+ACGT[2]+ACGT[3];}
else if (ACGTMAX==ACGT[1])
{putchar('C');Hamming+=ACGT[0]+ACGT[2]+ACGT[3];}
else if (ACGTMAX==ACGT[2])
{putchar('G');Hamming+=ACGT[0]+ACGT[1]+ACGT[3];}
else if (ACGTMAX==ACGT[3])
{putchar('T');Hamming+=ACGT[0]+ACGT[1]+ACGT[2];}
memset(ACGT,0,sizeof(ACGT));
}
printf("\n%d\n",Hamming);
}
system("pause");
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: