LA 3602 DNA Consensus String .

题目地址:http://vjudge.net/problem/UVALive-3602

第二次写,二周目

#include <bits/stdc++.h>
using namespace std;
#define REP(i,a,b) for(int i=a;i<=(b);++i)
char A[50+1][1000+3],Ans[1000+3];
std::map<char, int> idx;
const char DNA[]="ACGT";
int main(int argc, char const *argv[])
{
	int T; scanf("%d",&T);
	idx['A']=0;idx['C']=1;idx['G']=2;idx['T']=3;  //为了字典序
	while(T--){
		int n,m;  char ch;
		scanf("%d%d",&n,&m);
		REP(i,1,n) scanf("%s",A[i]+1);
		int cnt=0;
		REP(j,1,m){
			int B[4]={0};
			REP(i,1,n) B[idx[A[i][j]]]++;
			int sum=0,n=0;  
			REP(i,0,3) {
				sum+=B[i];
				if(n<B[i]) Ans[j]=DNA[i],n=B[i];
			}
			cnt+=sum-n;
		}
		Ans[m+1]='\0'; 
		printf("%s\n%d\n", Ans+1,cnt);
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值