zoj 3132 DNA Consensus String(字符串处理)

最新推荐文章于 2022-12-08 16:04:17 发布

溺水行舟

最新推荐文章于 2022-12-08 16:04:17 发布

阅读量1.2k

点赞数

CC 4.0 BY-SA版权

分类专栏：训练指南（刘汝佳）

本文链接：https://blog.youkuaiyun.com/liang5630/article/details/8211710

训练指南（刘汝佳）专栏收录该内容

25 篇文章

订阅专栏

本文介绍了一个C++程序，用于分析一组DNA序列数据。该程序通过计数每种碱基(A、C、G、T)出现的次数来确定最可能的序列，并计算所需变异数量以达到最优序列。适用于生物信息学中的序列比对问题。

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

#include <iostream>
#include<cstdio> 
using namespace std;
struct DNA
{
	int a,c,g,t;
}d[1005];
int main(int argc, char *argv[])
{
	int t,n,m,i,j,ans,max;
	char s[1005];
	cin>>t;
	while(t--)
	{   scanf("%d%d",&n,&m); 
		for(i=0;i<m;i++)
			d[i].a=d[i].c=d[i].g=d[i].t=0;
		for(i=0;i<n;i++)
		{scanf("%s",&s);
		 for(j=0;j<m;j++)
		 if(s[j]=='A') d[j].a++;
		 	else if(s[j]=='C') d[j].c++;
		 		else if(s[j]=='G') d[j].g++;
		 				else d[j].t++;
		}
		for(ans=i=0;i<m;i++)
		{   
			s[i]='A';max=d[i].a;
			if(d[i].c>max) s[i]='C',max=d[i].c;
			if(d[i].g>max) s[i]='G',max=d[i].g;
			if(d[i].t>max) s[i]='T',max=d[i].t;
			ans+=n-max;
		}
		s[m]=0;
		cout<<s<<endl<<ans<<endl;	
	}
	return 0;
}