Week 4
Week 4
Week 4
Programming Assignment 3:
Algorithmic Challenges
Revision: September 5, 2016
Introduction
Welcome to your third programming assignment of the Algorithms on Strings class! In this programming
assignment, you will be practicing implementing very efficient string algorithms.
Learning Outcomes
Upon completing this programming assignment you will be able to:
1. find all occurrences of a pattern in text;
Contents
1 Problem: Find All Occurrences of a Pattern in a String 3
4 Advanced Problem: Construct the Suffix Tree from the Suffix Array 10
1
5 General Instructions and Recommendations on Solving Algorithmic Problems 13
5.1 Reading the Problem Statement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
5.2 Designing an Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
5.3 Implementing Your Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
5.4 Compiling Your Program . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
5.5 Testing Your Program . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
5.6 Submitting Your Program to the Grading System . . . . . . . . . . . . . . . . . . . . . . . . . 15
5.7 Debugging and Stress Testing Your Program . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2
1 Problem: Find All Occurrences of a Pattern in a String
Problem Introduction
In this problem, we ask a simple question: how many times one string occurs as a substring of another?
Recall that different occurrences of a substring can overlap with each other. For example, ATA occurs three
times in CGATATATCCATAG.
Problem Description
Task. Find all occurrences of a pattern in a string.
Input Format. Strings 𝑃 𝑎𝑡𝑡𝑒𝑟𝑛 and 𝐺𝑒𝑛𝑜𝑚𝑒.
Constraints. 1 ≤ |𝑃 𝑎𝑡𝑡𝑒𝑟𝑛| ≤ 106 ; 1 ≤ |𝐺𝑒𝑛𝑜𝑚𝑒| ≤ 106 ; both strings are over A, C, G, T.
Output Format. All starting positions in 𝐺𝑒𝑛𝑜𝑚𝑒 where 𝑃 𝑎𝑡𝑡𝑒𝑟𝑛 appears as a substring (using 0-based
indexing as usual).
Time Limits.
language C C++ Java Python C# Haskell JavaScript Ruby Scala
time in seconds 1 1 12 4 1.5 2 4 4 24
Explanation:
The pattern is longer than the text and hence has no occurrences in the text.
Sample 2.
Input:
ATA
ATATA
Output:
02
Explanation:
The pattern appears at positions 1 and 3 (and these two occurrences overlap each other).
Sample 3.
Input:
ATAT
GATATATGCATATACTT
Output:
139
Explanation:
The pattern appears at positions 1, 3, and 9 in the text.
3
Starter Files
The starter solutions for this problem read the input data from the standard input, pass it to a blank
procedure, and then write the result to the standard output. You are supposed to implement your algorithm
in this blank procedure if you are using C++, Java, or Python3. For other programming languages, you need
to implement a solution from scratch. Filename: kmp
What To Do
To solve this problem, it is enough to implement carefully the corresponding algorithm covered in the lectures.
Need Help?
Ask a question or see the questions asked by other learners at this forum thread.
4
2 Problem: Construct the Suffix Array of a Long String
Problem Introduction
The goal in this problem is to construct the suffix array of a given string again, but this time for a longer
string. This will require you to implement an efficient algorithm.
Problem Description
Task. Construct the suffix array of a string.
Input Format. A string Text ending with a “$” symbol.
Constraints. 1 ≤ |Text(Text)| ≤ 2 · 105 ; except for the last symbol, Text contains symbols A, C, G, T only.
Output Format. SuffixArray(Text), that is, the list of starting positions of sorted suffixes separated by
spaces.
Time Limits.
language C C++ Java Python C# Haskell JavaScript Ruby Scala
time in seconds 1 1 4 30 1.5 2 30 30 8
5
Sample 3.
Input:
GAGAGAGA$
Output:
875316420
Explanation:
Sorted suffixes:
8 $
7 A$
5 AGA$
3 AGAGA$
1 AGAGAGA$
6 GA$
4 GAGA$
2 GAGAGA$
0 GAGAGAGA$
Sample 4.
Input:
AACGATAGCGGTAGA$
Output:
15 14 0 1 12 6 4 2 8 13 3 7 9 10 11 5
Explanation:
Sorted suffixes:
15 $
14 A$
0 AACGATAGCGGTAGA$
1 ACGATAGCGGTAGA$
12 AGA$
6 AGCGGTAGA$
4 ATAGCGGTAGA$
2 CGATAGCGGTAGA$
8 CGGTAGA$
13 GA$
3 GATAGCGGTAGA$
7 GCGGTAGA$
9 GGTAGA$
10 GTAGA$
11 TAGA$
5 TAGCGGTAGA$
Starter Files
The starter solutions for this problem read the input data from the standard input, pass it to a blank
procedure, and then write the result to the standard output. You are supposed to implement your algorithm
in this blank procedure if you are using C++, Java, or Python3. For other programming languages, you need
to implement a solution from scratch. Filename: suffix_array_long
What To Do
To solve this problem, it is enough to implement carefully the corresponding algorithm covered in the lectures.
6
Need Help?
Ask a question or see the questions asked by other learners at this forum thread.
7
3 Problem: Pattern Matching with the Suffix Array
Problem Introduction
In this problem, we will let you use the suffix array to solve the Multiple Pattern Matching Problem.
Problem Description
Task. Find all occurrences of a given collection of patterns in a string.
Input Format. The first line contains a string Text). The second line specifies an integer 𝑛. The last line
gives a collection of 𝑛 strings Patterns = {𝑝1 , . . . , 𝑝𝑛 } separated by spaces.
∑︀𝑛
Constraints. All strings contain symbols A, C, G, T only; 1 ≤ |Text| ≤ 105 ; 1 ≤ 𝑛 ≤ 104 ; 𝑖=1 |𝑝𝑖 | ≤ 105 .
Output Format. All starting positions (in any order) in Text where a pattern appears as a substring (using
0-based indexing as usual). If several patterns occur at the same position of the Text, still output this
position only once.
Time Limits.
language C C++ Java Python C# Haskell JavaScript Ruby Scala
time in seconds 1 1 4 12 1.5 2 12 12 8
Explanation:
There are no occurrences of the patterns in the text
Sample 3.
Input:
ATATATA
3
ATA C TATAT
Output:
4201
Explanation:
The pattern ATA appears at positions 0, 2, and 4, the pattern TATAT appears at position 1.
8
Starter Files
The starter solutions for this problem read the input data from the standard input, pass it to a blank
procedure, and then write the result to the standard output. You are supposed to implement your algorithm
in this blank procedure if you are using C++, Java, or Python3. For other programming languages, you need
to implement a solution from scratch. Filename: suffix_array_matching
What To Do
To solve this problem, it is enough to implement carefully the corresponding algorithm covered in the lectures.
Need Help?
Ask a question or see the questions asked by other learners at this forum thread.
9
4 Advanced Problem: Construct the Suffix Tree from the Suffix
Array
We strongly recommend you start solving advanced problems only when you are done with the basic problems
(for some advanced problems, algorithms are not covered in the video lectures and require additional ideas
to be solved; for some other advanced problems, algorithms are covered in the lectures, but implementing
them is a more challenging task than for other problems).
Problem Introduction
SuffixTree(Text) can be constructed in linear time from SuffixArray(Text) by using the longest common
prefix (LCP) array of Text, LCP(Text), which stores the length of the longest common prefix shared by
consecutive lexicographically ordered suffixes of Text. For example,
Problem Description
Task. Construct a suffix tree from the suffix array and LCP array of a string.
Input Format. The first line contains a string Text ending with a “$” symbol, the second line contains
SuffixArray(Text) as a list of |Text| integers separated by spaces, the last line contains LCP(Text) as
a list of |Text| − 1 integers separated by spaces.
Constraints. 1 ≤ |Text(Text)| ≤ 2 · 105 ; except for the last symbol, Text contains symbols A, C, G, T only.
Output Format. The output format in this problem differs from the output format in the problem “Suffix
Tree” from the Programming Assignment 2 and is somewhat tricky. It is because this problem is
harder: the input string can be longer, so it would take too long to output all the edge labels directly
and compare them with the correct ones, as their combined length can be Θ(|Text|2 ), which is too
much when the Text can be as long as 200 000 characters.
Output the 𝑇 𝑒𝑥𝑡 from the input on the first line. Then output all the edges of the suffix tree
in a specific order (see below), each on its own line. Output each edge as a pair of integers
(start, end), where start is the position in Text corresponding to the start of the edge label substring
in the Text and end is the position right after the end of the edge label in the Text. Note that start
must be a valid position in the Text, that is, 0 ≤ 𝑠𝑡𝑎𝑟𝑡 ≤ |Text| − 1, and end must be between 1 and
|Text| inclusive. Substring Text[start..end − 1] must be equal to the edge label of the corresponding
edge. For example, if Text = “ACACAA$” and the edge label is “CA”, you can output this edge either
as (1, 3) corresponding to Text[1..2] = “CA” or as (3, 5) corresponding to Text[3..4] = “CA” — both
variants will be accepted.
The order of the edges is important here — if you output all the correct edges in the wrong order, your
solution will not be accepted. However, you don’t need to construct this order yourself if you
write in C++, Java or Python3, because it is implemented for you in the starter files.
Output all the edges in the order of sorted suffixes: first, take the leaf of the suffix tree corresponding
to the smallest suffix of Text and output all the edges on the path from the root to this leaf. Then
take the leaf corresponding to the second smallest suffix of Text and output all the edges on the path
from the root to this leaf except for those edges which were printed before. Then take the leaf
corresponding to the third smallest suffix, fourth smallest suffix and so on. Print each edge only once
— as a part of the path corresponding to the smallest suffix of Text where this edge appears. This way,
you will only output 𝑂(|Text|) integers. See the examples below for clarification.
10
Time Limits.
language C C++ Java Python C# Haskell JavaScript Ruby Scala
time in seconds 2 2 20 10 3 4 10 10 40
Sample 2.
Input:
AAA$
3210
012
Output:
AAA$
34
01
34
12
34
24
Explanation:
$ A
𝑖 𝑆𝐴[𝑖] 𝐿𝐶𝑃 [𝑖] suffix
0 3 0 $ 3
1 2 1 A$ $ A
2 1 2 AA$
3 0 AAA$ 2
$ A$
1 0
11
Sample 3.
Input:
GTAGT$
523041
00201
Output:
GTAGT$
56
26
35
56
26
45
56
26
Explanation:
Starter Files
The starter solutions for this problem read the input data from the standard input, pass it to the blank
procedure which is supposed to build the suffix tree, then print the edges of the resulting suffix tree in the
order specified in the output format description. You don’t have to construct the correct order of
edges for the output yourself if you use these starter files. You are only supposed to implement the
procedure to build the suffix tree given the string, its suffix array and LCP array if you are using C++, Java,
or Python3. For other programming languages, you need to implement a solution from scratch. Filename:
suffix_tree_from_array
What To Do
To solve this problem, it is enough to implement carefully the corresponding algorithm covered in the lectures.
Need Help?
Ask a question or see the questions asked by other learners at this forum thread.
12
5 General Instructions and Recommendations on Solving Algorith-
mic Problems
Your main goal in an algorithmic problem is to implement a program that solves a given computational
problem in just few seconds even on massive datasets. Your program should read a dataset from the
standard input and write an answer to the standard output.
Below we provide general instructions and recommendations on solving such problems. Before reading
them, go through readings and screencasts in the first module that show a step by step process of solving
two algorithmic problems: link.
13
same way on your machine and on the testing machine (note that a buggy program may behave differently
when compiled by different compilers, or even by the same compiler with different flags).
If your C/C++ compiler does not recognize -std=c++14 flag, try replacing it with -std=c++0x flag or
compiling without this flag at all (all starter solutions can be compiled without it). On Linux and
MacOS, you most probably have the required compiler. On Windows, you may use your favorite
compiler or install, e.g., cygwin.
∙ Python 2 (CPython 2.7). File extensions: .py2 or .py (a file ending in .py needs to have a first line
which is a comment containing “python2”). No flags:
python2
∙ Python 3 (CPython 3.4). File extensions: .py3 or .py (a file ending in .py needs to have a first line
which is a comment containing “python3”). No flags:
python3
14
5.5 Testing Your Program
When your program is ready, you start testing it. It makes sense to start with small datasets — for example,
sample tests provided in the problem description. Ensure that your program produces a correct result.
You then proceed to checking how long does it take your program to process a massive dataset. For
this, it makes sense to implement your algorithm as a function like solve(dataset) and then implement an
additional procedure generate() that produces a large dataset. For example, if an input to a problem is a
sequence of integers of length 1 ≤ 𝑛 ≤ 105 , then generate a sequence of length exactly 105 , pass it to your
solve() function, and ensure that the program outputs the result quickly.
Also, check the boundary values. Ensure that your program processes correctly sequences of size 𝑛 =
1, 2, 105 . If a sequence of integers from 0 to, say, 106 is given as an input, check how your program behaves
when it is given a sequence 0, 0, . . . , 0 or a sequence 106 , 106 , . . . , 106 . Check also on randomly generated
data. For each such test check that you program produces a correct result (or at least a reasonably looking
result).
In the end, we encourage you to stress test your program to make sure it passes in the system at the first
attempt. See the readings and screencasts from the first week to learn about testing and stress testing: link.
15
6 Frequently Asked Questions
6.1 I submit the program, but nothing happens. Why?
You need to create submission and upload the file with your solution in one of the programming languages C,
C++, Java, or Python (see Subsections 5.3 and 5.4). Make sure that after uploading the file with your solution
you press on the blue “Submit” button in the bottom. After that, the grading starts, and the submission
being graded is enclosed in an orange rectangle. After the testing is finished, the rectangle disappears, and
the results of the testing of all problems is shown to you.
6.2 I submit the solution only for one problem, but all the problems in the
assignment are graded. Why?
Each time you submit any solution, the last uploaded solution for each problem is tested. Don’t worry: this
doesn’t affect your score even if the submissions for the other problems are wrong. As soon as you pass the
sufficient number of problems in the assignment (see in the pdf with instructions), you pass the assignment.
After that, you can improve your result if you successfully pass more problems from the assignment. We
recommend working on one problem at a time, checking whether your solution for any given problem passes
in the system as soon as you are confident in it. However, it is better to test it first, please refer to the
reading about stress testing: link.
6.3 What are the possible grading outcomes, and how to read them?
Your solution may either pass or not. To pass, it must work without crashing and return the correct answers
on all the test cases we prepared for you, and do so under the time limit and memory limit constraints
specified in the problem statement. If your solution passes, you get the corresponding feedback "Good job!"
and get a point for the problem. If your solution fails, it can be because it crashes, returns wrong answer,
works for too long or uses too much memory for some test case. The feedback will contain the number of
the test case on which your solution fails and the total number of test cases in the system. The tests for the
problem are numbered from 1 to the total number of test cases for the problem, and the program is always
tested on all the tests in the order from the test number 1 to the test with the biggest number.
Here are the possible outcomes:
Good job! Hurrah! Your solution passed, and you get a point!
Wrong answer. Your solution has output incorrect answer for some test case. If it is a sample test case from
the problem statement, or if you are solving Programming Assignment 1, you will also see the input
data, the output of your program and the correct answer. Otherwise, you won’t know the input, the
output, and the correct answer. Check that you consider all the cases correctly, avoid integer overflow,
output the required white space, output the floating point numbers with the required precision, don’t
output anything in addition to what you are asked to output in the output specification of the problem
statement. See this reading on testing: link.
Time limit exceeded. Your solution worked longer than the allowed time limit for some test case. If it
is a sample test case from the problem statement, or if you are solving Programming Assignment 1,
you will also see the input data and the correct answer. Otherwise, you won’t know the input and the
correct answer. Check again that your algorithm has good enough running time estimate. Test your
program locally on the test of maximum size allowed by the problem statement and see how long it
works. Check that your program doesn’t wait for some input from the user which makes it to wait
forever. See this reading on testing: link.
Memory limit exceeded. Your solution used more than the allowed memory limit for some test case. If it
is a sample test case from the problem statement, or if you are solving Programming Assignment 1,
16
you will also see the input data and the correct answer. Otherwise, you won’t know the input and the
correct answer. Estimate the amount of memory that your program is going to use in the worst case
and check that it is less than the memory limit. Check that you don’t create too large arrays or data
structures. Check that you don’t create large arrays or lists or vectors consisting of empty arrays or
empty strings, since those in some cases still eat up memory. Test your program locally on the test of
maximum size allowed by the problem statement and look at its memory consumption in the system.
Cannot check answer. Perhaps output format is wrong. This happens when you output something
completely different than expected. For example, you are required to output word “Yes” or “No”, but
you output number 1 or 0, or vice versa. Or your program has empty output. Or your program outputs
not only the correct answer, but also some additional information (this is not allowed, so please follow
exactly the output format specified in the problem statement). Maybe your program doesn’t output
anything, because it crashes.
Unknown signal 6 (or 7, or 8, or 11, or some other). This happens when your program crashes.
It can be because of division by zero, accessing memory outside of the array bounds, using uninitial-
ized variables, too deep recursion that triggers stack overflow, sorting with contradictory comparator,
removing elements from an empty data structure, trying to allocate too much memory, and many other
reasons. Look at your code and think about all those possibilities. Make sure that you use the same
compilers and the same compiler options as we do. Try different testing techniques from this reading:
link.
Internal error: exception... Most probably, you submitted a compiled program instead of a source
code.
Grading failed. Something very wrong happened with the system. Contact Coursera for help or write in
the forums to let us know.
17
make one’s programs work, one must test them really well. Sometimes, the programs still don’t work although
you tried really hard to test them, and you need to be both skilled and creative to fix your bugs. Solutions
to algorithmic problems are one of the hardest to implement correctly. That’s why in this Specialization you
will gain this important experience which will be invaluable in the future when you write programs which
you really need to get right.
It is crucial for you to learn to test and fix your programs yourself. In the real life, often there will be
no or only partial information about the failure of your program or service. Still, you will have to reproduce
the failure to fix it (or just guess what it is, but that’s rare, and you will still need to reproduce the failure
to make sure you have really fixed it). When you solve algorithmic problems, it is very frequent to make
subtle mistakes. That’s why you should apply the testing techniques described in this reading to find the
failing test case and fix your program.
6.6 My solution does not pass the tests? May I post it in the forum and ask
for a help?
No, please do not post any solutions in the forum or anywhere on the web, even if a solution does not
pass the tests (as in this case you are still revealing parts of a correct solution). Recall the third item
of the Coursera Honor Code: “I will not make solutions to homework, quizzes, exams, projects, and other
assignments available to anyone else (except to the extent an assignment explicitly permits sharing solutions).
This includes both solutions written by me, as well as any solutions provided by the course staff or others”
(link).
6.7 My implementation always fails in the grader, though I already tested and
stress tested it a lot. Would not it be better if you give me a solution to
this problem or at least the test cases that you use? I will then be able to
fix my code and will learn how to avoid making mistakes. Otherwise, I do
not feel that I learn anything from solving this problem. I am just stuck.
First of all, you always learn from your mistakes.
The process of trying to invent new test cases that might fail your program and proving them wrong is
often enlightening. This thinking about the invariants which you expect your loops, ifs, etc. to keep and
proving them wrong (or right) makes you understand what happens inside your program and in the general
algorithm you’re studying much more.
Also, it is important to be able to find a bug in your implementation without knowing a test case and
without having a reference solution. Assume that you designed an application and an annoyed user reports
that it crashed. Most probably, the user will not tell you the exact sequence of operations that led to a
crash. Moreover, there will be no reference application. Hence, once again, it is important to be able to
locate a bug in your implementation yourself, without a magic oracle giving you either a test case that your
program fails or a reference solution. We encourage you to use programming assignments in this class as a
way of practicing this important skill.
If you have already tested a lot (considered all corner cases that you can imagine, constructed a set of
manual test cases, applied stress testing), but your program still fails and you are stuck, try to ask for help
on the forum. We encourage you to do this by first explaining what kind of corner cases you have already
considered (it may happen that when writing such a post you will realize that you missed some corner cases!)
and only then asking other learners to give you more ideas for tests cases.
18