首页 > 其他分享> 文章详细

题解【loj537】「LibreOJ NOIP Round #1」DNA 序列

2019-07-21 10:02:55 阅读：376 来源： 互联网

标签：DNA LibreOJ ++ 题解碱基样例序列 include

题目描述

\(NOIP\)复赛之前\(HSD\)桑进行了一项研究，发现人某条染色体上的一段\(DNA\)序列中连续的\(k\)个碱基组成的碱基序列与做题的 \(AC\) 率有关！于是他想研究一下这种关系。

现在给出一段 \(DNA\) 序列，请帮他求出这段 \(DNA\) 序列中所有连续\(k\)个碱基形成的碱基序列中，出现最多的一种的出现次数。

输入格式

两行，第一行为一段 \(DNA\) 序列，保证 \(DNA\) 序列合法，即只含有 \(A, G, C, T\) 四种碱基；
第二行为一个正整数\(k\)，意义与题目描述相同。

输出格式

一行，一个正整数，为题目描述中所求答案。

样例

样例输入 1

AAAAA
1

样例输出 1

样例解释 1

对于这段 \(DNA\) 序列，连续的\(1\)个碱基组成的碱基序列只有

A，共出现\(5\)次，所以答案为\(5\)。

样例输入 2

ACTCACTC
4

样例输出 2

样例解释 2

对于这段 \(DNA\) 序列，连续的\(4\)个碱基组成的碱基序列为：

\(ACTC, CTCA, TCAC\)与 \(CACT\)。其中 \(ACTC\)出现\(2\)次，其余均出现 \(1\) 次，所以出现最多的次数为\(2\)，即为答案。

数据范围与提示

记 \(DNA\) 序列长度为\(n\)。

本题共\(10\)组数据，只有输出与标准输出一致才可以获得该测试点的分数。

下面给出每组数据的范围和满足性质情况：

性质：给出的 \(DNA\) 碱基序列中每个碱基均相同。
对于所有数据均保证 \(k \leq n\)

题解

这题很容易想到用\(STL\)的\(map\)解决。

然而，经过本人的尝试，\(map\)在这题中只能得到\(80\)分的成绩，出题人会卡。

但是，\(unordered\)_\(map\)在此题中并不会被卡。

所以用\(unordered\)_\(map\)就可以啦\(QwQ\)。

（前提是评测开启\(C++11\)）

代码

#pragma GCC optimize(2)
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cctype>
#include <string>
#include <unordered_map>

using namespace std;

inline int gi()
{
    int f = 1, x = 0; char c = getchar();
    while (c < '0' || c > '9') { if (c == '-') f = -1; c = getchar();}
    while (c >= '0' && c <= '9') { x = x * 10 + c - '0'; c = getchar();}
    return f * x;
}

string s;
int len, k, n, m, ans, js, p[9];
unordered_map <string, int> ap;

int main()
{
    cin >> s;
    k = gi();
    len = s.size();
    bool fl = true;
    for (int i = 1; i < len; i++)
    {
        if (s[i] != s[i - 1]) fl = false; 
    }
    if (fl)
    {
        if (k == 1)
        {
            printf("%d\n", len);
            return 0;
        }
        else
        {
            printf("%d\n", len - k + 1);
            return 0;
        }
    }
    else if (k == 1)
    {
        for (int i = 0; i < len; i++)
        {
            if (s[i] == 'A') ++p[1];
            else if (s[i] == 'G') ++p[2];
            else if (s[i] == 'C') ++p[3];
            else ++p[4];
        }
        printf("%d\n", max(p[1], max(p[2], max(p[3], p[4]))));
        return 0;
    }
    else
    {
        for (int i = 0; i < s.size() - k + 1; i++) ++ap[s.substr(i, k)];
        for (auto it : ap) ans = max(ans, it.second);
        printf("%d\n", ans);
    }
    return 0;
}

标签：DNA,LibreOJ,++,题解,碱基,样例,序列,include
来源： https://www.cnblogs.com/xsl19/p/11220279.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9