DNA序列（DNA Consensus String）

最新推荐文章于 2022-01-24 22:36:38 发布

「已注销」

最新推荐文章于 2022-01-24 22:36:38 发布

阅读量670

点赞数

CC 4.0 BY-SA版权

分类专栏： c语言算法

本文链接：https://blog.youkuaiyun.com/qq_42740816/article/details/86520139

这篇博客探讨了如何找出DNA序列的共识字符串，并计算汉明距离。通过使用二维数组压缩方法，找出每列出现最频繁的字符，然后组合成共识字符串。汉明距离则是逐列计算与共识字符串不同的位置数量，最终求和。文中提到，当存在多种解决方案时，输出的将是字典序最小的。内容受到了其他博客的启发。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

/*DNA序列*/
#include<stdio.h>
#include<string.h>
int main ()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		int m,n;
		scanf("%d%d",&m,&n);
		getchar();
		char dna[55][1005];
		memset(dna,0,sizeof(dna));
		for(int i=0;i<m;i++)
		{
			for(int j=0;j<n;j++)
				dna[i][j]=getchar();
			getchar();
		}
	/*	printf("-----------------\n");
		for(int i=0;i<m;i++)
		{
			for(int j=0;j<n;j++)
			   printf("%c",dna[i][j]);
			printf("\n");
		}
		printf("-------------\n");
	*/	
		int acgt[55][1005]={0};
		for(int i=0;i<m;i++)
		{
			for(int j=0;j<n;j++)
			{
				switch(dna[i][j])
				{
					case'A':acgt[0][j]++;break;
					case'C':acgt[1][j]++;break;
					case'G':acgt[2][j]++;break;
					case'T':acgt[3][j]++;break;
				}
				
			}
		}
		char maxx[1005]={0};
		for(int j=0;j<n;j++)
		{
			
			int t=0;
			for(int i=0,max=acgt[j][0];i<4;i++)
			{
				if(max<acgt[i][j])
				{
					max=acgt[i][j];
					t=i;
				}
	//