Refactor the multiprecision code of gamma_imp_bernoulli().