Abstract:Denote by $\Sigma n^m$ the sum of the $m$-th powers of the first $n$ positive integers $1^m+2^m+\ldots +n^m$. Similarly let $\Sigma^r n^m$ be the $r$-fold sum of the $m$-th powers of the first $n$ positive integers, defined such that $\Sigma^0 n^{m}=n^m$, and then recursively by $\Sigma^{r+1} n^{m}=\Sigma^{r} 1^{m}+\Sigma^{r} 2^{m}+\ldots + \Sigma^{r} n^{m}$. During the early 17th-century, polynomial expressions for the sums $\Sigma^r n^m$ and their factorisation and polynomial basis representation properties were examined by Johann Faulhaber, who published some remarkable theorems relating to these $r$-fold sums in his Academia Algebrae (1631). In this paper we consider families of polynomials related to the Fibonacci and Lucas polynomials which naturally lend themselves to representing sums and differences of integer powers, as well as $\Sigma^r n^m$. Using summations over polynomial basis representations for monomials we generalise this sum for any polynomial $q(x)\in\mathbb{Q}[x]$, encountering some interesting coefficient families in the process. Other results include using our polynomial expressions to state some theorems of Fermat, where in particular, we obtain an explicit expression for the quotient in Fermat's Little Theorem. In the final two sections we examine these sums in the context of the three pairs of Stirling number varieties and binomial decompositions. We also derive an expression for $\Sigma n^m$ in terms of the Stirling numbers and Pochhammer symbols, leading to a seemingly new relation between the Stirling numbers and the Bernoulli numbers.

Faulhaber's formula, Bernoulli numbers and the equation $f(x)+x^k=f(x+1)$

A summation formula for power series using eulerian fractions*

A summation formula for power series using Eulerian fractions

The Generalization of Faulhaber's Formula to Sums of Arbitrary Complex Powers

Formulas for Bernoulli Numbers and Polynomials

On Fibonacci Polynomial Expressions for Sums of $m$th Powers, their implications for Faulhaber's Formula and some Theorems of Fermat

On Convolution Formulae Involving Bernoulli and Euler Polynomials

Summation Formulas for the Products of the Frobenius–Euler Polynomials

Some Sum Relations Involving Bernoulli and Euler Polynomials

Some New Formulas for the Products of the Frobenius-Euler Polynomials

On an exponential power sum

Some Explicit Formulas for the Bernoulli and Euler Numbers and Polynomials

Explicit Expressions for Iterates of Power Series

Sums of powers of integers via differentiation

Explicit Formulas Associated with Some Families of Generalized Bernoulli and Euler Polynomials

A Note on the Bernoulli and Euler Polynomials

Solution of the Basel problem using the Feynman integral trick

Explicit, recurrent, determinantal expressions of the $k$th power of formal power series and applications to the generalized Bernoulli numbers

Obtaining More Karatsuba-like Formulae over the Binary Field

A Convolution Formula for Bernoulli Polynomials.

Bernoulli-Fibonacci Polynomials