Fast simulation of reconstructed phylogenetic trees under