Summary: The problem of finding formulas for sums of powers of natural numbers has been of interest to mathematicians for many centuries. Among these is Faulhabers well-known formula expressing the power sums as polynomials whose coefficients involve Bernoulli numbers. In this paper we give an elementary proof that the sum of $p$-th powers of the first $n$ natural numbers can be expressed as a polynomial in $n$ of degree $p + 1$. We also prove a novel identity involving Bernoulli numbers and use it to show the symmetry of this polynomial.

11B68, 11B37

number theory, power sum, Bernoulli number