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

Chai Wah Wu

2024-03-02

Abstract:In modern usage the Bernoulli numbers and polynomials follow Euler's approach and are defined using generating functions. We consider the functional equation $f(x)+x^k=f(x+1)$ and show that a solution can be derived from Faulhaber's formula for the sum of powers. We show how these solutions provide a characterization of Bernoulli numbers and related results.

Number Theory,Combinatorics,History and Overview

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is to solve the functional equation \( f(x)+x^k = f(x + 1)\) through the Faulhaber formula and Bernoulli numbers, and use these results to characterize Bernoulli numbers. Specifically, the author explores how to use the Faulhaber formula (a formula for calculating the sum of powers) to derive solutions that satisfy the above - mentioned functional equation, and shows the relationships between these solutions and Bernoulli numbers, Bernoulli polynomials and their related identities.

### Specific description of the problem

The functional equation considered in the paper is:
\[ f(x)+x^k = f(x + 1)\tag{1} \]

### Main methods and conclusions

1. **Faulhaber formula**: The Faulhaber formula is a polynomial expression used to represent the sum of consecutive integer powers. For non - negative integers \( n\) and \( p\), the Faulhaber formula can be written as:
\[
F(n,p)=\sum_{i = 1}^{n}i^p=\frac{1}{p + 1}\sum_{j=0}^{p}\binom{p+1}{j}B_j^+n^{p+1-j}
\]
where \( B_j^+\) is the second - kind Bernoulli number.

2. **Solving the functional equation**: The author defines a polynomial \( f_k(x)\) to solve the above - mentioned functional equation:
\[
f_k(x)=\frac{1}{k+1}\sum_{j = 0}^{k}\binom{k+1}{j}B_j^-x^{k+1-j}
\]
In this way, the author proves that \( f_k(x)\) satisfies \( f_k(x + 1)=f_k(x)+x^k\).

3. **Characterization of Bernoulli numbers**: By solving the above - mentioned functional equation, the author shows how to extract Bernoulli numbers from the solutions. Specifically, the coefficients \( a_j\) of the polynomial \( p(x)\) with the smallest degree can be used to define the first \( k + 2\) Bernoulli numbers \( B_j^-\):
\[
B_j^-=(k + 1)a_{k+1-j}/\binom{k+1}{j}
\]

4. **Generalization to the first - kind Bernoulli numbers**: Similarly, the author also shows how to use the first - kind Bernoulli numbers \( B_j^+\) to solve another related functional equation:
\[
f(x)+(x + 1)^k=f(x + 1)
\]

5. **Appell sequence properties**: Since Bernoulli polynomials form an Appell sequence, \( f_k(x)\) also has similar properties and satisfies:
\[
f_k(x + y)=\sum_{i=0}^{k}\binom{k}{i}f_i(x)y^{k - i}
\]

### Summary

The main contribution of this paper is to provide a systematic method for solving a specific form of functional equations through the Faulhaber formula and Bernoulli numbers, and to reveal the profound connections between these solutions and Bernoulli numbers. This not only deepens our understanding of Bernoulli numbers, but also provides new tools and perspectives for studying the sum of powers and related mathematical problems.