Articulo de referencia

Número de Bernoulli

En matemáticas , los números de Bernoulli ''n'' "}},"i":0}}]}"> B n son una secuencia de números racionales que aparecen con frecuencia en el análisis . Los números de Bernoulli...

En matemáticas , los números de Bernoulli B n son una secuencia de números racionales que aparecen con frecuencia en el análisis . Los números de Bernoulli aparecen en (y pueden definirse mediante) las expansiones en serie de Taylor de las funciones tangente e hiperbólica , en la fórmula de Faulhaber para la suma de potencias m -ésimas de los primeros n enteros positivos , en la fórmula de Euler-Maclaurin y en expresiones para ciertos valores de la función zeta de Riemann .

Los valores de los primeros 20 números de Bernoulli se dan en la tabla adyacente. En la literatura se utilizan dos convenciones, denotadas aquí porBnorte{\displaystyle B_{n}^{-{}}}yBnorte+{\displaystyle B_{n}^{+{}}}; difieren solo para n = 1 , dondeB1=1/2{\displaystyle B_{1}^{-{}}=-1/2}yB1+=+1/2{\displaystyle B_{1}^{+{}}=+1/2}Para todo n impar > 1 , B n = 0. Para todo n par > 0 , B n es negativo si n es divisible por 4 y positivo en caso contrario. Los números de Bernoulli son valores especiales de los polinomios de Bernoulli.Bnorte(incógnita){\displaystyle B_{n}(x)}, conBnorte=Bnorte(0){\displaystyle B_{n}^{-{}}=B_{n}(0)}yBnorte+=Bnorte(1){\displaystyle B_{n}^{+}=B_{n}(1)}. [ 1 ]

Los números de Bernoulli fueron descubiertos casi al mismo tiempo por el matemático suizo Jacob Bernoulli , de quien reciben su nombre, e independientemente por el matemático japonés Seki Takakazu . El descubrimiento de Seki fue publicado póstumamente en 1712 [ 2 ] [ 3 ] [ 4 ] en su obra Katsuyō Sanpō ; el de Bernoulli, también póstumamente, en su Ars Conjectandi de 1713. La nota G de Ada Lovelace sobre la Máquina Analítica de 1842 describe un algoritmo para generar números de Bernoulli con la máquina de Babbage ; [ 5 ] se discute si Lovelace o Babbage desarrollaron el algoritmo . Como resultado, los números de Bernoulli tienen la distinción de ser el tema del primer programa informático complejo publicado .

Notación

El superíndice ± utilizado en este artículo distingue las dos convenciones de signos para los números de Bernoulli. Solo se ve afectado el término n = 1 :

En las fórmulas siguientes, se puede cambiar de una convención de signos a otra con la relaciónBnorte+=(1)norteBnorte{\displaystyle B_{n}^{+}=(-1)^{n}B_{n}^{-}}o, para enteros n = 2 o mayores, simplemente ignórelo.

Dado que B n = 0 para todo n impar > 1 , y muchas fórmulas solo involucran números de Bernoulli de índice par, algunos autores escriben " B n " en lugar de B 2 n  . Este artículo no sigue esa notación.

Historia

Historia temprana

Los números de Bernoulli tienen su origen en los inicios del cálculo de sumas de potencias enteras, que han sido de interés para los matemáticos desde la antigüedad.

Una página del Katsuyō Sanpō (1712) de Seki Takakazu , que tabula coeficientes binomiales y números de Bernoulli.

Se conocían métodos para calcular la suma de los primeros n enteros positivos, la suma de los cuadrados y de los cubos de los primeros n enteros positivos, pero no existían fórmulas propiamente dichas, sino solo descripciones verbales. Entre los grandes matemáticos de la antigüedad que abordaron este problema se encuentran Pitágoras (c. 572-497  a. C., Grecia), Arquímedes (287-212  a. C., Italia), Aryabhata (n. 476, India), Al-Karaji (m. 1019, Persia) e Ibn al-Haytham (965-1039, Irak).

Durante finales del siglo XVI y principios del XVII, los matemáticos lograron avances significativos. En Occidente, Thomas Harriot (1560-1621) de Inglaterra, Johann Faulhaber (1580-1635) de Alemania, Pierre de Fermat (1601-1665) y su colega matemático francés Blaise Pascal (1623-1662) desempeñaron papeles importantes.

Thomas Harriot parece haber sido el primero en derivar y escribir fórmulas para sumas de potencias utilizando notación simbólica, pero incluso él solo calculaba hasta la suma de las cuartas potencias. Johann Faulhaber proporcionó fórmulas para sumas de potencias hasta la decimoséptima potencia en su obra Academia Algebrae de 1631 , mucho más altas que las de cualquier otro autor anterior, pero no ofreció una fórmula general.

Blaise Pascal demostró en 1654 la identidad de Pascal que relaciona ( n +1) k +1 con las sumas de las potencias p de los primeros n enteros positivos para p = 0, 1, 2, ..., k .

El matemático suizo Jacob Bernoulli (1654–1705) fue el primero en darse cuenta de la existencia de una única secuencia de constantes B 0 , B 1 , B 2 ,... que proporciona una fórmula uniforme para todas las sumas de potencias. [ 9 ]

La alegría que experimentó Bernoulli al dar con el patrón necesario para calcular rápida y fácilmente los coeficientes de su fórmula para la suma de las potencias c -ésimas para cualquier entero positivo c se puede apreciar en su comentario. Escribió:

"Con la ayuda de esta tabla, me llevó menos de un cuarto de hora descubrir que la suma de las potencias décimas de los primeros 1000 números da como resultado 91.409.924.241.424.243.424.241.924.242.500."

El resultado de Bernoulli se publicó póstumamente en Ars Conjectandi en 1713. Seki Takakazu descubrió de forma independiente los números de Bernoulli y su resultado se publicó un año antes, también póstumamente, en 1712. [ 2 ] Sin embargo, Seki no presentó su método como una fórmula basada en una secuencia de constantes.

La fórmula de Bernoulli para sumas de potencias es la formulación más útil y generalizable hasta la fecha. Los coeficientes de la fórmula de Bernoulli se denominan ahora números de Bernoulli, siguiendo una sugerencia de Abraham de Moivre .

La fórmula de Bernoulli a veces se denomina fórmula de Faulhaber, en honor a Johann Faulhaber, quien halló métodos notables para calcular la suma de potencias, pero nunca enunció la fórmula de Bernoulli. Según Knuth [ 9 ], Carl Jacobi publicó por primera vez una demostración rigurosa de la fórmula de Faulhaber en 1834. [ 10 ] El estudio exhaustivo de Knuth sobre la fórmula de Faulhaber concluye:

"Faulhaber nunca descubrió los números de Bernoulli; es decir, nunca se dio cuenta de que una sola secuencia de constantes B 0 , B 1 , B 2 , ... proporcionaría una distribución uniforme
nortemetro=1metro+1(B0nortemetro+1(metro+11)B1nortemetro+(metro+12)B2nortemetro1+(1)metro(metro+1metro)Bmetronorte){\textstyle \sum n^{m}={\frac {1}{m+1}}\left(B_{0}n^{m+1}-{\binom {m+1}{1}}B_{1}n^{m}+{\binom {m+1}{2}}B_{2}n^{m-1}-\cdots +(-1)^{m}{\binom {m+1}{m}}B_{m}n\right)}
para todas las sumas de potencias. Nunca mencionó, por ejemplo, el hecho de que casi la mitad de los coeficientes resultaron ser cero después de haber convertido sus fórmulas para Σ n m de polinomios en N a polinomios en n ." [ 11 ]

En lo anterior, Knuth quería decirB1{\displaystyle B_{1}^{-}}; en su lugar, utiliceB1+{\displaystyle B_{1}^{+}}La fórmula evita la resta:

nortemetro=1metro+1(B0nortemetro+1+(metro+11)B1+nortemetro+(metro+12)B2nortemetro1++(metro+1metro)Bmetronorte).{\textstyle \sum n^{m}={\frac {1}{m+1}}\left(B_{0}n^{m+1}+{\binom {m+1}{1}}B_{1}^{+}n^{m}+{\binom {m+1}{2}}B_{2}n^{m-1}+\cdots +{\binom {m+1}{m}}B_{m}n\right).}

Reconstrucción de la "Summæ Potestatum"

... Atque sic porrò ad altiores gradatim potestates pergere, levique negotio sequentem adornare laterculum licet: Summæ Potestatumnorte=12nortenorte+12norte{\displaystyle \textstyle \int n={\frac {1}{2}}nn+{\frac {1}{2}}n}nortenorte=13norte3+12nortenorte+16norte{\displaystyle \textstyle \int nn={\frac {1}{3}}n^{3}+{\frac {1}{2}}nn+{\frac {1}{6}}n}norte3=14norte4+12norte3+14nortenorte{\displaystyle \textstyle \int n^{3}={\frac {1}{4}}n^{4}+{\frac {1}{2}}n^{3}+{\frac {1}{4}}nn}norte4=15norte5+12norte4+13norte3130norte{\displaystyle \textstyle \int n^{4}={\frac {1}{5}}n^{5}+{\frac {1}{2}}n^{4}+{\frac {1}{3}}n^{3}-{\frac {1}{30}}n}norte5=16norte6+12norte5+512norte4112nortenorte{\displaystyle \textstyle \int n^{5}={\frac {1}{6}}n^{6}+{\frac {1}{2}}n^{5}+{\frac {5}{12}}n^{4}-{\frac {1}{12}}nn}norte6=17norte7+12norte6+12norte516norte3+142norte{\displaystyle \textstyle \int n^{6}={\frac {1}{7}}n^{7}+{\frac {1}{2}}n^{6}+{\frac {1}{2}}n^{5}-{\frac {1}{6}}n^{3}+{\frac {1}{42}}n}norte7=18norte8+12norte7+712norte6724norte4+112nortenorte{\displaystyle \textstyle \int n^{7}={\frac {1}{8}}n^{8}+{\frac {1}{2}}n^{7}+{\frac {7}{12}}n^{6}-{\frac {7}{24}}n^{4}+{\frac {1}{12}}nn}norte8=19norte9+12norte8+23norte7715norte5+29norte3130norte{\displaystyle \textstyle \int n^{8}={\frac {1}{9}}n^{9}+{\frac {1}{2}}n^{8}+{\frac {2}{3}}n^{7}-{\frac {7}{15}}n^{5}+{\frac {2}{9}}n^{3}-{\frac {1}{30}}n}norte9=110norte10+12norte9+34norte8710norte6+12norte4112nortenorte{\displaystyle \textstyle \int n^{9}={\frac {1}{10}}n^{10}+{\frac {1}{2}}n^{9}+{\frac {3}{4}}n^{8}-{\frac {7}{10}}n^{6}+{\frac {1}{2}}n^{4}-{\frac {1}{12}}nn}norte10=111norte11+12norte10+56norte91norte7+1norte512norte3+566norte{\displaystyle \textstyle \int n^{10}={\frac {1}{11}}n^{11}+{\frac {1}{2}}n^{10}+{\frac {5}{6}}n^{9}-1n^{7}+1n^{5}-{\frac {1}{2}}n^{3}+{\frac {5}{66}}n}Quin imò qui legem progresivois inibi attentius inspexerit, eundem etiam continuare poterit absque his ratiociniorum ambagibus: Sumtâ enimdo{\displaystyle c}pro potestatis cujuslibet exponente, ajuste summa omniumnortedo{\displaystyle n^{c}}sunortedo=1do+1nortedo+1+12nortedo+do2Anortedo1+dodo1do2234Bnortedo3{\displaystyle \textstyle \int n^{c}={\frac {1}{c+1}}n^{c+1}+{\frac {1}{2}}n^{c}+{\frac {c}{2}}An^{c-1}+{\frac {c\cdot c-1\cdot c-2}{2\cdot 3\cdot 4}}Bn^{c-3}}+dodo1do2do3do423456donortedo5{\displaystyle \qquad +{\frac {c\cdot c-1\cdot c-2\cdot c-3\cdot c-4}{2\cdot 3\cdot 4\cdot 5\cdot 6}}Cn^{c-5}}+dodo1do2do3do4do5do62345678Dnortedo7{\displaystyle \qquad +{\frac {c\cdot c-1\cdot c-2\cdot c-3\cdot c-4\cdot c-5\cdot c-6}{2\cdot 3\cdot 4\cdot 5\cdot 6\cdot 7\cdot 8}}Dn^{c-7}\ldots }& ita deinceps, exponentem potestatis ipsiusnorte{\displaystyle n}continuè minuendo binario, quosque perveniatur adnorte{\displaystyle n}velocidadnortenorte{\displaystyle nn}. Literæ capitalesA{\displaystyle A},B{\displaystyle B},do{\displaystyle C},D{\displaystyle D}&do. ordine denotan coëfficientes ultimorum terminorum pronortenorte{\displaystyle \textstyle \int nn},norte4{\displaystyle \textstyle \int n^{4}},norte6{\displaystyle \textstyle \int n^{6}},norte8{\displaystyle \textstyle \int n^{8}}etc. nempeA=16,B=130,do=142,D=130.{\displaystyle \textstyle A={\frac {1}{6}},B=-{\frac {1}{30}},C={\frac {1}{42}},D=-{\frac {1}{30}}.} De Ars Conjectandi de Jacob Bernoulli , 1713 [ a ] [ 13 ]

Los números de Bernoulli OEIS : A164555  (n)/ OEIS : A027642  (n) fueron introducidos por Jacob Bernoulli en el libro Ars Conjectandi , publicado póstumamente en 1713. [ 13 ] La fórmula principal se puede ver en la segunda mitad del facsímil correspondiente . Los coeficientes constantes denotados A , B , C y D por Bernoulli se asignan a la notación que ahora prevalece como A = B 2 , B = B 4 , C = B 6 , D = B 8 . La expresión c · c −1· c −2· c −3 significa c ·( c −1)·( c −2)·( c −3) – los puntos pequeños se usan como símbolos de agrupación. Usando la terminología actual, estas expresiones son potencias factoriales descendentes c k . La notación factorial k ! El uso de ∫ como abreviatura de 1 × 2 × ... × k fue introducido mucho más tarde, en 1808, por Christian Kramp . El símbolo integral ∫ en el lado izquierdo se remonta a Gottfried Wilhelm Leibniz en 1675, quien lo utilizó como una letra S larga para "summa" (suma). [ 14 ] La letra n en el lado izquierdo no es un índice de sumatoria, sino que indica el límite superior del rango de sumatoria, que debe entenderse como 1, 2, ..., n . En resumen, para c positivo , hoy en día un matemático probablemente escribiría la fórmula de Bernoulli como:

k=1nortekdo=nortedo+1do+1+12nortedo+k=2doBkk¡dok1_nortedok+1.{\displaystyle \sum _{k=1}^{n}k^{c}={\frac {n^{c+1}}{c+1}}+{\frac {1}{2}}n^{c}+\sum _{k=2}^{c}{\frac {B_{k}}{k!}}c^{\underline {k-1}}n^{c-k+1}.}

Esta fórmula sugiere establecer B 1 = 1 / 2 al cambiar de la llamada enumeración "arcaica", que usa solo los índices pares 2, 4, 6... a la forma moderna. Lo más llamativo en este contexto es el hecho de que el factorial descendente c k −1 tiene para k = 0 el valor 1 / c + 1. [ 15 ] Por lo tanto , la fórmula de Bernoulli se puede escribir

k=1nortekdo=k=0doBkk¡dok1_nortedok+1{\displaystyle \sum _{k=1}^{n}k^{c}=\sum _{k=0}^{c}{\frac {B_{k}}{k!}}c^{\underline {k-1}}n^{c-k+1}}

si B 1 = 1/2 , recuperando el valor que Bernoulli le dio al coeficiente en esa posición.

La fórmula parak=1nortek9{\displaystyle \textstyle \sum _{k=1}^{n}k^{9}}En la página 97 del Ars Conjectandi de Bernoulli hay un error en el último término; debería ser320norte2{\displaystyle -{\tfrac {3}{20}}n^{2}}en lugar de112norte2{\displaystyle -{\tfrac {1}{12}}n^{2}}.

Definiciones

En los últimos 300 años se han encontrado numerosas caracterizaciones de los números de Bernoulli, y cada una podría utilizarse para introducirlos. Aquí solo se mencionan cuatro de las más útiles:

  • una ecuación recursiva,
  • una fórmula explícita,
  • una función generadora,
  • una expresión integral.

Para la demostración de la equivalencia de los cuatro enfoques, véase Ireland & Rosen (1990) o Conway & Guy (1996) .

Definición recursiva

Los números de Bernoulli obedecen las fórmulas de suma [ 1 ].

k=0metro(metro+1k)Bk=δmetro,0k=0metro(metro+1k)Bk+=metro+1{\displaystyle {\begin{aligned}\sum _{k=0}^{m}{\binom {m+1}{k}}B_{k}^{-{}}&=\delta _{m,0}\\\sum _{k=0}^{m}{\binom {m+1}{k}}B_{k}^{+{}}&=m+1\end{aligned}}}

dóndemetro=0,1,2...{\displaystyle m=0,1,2...}y δ denota la delta de Kronecker .

La primera de ellas a veces se escribe [ 16 ] como la fórmula (para m > 1) (B+1)metroBmetro=0,{\displaystyle (B+1)^{m}-B_{m}=0,} donde la potencia se expande formalmente utilizando el teorema del binomio yBk{\displaystyle B^{k}}es reemplazado porBk{\displaystyle B_{k}}.

Resolver paraBmetro{\displaystyle B_{m}^{\mp {}}}proporciona las fórmulas recursivas [ 17 ]

Bmetro=δmetro,0k=0metro1(metrok)Bkmetrok+1Bmetro+=1k=0metro1(metrok)Bk+metrok+1.{\displaystyle {\begin{aligned}B_{m}^{-{}}&=\delta _{m,0}-\sum _{k=0}^{m-1}{\binom {m}{k}}{\frac {B_{k}^{-{}}}{m-k+1}}\\B_{m}^{+}&=1-\sum _{k=0}^{m-1}{\binom {m}{k}}{\frac {B_{k}^{+}}{m-k+1}}.\end{aligned}}}

Definición explícita

En 1893, Louis Saalschütz enumeró un total de 38 fórmulas explícitas para los números de Bernoulli, [ 18 ] generalmente dando alguna referencia en la literatura anterior. Una de ellas es (parametro1{\displaystyle m\geq 1}):

Bmetro=k=0metro1k+1j=0k(kj)(1)jjmetroBmetro+=k=0metro1k+1j=0k(kj)(1)j(j+1)metro.{\displaystyle {\begin{aligned}B_{m}^{-}&=\sum _{k=0}^{m}{\frac {1}{k+1}}\sum _{j=0}^{k}{\binom {k}{j}}(-1)^{j}j^{m}\\B_{m}^{+}&=\sum _{k=0}^{m}{\frac {1}{k+1}}\sum _{j=0}^{k}{\binom {k}{j}}(-1)^{j}(j+1)^{m}.\end{aligned}}}

Función generadora

Las funciones generadoras exponenciales son

tmit1=t2(cotht21)=metro=0Bmetrotmetrometro¡tmitmit1=t1mit=t2(cotht2+1)=metro=0Bmetro+tmetrometro¡.{\displaystyle {\begin{alignedat}{3}{\frac {t}{e^{t}-1}}&={\frac {t}{2}}\left(\operatorname {coth} {\frac {t}{2}}-1\right)&&=\sum _{m=0}^{\infty }{\frac {B_{m}^{-{}}t^{m}}{m!}}\\{\frac {te^{t}}{e^{t}-1}}={\frac {t}{1-e^{-t}}}&={\frac {t}{2}}\left(\operatorname {coth} {\frac {t}{2}}+1\right)&&=\sum _{m=0}^{\infty }{\frac {B_{m}^{+}t^{m}}{m!}}.\end{alignedat}}}

donde la sustitución estt{\displaystyle t\to -t}. La diferencia aritmética entre las funciones generadoras paraBmetro+{\displaystyle B_{m}^{+}}yBmetro{\displaystyle B_{m}^{-}}es t .

La función generadora (ordinaria)

z1ψ1(z1)=metro=0Bmetro+zmetro{\displaystyle z^{-1}\psi _{1}(z^{-1})=\sum _{m=0}^{\infty }B_{m}^{+}z^{m}}

es una serie asintótica . Contiene la función trigamma ψ 1 .

Expresión integral

A partir de las funciones generadoras anteriores, se puede obtener la siguiente fórmula integral para los números de Bernoulli pares:

B2norte=4norte(1)norte+10t2norte1mi2πt1dt{\displaystyle B_{2n}=4n(-1)^{n+1}\int _{0}^{\infty }{\frac {t^{2n-1}}{e^{2\pi t}-1}}\mathrm {d} t}

Números de Bernoulli y la función zeta de Riemann

Los números de Bernoulli, utilizando 1/2 para B 1 , están relacionados con la función zeta de Riemann de los números reales negativos.
Valor absoluto de los números de Bernoulli de índice par y su relación con la función zeta de Riemann.

Los números de Bernoulli se pueden expresar en términos de la función zeta de Riemann :

Bnorte+=norteζ(1norte){\displaystyle B_{n}^{+}=-n\,\zeta (1-n)\quad }para n ≥ 1  .

Aquí el argumento de la función zeta es 0 o negativo. Comoζ(k){\displaystyle \zeta (k)}es cero para los enteros pares negativos (los ceros triviales ), si n>1 es impar,ζ(1norte){\displaystyle \zeta (1-n)}es cero.

Mediante la ecuación funcional zeta y la fórmula de reflexión gamma se puede obtener la siguiente relación: [ 19 ]

B2norte=(1)norte+12(2norte)¡(2π)2norteζ(2norte){\displaystyle B_{2n}={\frac {(-1)^{n+1}2(2n)!}{(2\pi )^{2n}}}\zeta (2n)\quad }para enteros n ≥ 1  .

Ahora el argumento de la función zeta es positivo.

De ζ 1 ( n ) y de la fórmula de Stirling se deduce que

|B2norte|4πnorte(norteπmi)2norte{\displaystyle |B_{2n}|\sim 4{\sqrt {\pi n}}\left({\frac {n}{\pi e}}\right)^{2n}\quad }para n   .

Cálculo eficiente de los números de Bernoulli

Diagrama que muestra un triángulo invertido formado por números. La fila superior consta de 1, 1/2, 1/3, 1/4, y así sucesivamente. Cada número de las filas inferiores tiene dos flechas que apuntan hacia él desde dos números de la fila superior, lo que indica que esos dos números se utilizan para calcular el número inferior. Los números del borde izquierdo del triángulo están rodeados por un círculo; estos son los números de Bernoulli.
Visualización del algoritmo de Akiyama-Tanigawa para el cálculo de los números de Bernoulli.

Akiyama y Tanigawa dan un simple "algoritmo triangular" (en relación con el triángulo de Pascal ) para calcular los números de Bernoulli. [ 20 ] Primero seab0,metro=1metro+1{\textstyle b_{0,m}={\frac {1}{m+1}}}(para m ≥ 0). Entonces, los términos sucesivos en el triángulo se pueden calcular con la relación de recurrencia.

bnorte+1,metro=(metro+1)(bnorte,metrobnorte,metro+1){\displaystyle b_{n+1,m}=(m+1)(b_{n,m}-b_{n,m+1})}

Los términosbnorte,0{\textstyle b_{n,0}}corresponden al n -ésimo número de Bernoulli B n . [ 21 ] [ 22 ]

Brent y Harvey presentan varios algoritmos para calcular los números de Bernoulli, incluyendo un algoritmo simple que es más rápido y utiliza menos espacio que el algoritmo de Akiyama-Tanigawa. Este algoritmo utiliza una recurrencia para calcular los números tangentes T n y aplica

Tnorte=(1)norte122norte(22norte1)B2norte2norte{\displaystyle T_{n}=(-1)^{n-1}2^{2n}(2^{2n}-1){\frac {B_{2n}}{2n}}}

para calcular los números de Bernoulli. Esto está relacionado con un método dado por Knuth y Buckholtz. [ 23 ] [ 24 ]

En algunas aplicaciones es útil poder calcular los números de Bernoulli B 0 a B p − 3 módulo p , donde p es un primo; por ejemplo, para comprobar si la conjetura de Vandiver se cumple para p , o incluso simplemente para determinar si p es un primo irregular . No es factible realizar tal cálculo utilizando las fórmulas recursivas anteriores, ya que se requerirían al menos (un múltiplo constante de) p 2 operaciones aritméticas. Afortunadamente, se han desarrollado métodos más rápidos [ 25 ] que requieren solo O ( p (log p ) 2 ) operaciones (véase la notación O grande ).

David Harvey [ 26 ] describe un algoritmo para calcular los números de Bernoulli calculando B n módulo p para muchos primos pequeños p , y luego reconstruyendo B n mediante el teorema chino del resto . Harvey escribe que la complejidad temporal asintótica de este algoritmo es O ( log( n ) ² + ε ) y afirma que esta implementación es significativamente más rápida que las implementaciones basadas en otros métodos. Usando esta implementación , Harvey calculó B n para n = 10⁸ . La implementación de Harvey se ha incluido en SageMath desde la versión 3.1. Antes de eso, Bernd Kellner [ 27 ] calculó B n con precisión completa para n = 10⁶ en diciembre de 2002 y Oleksandr Pavlyk [ 28 ] para n = 10⁷ con Mathematica en abril de 2008.

* Digits debe entenderse como el exponente de 10 cuando B n se escribe como un número real en notación científica normalizada .

Aplicaciones de los números de Bernoulli

Análisis asintótico

Podría decirse que la aplicación más importante de los números de Bernoulli en matemáticas es su uso en la fórmula de Euler-Maclaurin . Suponiendo que f es una función suficientemente diferenciable, la fórmula de Euler-Maclaurin se puede escribir como [ 29 ].

k=ab1F(k)=abF(incógnita)dincógnita+k=1metroBkk¡(F(k1)(b)F(k1)(a))+R(F,metro).{\displaystyle \sum _{k=a}^{b-1}f(k)=\int _{a}^{b}f(x)\,dx+\sum _{k=1}^{m}{\frac {B_{k}^{-}}{k!}}(f^{(k-1)}(b)-f^{(k-1)}(a))+R_{-}(f,m).}

Esta formulación asume la convención B 1 = − 1 / 2 . Usando la convención B + 1 = + 1 / 2 la fórmula se convierte en

k=a+1bF(k)=abF(incógnita)dincógnita+k=1metroBk+k¡(F(k1)(b)F(k1)(a))+R+(F,metro).{\displaystyle \sum _{k=a+1}^{b}f(k)=\int _{a}^{b}f(x)\,dx+\sum _{k=1}^{m}{\frac {B_{k}^{+}}{k!}}(f^{(k-1)}(b)-f^{(k-1)}(a))+R_{+}(f,m).}

AquíF(0)=F{\displaystyle f^{(0)}=f}(es decir, la derivada de orden cero deF{\displaystyle f}es soloF{\displaystyle f}). Además, dejemosF(1){\displaystyle f^{(-1)}}denota una antiderivada deF{\displaystyle f}. Por el teorema fundamental del cálculo ,

abF(incógnita)dincógnita=F(1)(b)F(1)(a).{\displaystyle \int _{a}^{b}f(x)\,dx=f^{(-1)}(b)-f^{(-1)}(a).}

Así, la última fórmula se puede simplificar aún más a la siguiente forma concisa de la fórmula de Euler-Maclaurina.

k=a+1bF(k)=k=0metroBkk¡(F(k1)(b)F(k1)(a))+R(F,metro).{\displaystyle \sum _{k=a+1}^{b}f(k)=\sum _{k=0}^{m}{\frac {B_{k}}{k!}}(f^{(k-1)}(b)-f^{(k-1)}(a))+R(f,m).}

Esta forma es, por ejemplo, la fuente de la importante expansión de Euler-Maclaurin de la función zeta.

ζ(s)=k=0metroBk+k¡sk1¯+R(s,metro)=B00¡s1¯+B1+1¡s0¯+B22¡s1¯++R(s,metro)=1s1+12+112s++R(s,metro).{\displaystyle {\begin{aligned}\zeta (s)&=\sum _{k=0}^{m}{\frac {B_{k}^{+}}{k!}}s^{\overline {k-1}}+R(s,m)\\&={\frac {B_{0}}{0!}}s^{\overline {-1}}+{\frac {B_{1}^{+}}{1!}}s^{\overline {0}}+{\frac {B_{2}}{2!}}s^{\overline {1}}+\cdots +R(s,m)\\&={\frac {1}{s-1}}+{\frac {1}{2}}+{\frac {1}{12}}s+\cdots +R(s,m).\end{aligned}}}

Aquí s k denota la potencia factorial creciente . [ 30 ]

Los números de Bernoulli también se utilizan con frecuencia en otros tipos de expansiones asintóticas . El siguiente ejemplo es la expansión asintótica clásica de tipo Poincaré de la función digamma ψ .

ψ(z)lnzk=1Bk+kzk{\displaystyle \psi (z)\sim \ln z-\sum _{k=1}^{\infty }{\frac {B_{k}^{+}}{kz^{k}}}}

Suma de potencias

Los números de Bernoulli ocupan un lugar destacado en la expresión en forma cerrada de la suma de las m -ésimas potencias de los primeros n enteros positivos. Para m , n ≥ 0 , definimos

Smetro(norte)=k=1nortekmetro=1metro+2metro++nortemetro.{\displaystyle S_{m}(n)=\sum _{k=1}^{n}k^{m}=1^{m}+2^{m}+\cdots +n^{m}.}

Esta expresión siempre se puede reescribir como un polinomio en n de grado m + 1. Los coeficientes de estos polinomios están relacionados con los números de Bernoulli mediante la fórmula de Bernoulli :

Smetro(norte)=1metro+1k=0metro(metro+1k)Bk+nortemetro+1k=metro¡k=0metroBk+nortemetro+1kk¡(metro+1k)¡,{\displaystyle S_{m}(n)={\frac {1}{m+1}}\sum _{k=0}^{m}{\binom {m+1}{k}}B_{k}^{+}n^{m+1-k}=m!\sum _{k=0}^{m}{\frac {B_{k}^{+}n^{m+1-k}}{k!(m+1-k)!}},}

donde ( m + 1 k ) denota el coeficiente binomial .

Por ejemplo, tomando m igual a 1 se obtienen los números triangulares 0, 1, 3, 6, ... OEIS : A000217  .

1+2++norte=12(B0norte2+2B1+norte1)=12(norte2+norte).{\displaystyle 1+2+\cdots +n={\frac {1}{2}}(B_{0}n^{2}+2B_{1}^{+}n^{1})={\tfrac {1}{2}}(n^{2}+n).}

Tomando m igual a 2 se obtienen los números piramidales cuadrados 0, 1, 5, 14, ... OEIS : A000330  .

12+22++norte2=13(B0norte3+3B1+norte2+3B2norte1)=13(norte3+32norte2+12norte).{\displaystyle 1^{2}+2^{2}+\cdots +n^{2}={\frac {1}{3}}(B_{0}n^{3}+3B_{1}^{+}n^{2}+3B_{2}n^{1})={\tfrac {1}{3}}\left(n^{3}+{\tfrac {3}{2}}n^{2}+{\tfrac {1}{2}}n\right).}

Algunos autores utilizan la convención alternativa para los números de Bernoulli y enuncian la fórmula de Bernoulli de esta manera:

Smetro(norte)=1metro+1k=0metro(1)k(metro+1k)Bknortemetro+1k.{\displaystyle S_{m}(n)={\frac {1}{m+1}}\sum _{k=0}^{m}(-1)^{k}{\binom {m+1}{k}}B_{k}^{-{}}n^{m+1-k}.}

La fórmula de Bernoulli a veces se denomina fórmula de Faulhaber, en honor a Johann Faulhaber, quien también encontró formas notables de calcular sumas de potencias .

La fórmula de Faulhaber fue generalizada por V. Guo y J. Zeng a un q -análogo . [ 31 ]

Serie Taylor

Los números de Bernoulli aparecen en el desarrollo en serie de Taylor de muchas funciones trigonométricas e hiperbólicas .

broncearseincógnita=1incógnitanorte=1(1)norte122norte(22norte1)B2norte(2norte)¡incógnita2norte1,|incógnita|<π2.cunaincógnita=1incógnitanorte=0(1)norteB2norte(2incógnita)2norte(2norte)¡,0<|incógnita|<π.tanhincógnita=1incógnitanorte=122norte(22norte1)B2norte(2norte)¡incógnita2norte1,|incógnita|<π2.cothincógnita=1incógnitanorte=0B2norte(2incógnita)2norte(2norte)¡,0<|incógnita|<π.{\displaystyle {\begin{aligned}\tan x&={\hphantom {1 \over x}}\sum _{n=1}^{\infty }{\frac {(-1)^{n-1}2^{2n}(2^{2n}-1)B_{2n}}{(2n)!}}\;x^{2n-1},&&\left|x\right|<{\frac {\pi }{2}}.\\\cot x&={1 \over x}\sum _{n=0}^{\infty }{\frac {(-1)^{n}B_{2n}(2x)^{2n}}{(2n)!}},&0<&|x|<\pi .\\\tanh x&={\hphantom {1 \over x}}\sum _{n=1}^{\infty }{\frac {2^{2n}(2^{2n}-1)B_{2n}}{(2n)!}}\;x^{2n-1},&&|x|<{\frac {\pi }{2}}.\\\coth x&={1 \over x}\sum _{n=0}^{\infty }{\frac {B_{2n}(2x)^{2n}}{(2n)!}},&0<&|x|<\pi .\end{aligned}}}

Serie Laurent

Los números de Bernoulli aparecen en la siguiente serie de Laurent : [ 32 ]

Función Digamma :ψ(z)=lnzk=1Bk+kzk{\displaystyle \psi (z)=\ln z-\sum _{k=1}^{\infty }{\frac {B_{k}^{+{}}}{kz^{k}}}}

Uso en topología

La fórmula de Kervaire-Milnor para el orden del grupo cíclico de clases de difeomorfismos de (4 n − 1) -esferas exóticas que delimitan variedades paralelizable involucra números de Bernoulli. Sea ES n el número de tales esferas exóticas para n ≥ 2 , entonces

ESnorte=(22norte224norte3)Numerador(B4norte4norte).{\displaystyle {\textit {ES}}_{n}=(2^{2n-2}-2^{4n-3})\operatorname {Numerator} \left({\frac {B_{4n}}{4n}}\right).}

El teorema de la signatura de Hirzebruch para el género L de una variedad cerrada orientada y lisa de dimensión 4 n también involucra números de Bernoulli.

Conexiones con números combinatorios

La conexión del número de Bernoulli con varios tipos de números combinatorios se basa en la teoría clásica de las diferencias finitas y en la interpretación combinatoria de los números de Bernoulli como un ejemplo de un principio combinatorio fundamental, el principio de inclusión-exclusión .

Conexión con los números de Worpitzky

La definición que se utilizará fue desarrollada por Julius Worpitzky en 1883. Además de la aritmética elemental, solo se emplean la función factorial n ! y la función potencia k m . Los números de Worpitzky sin signo se definen como

Wnorte,k=v=0k(1)v+k(v+1)norte(kv) .{\displaystyle W_{n,k}=\sum _{v=0}^{k}(-1)^{v+k}(v+1)^{n}{k \choose v}\ .}

También pueden expresarse mediante los números de Stirling de segunda especie.

Wnorte,k=k¡{norte+1k+1}.{\displaystyle W_{n,k}=k!\left\{{n+1 \atop k+1}\right\}.}

A continuación , se introduce un número de Bernoulli como una suma de inclusión-exclusión de números de Worpitzky ponderados por la secuencia armónica 1 , 1/2 , 1/3 , ...   

Bnorte=k=0norte(1)kWnorte,kk+1 = k=0norte1k+1v=0k(1)v(v+1)norte(kv) .{\displaystyle B_{n}=\sum _{k=0}^{n}(-1)^{k}{\frac {W_{n,k}}{k+1}}\ =\ \sum _{k=0}^{n}{\frac {1}{k+1}}\sum _{v=0}^{k}(-1)^{v}(v+1)^{n}{k \choose v}\ .}
B 0 = 1
B 1 = 11 / 2
B 2 = 1 − 3 / 2 + 2 / 3
B 3 = 1 − 7 / 2 + 12 / 36 / 4
B 4 = 1 − 15 / 2 + 50 / 360 / 4 + 24 / 5
B 5 = 1 − 31 / 2 + 180 / 3390 / 4 + 360 / 5120 / 6
B 6 = 1 − 63 / 2 + 602 / 32100 / 4 + 3360 / 52520 / 6 + 720 / 7

Esta representación tiene B + 1 = + 1 / 2 .

Consideremos la secuencia s n , n ≥ 0 . A partir de los números de Worpitzky OEIS : A028246  , OEIS : A163626  aplicados a s 0 , s 0 , s 1 , s 0 , s 1 , s 2 , s 0 , s 1 , s 2 , s 3 , ... es idéntico a la transformada de Akiyama-Tanigawa aplicada a s n (véase Conexión con los números de Stirling de primera especie ). Esto se puede ver en la tabla:

La primera fila representa s 0 , s 1 , s 2 , s 3 , s 4 .

Por lo tanto, para los segundos números fraccionarios de Euler OEIS : A198631  ( n ) / OEIS : A006519  ( n + 1 ):

E 0 = 1
E 1 = 11 / 2
mi 2 = 1 − 3 / 2 + 2 / 4
mi 3 = 1 − 7 / 2 + 12 / 46 / 8
mi 4 = 1 − 15 / 2 + 50 / 460 / 8 + 24 / 16
mi 5 = 1 − 31 / 2 + 180 / 4390 / 8 + 360 / 16120 / 32
mi 6 = 1 − 63 / 2 + 602 / 42100 / 8 + 3360 / 162520 / 32 + 720 / 64

Una segunda fórmula que representa los números de Bernoulli mediante los números de Worpitzky es para n ≥ 1

Bnorte=norte2norte+12k=0norte1(2)kWnorte1,k.{\displaystyle B_{n}={\frac {n}{2^{n+1}-2}}\sum _{k=0}^{n-1}(-2)^{-k}\,W_{n-1,k}.}

La segunda representación simplificada de Worpitzky de los segundos números de Bernoulli es:

OEIS : A164555  ( n + 1 ) / OEIS : A027642  ( n + 1 ) =n + 1 / 2 n + 2 − 2 × OEIS : A198631 ( n ) / OEIS : A006519 ( n + 1 )  

que vincula los segundos números de Bernoulli con los segundos números fraccionarios de Euler. El comienzo es:

1/2 , 1/6 , 0 , -1/30 ,0 , 1/42 , ... = ( 1/2 , 1/3 , 3/14 , 2/15 , 5/62 , 1/21 , ... ) × ( 1 , 1/2 , 0 , -1/4 , 0 , 1/2 , ... )

Los numeradores del primer paréntesis son OEIS : A111701  (véase Conexión con los números Stirling de primera clase ).

Relación con los números Stirling de segundo tipo

Los números de Stirling de segundo tipo , S ( k , m ) , tienen la propiedad de que [ 33 ] :

incógnitak=metro=0kincógnitametro_S(k,metro){\displaystyle x^{k}=\sum _{m=0}^{k}{x^{\underline {m}}}S(k,m)}

donde x m denota la función factorial descendente .

Los polinomios de Bernoulli B k ( x ) se pueden escribir como: [ 34 ]

Bk(incógnita)=kmetro=0k1(incógnitametro+1)S(k1,metro)metro¡+Bk{\displaystyle B_{k}(x)=k\sum _{m=0}^{k-1}{\binom {x}{m+1}}S(k-1,m)m!+B_{k}}

donde B k para k = 0, 1, 2,... son los números de Bernoulli.

La siguiente propiedad del coeficiente binomial :

(incógnitametro)=(incógnita+1metro+1)(incógnitametro+1){\displaystyle {\binom {x}{m}}={\binom {x+1}{m+1}}-{\binom {x}{m+1}}}

esto implica que

incógnitak=Bk+1(incógnita+1)Bk+1(incógnita)k+1.{\displaystyle x^{k}={\frac {B_{k+1}(x+1)-B_{k+1}(x)}{k+1}}.}

También se tiene lo siguiente para los polinomios de Bernoulli, [ 34 ]

Bk(incógnita)=norte=0k(knorte)Bnorteincógnitaknorte.{\displaystyle B_{k}(x)=\sum _{n=0}^{k}{\binom {k}{n}}B_{n}x^{k-n}.}

El coeficiente de x en ( x m + 1 ) es (−1) m / m + 1 .

El coeficiente de x en la primera expresión paraBk(incógnita){\displaystyle B_{k}(x)}eskmetro=0k1(1)metrometro+1S(k1,metro)metro¡{\displaystyle k\sum _{m=0}^{k-1}{\frac {(-1)^{m}}{m+1}}S(k-1,m)m!}mientras que en la segunda expresión eskBk1.{\displaystyle kB_{k-1}.}Reemplazark1{\displaystyle k-1}pork{\displaystyle k}Esto produce:

Bk=metro=0k(1)metrometro¡metro+1S(k,metro){\displaystyle B_{k}=\sum _{m=0}^{k}(-1)^{m}{\frac {m!}{m+1}}S(k,m)}

(resultando en B 1 = + 1 / 2 ) que es una fórmula explícita para los números de Bernoulli y puede usarse para demostrar el teorema de Von-Staudt Clausen . [ 35 ] [ 36 ] [ 37 ]

Conexión con los números Stirling de primera clase.

Las dos fórmulas principales que relacionan los números de Stirling sin signo de primera especie [ n m ] con los números de Bernoulli (con B 1 = + 1 / 2 ) son

1metro¡k=0metro(1)k[metro+1k+1]Bk=1metro+1,{\displaystyle {\frac {1}{m!}}\sum _{k=0}^{m}(-1)^{k}\left[{m+1 \atop k+1}\right]B_{k}={\frac {1}{m+1}},}

y la inversión de esta suma (para n ≥ 0 , m ≥ 0 )

1metro¡k=0metro(1)k[metro+1k+1]Bnorte+k=Anorte,metro.{\displaystyle {\frac {1}{m!}}\sum _{k=0}^{m}(-1)^{k}\left[{m+1 \atop k+1}\right]B_{n+k}=A_{n,m}.}

Aquí, los números A n , m son los números racionales de Akiyama-Tanigawa, los primeros de los cuales se muestran en la siguiente tabla.

Los números de Akiyama-Tanigawa satisfacen una relación de recurrencia simple que puede utilizarse para calcular iterativamente los números de Bernoulli. Véase OEIS : A051714  / OEIS : A051715  .

Una autosecuencia es una secuencia cuya transformada binomial inversa es igual a la secuencia con signo. Si la diagonal principal es ceros = OEIS : A000004  , la autosecuencia es de primera especie. Ejemplo: OEIS : A000045  , los números de Fibonacci. Si la diagonal principal es la primera diagonal superior multiplicada por 2, es de segunda especie. Ejemplo: OEIS : A164555  / OEIS : A027642  , los segundos números de Bernoulli (véase OEIS : A190339  ). La transformada de Akiyama-Tanigawa aplicada a 2 n = 1/ OEIS : A000079  conduce a OEIS : A198631  ( n ) / OEIS : A06519  ( n + 1). Por lo tanto:

Consulte OEIS : A209308  y OEIS : A227577  . OEIS : A198631  ( n ) / OEIS : A006519  ( n + 1 ) son los segundos números de Euler (fraccionales) y una autosecuencia de segundo tipo.

( OEIS : A164555 ( n + 2 ) / OEIS : A027642 ( n + 2 ) = 1 / 6 , 0, − 1 / 30 , 0, 1 / 42 , ... ) × ( 2 n + 3 − 2 / n + 2 = 3, 14 / 3 , 15 / 2 , 62 / 5 , 21, ... ) = OEIS : A198631 ( n + 1 ) / OEIS : A006519 ( n + 2 ) = 1 / 2 , 0, − 1 / 4 , 0, 1 / 2 , ... .    

También valioso para OEIS : A027641  / OEIS : A027642  (ver Conexión con los números de Worpitzky ).

Conexión con el triángulo de Pascal

Hay fórmulas que conectan el triángulo de Pascal con los números de Bernoulli [ b ].

Bnorte+=|Anorte|(norte+1)¡   {\displaystyle B_{n}^{+}={\frac {|A_{n}|}{(n+1)!}}~~~}con definiciones de filas y columnas[Anorte]i,k:={0si k>1+i(i+1k1)de lo contrario{\displaystyle [A_{n}]_{i,k}:={\begin{cases}0&{\text{if }}k>1+i\\{i+1 \choose k-1}&{\text{otherwise}}\end{cases}}}

dónde|Anorte|{\displaystyle |A_{n}|}es el determinante de una matriz de Hessenberg de n x n que forma parte del triángulo de Pascal.

Ejemplo:

B6+=|120000133000146400151010501615201561721353521|7¡=1205040=142{\displaystyle B_{6}^{+}={\frac {\begin{vmatrix}1&2&0&0&0&0\\1&3&3&0&0&0\\1&4&6&4&0&0\\1&5&10&10&5&0\\1&6&15&20&15&6\\1&7&21&35&35&21\end{vmatrix}}{7!}}={\frac {120}{5040}}={\frac {1}{42}}}[ 38 ]

Relación con los números eulerianos

Existen fórmulas que conectan los números eulerianos n m con los números de Bernoulli:

metro=0norte(1)metronortemetro=2norte+1(2norte+11)Bnorte+1norte+1,metro=0norte(1)metronortemetro(nortemetro)1=(norte+1)Bnorte.{\displaystyle {\begin{aligned}\sum _{m=0}^{n}(-1)^{m}\left\langle {n \atop m}\right\rangle &=2^{n+1}(2^{n+1}-1){\frac {B_{n+1}}{n+1}},\\\sum _{m=0}^{n}(-1)^{m}\left\langle {n \atop m}\right\rangle {\binom {n}{m}}^{-1}&=(n+1)B_{n}.\end{aligned}}}

Ambas fórmulas son válidas para n ≥ 0 si B 1 se establece en 1 / 2 . Si B 1 se establece en − 1 / 2 ⁠, son válidas solo para n ≥ 1 y n ≥ 2 respectivamente.

Una representación en árbol binario

Los polinomios de Stirling σ n ( x ) están relacionados con los números de Bernoulli mediante B n = n ! σ n (1) . SC Woon describió un algoritmo para calcular σ n (1) como un árbol binario: [ 39 ]

El algoritmo recursivo de Woon (para n ≥ 1 ) comienza asignando al nodo raíz N = [1,2] . Dado un nodo N = [ a 1 , a 2 , ..., a k ] del árbol, el hijo izquierdo del nodo es L ( N ) = [− a 1 , a 2 + 1, a 3 , ..., a k ] y el hijo derecho R ( N ) = [ a 1 , 2, a 2 , ..., a k ] . Un nodo N = [ a 1 , a 2 , ..., a k ] se escribe como ±[ a 2 , ..., a k ] en la parte inicial del árbol representado arriba, donde ± denota el signo de a 1 .

Dado un nodo N, el factorial de N se define como

norte¡=a1k=2longitud(norte)ak¡.{\displaystyle N!=a_{1}\prod _{k=2}^{\operatorname {length} (N)}a_{k}!.}

Restringido a los nodos N de un nivel de árbol fijo n, la suma de 1 / N ! es σ n (1) , por lo tanto

Bnorte= nivel de árbol nortenorte nodo denorte¡norte¡.{\displaystyle B_{n}=\sum _{\stackrel {N{\text{ node of}}}{{\text{ tree-level }}n}}{\frac {n!}{N!}}.}

Por ejemplo:

B 1 = 1!( 1 / 2! )
B 2 = 2!(− 1 / 3! + 1 / 2!2! )
B 3 = 3!( 1 / 4!1 / 2!3!1 / 3!2! + 1 / 2!2!2! )

Representación integral y continuación

El integral

b(s)=2misiπ/20sts1mi2πtdtt=s¡2s1ζ(s)πs(i)s=2s¡ζ(s)(2πi)s{\displaystyle b(s)=2e^{si\pi /2}\int _{0}^{\infty }{\frac {st^{s}}{1-e^{2\pi t}}}{\frac {dt}{t}}={\frac {s!}{2^{s-1}}}{\frac {\zeta (s)}{{}\pi ^{s}{}}}(-i)^{s}={\frac {2s!\zeta (s)}{(2\pi i)^{s}}}}

tiene como valores especiales b (2 n ) = B 2 n para n > 0 .

Por ejemplo, b (3) = 3 / 2 ζ (3) π −3 i y b (5) = − 15 / 2 ζ (5) π −5 i . Aquí, ζ es la función zeta de Riemann e i es la unidad imaginaria . Leonhard Euler ( Opera Omnia , Ser. 1, Vol. 10, p.  351) consideró estos números y calculó

pag=32π3(1+123+133+)=0,0581522q=152π5(1+125+135+)=0,0254132{\displaystyle {\begin{aligned}p&={\frac {3}{2\pi ^{3}}}\left(1+{\frac {1}{2^{3}}}+{\frac {1}{3^{3}}}+\cdots \right)=0.0581522\ldots \\q&={\frac {15}{2\pi ^{5}}}\left(1+{\frac {1}{2^{5}}}+{\frac {1}{3^{5}}}+\cdots \right)=0.0254132\ldots \end{aligned}}}

Otra representación integral similar es

b(s)=misiπ/22s10stssinhπtdtt=2misiπ/22s10miπtsts1mi2πtdtt.{\displaystyle b(s)=-{\frac {e^{si\pi /2}}{2^{s}-1}}\int _{0}^{\infty }{\frac {st^{s}}{\sinh \pi t}}{\frac {dt}{t}}={\frac {2e^{si\pi /2}}{2^{s}-1}}\int _{0}^{\infty }{\frac {e^{\pi t}st^{s}}{1-e^{2\pi t}}}{\frac {dt}{t}}.}

La relación con los números de Euler y π

Los números de Euler son una secuencia de enteros íntimamente relacionados con los números de Bernoulli. Al comparar las expansiones asintóticas de los números de Bernoulli y de Euler, se observa que los números de Euler E²ⁿ son aproximadamente 2 / π ( 4²ⁿ - 2²ⁿ ) veces mayores que los números de Bernoulli B²ⁿ . En consecuencia :

π2(22norte42norte)B2nortemi2norte.{\displaystyle \pi \sim 2(2^{2n}-4^{2n}){\frac {B_{2n}}{E_{2n}}}.}

Esta ecuación asintótica revela que π se encuentra en la raíz común de los números de Bernoulli y Euler. De hecho, π podría calcularse a partir de estas aproximaciones racionales.

Los números de Bernoulli se pueden expresar mediante los números de Euler y viceversa. Dado que, para n impar , B n = E n = 0 (con la excepción de B 1 ), basta con considerar el caso en que n es par.

Bnorte=k=0norte1(norte1k)norte4norte2nortemiknorte=2,4,6,minorte=k=1norte(nortek1)2k4kkBknorte=2,4,6,{\displaystyle {\begin{aligned}B_{n}&=\sum _{k=0}^{n-1}{\binom {n-1}{k}}{\frac {n}{4^{n}-2^{n}}}E_{k}&n&=2,4,6,\ldots \\[6pt]E_{n}&=\sum _{k=1}^{n}{\binom {n}{k-1}}{\frac {2^{k}-4^{k}}{k}}B_{k}&n&=2,4,6,\ldots \end{aligned}}}

Estas fórmulas de conversión expresan una conexión entre los números de Bernoulli y los de Euler. Pero, más importante aún, existe una raíz aritmética profunda común a ambos tipos de números, que puede expresarse a través de una secuencia de números más fundamental, también estrechamente ligada a π . Estos números se definen para n ≥ 1 como [ 40 ] [ 41 ].

Snorte=2(2π)nortek=0(1)knorte(2k+1)norte=2(2π)nortelímiteKk=KK(4k+1)norte.{\displaystyle S_{n}=2\left({\frac {2}{\pi }}\right)^{n}\sum _{k=0}^{\infty }{\frac {(-1)^{kn}}{(2k+1)^{n}}}=2\left({\frac {2}{\pi }}\right)^{n}\lim _{K\to \infty }\sum _{k=-K}^{K}(4k+1)^{-n}.}

La magia de estos números reside en el hecho de que resultan ser números racionales. Esto fue demostrado por primera vez por Leonhard Euler en un artículo fundamental, De summis serierum reciprocarum (Sobre las sumas de series de recíprocos), y ha fascinado a los matemáticos desde entonces. [ 42 ] Los primeros de estos números son:

Snorte=1,1,12,13,524,215,61720,17315,2778064,622835,{\displaystyle S_{n}=1,1,{\frac {1}{2}},{\frac {1}{3}},{\frac {5}{24}},{\frac {2}{15}},{\frac {61}{720}},{\frac {17}{315}},{\frac {277}{8064}},{\frac {62}{2835}},\ldots }( OEIS : A099612  / OEIS : A099617  )

Estos son los coeficientes en la expansión de sec x + tan x .

Los números de Bernoulli y los números de Euler pueden entenderse como vistas especiales de estos números, seleccionados de la secuencia S n y escalados para su uso en aplicaciones especiales.

Bnorte=(1)norte2[norte incluso]norte¡2norte4norteSnorte ,norte=2,3,minorte=(1)norte2[norte incluso]norte¡Snorte+1norte=0,1,{\displaystyle {\begin{aligned}B_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]{\frac {n!}{2^{n}-4^{n}}}\,S_{n}\ ,&n&=2,3,\ldots \\E_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]n!\,S_{n+1}&n&=0,1,\ldots \end{aligned}}}

La expresión [ n par] tiene el valor 1 si n es par y 0 en caso contrario ( corchete de Iverson ).

Estas identidades muestran que el cociente de los números de Bernoulli y Euler al comienzo de esta sección es solo el caso especial de R n = 2 S n / S n + 1 cuando n es par. Los R n son aproximaciones racionales a π y dos términos sucesivos siempre encierran el valor verdadero de π . Comenzando con n = 1, la secuencia comienza ( OEIS : A132049  / OEIS : A132050  ):

2,4,3,165,258,19261,427136,43521385,124653968,15872050521,π.{\displaystyle 2,4,3,{\frac {16}{5}},{\frac {25}{8}},{\frac {192}{61}},{\frac {427}{136}},{\frac {4352}{1385}},{\frac {12465}{3968}},{\frac {158720}{50521}},\ldots \quad \longrightarrow \pi .}

Estos números racionales también aparecen en el último párrafo del artículo de Euler citado anteriormente.

Consideremos la transformación de Akiyama-Tanigawa para la secuencia OEIS : A046978  ( n + 2 ) / OEIS : A016116  ( n + 1 ):

A partir del segundo, los numeradores de la primera columna son los denominadores de la fórmula de Euler. La primera columna es − 1 / 2 × OEIS : A163982 . 

Una perspectiva algorítmica: el triángulo de Seidel

La secuencia S n tiene otra propiedad inesperada pero importante: los denominadores de S n +1 dividen al factorial n ! . En otras palabras: los números T n  = S n + 1 n !  , a veces llamados números en zigzag de Euler , son enteros.

Tnorte=1,1,1,2,5,16,61,272,1385,7936,50521,353792,norte=0,1,2,3,{\displaystyle T_{n}=1,\,1,\,1,\,2,\,5,\,16,\,61,\,272,\,1385,\,7936,\,50521,\,353792,\ldots \quad n=0,1,2,3,\ldots }( OEIS : A000111  ). Véase ( OEIS : A253671  ).

Su función generadora exponencial es la suma de las funciones secante y tangente .

norte=0Tnorteincógnitanortenorte¡=broncearse(π4+incógnita2)=segundoincógnita+broncearseincógnita{\displaystyle \sum _{n=0}^{\infty }T_{n}{\frac {x^{n}}{n!}}=\tan \left({\frac {\pi }{4}}+{\frac {x}{2}}\right)=\sec x+\tan x}.

Así, las representaciones anteriores de los números de Bernoulli y Euler pueden reescribirse en términos de esta secuencia como

Bnorte=(1)norte2[norte incluso]norte2norte4norteTnorte1 norte2minorte=(1)norte2[norte incluso]Tnortenorte0{\displaystyle {\begin{aligned}B_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]{\frac {n}{2^{n}-4^{n}}}\,T_{n-1}\ &n&\geq 2\\E_{n}&=(-1)^{\left\lfloor {\frac {n}{2}}\right\rfloor }[n{\text{ even}}]T_{n}&n&\geq 0\end{aligned}}}

Estas identidades facilitan el cálculo de los números de Bernoulli y Euler: los números de Euler E 2 n se obtienen inmediatamente a partir de T 2 n y los números de Bernoulli B 2 n son fracciones que se obtienen a partir de T 2 n - 1 mediante un sencillo cambio de términos, evitando la aritmética racional.

Lo que queda es encontrar una manera conveniente de calcular los números T n . Sin embargo, ya en 1877 Philipp Ludwig von Seidel publicó un ingenioso algoritmo que simplifica el cálculo de T n . [ 43 ]

1112212455161614105{\displaystyle {\begin{array}{crrrcc}{}&{}&{\color {red}1}&{}&{}&{}\\{}&{\rightarrow }&{\color {blue}1}&{\color {red}1}&{}\\{}&{\color {red}2}&{\color {blue}2}&{\color {blue}1}&{\leftarrow }\\{\rightarrow }&{\color {blue}2}&{\color {blue}4}&{\color {blue}5}&{\color {red}5}\\{\color {red}16}&{\color {blue}16}&{\color {blue}14}&{\color {blue}10}&{\color {blue}5}&{\leftarrow }\end{array}}}
El algoritmo de Seidel para T n
  1. Comience colocando 1 en la fila 0 y sea k el número de la fila que se está llenando actualmente.
  2. Si k es impar, entonces coloca el número del extremo izquierdo de la fila k − 1 en la primera posición de la fila k , y llena la fila de izquierda a derecha, de modo que cada entrada sea la suma del número de la izquierda y el número de la derecha.
  3. Al final de la fila, duplique el último número.
  4. Si k es par, proceda de forma similar en la otra dirección.

El algoritmo de Seidel es de hecho mucho más general (véase la exposición de Dominique Dumont [ 44 ] ) y fue redescubierto varias veces posteriormente.

De forma similar al enfoque de Seidel, DE Knuth y TJ Buckholtz dieron una ecuación de recurrencia para los números T 2 n y recomendaron este método para calcular B 2 n y E 2 n 'en computadoras electrónicas usando solo operaciones simples con números enteros'. [ 24 ]

VI Arnold [ 45 ] redescubrió el algoritmo de Seidel y más tarde Millar, Sloane y Young popularizaron el algoritmo de Seidel bajo el nombre de transformada bustrophedon .

Forma triangular:

Solo OEIS : A000657  , con un 1, y OEIS : A214267  , con dos 1, están en el OEIS.

Distribución con un 1 suplementario y un 0 en las siguientes filas:

Este es OEIS : A239005  , una versión firmada de OEIS : A008280  . La diagonal principal es OEIS : A122045  . La diagonal principal es OEIS : A155585  . La columna central es OEIS : A099023  . Sumas de fila: 1, 1, −2, −5, 16, 61.... Ver OEIS : A163747  . Ver la matriz que comienza con 1, 1, 0, −2, 0, 16, 0 a continuación.

El algoritmo de Akiyama-Tanigawa aplicado a OEIS : A046978  ( n + 1 )/ OEIS : A016116  ( n ) produce:

1. La primera columna es OEIS : A122045  . Su transformación binomial da como resultado:

La primera fila de esta matriz es OEIS : A155585 .  Los valores absolutos de las antidiagonales crecientes son OEIS : A008280  . La suma de las antidiagonales es −OEIS : A163747  ( n + 1 ).

2. La segunda columna es 1 1 −1 −5 5 61 −61 −1385 1385... . Su transformación binomial produce:

La primera fila de esta matriz es 1 2 2 −4 −16 32 272 544 −7936 15872 353792 −707584... . Los valores absolutos de la segunda bisección son el doble de los valores absolutos de la primera bisección.

Considere el algoritmo de Akiyama-Tanigawa aplicado a OEIS : A046978  ( n )/( OEIS : A158780  ( n + 1 ) =  abs( OEIS : A117575  ( n )) + 1 = 1 , 2, 2, 3/2 , 1 , 3/4 , 3 / 4 , 7 / 8 , 1, 17 / 16 , 17 / 16 , 33 / 32 ... .

La primera columna cuyos valores absolutos son OEIS : A000111  podría ser el numerador de una función trigonométrica.

OEIS : A163747  es una autosecuencia de primer tipo (la diagonal principal es OEIS : A000004  ). La matriz correspondiente es:

Las dos primeras diagonales superiores son −1 3 −24 402... = (−1) n + 1  × OEIS : A002832 . La suma de las antidiagonales es 0 −2 0 10... = 2 × OEIS : A122045 ( n + 1).       

OEIS : A163982  es una autosecuencia de segundo tipo, como por ejemplo OEIS : A164555  / OEIS : A027642  . Por lo tanto, el array:

La diagonal principal, aquí 2 −2 8 −92... , es el doble de la primera superior, aquí OEIS : A099023  . La suma de las antidiagonales es 2 0 −4 0... = 2  × OEIS : A155585 ( n + 1). OEIS : A163747OEIS : A163982 = 2 × OEIS : A122045 .          

Una perspectiva combinatoria: permutaciones alternas

Alrededor de 1880, tres años después de la publicación del algoritmo de Seidel, Désiré André demostró un resultado ahora clásico del análisis combinatorio. [ 46 ] [ 47 ] Al observar los primeros términos del desarrollo de Taylor de las funciones trigonométricas tan x y sec x , André hizo un descubrimiento.

broncearseincógnita=incógnita+2incógnita33¡+16incógnita55¡+272incógnita77¡+7936incógnita99¡+segundoincógnita=1+incógnita22¡+5incógnita44¡+61incógnita66¡+1385incógnita88¡+50521incógnita1010¡+{\displaystyle {\begin{aligned}\tan x&=x+{\frac {2x^{3}}{3!}}+{\frac {16x^{5}}{5!}}+{\frac {272x^{7}}{7!}}+{\frac {7936x^{9}}{9!}}+\cdots \\[6pt]\sec x&=1+{\frac {x^{2}}{2!}}+{\frac {5x^{4}}{4!}}+{\frac {61x^{6}}{6!}}+{\frac {1385x^{8}}{8!}}+{\frac {50521x^{10}}{10!}}+\cdots \end{aligned}}}

Los coeficientes son los números de Euler en zigzag de índice impar y par, respectivamente. En consecuencia, la expansión ordinaria de tan x + sec x tiene como coeficientes los números racionales S n .

broncearseincógnita+segundoincógnita=1+incógnita+12incógnita2+13incógnita3+524incógnita4+215incógnita5+61720incógnita6+{\displaystyle \tan x+\sec x=1+x+{\tfrac {1}{2}}x^{2}+{\tfrac {1}{3}}x^{3}+{\tfrac {5}{24}}x^{4}+{\tfrac {2}{15}}x^{5}+{\tfrac {61}{720}}x^{6}+\cdots }

Posteriormente, André logró demostrar, mediante un argumento de recurrencia, que las permutaciones alternas de tamaño impar se enumeran mediante los números de Euler de índice impar (también llamados números tangentes) y las permutaciones alternas de tamaño par mediante los números de Euler de índice par (también llamados números secantes).

La media aritmética del primer y segundo número de Bernoulli son los números de Bernoulli asociados: B 0 = 1 , B 1 = 0 , B 2 = 1 / 6 , B 3 = 0 , B 4 = − 1 / 30 , OEIS : A176327  / OEIS : A027642  . A través de la segunda fila de su transformada inversa de Akiyama-Tanigawa OEIS : A177427  , conducen a la serie de Balmer OEIS : A061037  / OEIS : A061038  .

El algoritmo de Akiyama-Tanigawa aplicado a OEIS : A060819  ( n + 4 ) / OEIS : A145979  ( n ) conduce a los números de Bernoulli OEIS : A027641  / OEIS : A027642  , OEIS : A164555  / OEIS : A027642  , o OEIS : A176327  OEIS : A176289  sin B 1 , denominados números de Bernoulli intrínsecos B i ( n ) .

Por lo tanto, otro vínculo entre los números intrínsecos de Bernoulli y la serie de Balmer a través de OEIS : A145979  ( n ).

OEIS : A145979  ( n − 2 ) = 0, 2, 1, 6,... es una permutación de los números no negativos.

Los términos de la primera fila son f(n) = 1/2 + 1/ n + 2. 2 , f ( n ) es una autosecuencia de segundo tipo. 3/2, f(n) conduce por su transformada binomial inversa a 3/2 −1/2 1/3 −1/4 1/5 ... = 1/2 + log 2.

Consideremos g(n) = 1/2 – 1 / (n+2) = 0, 1/6, 1/4, 3/10, 1/3. Las transformadas de Akiyama-Tanagiwa dan como resultado:

0, g(n), es una autosecuencia de segundo tipo.

Euler OEIS : A198631  ( n ) / OEIS : A006519  ( n + 1 ) sin el segundo término ( 1 / 2 ) son los números intrínsecos fraccionarios de Euler E i ( n ) = 1, 0, − 1 / 4 , 0, 1 / 2 , 0, − 17 / 8 , 0, ... La transformada de Akiyama correspondiente es:

La primera línea es Eu ( n ) . Eu ( n ) precedido por un cero es una autosecuencia de primera especie. Está vinculado a los números de Oresme. Los numeradores de la segunda línea son OEIS : A069834  precedido por 0. La tabla de diferencias es:

Propiedades aritméticas de los números de Bernoulli

Los números de Bernoulli pueden expresarse en términos de la función zeta de Riemann como B n = − (1 − n ) para enteros n ≥ 0 siempre que para n = 0 la expresión (1 − n ) se entienda como el valor límite y se utilice la convención B 1 = 1 / 2 . Esto los relaciona íntimamente con los valores de la función zeta en enteros negativos. Como tal, cabría esperar que tuvieran, y de hecho tienen, profundas propiedades aritméticas. Por ejemplo, la conjetura de Agoh-Giuga postula que p es un número primo si y solo si pB p − 1 es congruente con −1 módulo p . Las propiedades de divisibilidad de los números de Bernoulli están relacionadas con los grupos de clases ideales de cuerpos ciclotómicos mediante un teorema de Kummer y su fortalecimiento en el teorema de Herbrand-Ribet , y con los números de clases de cuerpos cuadráticos reales mediante Ankeny-Artin-Chowla .

Los teoremas de Kummer

Los números de Bernoulli están relacionados con el Último Teorema de Fermat (LTF) mediante el teorema de Kummer , [ 48 ] que dice:

Si el primo impar p no divide a ninguno de los numeradores de los números de Bernoulli B 2 , B 4 , ..., B p − 3 entonces x p + y p + z p = 0 no tiene soluciones en enteros distintos de cero.

Los números primos con esta propiedad se denominan primos regulares . Otro resultado clásico de Kummer son las siguientes congruencias . [ 49 ]

Sea p un número primo impar y b un número par tal que p  1 no divide a b . Entonces, para cualquier entero no negativo k
Bk(pag1)+bk(pag1)+bBbb(modpag).{\displaystyle {\frac {B_{k(p-1)+b}}{k(p-1)+b}}\equiv {\frac {B_{b}}{b}}{\pmod {p}}.}

Una generalización de estas congruencias se conoce como continuidad p -ádica.

continuidad p -ádica

Si b , m y n son enteros positivos tales que m y n no son divisibles por p − 1 y mn (mod p b − 1 ( p − 1)) , entonces

(1pagmetro1)Bmetrometro(1pagnorte1)Bnortenorte(modpagb).{\displaystyle (1-p^{m-1}){\frac {B_{m}}{m}}\equiv (1-p^{n-1}){\frac {B_{n}}{n}}{\pmod {p^{b}}}.}

Dado que B n = − (1 − n ) , esto también se puede escribir

(1pag)ζ()(1pagv)ζ(v)(modpagb),{\displaystyle \left(1-p^{-u}\right)\zeta (u)\equiv \left(1-p^{-v}\right)\zeta (v){\pmod {p^{b}}},}

donde u = 1 − m y v = 1 − n , de modo que u y v son no positivos y no congruentes con 1 módulo p − 1 . Esto nos dice que la función zeta de Riemann, con 1 − p s tomado de la fórmula del producto de Euler, es continua en los números p -ádicos sobre enteros negativos impares congruentes módulo p − 1 a un a ≢ 1 mod ( p − 1) particular , y por lo tanto puede extenderse a una función continua ζ p ( s ) para todos los enteros p -ádicosZpag,{\displaystyle \mathbb {Z} _{p},}la función zeta p -ádica .

Las congruencias de Ramanujan

Las siguientes relaciones, debidas a Ramanujan , proporcionan un método para calcular los números de Bernoulli que es más eficiente que el dado por su definición recursiva original:

(metro+3metro)Bmetro={metro+33j=1metro6(metro+3metro6j)Bmetro6j,si metro0(mod6);metro+33j=1metro26(metro+3metro6j)Bmetro6j,si metro2(mod6);metro+36j=1metro46(metro+3metro6j)Bmetro6j,si metro4(mod6).{\displaystyle {\binom {m+3}{m}}B_{m}={\begin{cases}{\frac {m+3}{3}}-\sum \limits _{j=1}^{\frac {m}{6}}{\binom {m+3}{m-6j}}B_{m-6j},&{\text{if }}m\equiv 0{\pmod {6}};\\{\frac {m+3}{3}}-\sum \limits _{j=1}^{\frac {m-2}{6}}{\binom {m+3}{m-6j}}B_{m-6j},&{\text{if }}m\equiv 2{\pmod {6}};\\-{\frac {m+3}{6}}-\sum \limits _{j=1}^{\frac {m-4}{6}}{\binom {m+3}{m-6j}}B_{m-6j},&{\text{if }}m\equiv 4{\pmod {6}}.\end{cases}}}

Teorema de Von Staudt-Clausen

El teorema de von Staudt-Clausen fue propuesto por Karl Georg Christian von Staudt [ 50 ] y Thomas Clausen [ 51 ] de forma independiente en 1840. El teorema establece que para cada n > 0 ,

B2norte+(pag1)2norte1pag{\displaystyle B_{2n}+\sum _{(p-1)\,\mid \,2n}{\frac {1}{p}}}

es un número entero. La suma se extiende sobre todos los primos p para los cuales p − 1 divide a 2 n .

Una consecuencia de esto es que el denominador de B 2 n viene dado por el producto de todos los primos p para los que p − 1 divide a 2 n . En particular, estos denominadores no son cuadrados y son divisibles por  6.

¿Por qué desaparecen los números impares de Bernoulli?

La suma

φk(norte)=i=0norteiknortek2{\displaystyle \varphi _{k}(n)=\sum _{i=0}^{n}i^{k}-{\frac {n^{k}}{2}}}

puede evaluarse para valores negativos del índice n . Al hacerlo, se demostrará que es una función impar para valores pares de k , lo que implica que la suma solo tiene términos de índice impar. Esto y la fórmula para la suma de Bernoulli implican que B 2 k + 1 − m es 0 para m par y 2 k + 1 − m > 1 ; y que el término para B 1 se cancela con la resta. El teorema de von Staudt-Clausen combinado con la representación de Worpitzky también da una respuesta combinatoria a esta pregunta (válida para n > 1).

Del teorema de von Staudt-Clausen se sabe que para n impar > 1 el número 2 B n es un entero. Esto parece trivial si se sabe de antemano que el entero en cuestión es cero. Sin embargo, al aplicar la representación de Worpitzky se obtiene

2Bnorte=metro=0norte(1)metro2metro+1metro¡{norte+1metro+1}=0(norte>1 es extraño){\displaystyle 2B_{n}=\sum _{m=0}^{n}(-1)^{m}{\frac {2}{m+1}}m!\left\{{n+1 \atop m+1}\right\}=0\quad (n>1{\text{ is odd}})}

como una suma de enteros , lo cual no es trivial. Aquí sale a la luz un hecho combinatorio que explica la desaparición de los números de Bernoulli en índices impares. Sea S n , m el número de aplicaciones sobreyectivas de {1, 2, ..., n } a {1, 2, ..., m }, entonces S n , m = m ! { n m } . La última ecuación solo puede cumplirse si

extraño metro=1norte12metro2Snorte,metro=incluso metro=2norte2metro2Snorte,metro(norte>2 es incluso).{\displaystyle \sum _{{\text{odd }}m=1}^{n-1}{\frac {2}{m^{2}}}S_{n,m}=\sum _{{\text{even }}m=2}^{n}{\frac {2}{m^{2}}}S_{n,m}\quad (n>2{\text{ is even}}).}

Esta ecuación se puede demostrar por inducción. Los dos primeros ejemplos de esta ecuación son:

n = 4: 2 + 8 = 7 + 3 ,
n = 6: 2 + 120 + 144 = 31 + 195 + 40 .

Por lo tanto, los números de Bernoulli se anulan en índices impares porque algunas identidades combinatorias no obvias están incorporadas en los números de Bernoulli.

Una reformulación de la hipótesis de Riemann

La conexión entre los números de Bernoulli y la función zeta de Riemann es lo suficientemente fuerte como para proporcionar una formulación alternativa de la hipótesis de Riemann (HR) que utiliza solo los números de Bernoulli. De hecho, Marcel Riesz demostró que la HR es equivalente a la siguiente afirmación: [ 52 ]

Para cada ε > 1 / 4 existe una constante C ε > 0 (que depende de ε ) tal que | R ( x ) | < C ε x ε cuando x → ∞ .

Aquí R ( x ) es la función de Riesz.

R(incógnita)=2k=1kk¯incógnitak(2π)2k(B2k2k)=2k=1kk¯incógnitak(2π)2kβ2k.{\displaystyle R(x)=2\sum _{k=1}^{\infty }{\frac {k^{\overline {k}}x^{k}}{(2\pi )^{2k}\left({\frac {B_{2k}}{2k}}\right)}}=2\sum _{k=1}^{\infty }{\frac {k^{\overline {k}}x^{k}}{(2\pi )^{2k}\beta _{2k}}}.}

n k denota la potencia factorial ascendente en la notación de DE Knuth . Los números β n = B n / n aparecen con frecuencia en el estudio de la función zeta y son significativos porque β n es un p -entero para primos p donde p − 1 no divide a n . Los β n se denominan números de Bernoulli divididos .

Números de Bernoulli generalizados

Los números de Bernoulli generalizados son ciertos números algebraicos , definidos de forma similar a los números de Bernoulli, que están relacionados con valores especiales de las funciones L de Dirichlet del mismo modo que los números de Bernoulli están relacionados con valores especiales de la función zeta de Riemann.

Sea χ un carácter de Dirichlet módulo f . Los números de Bernoulli generalizados asociados a χ se definen por

a=1Fχ(a)tmiatmiFt1=k=0Bk,χtkk¡.{\displaystyle \sum _{a=1}^{f}\chi (a){\frac {te^{at}}{e^{ft}-1}}=\sum _{k=0}^{\infty }B_{k,\chi }{\frac {t^{k}}{k!}}.}

Aparte de la excepcional B 1,1 = 1 / 2 , tenemos, para cualquier carácter de Dirichlet χ , que B k , χ = 0 si χ (−1) ≠ (−1) k .

Generalizando la relación entre los números de Bernoulli y los valores de la función zeta de Riemann en enteros no positivos, se tiene que para todos los enteros k ≥ 1 :

L(1k,χ)=Bk,χk,{\displaystyle L(1-k,\chi )=-{\frac {B_{k,\chi }}{k}},}

donde L ( s , χ ) es la función L de Dirichlet de χ . [ 53 ]

Número de Eisenstein-Kronecker

Los números de Eisenstein-Kronecker son un análogo de los números de Bernoulli generalizados para cuerpos cuadráticos imaginarios . [ 54 ] [ 55 ] Están relacionados con los valores L críticos de los caracteres de Hecke . [ 55 ]

Apéndice

Identidades diversas

  • El cálculo umbral proporciona una forma compacta de la fórmula de Bernoulli, mediante el uso de un símbolo abstracto B ,
    Smetro(norte)=1metro+1((B+norte)metro+1Bmetro+1),{\displaystyle S_{m}(n)={\frac {1}{m+1}}((\mathbf {B} +n)^{m+1}-B_{m+1}),}

    donde el símbolo B k que aparece durante la expansión binomial del término entre paréntesis debe ser reemplazado por el número de Bernoulli B k (y B 1 = + 1 / 2 ). De manera más sugerente y mnemotécnica, esto puede escribirse como

    Smetro(norte)=0norte(B+incógnita)metrodincógnita.{\displaystyle S_{m}(n)=\int _{0}^{n}(\mathbf {B} +x)^{m}\,dx.}

    Otras identidades de Bernoulli se pueden escribir de forma compacta con este símbolo. Por ejemplo,

    (12B)metro=(22metro)Bmetro.{\displaystyle (1-2\mathbf {B} )^{m}=(2-2^{m})B_{m}.}
  • Sea n un número par y no negativo. Entonces
    ζ(norte)=(1)norte21Bnorte(2π)norte2(norte¡).{\displaystyle \zeta (n)={\frac {(-1)^{{\frac {n}{2}}-1}B_{n}(2\pi )^{n}}{2(n!)}}.}
  • El n -ésimo cumulante de la distribución de probabilidad uniforme en el intervalo [−1, 0] es B n / n . 
  • Para n ≥ 1 , B n viene dado por el determinante de una matriz de ( n + 1) × ( n + 1) . [ 56 ] Es decir,
    Bnorte=norte¡|100112¡1001norte¡1(norte1)¡101(norte+1)¡1norte¡12¡0|.{\displaystyle B_{n}=n!{\begin{vmatrix}1&0&\cdots &0&1\\{\frac {1}{2!}}&1&\cdots &0&0\\\vdots &\vdots &\ddots &\vdots &\vdots \\{\frac {1}{n!}}&{\frac {1}{(n-1)!}}&\cdots &1&0\\{\frac {1}{(n+1)!}}&{\frac {1}{n!}}&\cdots &{\frac {1}{2!}}&0\end{vmatrix}}.}
    Entonces, el determinante es σ n (1) , el polinomio de Stirling en x = 1 .
  • El número de Bernoulli par B 2 n viene dado por el determinante de una matriz ( n + 1) × ( n + 1) . [ 56 ] Es decir,
    B2norte=(2norte)¡22norte2|1000113¡100015¡13¡1001(2norte+1)¡1(2norte1)¡1(2norte3)¡13¡0|.{\displaystyle B_{2n}=-{\frac {(2n)!}{2^{2n}-2}}{\begin{vmatrix}1&0&0&\cdots &0&1\\{\frac {1}{3!}}&1&0&\cdots &0&0\\{\frac {1}{5!}}&{\frac {1}{3!}}&1&\cdots &0&0\\\vdots &\vdots &\vdots &\ddots &\vdots &\vdots \\{\frac {1}{(2n+1)!}}&{\frac {1}{(2n-1)!}}&{\frac {1}{(2n-3)!}}&\cdots &{\frac {1}{3!}}&0\end{vmatrix}}.}
  • Sea n ≥ 1. Entonces [ 57 ]
    1nortek=1norte(nortek)BkBnortek+Bnorte1=Bnorte.{\displaystyle {\frac {1}{n}}\sum _{k=1}^{n}{\binom {n}{k}}B_{k}B_{n-k}+B_{n-1}=-B_{n}.}
  • Sea n ≥ 1. Entonces [ 58 ]
    k=0norte(norte+1k)(norte+k+1)Bnorte+k=0.{\displaystyle \sum _{k=0}^{n}{\binom {n+1}{k}}(n+k+1)B_{n+k}=0.}
  • Sea n ≥ 0. Entonces ( Leopold Kronecker 1883)
    Bnorte=k=1norte+1(1)kk(norte+1k)j=1kjnorte.{\displaystyle B_{n}=-\sum _{k=1}^{n+1}{\frac {(-1)^{k}}{k}}{\binom {n+1}{k}}\sum _{j=1}^{k}j^{n}.}
  • Sea n ≥ 1 y m ≥ 1. Entonces [ 59 ]
    (1)metror=0metro(metror)Bnorte+r=(1)nortes=0norte(nortes)Bmetro+s.{\displaystyle (-1)^{m}\sum _{r=0}^{m}{\binom {m}{r}}B_{n+r}=(-1)^{n}\sum _{s=0}^{n}{\binom {n}{s}}B_{m+s}.}
  • Sea n ≥ 4 y H n el n -ésimo número armónico , es decir,
    Hnorte=k=1nortek1.{\displaystyle H_{n}=\sum _{k=1}^{n}k^{-1}.}

    Entonces (H. Miki 1978)

    norte2k=2norte2BnorteknortekBkkk=2norte2(nortek)BnorteknortekBk=HnorteBnorte.{\displaystyle {\frac {n}{2}}\sum _{k=2}^{n-2}{\frac {B_{n-k}}{n-k}}{\frac {B_{k}}{k}}-\sum _{k=2}^{n-2}{\binom {n}{k}}{\frac {B_{n-k}}{n-k}}B_{k}=H_{n}B_{n}.}
  • Sea n ≥ 4 . Entonces ( Yuri Matiyasevich 1997)
    (norte+2)k=2norte2BkBnortek2l=2norte2(norte+2l)BlBnortel=norte(norte+1)Bnorte.{\displaystyle (n+2)\sum _{k=2}^{n-2}B_{k}B_{n-k}-2\sum _{l=2}^{n-2}{\binom {n+2}{l}}B_{l}B_{n-l}=n(n+1)B_{n}.}
  • La identidad de Faber– PandharipandeZagier –Gessel establece que para n ≥ 1 ,
    norte2(Bnorte1(incógnita)+k=1norte1Bk(incógnita)kBnortek(incógnita)nortek)k=0norte1(nortek)BnorteknortekBk(incógnita)=Hnorte1Bnorte(incógnita).{\displaystyle {\frac {n}{2}}\left(B_{n-1}(x)+\sum _{k=1}^{n-1}{\frac {B_{k}(x)}{k}}{\frac {B_{n-k}(x)}{n-k}}\right)-\sum _{k=0}^{n-1}{\binom {n}{k}}{\frac {B_{n-k}}{n-k}}B_{k}(x)=H_{n-1}B_{n}(x).}
    Elegir x = 0 o x = 1 da como resultado la identidad numérica de Bernoulli.
  • Para n ≥ 0 si B 1 = 1 / 2 , o n ≥ 1 si B 1 = − 1 / 2 ,
    k=0norte(nortek)Bknortek+2=Bnorte+1norte+1.{\displaystyle \sum _{k=0}^{n}{\binom {n}{k}}{\frac {B_{k}}{n-k+2}}={\frac {B_{n+1}}{n+1}}.}
  • Sea n ≥ 0. Entonces
    1+k=0norte(nortek)2nortek+1nortek+1Bk(1)=2norte{\displaystyle -1+\sum _{k=0}^{n}{\binom {n}{k}}{\frac {2^{n-k+1}}{n-k+1}}B_{k}(1)=2^{n}}

    y

    1+k=0norte(nortek)2nortek+1nortek+1Bk(0)=δnorte,0.{\displaystyle -1+\sum _{k=0}^{n}{\binom {n}{k}}{\frac {2^{n-k+1}}{n-k+1}}B_{k}(0)=\delta _{n,0}.}
  • Una relación de reciprocidad de M.  B.  Gelfand [ 60 ] establece que
    (1)metro+1j=0k(kj)Bmetro+1+jmetro+1+j+(1)k+1j=0metro(metroj)Bk+1+jk+1+j=k¡metro¡(k+metro+1)¡.{\displaystyle (-1)^{m+1}\sum _{j=0}^{k}{\binom {k}{j}}{\frac {B_{m+1+j}}{m+1+j}}+(-1)^{k+1}\sum _{j=0}^{m}{\binom {m}{j}}{\frac {B_{k+1+j}}{k+1+j}}={\frac {k!m!}{(k+m+1)!}}.}

Véase también

Notas

  1. Traducción del texto: [ 12 ] " ... Así podemos, paso a paso, alcanzar potencias cada vez mayores y, con poco esfuerzo, formar la siguiente tabla: Sumas de potenciasnorte=k=1nortek=12norte2+12norte{\displaystyle \textstyle \int n=\sum _{k=1}^{n}k={\frac {1}{2}}n^{2}+{\frac {1}{2}}n}norte10=k=1nortek10=111norte11+12norte10+56norte91norte7+1norte512norte3+566norte{\displaystyle \textstyle \int n^{10}=\sum _{k=1}^{n}k^{10}={\frac {1}{11}}n^{11}+{\frac {1}{2}}n^{10}+{\frac {5}{6}}n^{9}-1n^{7}+1n^{5}-{\frac {1}{2}}n^{3}+{\frac {5}{66}}n}En efecto, si uno ha examinado diligentemente la ley de progresión aritmética allí, también podrá continuar lo mismo sin estos cálculos tortuosos: Porque sido{\displaystyle \textstyle c}se toma como exponente de cualquier potencia, la suma de todosnortedo{\displaystyle \textstyle n^{c}}se produce o nortedo=k=1nortekdo=1do+1nortedo+1+12nortedo+do2Anortedo1+do(do1)(do2)234Bnortedo3+do(do1)(do2)(do3)(do4)23456donortedo5+do(do1)(do2)(do3)(do4)(do5)(do6)2345678Dnortedo7+{\displaystyle \textstyle \int n^{c}=\sum _{k=1}^{n}k^{c}={\frac {1}{c+1}}n^{c+1}+{\frac {1}{2}}n^{c}+{\frac {c}{2}}An^{c-1}+{\frac {c(c-1)(c-2)}{2\cdot 3\cdot 4}}Bn^{c-3}+{\frac {c(c-1)(c-2)(c-3)(c-4)}{2\cdot 3\cdot 4\cdot 5\cdot 6}}Cn^{c-5}+{\frac {c(c-1)(c-2)(c-3)(c-4)(c-5)(c-6)}{2\cdot 3\cdot 4\cdot 5\cdot 6\cdot 7\cdot 8}}Dn^{c-7}+\cdots }y así sucesivamente, el exponente de su podernorte{\displaystyle n}disminuyendo continuamente en 2 hasta que llega anorte{\displaystyle n}onorte2{\displaystyle n^{2}}Las letras mayúsculasA,B,do,D,{\displaystyle \textstyle A,B,C,D,}etc. denotan en orden los coeficientes de los últimos términos paranorte2,norte4,norte6,norte8{\displaystyle \textstyle \int n^{2},\int n^{4},\int n^{6},\int n^{8}}, etc., es decir A=16,B=130,do=142,D=130{\displaystyle \textstyle A={\frac {1}{6}},B=-{\frac {1}{30}},C={\frac {1}{42}},D=-{\frac {1}{30}}}.
  2. Esta fórmula fue descubierta (o quizás redescubierta) por Giorgio Pietrocola. Su demostración está disponible en italiano ( Pietrocola 2008 ).

Referencias

  1. 1 2 3 Weisstein, Eric W. , "Número de Bernoulli" , MathWorld
  2. 1 2 Selin, Helaine , ed. (1997), Enciclopedia de la historia de la ciencia, la tecnología y la medicina en culturas no occidentales , Springer, pág. 819 (pág. 891), Bibcode : 2008ehst.book.....S , ISBN 0-7923-4066-3
  3. Smith, David Eugene; Mikami, Yoshio (1914), Historia de las matemáticas japonesas , editorial Open Court, pág. 108 ; reimpreso , Dover Publications, 2005, ISBN 9780486434827
  4. Kitagawa, Tomoko L. (23 de julio de 2021), "El origen de los números de Bernoulli: las matemáticas en Basilea y Edo a principios del siglo XVIII", The Mathematical Intelligencer , 44 : 46–56 , doi : 10.1007/s00283-021-10072-y , ISSN 0343-6993 
  5. Menabrea, LF (1842), «Esbozo de la máquina analítica inventada por Charles Babbage, con notas sobre la memoria de la traductora Ada Augusta, condesa de Lovelace» , Bibliothèque Universelle de Genève , 82 , véase la nota G.
  6. Arfken (1970) , pág. 278.
  7. Donald Knuth (2022), Noticias recientes (2022): Matemáticas concretas y Bernoulli .
    Pero el año pasado examiné detenidamente el manifiesto de Bernoulli de Peter Luschny, donde ofrece más de una docena de buenas razones por las que el valor de $B_1$ debería ser en realidad más un medio. Explica que algunos matemáticos de principios del siglo XX habían cambiado unilateralmente las convenciones, porque algunas de sus fórmulas resultaban un poco más elegantes al usar el valor negativo. Fue su decisión bienintencionada, pero en última instancia desacertada, la que condujo a lo que me enseñaron en la década de 1950. […] A estas alturas, lamentablemente, se han escrito cientos de libros que utilizan la convención de "menos un medio". Peor aún, todos los principales sistemas de software para matemáticas simbólicas tienen esa aberración del siglo XX profundamente arraigada. Sin embargo, Luschny me convenció de que todos hemos estado equivocados y que ya es hora de volver a la definición correcta antes de que la situación empeore aún más.
  8. Peter Luschny (2013), El Manifiesto de Bernoulli
  9. 1 2 Knuth (1993) .
  10. Jacobi, CGJ (1834), "De usu legitimo formulas summatoriae Maclaurinianae" , Journal für die reine und angewandte Mathematik , 12 : 263– 272
  11. Knuth (1993) , pág. 14.
  12. Smith, David Eugene (1929). "Jacques (I) Bernoulli: Sobre los 'números de Bernoulli'"" . Un libro de referencia en matemáticas . Nueva York: McGraw-Hill. págs. 85–90 . 
  13. ^ Bernoulli , Jacob (1713). Ars Conjectandi (en latín). Basilea: Impensis Thurnisiorum, Fratrum. págs. 97– 98. doi : 10.5479/sil.262971.39088000323931 . 
  14. Miller, Jeff (23 de junio de 2017), Primeros usos de los símbolos del cálculo , consultado el 3 de enero de 2026
  15. ^ Graham, Knuth y Patashnik (1989) , sección 2.51.
  16. Jordan (1950) pág. 233
  17. Ireland y Rosen (1990) pág. 229
  18. Saalschütz, Louis (1893), Vorlesungen über die Bernoullischen Zahlen, ihren Zusammenhang mit den Secanten-Coficienten und ihre wichtigeren Anwendungen , Berlín: Julius Springer.
  19. Arfken (1970) , pág. 279.
  20. Akiyama, Shigeki; Tanigawa, Yoshio (1 de diciembre de 2001). "Múltiples valores Zeta en números enteros no positivos" (PDF) . El diario Ramanujan . 5 (4): 327– 351. doi : 10.1023/A:1013981102941 . Consultado el 11 de febrero de 2026 .
  21. ^ Kaneko, Masanobu (12 de diciembre de 2000). "El algoritmo de Akiyama-Tanigawa para números de Bernoulli" (PDF) . Diario de secuencias enteras . 3 (1) . Consultado el 11 de febrero de 2026 .
  22. Kawasaki, Naho; Ohno, Yasuo (12 de junio de 2023). "El algoritmo triangular para polinomios de Bernoulli" (PDF) . Enteros . 23. doi : 10.5281/zenodo.8028914 . Recuperado el 11 de febrero de 2026 .
  23. Brent, Richard P.; Harvey, David (2013). "Cálculo rápido de números de Bernoulli, tangentes y secantes" . Matemáticas Computacionales y Analíticas . 50 : 127–142 . arXiv : 1108.0286 . doi : 10.1007/978-1-4614-7621-4_8 . Consultado el 11 de febrero de 2026 .
  24. 1 2 Knuth, DE ; Buckholtz, TJ (1967), "Cálculo de números tangentes, de Euler y de Bernoulli", Mathematics of Computation , 21 (100), American Mathematical Society: 663– 688, doi : 10.2307/2005010 , JSTOR 2005010 
  25. Buhler, J.; Crandall, R.; Ernvall, R.; Metsankyla, T.; Shokrollahi, M. (2001), "Números primos irregulares e invariantes ciclotómicos hasta 12 millones", Journal of Symbolic Computation , 31 ( 1–2 ): 89–96 , doi : 10.1006/jsco.1999.1011
  26. Harvey, David (2010), "Un algoritmo multimodular para calcular números de Bernoulli", Math. Comput. , 79 (272): 2361– 2370, arXiv : 0807.1347 , doi : 10.1090/S0025-5718-2010-02367-1 , S2CID 11329343 , Zbl 1215.11016  
  27. Kellner, Bernd (2002), Program Calcbn – A program for calculate Bernoulli numbers.
  28. Pavlyk, Oleksandr (29 de abril de 2008), "Hoy batimos el récord de Bernoulli: del motor analítico a Mathematica" , Wolfram News.
  29. ^ Graham, Knuth y Patashnik (1989) , 9,67.
  30. ^ Graham, Knuth y Patashnik (1989) , 2,44, 2,52.
  31. Guo, Victor JW; Zeng, Jiang (30 de agosto de 2005), "Un q-análogo de la fórmula de Faulhaber para sumas de potencias", The Electronic Journal of Combinatorics , 11 (2) R19, arXiv : math/0501441 , Bibcode : 2005math......1441G , doi : 10.37236/1876 , S2CID 10467873 
  32. Arfken (1970) , pág. 463.
  33. Comtet, L. (1974). Combinatoria avanzada. El arte de las expansiones finitas e infinitas (Edición revisada y ampliada ). Dordrecht-Boston: D. Reidel Publ. 
  34. 1 2 Rademacher, H. (1973), Teoría analítica de números , Nueva York: Springer-Verlag.
  35. Boole, G. (1880), Tratado del cálculo de diferencias finitas (3.ª ed.), Londres: Macmillan .
  36. Gould, Henry W. (1972), "Fórmulas explícitas para los números de Bernoulli", Amer. Math. Monthly , 79 (1): 44– 51, doi : 10.2307/2978125 , JSTOR 2978125 
  37. Apostol, Tom M. (2010), Introducción a la teoría analítica de números , Springer-Verlag, pág. 197 
  38. Pietrocola, Giorgio (2017), Sobre polinomios para el cálculo de sumas de potencias de enteros sucesivos y números de Bernoulli deducidos del triángulo de Pascal , Academia.edu.
  39. Woon, SC (1997), "Un árbol para generar números de Bernoulli", Math. Mag. , 70 (1): 51– 56, doi : 10.2307/2691054 , JSTOR 2691054 
  40. Stanley, Richard P. (2010), "Un estudio de permutaciones alternas", Combinatoria y grafos , Matemáticas contemporáneas, vol. 531, Providence, RI: American Mathematical Society, pp. 165–196 , arXiv : 0912.4240 , doi : 10.1090/conm/531/10466 , ISBN   978-0-8218-4865-4, MR 2757798 , S2CID 14619581  
  41. Elkies, ND (2003), "Sobre las sumas Sum_(k=-infinito...infinito) (4k+1)^(-n)", Amer. Math. Monthly , 110 (7): 561– 573, arXiv : math.CA/0101168 , doi : 10.2307/3647742 , JSTOR 3647742 
  42. ^ Euler, Leonhard (1735), "De summis serierum reciprocarum", Opera Omnia , I.14, E 41: 73– 86, arXiv : math/0506415 , Bibcode : 2005math......6415E
  43. ^ Seidel, L. (1877), "Über eine einfache Entstehungsweise der Bernoullischen Zahlen und einiger verwandten Reihen", Sitzungsber. Mascar. Akád. , 4 : 157-187
  44. ^ Dumont, D. (1981), "Matrices d'Euler-Seidel" , Séminaire Lotharingien de Combinatoire , B05c
  45. Arnold, VI (1991), "Números ascendentes y descendentes de Bernoulli-Euler asociados con singularidades de funciones, su combinatoria y aritmética", Duke Math. J. , 63 (2): 537– 555, doi : 10.1215/s0012-7094-91-06323-4
  46. ^ André, D. (1879), "Développements de sec x et tan x", CR Acad. Ciencia. , 88 : 965-967
  47. ^ André, D. (1881), "Mémoire sur les permutations alternées", Journal de Mathématiques Pures et Appliquées , 7 : 167– 184
  48. ^ Kummer, EE (1850), "Allgemeiner Beweis des Fermat'schen Satzes, dass die Gleichung x λ + y λ = z λ durch ganze Zahlen unlösbar ist, für alle diejenigen Potenz-Exponenten λ, welche ungerade Primzahlen sind und in den Zählern der ersten (λ-3)/2 Bernoulli'schen Zahlen als Factoren nicht vorkommen" , J. Reine Angew. Matemáticas. , 40 : 131-138
  49. ^ Kummer, EE (1851), "Über eine allgemeine Eigenschaft der racionalen Entwicklungscoficienten einer bestimmten Gattung analytischer Functionen" , J. Reine Angew. Matemáticas. , 1851 ( 41): 368-372
  50. von Staudt, KG cap. (1840), "Beweis eines Lehrsatzes, die Bernoullischen Zahlen betreffend", Journal für die reine und angewandte Mathematik , 21 : 372– 374
  51. ^ Clausen, Thomas (1840), "Lehrsatz aus einer Abhandlung über die Bernoullischen Zahlen", Astron. Nachr. , 17 (22): 351– 352, doi : 10.1002/asna.18400172205
  52. ^ Riesz, M. (1916), "Sur l'hypothèse de Riemann", Acta Mathematica , 40 : 185– 90, doi : 10.1007/BF02418544
  53. Neukirch, Jürgen (1999), Algebraische Zahlentheorie , Grundlehren der mathematischen Wissenschaften , vol. 322, Berlín: Springer-Verlag , ISBN  978-3-540-65399-8, MR 1697859 , Zbl 0956.11021  §VII.2.
  54. Charollois, Pierre; Sczech, Robert (2016), "Funciones elípticas según Eisenstein y Kronecker: una actualización", EMS Newsletter , 2016–9 (101): 8–14 , doi : 10.4171/NEWS/101/4 , ISSN 1027-488X , S2CID 54504376  
  55. ^ Bannai , Kenichi; Kobayashi, Shinichi (2010), "Funciones theta algebraicas y la interpolación p-ádica de números de Eisenstein-Kronecker" , Duke Mathematical Journal , 153 (2), arXiv : math/0610163 , doi : 10.1215/00127094-2010-024 , ISSN 0012-7094 , S2CID 9262012  
  56. 1 2 Malenfant, Jerome (2011), "Expresiones finitas de forma cerrada para la función de partición y para los números de Euler, Bernoulli y Stirling", arXiv : 1103.1585 [ math.NT ]
  57. ^ Euler, E41, Inventio summae cuiusque seriei ex dato termino generali
  58. von Ettingshausen, A. (1827), Vorlesungen über die höhere Mathematik , vol. 1, Viena: Carl Gerold 
  59. Carlitz, L. (1968), "Números de Bernoulli", Fibonacci Quarterly , 6 (3): 71– 85, doi : 10.1080/00150517.1968.12431229
  60. Agoh, Takashi; Dilcher, Karl (2008), "Relaciones de reciprocidad para números de Bernoulli", American Mathematical Monthly , 115 (3): 237–244 , doi : 10.1080/00029890.2008.11920520 , JSTOR 27642447 , S2CID 43614118  

Bibliografía

  • Abramowitz, M.; Stegun, IA (1972), "§23.1: Polinomios de Bernoulli y Euler y la fórmula de Euler-Maclaurin", Manual de funciones matemáticas con fórmulas, gráficos y tablas matemáticas ( 9.ª edición  ), Nueva York: Dover Publications, págs. 804–806 .
  • Arfken, George (1970), Métodos matemáticos para físicos (2.ª  ed.), Academic Press, Bibcode : 1970mmp..book.....A , ISBN 978-0120598519
  • Arlettaz, D. (1998), "Die Bernoulli-Zahlen: eine Beziehung zwischen Topologie und Gruppentheorie", Math. Semestre , 45 : 61– 75, doi : 10.1007/s005910050037 , S2CID 121753654 .
  • Ayoub, A. (1981), "Euler y la función zeta", Amer. Math. Monthly , 74 (2): 1067– 1086, doi : 10.2307/2319041 , JSTOR 2319041 .
  • Conway, John ; Guy, Richard (1996), El libro de los números , Springer-Verlag.
  • Dilcher, K.; Skula, L.; Slavutskii, I. Sh. (1991), "Números de Bernoulli. Bibliografía (1713–1990)" , Queen's Papers in Pure and Applied Mathematics (87), Kingston, Ontario.
  • Dumont, D.; Viennot, G. (1980), "Una interpretación combinatoria de la generación de Seidel de los números de Genocchi", Ann. Discrete Math. , Annals of Discrete Mathematics, vol.  6, pp. 77–87 , doi : 10.1016/S0167-5060(08)70696-4 , ISBN  978-0-444-86048-4.
  • Entringer, RC (1966), "Una interpretación combinatoria de los números de Euler y Bernoulli", Nieuw. Arco. V. Wiskunde , 14 : 241-6.
  • Fee, G.; Plouffe, S. (2007), "Un algoritmo eficiente para el cálculo de números de Bernoulli", arXiv : math/0702300.
  • Graham, R.; Knuth , DE ; Patashnik, O. (1989), Matemáticas concretas (2.ª  ed.), Addison-Wesley, ISBN 0-201-55802-5
  • Ireland, Kenneth; Rosen, Michael (1990), Introducción clásica a la teoría moderna de números (2.ª  ed.), Springer-Verlag, ISBN 0-387-97329-X
  • Jordan, Charles (1950), Cálculo de diferencias finitas , Nueva York: Chelsea Publ. Co..
  • Kaneko, M. (2000), "El algoritmo de Akiyama-Tanigawa para números de Bernoulli" , Journal of Integer Sequences , 12 : 29, Bibcode : 2000JIntS...3...29K.
  • Knuth, DE (1993), "Johann Faulhaber y las sumas de potencias", Matemáticas de la computación , 61 (203), Sociedad Matemática Americana: 277–294 , arXiv : math/9207222 , doi : 10.2307/2152953 , JSTOR 2152953 
  • Kouba, Omran (2016). "Apuntes de clase, polinomios de Bernoulli y aplicaciones". arXiv : 1309.7560v2 [ math.CA ].
  • Luschny, Peter (2007), Una inclusión de los números de Bernoulli.
  • Luschny, Peter (8 de octubre de 2011), "TheLostBernoulliNumbers" , OeisWiki , consultado el 11 de mayo de 2019..
  • El Proyecto de Genealogía Matemática , Fargo: Departamento de Matemáticas, Universidad Estatal de Dakota del Norte, s.f., archivado del original el 10 de mayo de 2019 , recuperado el 11 de mayo de 2019..
  • Milnor, John W.; Stasheff , James D. (1974), "Apéndice B: Números de Bernoulli", Clases características , Annals of Mathematics Studies , vol.  76, Princeton University Press y University of Tokyo Press, pp. 281–287 , ISBN  978-0-691-08122-9, MR 0440554 .
  • Pietrocola, Giorgio (31 de octubre de 2008), "Esplorando un antico sentiero: teoremi sulla somma di potenze di interi Successivi (Corollario 2b)" , Maecla (en italiano) , consultado el 8 de abril de 2017.
  • Slavutskii, Ilya Sh. (1995), "Staudt y las propiedades aritméticas de los números de Bernoulli", Historia Scientiarum , 2 : 69– 74.
  • von Staudt, KG Cap. (1845), "De numeris Bernoullianis, commentationem alteram", Erlangen.
  • Sun, Zhi-Wei (2005–2006), Algunos resultados curiosos sobre los polinomios de Bernoulli y Euler , archivado del original el 31 de octubre de 2001..
  • Woon, SC (1998), "Generalización de una relación entre la función zeta de Riemann y los números de Bernoulli", arXiv : math.NT/9812143.
  • Worpitzky, J. (1883), "Studien über die Bernoullischen und Eulerschen Zahlen" , Journal für die reine und angewandte Mathematik , 94 : 203– 232.
  • "Números de Bernoulli" , Enciclopedia de Matemáticas , EMS Press, 2001 [1994]
  • Los primeros 498 números de Bernoulli del Proyecto Gutenberg.
  • Un algoritmo multimodular para el cálculo de números de Bernoulli
  • La página del número de Bernoulli
  • Programas de numeración de Bernoulli en LiteratePrograms
  • P. Luschny, El cálculo de números primos irregulares
  • P. Luschny, El cálculo y la asintótica de los números de Bernoulli
  • Gottfried Helms, Números de Bernoulli en el contexto de la matriz de Pascal (binomial) (PDF) , archivado (PDF) del original el 9 de octubre de 2022
  • Gottfried Helms, suma de potencias iguales en el contexto de la matriz de Pascal/Bernoulli (PDF) , archivado (PDF) del original el 9 de octubre de 2022
  • Gottfried Helms, Algunas propiedades especiales, sumas de números de Bernoulli y números relacionados (PDF) , archivado (PDF) del original el 9 de octubre de 2022