UVa1368 DNA Consensus string AC

最新推荐文章于 2023-12-02 15:58:17 发布

qicayaya

最新推荐文章于 2023-12-02 15:58:17 发布

阅读量268

点赞数

本文链接：https://blog.csdn.net/qicayaya/article/details/72591092

版权

#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;
}