View source for BK-tree

{{technical|date=March 2019}}
A '''BK-tree''' is a [[metric trees|metric tree]] suggested by Walter Austin Burkhard and [[Robert M. Keller]]{{ref|BK73}} specifically adapted to discrete [[metric space]]s.
For simplicity, consider '''integer''' discrete metric <math>d(x,y)</math>. Then, BK-tree is defined in the following way. An arbitrary element ''a'' is selected as root node.  The root node may have zero or more subtrees.  The ''k-th'' subtree is recursively built of all elements ''b'' such that <math>d(a,b) = k</math>. BK-trees can be used for [[approximate string matching]] in a dictionary.{{ref|BN98}}{{example needed|date=March 2019}}

== Example==

[[File:Bk tree.svg|thumb|right|An example of BK-tree]]

This picture depicts the BK-tree for the set <math>W</math> of words {"book", "books", "cake", "boo", "boon", "cook", "cake", "cape", "cart"} obtained by using the [[Levenshtein distance]]
* each node <math>u</math> is labeled by a string of <math>w_u \in W</math>;
* each [[Directed graph|arc]] <math>(u,v)</math> is labeled by <math>d_{uv} = d(w_u,w_v)</math> where <math>w_u</math> denotes the word assigned to <math>u</math>.

The BK-tree is built so that:
* for all node <math>u</math> of the BK-tree, the weight assigned to its egress arcs are distinct;
* for all arc <math>e=(u,v)</math> labeled by <math>k</math>, each descendant <math>v'</math> of <math>v</math> satisfies the following equation: <math>d(w_u, w_{v'}) = k</math>:
** ''Example 1:'' Consider the arc from "book" to "books". The distance between "book" and any word in {"books", "boo", "boon", "cook"} is equal to 1;
** ''Example 2:'' Consider the arc from "books" to "boo".  The distance between "books" and any word in {"boo", "boon", "cook"} is equal to 2.

== Insertion ==

The insertion primitive is used to populate a BK-tree <math>t</math> according to a discrete metric <math>d</math>.

'''Input:'''
* <math>t</math>: the BK-tree;
** <math>d_{uv}</math> denotes the weight assigned to an arc <math>(u, v)</math>;
** <math>w_u</math> denotes word assigned to a node <math>u</math>;
* <math>d</math>: the discrete metric used by <math>t</math>  (e.g. the [[Levenshtein distance]]);
* <math>w</math>: the element to be inserted into <math>t</math>;

'''Output:'''
* The node of <math>t</math> corresponding to <math>w</math>

'''Algorithm:'''
* '''If''' the <math>t</math> is empty:
** Create a root node <math>r</math> in <math>t</math>
** <math>w_r \leftarrow w</math>
** '''Return''' <math>r</math>
* Set <math>u</math> to the root of <math>t</math>
* '''While''' <math>u</math> exists:
** <math>k \leftarrow d(w_u, w)</math>
** '''If''' <math>k = 0</math>:
*** '''Return''' <math>u</math>
** Find <math>v</math> the child of <math>u</math> such that <math>d_{uv} = k</math>
** '''If '''<math>v</math> is not found:
*** Create the node <math>v</math>
*** <math>w_v \leftarrow w</math>
*** Create the arc <math>(u, v)</math>
*** <math>d_{uv} \leftarrow k</math>
*** '''Return''' <math>v</math>
** <math>u \leftarrow v</math>

== Lookup ==

Given a searched element <math>w</math>, the lookup primitive traverses the BK-tree to find the closest element of <math>w</math>. The key idea is to restrict the exploration of <math>t</math> to nodes that can only improve the best candidate found so far by taking advantage of the BK-tree organization and of the triangle inequality (cut-off criterion).

'''Input:'''
* <math>t</math>: the BK-tree;
* <math>d</math>: the corresponding discrete metric (e.g. the [[Levenshtein distance]]);
* <math>w</math>: the searched element;
* <math>d_{max}</math>: the maximum distance allowed between the best match and <math>w</math>, defaults to <math>+\infty</math>;

'''Output:'''
* <math>w_{best}</math>: the closest element to <math>w</math> stored in <math>t</math> and according to <math>d</math> or <math>\perp</math> if not found;

'''Algorithm:'''
* '''If''' <math>t</math> is empty:
** '''Return''' <math>\perp</math>
* Create <math>S</math> a set of nodes to process, and insert the root of <math>t</math> into <math>S</math>.
* <math>(w_{best}, d_{best}) \leftarrow (\perp, d_{max})</math>
* '''While''' <math>S \ne \emptyset</math>:
** Pop an arbitrary node <math>u</math> from <math>S</math>
** <math>d_u \leftarrow d(w, w_u)</math>
** '''If''' <math>d_u < d_{best}</math>:
*** <math>(w_{best}, d_{best}) \leftarrow (w_u, d_u)</math>
** '''For each''' egress-arc <math>(u, v)</math>:
*** '''If''' <math>|d_{uv} - d_u| < d_{best}</math>: (cut-off criterion)
**** Insert <math>v</math> into <math>S</math>.
* '''Return''' <math>w_{best}</math>

== Example of the lookup algorithm ==

Consider the example 8-node B-K Tree shown above and set <math>w=</math>"cool". <math>S</math> is initialized to contain the root of the tree, which is subsequently popped as the first value of <math>u</math> with <math>w_u</math>="book". Further <math>d_u=2</math> since the distance from "book" to "cool" is 2, and <math>d_{best}=2</math> as this is the best (i.e. smallest) distance found thus far. Next each outgoing arc from the root is considered in turn: the arc from "book" to "books" has weight 1, and since <math>|1-2|=1</math> is less than <math>d_{best}=2</math>, the node containing "books" is inserted into <math>S</math> for further processing. 
The next arc, from "book" to "cake," has weight 4, and since <math>|4-2|=2</math> is not less than <math>d_{best}=2</math>, the node containing "cake" is not inserted into <math>S</math>. Therefore, the subtree rooted at "cake" will be pruned from the search, as the word closest to "cool" cannot appear in that subtree.  To see why this pruning is correct, notice that a candidate word <math>c</math> appearing in "cake"s subtree having distance less than 2 to "cool" would violate the triangle inequality: the triangle inequality requires that for this set of three numbers (as sides of a triangle), no two can sum to less than the third, but here the distance from "cool" to "book" (which is 2) plus the distance from "cool" to <math>c</math> (which is less than 2) cannot reach or exceed the distance from "book" to "cake" (which is 4).  Therefore, it is safe to disregard the entire subtree rooted at "cake".

Next the node containing "books" is popped from <math>S</math> and now <math>d_u=3</math>, the distance from "cool" to "books."  As <math>d_u > d_{best}</math>, <math>d_{best}</math> remains set at 2 and the single outgoing arc from the node containing "books" is considered.  Next, the node containing "boo" is popped from <math>S</math> and <math>d_u=2</math>, the distance from "cool" to "boo."  This again does not improve upon <math>d_{best} = 2</math>.  Each outgoing arc from "boo" is now considered; the arc from "boo" to "boon" has weight 1, and since <math>|2-1|=1 < d_{best}=2</math>, "boon" is added to <math>S</math>.  Similarly, since <math>|2-2|=0 < d_{best}</math>, "cook" is also added to <math>S</math>.

Finally each of the two last elements in <math>S</math> are considered in arbitrary order: suppose the node containing "cook" is popped first, improving <math>d_{best}</math> to distance 1, then the node containing "boon" is popped last, which has distance 2 from "cool" and therefore does not improve the best result.  Finally, "cook" is returned as the answer <math>w_{best}</math> with <math>d_{best}=1</math>.

== See also ==

* [[Levenshtein distance]] – the distance metric commonly used when building a BK-tree
* [[Damerau–Levenshtein distance]] – a modified form of Levenshtein distance that allows transpositions

== References ==
* {{note|BK73}} [http://doi.acm.org/10.1145/362003.362025 W. Burkhard and R. Keller. Some approaches to best-match file searching, CACM, 1973]
* {{note|BCMW04}} R. Baeza-Yates, W. Cunto, U. Manber, and S. Wu. Proximity matching using fixed queries trees. In M. Crochemore and D. Gusfield, editors, 5th Combinatorial   Pattern Matching, LNCS 807, pages 198–212, Asilomar, CA, June 1994.
* {{note|BN98}}  [https://web.archive.org/web/20110825210139/http://reference.kfupm.edu.sa/content/f/a/fast_approximate_string_matching_in_a_di_269284.pdf Ricardo Baeza-Yates and Gonzalo Navarro. Fast Approximate String Matching in a Dictionary. Proc. SPIRE'98]

== External links ==
* A BK-tree implementation in [https://github.com/vy/bk-tree Common Lisp] with test results and performance graphs.
* An explanation of BK-Trees and their relationship to metric spaces [http://blog.notdot.net/2007/4/Damn-Cool-Algorithms-Part-1-BK-Trees]
* An explanation of BK-Trees with an implementation in C# [http://nullwords.wordpress.com/2013/03/13/the-bk-tree-a-data-structure-for-spell-checking/]
* A BK-tree implementation in [[Lua (programming language)|Lua]] [https://profan.github.io/lua-bk-tree/]
* A BK-tree implementation in [[Python (programming language)|Python]] [https://github.com/benhoyt/pybktree]

{{CS-Trees}}

[[Category:Trees (data structures)]]