Aho-Corasick algorithm
Let there be a set of strings with the total length (m) (sum of all lengths).
The Aho-Corasick algorithm constructs a data structure similar to a trie with some additional links, and then constructs a finite state machine (automaton) in (O(m k)) time, where (k) is the size of the used alphabet.
The algorithm was proposed by Alfred Aho and Margaret Corasick in 1975.
Construction of the trie
Formally a trie is a rooted tree, where each edge of the tree is labeled by some letter.
All outgoing edge from one vertex mush have different labels.
Consider any path in the trie from the root to any vertex.
If we write out the labels of all edges on the path, we get a string that corresponds to this path.
For any vertex in the trie we will associate the string from the root to the vertex.
Each vertex will also have a flag ( ext{leaf}) which will be true, if any string from the given set corresponds to this vertex.
Accordingly to build a trie for a set of strings means to build a trie such that each ( ext{leaf}) vertex will correspond to one string from the set, and conversely that each string of the set corresponds to one ( ext{leaf}) vertex.
We now describe how to construct a trie for a given set of strings in linear time with respect to their total length.
We introduce a structure for the vertices of the tree.
const int K = 26;
struct Vertex {
int next[K];
bool leaf = false;
Vertex() {
fill(begin(next), end(next), -1);
}
};
vector<Vertex> trie(1);
Here we store the trie as an array of ( ext{Vertex}).
Each ( ext{Vertex}) contains the flag ( ext{leaf}), and the edges in the form of ans array ( ext{next}[]), where ( ext{next}[i]) is the index to the vertex that we reach by following the character (i), or (-1), if there is no such edge.
Initially the trie consists of only one vertex - the root - with the index (0).
Now we implement a function that will add a string (s) to the trie.
The implementation is extremely simple:
we start at the root node, and as long as there are edges corresponding to the characters of (s) we follow them.
If there is no edge for one character, we simply generate a new vertex and connect it via an edge.
At the end of the process we mark the last vertex with flag ( ext{leaf}).
void add_string(string const& s) {
int v = 0;
for (char ch : s) {
int c = ch - 'a';
if (trie[v].next[c] == -1) {
trie[v].next[c] = trie.size();
trie.emplace_back();
}
v = trie[v].next[c];
}
trie[v].leaf = true;
}
The implementation obviously runs in linear time.
And since every vertex store (k) links, it will use (O(m k)) memory.
It is possible to decrease the memory consumption to (O(m)) by using a map instead of an array in each vertex.
However this will increase the complexity to (O(n log k)).
Construction of an automaton
Suppose we have built a trie for the given set of strings.
Now let's look at it from a different side.
If we look at any vertex.
The string that corresponds to it is a prefix of one or more strings in the set, thus each vertex of the trie can be interpreted as a position in one or more strings from the set.
In fact the trie vertices can be interpreted as states in a finite deterministic automaton.
From any state we can transition - using some input letter - to other states, i.e. to another position in the set of strings.
For example, if there is only one string in the trie (abc), and we are standing at vertex (2) (which corresponds to the string (ab)), then using the letter (c) we can transition to the state (3).
Thus we can understand the edges of the trie as transitions in an automaton according to the corresponding letter.
However for an automaton we cannot restrict the possible transitions for each state.
If we try to perform a transition using a letter, and there is no corresponding edge in the trie, then we nevertheless must go into some state.
More strictly, let us be in a state (p) corresponding to the string (t), and we want to transition to a different state with the character (c).
If the is an edge labeled with this letter (c), then we can simply go over this edge, and get the vertex corresponding to (t + c).
If there is no such edge, then we must find the state corresponding to the longest proper suffix of the string (t) (the longest available in the trie), and try to perform a transition via (c) from there.
For example let the trie be constructed by the strings (ab) and (bc), and we are currently at the vertex corresponding to (ab), which is a ( ext{leaf}).
For a transition with the letter (c), we are forced to go to the state corresponding to the string (b), and from there follow the edge with the letter (c).
A suffix link for a vertex (p) is a edge that points to the longest proper suffix of the string corresponding to the vertex (p).
The only special case is the root of the trie, the suffix link will point to itself.
Now we can reformulate the statement about the transitions in the automaton like this:
while from the current vertex of the trie there is no transition using the current letter (or until we reach the root), we follow the suffix link.
Thus we reduced the problem of constructing an automaton to the problem of finding suffix links for all vertices of the trie.
However we will build these suffix links, oddly enough, using the transitions constructed in the automaton.
Note that if we want to find a suffix link for some vertex (v), then we can go to the ancestor (p) of the current vertex (let (c) be the letter of the edge from (p) to (v)), then follow its suffix link, and perform from there the transition with the letter (c).
Thus the problem of finding the transitions has been reduced to the problem of finding suffix links, and the problem of finding suffix links has been reduced to the problem of finding a suffix link and a transition, but for vertices closer to the root.
So we have a recursive dependence that we can resolve in linear time.
Let's move to the implementation.
Note that we now will store the ancestor (p) and the character (pch) of the edge from (p) to (v) for each vertex (v).
Also at each vertex we will store the suffix link ( ext{link}) (or (-1) if it hasn't been calculated yet), and in the array ( ext{go}[k]) the transitions in the machine for each symbol (again (-1) if it hasn't been calculated yet).
const int K = 26;
struct Vertex {
int next[K];
bool leaf = false;
int p = -1;
char pch;
int link = -1;
int go[K];
Vertex(int p=-1, char ch='$') : p(p), pch(ch) {
fill(begin(next), end(next), -1);
fill(begin(go), end(go), -1);
}
};
vector<Vertex> t(1);
void add_string(string const& s) {
int v = 0;
for (char ch : s) {
int c = ch - 'a';
if (t[v].next[c] == -1) {
t[v].next[c] = t.size();
t.emplace_back(v, ch);
}
v = t[v].next[c];
}
t[v].leaf = true;
}
int go(int v, char ch);
int get_link(int v) {
if (t[v].link == -1) {
if (v == 0 || t[v].p == 0)
t[v].link = 0;
else
t[v].link = go(get_link(t[v].p), t[v].pch);
}
return t[v].link;
}
int go(int v, char ch) {
int c = ch - 'a';
if (t[v].go[c] == -1) {
if (t[v].next[c] != -1)
t[v].go[c] = t[v].next[c];
else
t[v].go[c] = v == 0 ? 0 : go(get_link(v), ch);
}
return t[v].go[c];
}
It is easy to see, that due to the memorization of the found suffix links and transitions the total time for finding all suffix links and transitions will be linear.
Applications
Find all strings from a given set in a text
Given a set of strings and a text.
We have to print all occurrences of all strings from the set in the given text in (O( ext{len} + ext{ans})), where ( ext{len}) is the length of the text and ( ext{ans}) is the size of the answer.
We construct an automaton for this set of strings.
We will now process the text letter by letter, transitioning during the different states.
Initially we are at the root of the trie.
If we are at any time at state (v), and the next letter is (c), then we transition to the next state with ( ext{go}(v, c)), thereby either increasing the length of the current match substring by (1), or decreasing it by following a suffix link.
How can we find out for a state (v), if there are any matches with strings for the set?
First, it is clear that if we stand on a ( ext{leaf}) vertex, then the string corresponding to the vertex ends at this position in the text.
However this is by no means the only possible case of achieving a match:
if we can reach one or more ( ext{leaf}) vertices by moving along the suffix links, then there will be also a match corresponding to each found ( ext{leaf}) vertex.
A simple example demonstrating this situation can be creating using the set of strings (\{dabce, abc, bc\}) and the text (dabc).
Thus if we store in each ( ext{leaf}) vertex the index of the string corresponding to it (or the list of indices if duplicate strings appear in the set), then we can find in (O(n)) time the indices of all strings which match the current state, by simply following the suffix links from the current vertex to the root.
However this is not the most efficient solution, since this gives us (O(n ~ ext{len})) complexity in total.
However this can be optimized by computing and storing the nearest ( ext{leaf}) vertex that is reachable using suffix links (this is sometimes called the exit link).
This value we can compute lazily in linear time.
Thus for each vertex we can advance in (O(1)) time to the next marked vertex in the suffix link path, i.e. to the next match.
Thus for each match we spend (O(1)) time, and therefore we reach the complexity (O( ext{len} + ext{ans})).
If you only want to count the occurrences and not find the indices themselves, you can calculate the number of marked vertices in the suffix link path for each vertex (v).
This can be calculated in (O(n)) time in total.
Thus we can sum up all matches in (O( ext{len})).
Finding the lexicographical smallest string of a given length that doesn't match any given strings
A set of strings and a length (L) is given.
We have to find a string of length (L), which does not contain any of the string, and derive the lexicographical smallest of such strings.
We can construct the automaton for the set of strings.
Let's remember, that the vertices from which we can reach a ( ext{leaf}) vertex are the states, at which we have a match with a string from the set.
Since in this task we have to avoid matches, we are not allowed to enter such states.
On the other hand we can enter all other vertices.
Thus we delete all "bad" vertices from the machine, and in the remaining graph of the automaton we find the lexicographical smallest path of length (L).
This task can be solved in (O(L)) for example by depth first search.
Finding the shortest string containing all given strings
Here we use the same ideas.
For each vertex we store a mask that denotes the strings which match at this state.
Then the problem can be reformulated as follows:
initially being in the state ((v = ext{root},~ ext{mask} = 0)), we want to reach the state ((v,~ ext{mask} = 2^n - 1)), where (n) is the number of strings in the set.
When we transition from one state to another using a letter, we update the mask accordingly.
By running a breath first search we can find a path to the state ((v,~ ext{mask} = 2^n - 1)) with the smallest length.
Finding the lexicographical smallest string of length (L) containing (k) strings
As in the previous problem, we calculate for each vertex the number of matches that correspond to it (that is the number of marked vertices reachable using suffix links).
We reformulate the problem: the current state is determined by a triple of numbers ((v,~ ext{len},~ ext{cnt})), and we want to reach from the state (( ext{root},~ 0,~ 0)) the state ((v,~ L,~ k)), where (v) can be any vertex.
Thus we can find such a path using depth first search (and if the search looks at the edges in their natural order, then the found path will automatically be the lexicographical smallest).