File: string2.cpp

package info (click to toggle)
seqan2 2.5.2-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 228,748 kB
  • sloc: cpp: 257,602; ansic: 91,967; python: 8,326; sh: 1,056; xml: 570; makefile: 229; awk: 51; javascript: 21
file content (42 lines) | stat: -rw-r--r-- 1,185 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
#include <seqan/sequence.h>
#include <seqan/stream.h>

using namespace seqan2;
int main()
{
    // Creating text
    String<char> text = "to be";
    std::cout << text << std::endl;
    appendValue(text, ' ');
    std::cout << "Last sign is whitespace? " << endsWith(text, " ") << std::endl;
    // Erasing whitespaces in text
    eraseBack(text);
    erase(text, 2);
    // Appending another string
    append(text, "ornottobe");
    std::cout << text << std::endl;

    // Pattern
    String<char> pattern = "be";
    // Number of consecutive matching characters per position
    String<int> score;
    resize(score, length(text) - length(pattern) + 1);

    // Brute force pattern matching for every position
    for (unsigned i = 0; i < length(text) - length(pattern) + 1; ++i)
    {
        int localScore = 0;
        for (unsigned j = 0; j < length(pattern); ++j)
            if (text[i + j] == pattern[j])
                ++localScore;
        score[i] = localScore;
    }

    std::cout << "hit at ";
    for (unsigned i = 0; i < length(score); ++i)
        if (score[i] == (int)length(pattern))
            std::cout << i << " ";
    std::cout << std::endl;

    return 0;
}