small speedup, use boost::math

use boost::math for fact(), binomial_coeficient()