Uva-1368-DNA Consensus String

150 篇文章 0 订阅

这个题要求求出一个DNA序列使得其到各序列的距离最小,如果答案有多组找字典序最小的。

依次找到最优解就好。

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int maxn=10e4+10;
char map[60][maxn];
int main()
{
    int cas;
    scanf("%d",&cas);
    while(cas--)
    {
	int n,m;
	scanf("%d%d",&n,&m);
	for(int i=0;i<n;i++)
	    scanf("%s",map[i]);
	char str[maxn];
	int ans=0;
	for(int i=0;i<m;i++)
	{
	    int ita=0,itb=0,itc=0,itd=0;
	    for(int j=0;j<n;j++)
	    {
		switch(map[j][i])
		{
		    case 'A': ita++;break;
		    case 'C': itb++;break;
		    case 'G': itc++;break;
		    case 'T': itd++;break;
		}
	    }
	    if(ita>=itb&&ita>=itc&&ita>=itd)
	    {	str[i]='A';ans+=n-ita;}
	    else if(itb>=ita&&itb>=itc&&itb>=itd)
	    {	str[i]='C';ans+=n-itb;}
	    else if(itc>=ita&&itc>=itb&&itc>=itd)
	    {	str[i]='G';ans+=n-itc;}
	    else
	    {	str[i]='T';ans+=n-itd;}
	}
	str[m]='\0';
	printf("%s\n%d\n",str,ans);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值