UVa 1368 DNA Consensus Strings

这题也是题目吓人,写起来不难,文字一大长串,前面都是没用的话,笑话,在中国DNA早就学烂了,题目就是要找一列中重复次数最多的那个字母,如果重复次数相同,就选字典序最小的。直接用map一直都是WA,检查了好久,都没问题,实在忍不住,就重新来过,用了死办法,一下就A了,玄学。
代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<map>
using namespace std;
char a[55][1010];
map<char, int>q;
map<int, char>p;
int nub[5][1010];
int n, m;
int Number(char s[]) {
	int cnt = 0;
	for (int i = 1; i <= m; i++)
		for (int j = 1; j <= n; j++)
			if (s[j] != a[i][j]) cnt++;
	return cnt;
}
int main() {
	int t;
	cin >> t;
	q['A'] = 1; q['C'] = 2; q['G'] = 3; q['T'] = 4;
	p[1] = 'A'; p[2] = 'C'; p[3] = 'G'; p[4] = 'T';
	while (t--) {
		char ans[1010];
		cin >> m >> n;
		memset(nub, 0, sizeof(nub));
		for (int i = 1; i <= m; i++) {
			getchar();
			for (int j = 1; j <= n; j++) {
				scanf("%c", &a[i][j]);
				nub[q[a[i][j]]][j]++;
			}
		}
		for (int j = 1; j <= n; j++) {
			int maxn = 0;
			for (int i = 1; i <= 4; i++) {
				if (maxn < nub[i][j]) {
					ans[j] = p[i];
					maxn = nub[i][j];
				}
			}
		}
		int cnt = Number(ans);
		for (int i = 1; i <= n; i++)
			printf("%c", ans[i]);
		printf("\n%d\n", cnt);
	}
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值