sicily--1035. DNA matching

  1. 水题
  2. 用一个map 容器储存系统所给的“核苷酸链”,以及其出现的次数; 
  3. 遍历map 容器, 要配对成功,则必须在系统给出的“核苷酸链”即我们的map 中有存在, 而可以组成最多的链数就要看这互补的两条链出现的次数最少的那一方
  4. 因为在遍历map时重复计算了同一条“双螺旋核苷酸链”, 顾最后结果需要除以2

#include<iostream>
#include<map>
#include<string>
using namespace std;

int main()
{
	int caseNum;
	cin >> caseNum;
	while(caseNum--)
	{
		map<string, int> DNA;//string 为核苷酸序列, int 为出现的次数
		map<string, int>::iterator it;

		int dnaNum;
		cin >> dnaNum;
		while(dnaNum--)
		{
			string temp;
			cin >> temp;
			it = DNA.find(temp);//检查是否已经存在过了
			if(it == DNA.end())//第一次出现
			{
				DNA[temp] = 1;
			}
			else//已经出现过
			{
				(*it).second++;//次数加一
			}
		}

		int sum = 0;
		for(it = DNA.begin(); it != DNA.end(); it++)
		{
			map<string, int>::iterator itTemp;
			string goal = "";
			for(int i = 0; i < (*it).first.length(); i++)
			{
				switch((*it).first[i])//找出与其互补的链
				{
				case 'A':
					goal.push_back('T');
					break;
				case 'T':
					goal.push_back('A');
					break;
				case 'G':
					goal.push_back('C');
					break;
				case 'C':
					goal.push_back('G');
					break;
				}
			}

			itTemp = DNA.find(goal);
			if(itTemp != DNA.end())//能找到这样的链
			{
				if((*it).second >= (*itTemp).second)//加最小的值
				{
					sum = sum + (*itTemp).second;
				}
				else
				{
					sum = sum + (*it).second;
				}
				//DNA.erase(it);//防止再次计算
				//DNA.erase(itTemp);
			}
			else//找不到
			{
				continue;
			}
		}
		cout << sum / 2 << endl;//因为重复计算了两次
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值