Mi profesor y yo derivamos el $n$ número de Bernoulli como la recursión $$B_n=\frac{-1}{n+1}\sum_{k=0}^{n-1}\binom{n+1}{k}B_k.$$ Más tarde, en un artículo que estaba leyendo, hay una identidad similar atribuida a Euler, donde $$-(n+1)B_n=\sum_{k=2}^{n-2}\binom{n}{k}B_kB_{n-k}$$
Ciertamente, reindexar la primera ecuación ayudaría a conseguir $$-(n-1)B_n=\sum_{k=0}^{n-2}\binom{n}{k}B_k$$
Pero ahora el problema radica en la multiplicación de cada sumando por el $(n-k)$ número de Bernoulli y esto, junto con el índice que comienza en 2 da la identidad. No consigo ver las matemáticas en esto. El documento que estoy leyendo dice que $b(x)=\frac{x}{e^x-1}=\sum_{k=0}^\infty\frac{B_nx^n}{n!}$ da el resultado, pero no veo cómo. ¿Hay alguna prueba en línea que pueda leer para ayudar a aclarar?
EDIT: El documento que estoy leyendo en realidad afirma que la identidad de Euler es una consecuencia de $$b(x)^2=(1-x)b(x)-xb'(x)$$ Aquí es donde realmente veo la confusión. He demostrado la primera identidad de forma independiente utilizando las ideas proporcionadas en la respuesta de abajo, pero al no estar familiarizado con las series de Laurent estoy teniendo problemas. Espero que, con la nueva identidad de abajo que alguien puede ofrecer un poco más de conocimiento. Gracias.