DNA匹配问题

描述

DNA（脱氧核糖核酸）为人体中存储基因信息的媒介。核苷酸是一种由亚基构成，排列起来聚合物链。DNA聚合物链也称DNA线。
根据化学成份或获取方式，DNA中的核苷酸基因可以分为四种：adenine，guanine, cytosine and thymine，我们分别以A，G，C和T表示。单个核苷酸通过端到端的化学行为链接成DNA，简单地说，我们可以作用一系列A，T，C和G构成的串代表一个链，如ATTCGAC，但我们应该注意核苷酸组成的串有方向性，所以我们认为ATTCGAC和CAGCTTA代表两种不同的链。
基因通常不会以自由单链的形式存在。在适当条件下单链将成对出现，并绕着对方，形成了著名的双螺旋结构。这种配对出现于As与Ts，Gs与Ts，因为相互吸引，称为氢键。因此称A/T和G/C为互补碱基对。
Molecular生物实验室是进行DNA实验的，其一项重要的工作是配比两个碱基对线，制成DNA双链。我们要求两个单碱基对线的长度必须相同，同时要求两个单链中相同位置的核苷酸基因为互补碱基对。例如ATTCGAC和TAAGCTG是互补，但CAGCTTA和TAAGCTG不是，ATTCGAC和GTAAGCT也不是。
作为生物研究的助手的你，你的老板交给你一个任务：给定n个单链，找出可以配对双链的最大数目（当然，每个单链最多只能使用一次）。如果n很小时，可以很容易用笔算出，但当n很大时，就不行了。现在要求使用编写程序实现。

输入

输入包括多个测试用例，第一行是一个正整数T(T<=20)，表示测试用例个数。在每一个测试用例中，第一行是一个正整数n(n<=10000)，表示单链的个数。接下来的n行，每行代表一个单链。每一个单链的长度不超过1000

输出

对每一个测试用例，用一行输出一个整数，即可以配成的最大双链的个数。

样例输入

 

2
3
ATCG
TAGC
TAGG
2
AATT
ATTA

 

样例输出

 

1
0

 

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

char a[10001][1010];

int main()
{
    int n, m, i, j;
    bool use[100010];
    cin>>n;
    while(n--)
    {
        memset(use,true,sizeof use);
        cin>>m;
        int sum = 0;getchar();
        for(i= 0; i < m; i++)
            gets(a[i]);
        for(i = 0; i < m-1; i++)
        {
            for(j = 0; a[i][j]!='\0'; j++)
            {
                if(a[i][j] == 'T')
                    a[i][j] = 'A';
                else if(a[i][j] == 'A')
                    a[i][j] = 'T';
                else if(a[i][j] == 'C')
                    a[i][j] = 'G';
                else if(a[i][j] == 'G')
                    a[i][j] = 'C';
            }
            for(j = i+1; j < m; j++)
            {
                if(strcmp(a[i],a[j])==0 && use[i] && use[j])
                {
                    sum++; 
                    use[i] = false;
                    use[j] = false;
                }
            }
        }
        cout<<sum<<endl;
    }
    return 0;
}

 

 

转载于:https://www.cnblogs.com/jianqian/p/4975506.html