An algorithm for the Faulhaber polynomials

Let $S_p(n)$ denote the sum of $p$th powers of the first $n$ positive integers $1^p+2^p+\cdots +n^p$. In this paper, first we express $S_p(n)$ in the so-called Faulhaber form, namely, as an even or odd polynomial in $(n + 1/2)$,accordingas$p$isoddoreven\mathrm{.}Then,usingtherelation$S_p(n) - S_p(n-1) = n^p$, we derive a recursive formula for the associated Faulhaber coefficients. Applying Cramer’s rule to the corresponding system of equations, we obtain an explicit determinant formula for the said coefficients. Furthermore, we show how to convert the (even or odd) Faulhaber polynomials in $(n+1\mathrm{/}2)$ into polynomials in $S_1(n)$ for any arbitrary $p$, and vice versa.