separate computation of Bernoulli coefficients and Faulhaber polynomials