In this note, we first derive a recursive formula for the sum of powers S k (n) = 1 k + 2 k + · · · + n k , with k and n non-negative integers. We then apply it to establish, via Cramer's rule, an explicit determinant formula for S k (n) involving the Bernoulli numbers and the binomial (n + 1 2 ). Evaluating the determinant gives us directly S k (n) in the form of the so-called Faulhaber polynomial, namely as a sum of even or odd powers of (n + 1 2 ). Furthermore, a connection with Hessenberg matrices is shown.