En Matemáticas, la fórmula de Faulhaber, en honor de Johann Faulhaber, expresa la suma de las potencias de los primeros n números naturales
![{\displaystyle \sum _{k=1}^{n}k^{p}=1^{p}+2^{p}+3^{p}+\cdots +n^{p}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cdfa28bc350e73f808fc51da16d427df1a45fd28)
como un polinomio en n de grado
cuyos coeficientes se construyen a partir de los números de Bernoulli:
.
La fórmula es la siguiente:
![{\displaystyle \sum _{k=1}^{n}k^{p}={1 \over p+1}\sum _{j=0}^{p}{p+1 \choose j}B_{j}n^{p+1-j}\qquad \left({\mbox{con }}B_{1}=+{1 \over 2}{\mbox{ en vez de }}-{1 \over 2}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8dac92eca0ca9fbd80c2524d3810ac877c384cb3)
Faulhaber no conoció nunca esta fórmula general; lo que sí conoció fueron al menos los primeros 17 casos y el hecho de que, si el exponente es impar, entonces la suma es una función polinomial de la suma en el caso especial en el que el exponente sea 1. También hizo algunas generalizaciones (véase Knuth).
La demostración de la fórmula de Faulhaber se puede encontrar en The Book of Numbers de John Horton Conway y Richard Guy.
Los primeros casos[editar]
![{\displaystyle 1+2+3+\cdots +n={n(n+1) \over 2}={n^{2}+n \over 2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/70f43d8f46e7fc6f4e151b5ffc23055cdc8b694c)
![{\displaystyle 1^{2}+2^{2}+3^{2}+\cdots +n^{2}={n(n+1)(2n+1) \over 6}={2n^{3}+3n^{2}+n \over 6}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/62cf67a599782333c4a1b223136a367bcb61e5e7)
![{\displaystyle 1^{3}+2^{3}+3^{3}+\cdots +n^{3}=\left({n^{2}+n \over 2}\right)^{2}={n^{4}+2n^{3}+n^{2} \over 4}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/20c7cac2bf61a8002eb569870dd9488bbc5ae7bf)
![{\displaystyle 1^{4}+2^{4}+3^{4}+\cdots +n^{4}={6n^{5}+15n^{4}+10n^{3}-n \over 30}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0cd150a4e85c59a5ac13ef2546c8de444b0231eb)
![{\displaystyle 1^{5}+2^{5}+3^{5}+\cdots +n^{5}={2n^{6}+6n^{5}+5n^{4}-n^{2} \over 12}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/635d7e9259c3534ad7e1ff923a0910392341e5c5)
![{\displaystyle 1^{6}+2^{6}+3^{6}+\cdots +n^{6}={6n^{7}+21n^{6}+21n^{5}-7n^{3}+n \over 42}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/63045680aa78ab75d7ffe7474eaab26fed381ada)
Forma alternativa[editar]
Si el índice de suma de la serie va desde 1 hasta
en vez desde 1 hasta n, estas fórmulas se modifican de tal manera que el único cambio es que tomamos
en vez de +1/2 (es decir, en este caso en la fórmula sólo intervienen números de Bernoulli); así, el segundo término de mayor orden en todos los resultados anteriores cambia el signo de suma por el de diferencia.
Relación con los polinomios de Bernoulli[editar]
La fórmula de Faulhaber se puede escribir en función de los polinomios de Bernoulli así:
![{\displaystyle \sum _{k=0}^{n}k^{p}={\frac {B_{p+1}(n+1)-B_{p+1}(0)}{p+1}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a53f8c0bec4c7c36e6b1928d9e41a6e102b4e3ec)
Forma Umbral[editar]
En el cálculo umbral clásico, se trata formalmente a los índices j en una secuencia
como si estos fueran exponentes. Haciendo esto, y siempre considerando la variante
[1], podemos aplicar el teorema del binomio y obtener:
![{\displaystyle \sum _{k=1}^{n}k^{p}={1 \over p+1}\sum _{j=0}^{p}{p+1 \choose j}B_{j}n^{p+1-j}={1 \over p+1}\sum _{j=0}^{p}{p+1 \choose j}B^{j}n^{p+1-j}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a6d340b3d93eb7c703f937f31c2c370db06ec2b)
![{\displaystyle ={(B+n)^{p+1}-B^{p+1} \over p+1}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eb309ef91420ee468b9f2f6ced895e5fb7a10d2e)
En el cálculo umbral moderno, se construye el funcional lineal T en el espacio vectorial de polinomios en una variable b dada por:
![{\displaystyle T(b^{j})=B_{j}.\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d11005cc296d998ce8f1b1630325e7462fb3178d)
Entonces se obtiene
![{\displaystyle \sum _{k=1}^{n}k^{p}={1 \over p+1}\sum _{j=0}^{p}{p+1 \choose j}B_{j}n^{p+1-j}={1 \over p+1}\sum _{j=0}^{p}{p+1 \choose j}T(b^{j})n^{p+1-j}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eeaa8fdb9183572838a77c92bfcbb81bd78d0676)
![{\displaystyle ={1 \over p+1}T\left(\sum _{j=0}^{p}{p+1 \choose j}b^{j}n^{p+1-j}\right)=T\left({(b+n)^{p+1}-b^{p+1} \over p+1}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a4fe4fafd6e1cc82f1812369dde0201fd963049e)
Polinomios de Faulhaber[editar]
Faulhaber observó que, si p es impar, entonces
![{\displaystyle 1^{p}+2^{p}+3^{p}+\cdots +n^{p}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/58eb615a9a5bfadc20dac44757e89f3d5e0eeb6f)
es un polinomio en a, donde a es la suma de los n primeros naturales:
![{\displaystyle a=1+2+3+\cdots +n.\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9506d1403f677b48a87df8d72d536a1e1d456f20)
En particular se tiene:
![{\displaystyle 1^{3}+2^{3}+3^{3}+\cdots +n^{3}=a^{2}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f2980460457210200b1729adc99b63c4229cda95)
![{\displaystyle 1^{5}+2^{5}+3^{5}+\cdots +n^{5}={4a^{3}-a^{2} \over 3}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e3e2cfbf8830769853f626d48f3d1bc55ffcbe8)
![{\displaystyle 1^{7}+2^{7}+3^{7}+\cdots +n^{7}={12a^{4}-8a^{3}+2a^{2} \over 6}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cc09a591f0f54222e643ff480ec84a1887ad8e95)
![{\displaystyle 1^{9}+2^{9}+3^{9}+\cdots +n^{9}={16a^{5}-20a^{4}+12a^{3}-3a^{2} \over 5}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9691b106ed82848098cae4987092f7ddda4cdb48)
![{\displaystyle 1^{11}+2^{11}+3^{11}+\cdots +n^{11}={32a^{6}-64a^{5}+68a^{4}-40a^{3}+5a^{2} \over 6}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4d8e66e23f5989ab4513e17051f61720c6e7a611)
La primera de estas identidades es el teorema de Nicomachus. Algunos autores llaman a los polinomios de la derecha de estas identidades "polinomios de Faulhaber en a".
Método por Recurrencia Integral[editar]
Las sumas de las m-ésimas potencias de los primeros números naturales, están dadas por polinomios de grado m+1, y las fórmulas para potencias mayores se deducen y demuestran en [2], esto es la siguiente recurrencia integral:
![{\displaystyle {\frac {S_{m+1}(n)}{m+1}}=\int _{0}^{n}S_{m}(t)dt-n\int _{-1}^{0}S_{m}(t)dt,\quad {\text{ para todo }}m\geq 0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/685ba9aed3f5cce6cd892e87ff57e8b621516935)
Aplicando la recurrencia anterior, tras algunos cálculos, se obtienen los polinomios:
![{\displaystyle S_{0}(n)=n,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/decebed523538f0c941ea1ac87ee73918c2298c7)
![{\displaystyle S_{1}(n)={\cfrac {n^{2}}{2}}+{\cfrac {n}{2}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9c171045d5b2ad0cd7ba521f84f98273b4dd60f0)
![{\displaystyle S_{2}(n)={\cfrac {n^{3}}{3}}+{\cfrac {n^{2}}{2}}+{\cfrac {n}{6}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/28f568eab7d52d99ca7435b56c9f7d1518133b97)
![{\displaystyle S_{3}(n)={\cfrac {n^{4}}{4}}+{\cfrac {n^{3}}{2}}+{\cfrac {n^{2}}{4}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/502f4cd9fabfb8328b3fd8af60ede7906f317883)
![{\displaystyle S_{4}(n)={\cfrac {n^{5}}{5}}+{\cfrac {n^{4}}{2}}+{\cfrac {n^{3}}{3}}-{\cfrac {n}{30}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fb912fff3d61421aae392b8383172569f0bb11a5)
![{\displaystyle S_{5}(n)={\cfrac {n^{6}}{6}}+{\cfrac {n^{5}}{2}}+{\cfrac {5n^{4}}{12}}-{\cfrac {n^{2}}{12}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67d1efd87a8f906543ce49dd1ebbd3238a37e17f)
Referencias[editar]
- ↑ Conway,Guy 1996
- ↑ Hurtado Benavides, Miguel Ángel. (2020). De las sumas de potencias a las sucesiones de Appell y su caracterización a través de funcionales. [Tesis de maestría]. Universidad Sergio Arboleda.
- The Book of Numbers, John H. Conway, Richard Guy, Spring, 1998, ISBN 0-387-97993-X, page 107
- CRC Concise Encyclopedia of Mathematics, Eric Weisstein, Chapman & Hall/CRC, 2003, ISBN 1-58488-347-2, page 2331
- "Johann Faulhaber and Sums of Powers" por Donald Knuth
- MathWorld: urlname: FaulhabersFormula. Faulhaber's formula
- "Darinnen die miraculosische Inventiones zu den höchsten Cossen weiters continuirt und profitiert werden", Academia Algebrae, Johann Faulhaber, Augpurg, bey Johann Ulrich Schöigs, 1631.
- Schumacher, Raphael (2016). «An Extended Version of Faulhaber’s Formula». Journal of Integer Sequences 19.