In this paper, we propose a novel exact algorithm for structure discovery in Bayesian networks
of a moderate size (say, 25 variables or less). In fact, we consider two versions of the algorithm: one
for computing posterior probabilities of structural features; another for finding an optimal structure.
We also present a rigorous complexity analysis of the algorithm. This work is motivated by three
serial observations summarized below.