Número de Bernoulli - Bernoulli number

Números de Bernoulli B±
n
norte fracción decimal
0 1 +1.000000000
1 ± 1/2 ± 0.500000000
2 1/6 +0.166666666
3 0 +0,000000000
4 -1/30 −0,033333333
5 0 +0,000000000
6 1/42 +0.023809523
7 0 +0,000000000
8 -1/30 −0,033333333
9 0 +0,000000000
10 5/66 +0.075757575
11 0 +0,000000000
12 -691/2730 −0,253113553
13 0 +0,000000000
14 7/6 +1.166666666
15 0 +0,000000000
dieciséis -3617/510 −7,092156862
17 0 +0,000000000
18 43867/798 +54.97117794
19 0 +0,000000000
20 -174611/330 −529.1242424

En matemáticas , los números de Bernoulli B n son una secuencia de números racionales que ocurren con frecuencia en la teoría de números . Los números de Bernoulli aparecen en (y pueden definirse mediante) las expansiones de la serie de Taylor de las funciones tangente e hiperbólica tangente , en la fórmula de Faulhaber para la suma de m -ésimas potencias 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, indicadas aquí por y ; difieren solo para n = 1 , donde y . Para cada n impar > 1 , B n = 0 . Para cada 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 , con y .

Los números de Bernoulli fueron descubiertos casi al mismo tiempo por el matemático suizo Jacob Bernoulli , de quien fueron nombrados, e independientemente por el matemático japonés Seki Takakazu . El descubrimiento de Seki fue publicado póstumamente en 1712 en su obra Katsuyō Sanpō ; 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 . 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 el término n = 1 se ve afectado:

  • B-
    n
    con B-
    1
    = -1/2
    ( OEISA027641 / OEISA027642 ) es la convención de signos prescrita por el NIST y la mayoría de los libros de texto modernos.
  • B+
    n
    con B+
    1
    = +1/2
    ( OEISA164555 / OEISA027642 ) se utiliza a veces en la literatura más antigua.

En las fórmulas siguientes, uno puede cambiar de una convención de signos a otra con la relación , o para un número entero n = 2 o mayor, simplemente ignórelo.

Dado que B n = 0 para todos los n impares > 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 sus raíces en la historia temprana 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 de Katsuyō Sanpō de Seki Takakazu (1712), tabulando 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' reales, solo descripciones dadas enteramente en palabras. Entre los grandes matemáticos de la antigüedad que consideraron este problema se encontraban Pitágoras (c. 572-497 a. C., Grecia), Arquímedes (287-212 a. C., Italia), Aryabhata (n. 476, India), Abu Bakr al-Karaji (m. 1019, Persia) y Abu Ali al-Hasan ibn al-Hasan ibn al-Haytham (965-1039, Iraq).

Durante finales del siglo XVI y principios del XVII, los matemáticos lograron avances significativos. En el oeste, 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) jugaron papeles importantes.

Thomas Harriot parece haber sido el primero en derivar y escribir fórmulas para sumas de potencias usando notación simbólica, pero incluso él calculó solo hasta la suma de las cuartas potencias. Johann Faulhaber dio fórmulas para sumas de poderes hasta el poder 17 en su Academia Algebrae de 1631 , mucho más alto que nadie antes que él, pero no dio una fórmula general.

Blaise Pascal en 1654 probó la identidad de Pascal relacionando las sumas de las p- ésimas potencias de los primeros n enteros positivos para p = 0, 1, 2, ..., k .

El matemático suizo Jakob 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 proporcionan una fórmula uniforme para todas las sumas de potencias.

La alegría que experimentó Bernoulli cuando encontró el patrón necesario para calcular rápida y fácilmente los coeficientes de su fórmula para la suma de las potencias c para cualquier entero positivo c se puede ver en su comentario. El escribio:

"Con la ayuda de esta tabla, me tomó menos de medio cuarto de hora encontrar que las décimas potencias de los primeros 1000 números que se suman dan la suma 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. 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 poderes es la formulación más útil y generalizable hasta la fecha. Los coeficientes en la fórmula de Bernoulli ahora se llaman números de Bernoulli, siguiendo una sugerencia de Abraham de Moivre .

La fórmula de Bernoulli a veces se llama fórmula de Faulhaber en honor a Johann Faulhaber, quien encontró formas notables de calcular la suma de potencias, pero nunca declaró la fórmula de Bernoulli. Según Knuth, Carl Jacobi publicó por primera vez una prueba rigurosa de la fórmula de Faulhaber en 1834. El estudio en profundidad de Knuth de la fórmula de Faulhaber concluye (la notación no estándar en el LHS se explica más adelante):

"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
o
para todas las sumas de poderes. 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 ".

Reconstrucción de "Summae Potestatum"

"Summae Potestatum" de Jakob Bernoulli, 1713

Los números de Bernoulli OEISA164555 (n) / OEISA027642 (n) fueron introducidos por Jakob Bernoulli en el libro Ars Conjectandi publicado póstumamente en 1713 página 97. La fórmula principal puede verse en la segunda mitad del facsímil correspondiente. Los coeficientes constantes denotados por 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 utilizan como símbolos de agrupación. Usando la terminología actual, estas expresiones son potencias factoriales decrecientes c k . La notación factorial k ! como atajo para 1 × 2 × ... × k no se introdujo hasta 100 años después. El símbolo integral en el lado izquierdo se remonta a Gottfried Wilhelm Leibniz en 1675, quien lo usó como una letra S larga para "summa" (suma). La letra n en el lado izquierdo no es un índice de suma, pero da el límite superior del rango de suma que debe entenderse como 1, 2, ..., n . Poniendo las cosas juntas, para c positivo , hoy en día es probable que un matemático escriba la fórmula de Bernoulli como:

Esta fórmula sugiere establecer B 1 =1/2al cambiar de la denominada enumeración 'arcaica' que usa solo los índices pares 2, 4, 6 ... a la forma moderna (más sobre diferentes convenciones en el siguiente párrafo). Lo más sorprendente en este contexto es el hecho de que el factorial descendente c k −1 tiene para k = 0 el valor1/c + 1. Así, la fórmula de Bernoulli se puede escribir

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

La fórmula para en la primera mitad contiene un error en el último término; debería ser en lugar de .

Definiciones

Se han encontrado muchas caracterizaciones de los números de Bernoulli en los últimos 300 años, y cada una podría usarse para introducir estos números. Aquí solo se mencionan tres de los más útiles:

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

Para la prueba de la equivalencia de los tres enfoques.

Definición recursiva

Los números de Bernoulli obedecen a las fórmulas de suma

donde y δ denota el delta de Kronecker . Resolver para da las fórmulas recursivas

Definición explícita

En 1893 Louis Saalschütz enumeró un total de 38 fórmulas explícitas para los números de Bernoulli, por lo general dando alguna referencia en la literatura más antigua. Uno de ellos es:

Función generadora

Las funciones generadoras exponenciales son

donde está la sustitución .

La función generadora (ordinaria)

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

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

Los números de Bernoulli dados por la función zeta de Riemann.

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

B+
n
= - (1 - n )
          para n ≥ 1  .

Aquí el argumento de la función zeta es 0 o negativo.

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

para n ≥ 1  .

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

Luego se sigue de ζ → 1 ( n → ∞ ) y la fórmula de Stirling que

para n → ∞  .

Cálculo eficiente de números de Bernoulli

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 número primo; por ejemplo, para probar si la conjetura de Vandiver es válida para p , o incluso simplemente para determinar si p es un número primo irregular . No es factible realizar dicho 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 que requieren solo operaciones O ( p (log p ) 2 ) (vea la notación O grande ).

David Harvey describe un algoritmo para calcular 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 ( n 2 log ( n ) 2 + ε ) 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 8 . La implementación de Harvey se ha incluido en SageMath desde la versión 3.1. Antes de eso, Bernd Kellner calculó B n con precisión total para n = 10 6 en diciembre de 2002 y Oleksandr Pavlyk para n = 10 7 con Mathematica en abril de 2008.

Computadora Año norte Dígitos *
J. Bernoulli ~ 1689 10 1
L. Euler 1748 30 8
JC Adams 1878 62 36
DE Knuth, TJ Buckholtz 1967 1 672 3 330
G. Fee, S. Plouffe 1996 10 000 27 677
G. Fee, S. Plouffe 1996 100 000 376 755
BC Kellner 2002 1 000 000 4 767 529
O. Pavlyk 2008 10 000 000 57 675 260
D. Harvey 2008 100 000 000 676 752 569
* Dígitos 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 diferenciable con suficiente frecuencia, la fórmula de Euler-Maclaurin se puede escribir como

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

Aquí (es decir, la derivada de orden cero de is just ). Además, denotemos una antiderivada de . Por el teorema fundamental del cálculo ,

Por lo tanto, la última fórmula se puede simplificar aún más a la siguiente forma sucinta de la fórmula de Euler-Maclaurin

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

Aquí s k denota la potencia factorial ascendente .

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 ψ .

Suma de poderes

Números de Bernoulli ocupan un lugar destacado en la forma cerrada expresión de la suma de la m º potencias de los primeros n enteros positivos. Para m , n ≥ 0 definir

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 :

donde (m + 1
k
)
denota elcoeficiente binomial.

Por ejemplo, tomar m como 1 da los números triangulares 0, 1, 3, 6, ... OEISA000217 .

Tomar m como 2 da los números piramidales cuadrados 0, 1, 5, 14, ... OEISA000330 .

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

La fórmula de Bernoulli a veces se llama 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 .

Serie de taylor

Los números de Bernoulli aparecen en la expansión de la serie de Taylor de muchas funciones trigonométricas y funciones hiperbólicas .

Tangente
Cotangente
Tangente hiperbólica
Cotangente hiperbólica

Serie Laurent

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

Función Digamma :

Uso en topología

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

El teorema de la firma de Hirzebruch para el género L de una variedad cerrada de orientación suave 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 ejemplo de un principio combinatorio fundamental, el principio de inclusión-exclusión .

Conexión con números Worpitzky

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

También se pueden expresar mediante los números de Stirling del segundo tipo.

Luego 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/21/3...

B 0 = 1
B 1 = 1 -1/2
B 2 = 1 -3/2 + 2/3
B 3 = 1 -7/2 + 12/3 - 6/4
B 4 = 1 -15/2 + 50/3 - 60/4 + 24/5
B 5 = 1 -31/2 + 180/3 - 390/4 + 360/5 - 120/6
B 6 = 1 -63/2 + 602/3 - 2100/4 + 3360/5 - 2520/6 + 720/7

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

Considere la secuencia s n , n ≥ 0 . De los números de Worpitzky OEISA028246 , OEISA163626 aplicado a s 0 , s 0 , s 1 , s 0 , s 1 , s 2 , s 0 , s 1 , s 2 , s 3 , ... es idéntico al Akiyama –Transformada de Tanigawa aplicada a s n (ver Conexión con números de Stirling del primer tipo ). Esto se puede ver a través de la tabla:

Identidad de
la representación de Worpitzky y la transformación de Akiyama-Tanigawa
1 0 1 0 0 1 0 0 0 1 0 0 0 0 1
1 −1 0 2 −2 0 0 3 −3 0 0 0 4 −4
1 −3 2 0 4 −10 6 0 0 9 −21 12
1 −7 12 −6 0 8 −38 54 −24
1 −15 50 −60 24

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

Por lo tanto, para los segundos números de Euler fraccionarios OEISA198631 ( n ) / OEISA006519 ( n + 1 ):

E 0 = 1
E 1 = 1 -1/2
E 2 = 1 -3/2 + 2/4
E 3 = 1 -7/2 + 12/4 - 6/8
E 4 = 1 -15/2 + 50/4 - 60/8 + 24/dieciséis
E 5 = 1 -31/2 + 180/4 - 390/8 + 360/dieciséis - 120/32
E 6 = 1 -63/2 + 602/4 - 2100/8 + 3360/dieciséis - 2520/32 + 720/64

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

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

OEISA164555 ( n + 1 ) / OEISA027642 ( n + 1 ) =n + 1/2 n + 2 - 2× OEISA198631 ( n ) / OEISA006519 ( n + 1 )

que vincula los segundos números de Bernoulli con los segundos números de Euler fraccionarios. 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 OEISA111701 (consulte Conexión con números de Stirling del primer tipo ).

Conexión con números de Stirling del segundo tipo

Si S ( k , m ) denota números de Stirling del segundo tipo, entonces uno tiene:

donde j m denota el factorial descendente .

Si se definen los polinomios de Bernoulli B k ( j ) como:

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

Luego, después de la siguiente propiedad del coeficiente binomial :

uno tiene,

Uno también tiene lo siguiente para los polinomios de Bernoulli,

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

Comparando el coeficiente de j en las dos expresiones de polinomios de Bernoulli, uno tiene:

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

Conexión con números de Stirling del primer tipo

Las dos fórmulas principales que relacionan los números de Stirling sin signo del primer tipo [n
m
]
a los números de Bernoulli (con B 1 = +1/2) están

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

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

Número de Akiyama – Tanigawa
metro
norte
0 1 2 3 4
0 1 1/2 1/3 1/4 1/5
1 1/2 1/3 1/4 1/5 ...
2 1/6 1/6 3/20 ... ...
3 0 1/30 ... ... ...
4 -1/30 ... ... ... ...

Los números de Akiyama-Tanigawa satisfacen una relación de recurrencia simple que se puede aprovechar para calcular iterativamente los números de Bernoulli. Esto conduce al algoritmo que se muestra en la sección 'descripción algorítmica' anterior. Consulte OEISA051714 / OEISA051715 .

Una autosecuencia es una secuencia que tiene su transformada binomial inversa igual a la secuencia con signo. Si la diagonal principal es ceros = OEISA000004 , la autosecuencia es del primer tipo. Ejemplo: OEISA000045 , los números de Fibonacci. Si la diagonal principal es la primera diagonal superior multiplicada por 2, es del segundo tipo. Ejemplo: OEISA164555 / OEISA027642 , los segundos números de Bernoulli (ver OEISA190339 ). La transformada Akiyama – Tanigawa aplicada a 2 - n = 1 / OEISA000079 conduce a OEISA198631 ( n ) / OEISA06519 ( n + 1). Por eso:

Transformada de Akiyama-Tanigawa para los segundos números de Euler
metro
norte
0 1 2 3 4
0 1 1/2 1/4 1/8 1/dieciséis
1 1/2 1/2 3/8 1/4 ...
2 0 1/4 3/8 ... ...
3 -1/4 -1/4 ... ... ...
4 0 ... ... ... ...

Consulte OEISA209308 y OEISA227577 . OEISA198631 ( n ) / OEISA006519 ( n + 1 ) son los segundos números de Euler (fraccionarios) y una autosecuencia del segundo tipo.

(OEISA164555 ( n + 2 )/OEISA027642 ( n + 2 ) = 1/6, 0, -1/30, 0, 1/42, ... ) × (2 n + 3 - 2/n + 2= 3,14/3, 15/2, 62/5, 21, ... ) =OEISA198631 ( n + 1 )/OEISA006519 ( n + 2 ) = 1/2, 0, -1/4, 0, 1/2, ... .

También valioso para OEISA027641 / OEISA027642 (consulte Conexión con números 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

donde es el determinante de una matriz de Hessenberg n-por-n parte del triángulo de Pascal cuyos elementos son:

Ejemplo:

Conexión con números eulerianos

Hay fórmulas que conectan números eulerianos n
m
⟩ A
números de Bernoulli:

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

Una representación de árbol binario

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

SCWoonTree.png

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 con ± denota el signo de un 1 .

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

Restringido a los nodos N de un nivel de árbol fijo n la suma de1/¡ N !es σ n (1) , entonces

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

La integral

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 yo . 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ó

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

Los números de Euler son una secuencia de números enteros íntimamente conectados con los números de Bernoulli. La comparación de las expansiones asintóticas de los números de Bernoulli y Euler muestra que los números de Euler E 2 n tienen una magnitud aproximada2/π(4 2 n - 2 2 n ) veces mayor que los números de Bernoulli B 2 n . En consecuencia:

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 el que n es par.

Estas fórmulas de conversión expresan una conexión entre los números de Bernoulli y Euler. Pero lo que es más importante, existe una raíz aritmética profunda común a ambos tipos de números, que se puede expresar 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

y S 1 = 1 por convención. La magia de estos números radica en el hecho de que resultan ser números racionales. Esto fue probado por primera vez por Leonhard Euler en un artículo histórico De summis serierum reciprocarum (Sobre las sumas de series de recíprocos) y desde entonces ha fascinado a los matemáticos. Los primeros de estos números son

( OEISA099612 / OEISA099617 )

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

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

La expresión [ n par] tiene el valor 1 si n es par y 0 en caso contrario ( corchete 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 + 1cuando n es par. Los R n son aproximaciones racionales de π y dos términos sucesivos siempre encierran el verdadero valor de π . A partir de n = 1 comienza la secuencia ( OEISA132049 / OEISA132050 ):

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

Considere la transformada de Akiyama-Tanigawa para la secuencia OEISA046978 ( n + 2 ) / OEISA016116 ( n + 1 ):

0 1 1/2 0 -1/4 -1/4 -1/8 0
1 1/2 1 3/4 0 -5/8 -3/4
2 -1/2 1/2 9/4 5/2 5/8
3 −1 -7/2 -3/4 15/2
4 5/2 -11/2 -99/4
5 8 77/2
6 -61/2

A partir de la segunda, los numeradores de la primera columna son los denominadores de la fórmula de Euler. La primera columna es:1/2× OEISA163982 .

Una visión algorítmica: el triángulo de Seidel

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

( OEISA000111 ). Ver ( OEISA253671 ).

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

Estas identidades facilitan el cálculo de los números de Bernoulli y Euler: los números de Euler E n vienen dados inmediatamente por T 2 n + 1 y los números de Bernoulli B 2 n se obtienen a partir de T 2 n mediante algunos cambios fáciles, evitando la aritmética racional.

Lo que queda es encontrar una forma 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 .

Algoritmo de Seidel para T n
  1. Comience poniendo 1 en la fila 0 y deje que k denote el número de la fila que se está llenando actualmente
  2. Si k es impar, coloque el número en el extremo izquierdo de la fila k - 1 en la primera posición de la fila k , y complete la fila de izquierda a derecha, siendo cada entrada la suma del número al a la izquierda y el número en la parte superior
  3. Al final de la fila, duplique el último número.
  4. Si k es par, proceda de manera similar en la otra dirección.

El algoritmo de Seidel es de hecho mucho más general (ver la exposición de Dominique Dumont) y fue redescubierto varias veces a partir de entonces.

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'.

VI Arnold redescubrió el algoritmo de Seidel y más tarde Millar, Sloane y Young popularizaron el algoritmo de Seidel con el nombre de transformada boustrophedon .

Forma triangular:

1
1 1
2 2 1
2 4 5 5
dieciséis dieciséis 14 10 5
dieciséis 32 46 56 61 61
272 272 256 224 178 122 61

Solo OEISA000657 , con un 1, y OEISA214267 , con dos 1, están en la OEIS.

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

1
0 1
−1 −1 0
0 −1 −2 −2
5 5 4 2 0
0 5 10 14 dieciséis dieciséis
−61 −61 −56 −46 −32 −16 0

Esta es OEISA239005 , una versión firmada de OEISA008280 . El andiagonal principal es OEISA122045 . La diagonal principal es OEISA155585 . La columna central es OEISA099023 . Sumas de filas: 1, 1, −2, −5, 16, 61 .... Ver OEISA163747 . Vea la matriz que comienza con 1, 1, 0, −2, 0, 16, 0 a continuación.

El algoritmo Akiyama – Tanigawa aplicado a OEISA046978 ( n + 1 ) / OEISA016116 ( n ) produce:

1 1 1/2 0 -1/4 -1/4 -1/8
0 1 3/2 1 0 -3/4
−1 −1 3/2 4 15/4
0 −5 -15/2 1
5 5 -51/2
0 61
−61

1. La primera columna es OEISA122045 . Su transformada binomial conduce a:

1 1 0 −2 0 dieciséis 0
0 −1 −2 2 dieciséis −16
−1 −1 4 14 −32
0 5 10 −46
5 5 −56
0 −61
−61

La primera fila de esta matriz es OEISA155585 . Los valores absolutos de los antidiagonales crecientes son OEISA008280 . La suma de los antidiagonales es - OEISA163747 ( n + 1 ).

2. La segunda columna es 1 1 -1 -5 5 61 -61 -1385 1385 ... . Su transformada binomial produce:

1 2 2 −4 −16 32 272
1 0 −6 −12 48 240
−1 −6 −6 60 192
−5 0 66 32
5 66 66
61 0
−61

La primera fila de esta matriz es 1 2 2 -4 -16 32 272 544 -7936 15 872 353 792 -707 584 ... . Los valores absolutos de la segunda bisección son el doble de los valores absolutos de la primera bisección.

Considere el algoritmo Akiyama-Tanigawa aplicado a OEISA046978 ( n ) / ( OEISA158780 ( n + 1 ) = abs ( OEISA117575 ( n )) + 1 = 1, 2, 2,3/2, 1, 3/4, 3/4, 7/8, 1, 17/dieciséis, 17/dieciséis, 33/32... .

1 2 2 3/2 1 3/4 3/4
−1 0 3/2 2 5/4 0
−1 −3 -3/2 3 25/4
2 −3 -27/2 −13
5 21 -3/2
−16 45
−61

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

OEISA163747 es una secuencia automática del primer tipo (la diagonal principal es OEISA000004 ). La matriz correspondiente es:

0 −1 −1 2 5 −16 −61
−1 0 3 3 −21 −45
1 3 0 −24 −24
2 −3 −24 0
−5 −21 24
−16 45
−61

Las dos primeras diagonales superiores son −1 3 −24 402 ... = (−1) n + 1  ×  OEISA002832 . La suma de los antidiagonales es 0 −2 0 10 ... = 2 ×  OEISA122045 ( n  + 1).

- OEISA163982 es una secuencia automática del segundo tipo, como por ejemplo OEISA164555 / OEISA027642 . De ahí la matriz:

2 1 −1 −2 5 dieciséis −61
−1 −2 −1 7 11 −77
−1 1 8 4 −88
2 7 −4 −92
5 −11 −88
−16 −77
−61

La diagonal principal, aquí 2 −2 8 −92 ... , es el doble de la primera superior, aquí OEISA099023 . La suma de los antidiagonales es 2 0 −4 0 ... = 2 ×  OEISA155585 ( n + 1). OEISA163747  -  OEISA163982 = 2 ×  OEISA122045 .

Una visión combinatoria: permutaciones alternas

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

Los coeficientes son los números de Euler de índices pares e impares, respectivamente. En consecuencia, la expansión ordinaria de tan x + sec x tiene como coeficientes los números racionales S n .

Luego, André tuvo éxito mediante un argumento de recurrencia para mostrar que las permutaciones alternas de tamaño impar se enumeran por los números de Euler de índice impar (también llamados números de tangente) y las permutaciones alternas de tamaño par por los números de Euler de índice par (también llamado números secantes).

Secuencias relacionadas

La media aritmética del primero y segundo números 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, OEISA176327 / OEISA027642 . A través de la segunda fila de su transformación inversa Akiyama – Tanigawa OEISA177427 , conducen a la serie Balmer OEISA061037 / OEISA061038 .

El algoritmo Akiyama-Tanigawa aplica a OEISA060819 ( n + 4 ) / OEISA145979 ( n ) conduce a los números de Bernoulli OEISA027641 / OEISA027642 , OEISA164555 / OEISA027642 , o OEISA176327 OEISA176289 sin B 1 , denominados números intrínsecos de Bernoulli B i ( n ) .

1 5/6 3/4 7/10 2/3
1/6 1/6 3/20 2/15 5/42
0 1/30 1/20 2/35 5/84
-1/30 -1/30 -3/140 -1/105 0
0 -1/42 -1/28 -4/105 -1/28

De ahí otro vínculo entre los números intrínsecos de Bernoulli y la serie de Balmer a través de OEISA145979 ( n ).

OEISA145979 ( 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 del 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.

Considere g (n) = 1/2 - 1 / (n + 2) = 0, 1/6, 1/4, 3/10, 1/3. Las transformaciones de Akiyama-Tanagiwa dan:

0 1/6 1/4 3/10 1/3 5/14 ...
-1/6 -1/6 -3/20 -2/15 -5/42 -3/28 ...
0 -1/30 -1/20 -2/35 -5/84 -5/84 ...
1/30 1/30 3/140 1/105 0 -1/140 ...

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

Euler OEISA198631 ( n ) / OEISA006519 ( n + 1 ) sin el segundo término (1/2) son los números de Euler intrínsecos fraccionarios E i ( n ) = 1, 0, -1/4, 0, 1/2, 0, -17/8, 0, ... La transformación de Akiyama correspondiente es:

1 1 7/8 3/4 21/32
0 1/4 3/8 3/8 5/dieciséis
-1/4 -1/4 0 1/4 25/64
0 -1/2 -3/4 -9/dieciséis -5/32
1/2 1/2 -9/dieciséis -13/8 -125/64

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

0 1 1 7/8 3/4 21/32 19/32
1 0 -1/8 -1/8 -3/32 -1/dieciséis -5/128
−1 -1/8 0 1/32 1/32 3/128 1/64

Propiedades aritméticas de los números de Bernoulli

Los números de Bernoulli se pueden expresar en términos de la función zeta de Riemann como B n = - (1 - n ) para números enteros n ≥ 0 siempre que para n = 0 la expresión - (1 - n ) se entienda como el valor límite y la convención B 1 =1/2se utiliza. Esto los relaciona íntimamente con los valores de la función zeta en enteros negativos. Como tales, se podría esperar que tuvieran y tengan propiedades aritméticas profundas. 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 campos ciclotómicos mediante un teorema de Kummer y su fortalecimiento en el teorema de Herbrand-Ribet , y con los números de clase de campos cuadráticos reales de Ankeny-Artin-Chowla .

Los teoremas de Kummer

Los números de Bernoulli están relacionados con el último teorema de Fermat (FLT) por el teorema de Kummer , que dice:

Si el primo impar p no divide 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 números 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 .

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

Una generalización de estas congruencias se conoce con el nombre de continuidad p -ádica.

p -continuidad ádica

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

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

donde u = 1 - m y v = 1 - n , de manera que u y v son no positiva y no congruente con 1 módulo p - 1 . Esto nos dice que la función zeta de Riemann, con 1 - p - s sacado de la fórmula del producto de Euler, es continua en los números p -ádicos en números enteros negativos impares congruente módulo p - 1 a un particular a ≢ 1 mod ( p - 1) , y así se puede extender a una función continua ζ p ( s ) para todas las p números enteros -adic la p función zeta -adic .

Congruencias de Ramanujan

Las siguientes relaciones, debido 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:

Teorema de Von Staudt-Clausen

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

es un entero. La suma se extiende a todos los primos p para los cuales p - 1 divide 2 n .

Una consecuencia de esto es que el denominador de B 2 n está dado por el producto de todos los primos p para los cuales p - 1 divide 2 n . En particular, estos denominadores son cuadrados libres y divisibles por 6.

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

La suma

se puede evaluar para valores negativos del índice n . Hacerlo mostrará 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 es cancelado por 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).

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

como una suma de números enteros , lo cual no es trivial. Aquí surge un hecho combinatorio que explica la desaparición de los números de Bernoulli en un índice impar. Sea S n , m el número de mapas sobreyectivos de {1, 2, ..., n } a {1, 2, ..., m }, entonces S n , m = m ! {n
m
}
. La última ecuación solo puede sostenerse si

Esta ecuación se puede probar 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 desaparecen con un índice impar 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 (RH) que usa solo el número de Bernoulli. De hecho, Marcel Riesz demostró que la RH es equivalente a la siguiente afirmación:

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

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

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

Números de Bernoulli generalizados

Los números de Bernoulli generalizadas son ciertos números algebraicos , que se define de manera similar a los números de Bernoulli, que están relacionados con los valores especiales de Dirichlet L -Funciones de la misma manera que los números de Bernoulli están relacionados con los valores especiales de la función zeta de Riemann.

Sea χ un carácter de Dirichlet módulo f . Los números de Bernoulli generalizados adjuntos a χ están definidos por

Aparte del 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 para todos los enteros k ≥ 1 :

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

Apéndice

Identidades variadas

  • El cálculo de umbrales da una forma compacta de la fórmula de Bernoulli usando un símbolo abstracto B :

    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 una integral definida:

    Muchas otras identidades de Bernoulli se pueden escribir de forma compacta con este símbolo, p. Ej.

  • Sea n no negativo e incluso
  • El n- ésimo acumulativo de la distribución de probabilidad uniforme en el intervalo [−1, 0] esB n/norte.
  • Deja n ? =1/n !y n ≥ 1 . Entonces B n es el siguiente determinante ( n + 1) × ( n + 1) :
    Por tanto, el determinante es σ n (1) , el polinomio de Stirling en x = 1 .
  • Para los números pares de Bernoulli, B 2 p viene dado por el determinante ( p + 1) × ( p + 1) :
  • Sea n ≥ 1 . Entonces ( Leonhard Euler )
  • Sea n ≥ 1 . Luego
  • Sea n ≥ 0 . Entonces ( Leopold Kronecker 1883)
  • Sean n ≥ 1 y m ≥ 1 . Luego
  • Sea n ≥ 4 y
    el número armónico . Entonces (H. Miki 1978)
  • Sea n ≥ 4 . Yuri Matiyasevich encontrado (1997)
  • Faber– Pandharipande - Zagier –Gessel identidad : para n ≥ 1 ,
    La elección de x = 0 o x = 1 da como resultado la identidad número Bernoulli en una u otra convención.
  • La siguiente fórmula es verdadera para n ≥ 0 si B 1 = B 1 (1) =1/2, pero solo para n ≥ 1 si B 1 = B 1 (0) = -1/2.
  • Sea n ≥ 0 . Luego
    y
  • Una relación de reciprocidad de M. B. Gelfand:

Ver también

Notas

Referencias

Notas al pie

enlaces externos