Articulo de referencia

Método de Egorychev

El método de Egorychev es un conjunto de técnicas introducidas por Georgy Egorychev para hallar identidades entre sumas de coeficientes binomiales , números de Stirling , número...

El método de Egorychev es un conjunto de técnicas introducidas por Georgy Egorychev para hallar identidades entre sumas de coeficientes binomiales , números de Stirling , números de Bernoulli , números armónicos , números de Catalan y otros números combinatorios. El método se basa en dos observaciones. Primero, muchas identidades pueden demostrarse extrayendo los coeficientes de las funciones generadoras . Segundo, muchas funciones generadoras son series de potencias convergentes , y la extracción de coeficientes puede realizarse utilizando el teorema de los residuos de Cauchy (generalmente integrando sobre un pequeño contorno circular que encierra el origen). La identidad buscada puede hallarse mediante manipulaciones de integrales. Algunas de estas manipulaciones no son claras desde la perspectiva de la función generadora. Por ejemplo, el integrando suele ser una función racional , y la suma de los residuos de una función racional es cero, lo que produce una nueva expresión para la suma original. El residuo en el infinito es particularmente importante en estas consideraciones. Si durante la suma aparece una serie que no es finita, los contornos deben elegirse de manera que la serie converja. Algunas de las integrales empleadas por el método de Egorychev son:

  • Primera integral del coeficiente binomial
(nortek)=rmisz(1+z)nortezk+1=12πi|z|=ρ(1+z)nortezk+1dz{\displaystyle {n \choose k}={\underset {z}{\mathrm {res} }}\;{\frac {(1+z)^{n}}{z^{k+1}}}={\frac {1}{2\pi i}}\int _{|z|=\rho }{\frac {(1+z)^{n}}{z^{k+1}}}\;dz}

dónde0<ρ<{\displaystyle 0<\rho <\infty }

  • Integral del segundo coeficiente binomial
(nortek)=rmisz1(1z)k+1znortek+1=12πi|z|=ρ1(1z)k+1znortek+1dz{\displaystyle {n \choose k}={\underset {z}{\mathrm {res} }}\;{\frac {1}{(1-z)^{k+1}z^{n-k+1}}}={\frac {1}{2\pi i}}\int _{|z|=\rho }{\frac {1}{(1-z)^{k+1}z^{n-k+1}}}\;dz}

dónde0<ρ<1{\displaystyle 0<\rho <1}

nortek=k¡rmiszexp(nortez)zk+1=k¡2πi|z|=ρexp(nortez)zk+1dz:{\displaystyle n^{k}=k!\;{\underset {z}{\mathrm {res} }}\;{\frac {\exp(nz)}{z^{k+1}}}={\frac {k!}{2\pi i}}\int _{|z|=\rho }{\frac {\exp(nz)}{z^{k+1}}}\;dz:}

dónde0<ρ<{\displaystyle 0<\rho <\infty }

[[knorte]]=rmiszzkznorte+111z=12πi|z|=ρzkznorte+111zdz{\displaystyle [[k\leq n]]={\underset {z}{\mathrm {res} }}\;{\frac {z^{k}}{z^{n+1}}}{\frac {1}{1-z}}={\frac {1}{2\pi i}}\int _{|z|=\rho }{\frac {z^{k}}{z^{n+1}}}{\frac {1}{1-z}}\;dz}

dónde0<ρ<1{\displaystyle 0<\rho <1}

[nortek]=norte¡k¡rmisz1znorte+1(registro11z)k=norte¡k¡12πi|z|=ρ1znorte+1(registro11z)kdz{\displaystyle \left[{n \atop k}\right]={\frac {n!}{k!}}\;{\underset {z}{\mathrm {res} }}\;{\frac {1}{z^{n+1}}}\left(\log {\frac {1}{1-z}}\right)^{k}={\frac {n!}{k!}}{\frac {1}{2\pi i}}\int _{|z|=\rho }{\frac {1}{z^{n+1}}}\left(\log {\frac {1}{1-z}}\right)^{k}\;dz}

dónde0<ρ<1{\displaystyle 0<\rho <1}

{nortek}=norte¡k¡rmisz(exp(z)1)kznorte+1=norte¡k¡12πi|z|=ρ(exp(z)1)kznorte+1dz{\displaystyle \left\{{n \atop k}\right\}={\frac {n!}{k!}}\;{\underset {z}{\mathrm {res} }}\;{\frac {(\exp(z)-1)^{k}}{z^{n+1}}}={\frac {n!}{k!}}{\frac {1}{2\pi i}}\int _{|z|=\rho }{\frac {(\exp(z)-1)^{k}}{z^{n+1}}}\;dz}

dónde0<ρ<.{\displaystyle 0<\rho <\infty .}

Ejemplo I

Supongamos que buscamos evaluar

Sj(norte)=k=0norte(1)k(nortek)(norte+kk)(kj){\displaystyle S_{j}(n)=\sum _{k=0}^{n}(-1)^{k}{n \choose k}{n+k \choose k}{k \choose j}}

que se afirma que es  :(1)norte(nortej)(norte+jj).{\displaystyle (-1)^{n}{n \choose j}{n+j \choose j}.}

Introducir  :(norte+kk)=12πi|z|=ε(1+z)norte+kzk+1dz{\displaystyle {n+k \choose k}={\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n+k}}{z^{k+1}}}\;dz}

y  :(kj)=12πi|w|=γ(1+w)kwj+1dw.{\displaystyle {k \choose j}={\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {(1+w)^{k}}{w^{j+1}}}\;dw.}

Esto da como resultado la suma  :

12πi|z|=ε(1+z)nortez12πi|w|=γ1wj+1k=0norte(1)k(nortek)(1+z)k(1+w)kzkdwdz=12πi|z|=ε(1+z)nortez12πi|w|=γ1wj+1(1(1+w)(1+z)z)nortedwdz=12πi|z|=ε(1+z)norteznorte+112πi|w|=γ1wj+1(1wwz)nortedwdz=(1)norte2πi|z|=ε(1+z)norteznorte+112πi|w|=γ1wj+1(1+w+wz)nortedwdz.{\displaystyle {\begin{aligned}&{\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n}}{z}}{\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {1}{w^{j+1}}}\sum _{k=0}^{n}(-1)^{k}{n \choose k}{\frac {(1+z)^{k}(1+w)^{k}}{z^{k}}}\;dw\;dz\\[6pt]={}&{\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n}}{z}}{\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {1}{w^{j+1}}}\left(1-{\frac {(1+w)(1+z)}{z}}\right)^{n}\;dw\;dz\\[6pt]={}&{\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n}}{z^{n+1}}}{\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {1}{w^{j+1}}}(-1-w-wz)^{n}\;dw\;dz\\[6pt]={}&{\frac {(-1)^{n}}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n}}{z^{n+1}}}{\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {1}{w^{j+1}}}(1+w+wz)^{n}\;dw\;dz.\end{alineado}}}

Esto es

(1)norte2πi|z|=ε(1+z)norteznorte+112πi|w|=γ1wj+1q=0norte(norteq)wq(1+z)qdwdz.{\displaystyle {\frac {(-1)^{n}}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n}}{z^{n+1}}}{\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {1}{w^{j+1}}}\sum _{q=0}^{n}{n \choose q}w^{q}(1+z)^{q}\;dw\;dz.}

Extrayendo el residuo enw=0{\displaystyle w=0}obtenemos

(1)norte2πi|z|=ε(1+z)norteznorte+1(nortej)(1+z)jdz=(nortej)(1)norte2πi|z|=ε(1+z)norte+jznorte+1dz=(1)norte(nortej)(norte+jnorte){\displaystyle {\begin{aligned}&{\frac {(-1)^{n}}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n}}{z^{n+1}}}{n \choose j}(1+z)^{j}\;dz\\[6pt]={}&{n \choose j}{\frac {(-1)^{n}}{2\pi i}}\int _{|z|=\varepsilon }{\frac {(1+z)^{n+j}}{z^{n+1}}}\;dz\\[6pt]={}&(-1)^{n}{n \choose j}{n+j \choose n}\end{aligned}}}

probando así la afirmación. No hay problemas de convergencia aquí ya que las sumas involucradas son finitas y connorte+k{\displaystyle n+k}yk{\displaystyle k} Al no ser negativo, podemos elegir cualquier valor finito distinto de cero para ε{\displaystyle \varepsilon }yγ{\displaystyle \gamma }.

Ejemplo II

Supongamos que buscamos evaluark=1nortek(2nortenorte+k).{\displaystyle \sum _{k=1}^{n}k{2n \choose n+k}.}

Introducir

(2nortenorte+k)=12πi|z|=ε1znortek+11(1z)norte+k+1dz.{\displaystyle {2n \choose n+k}={\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {1}{z^{n-k+1}}}{\frac {1}{(1-z)^{n+k+1}}}\;dz.}

Observe que esto es cero cuandok>norte{\displaystyle k>n}para que podamos extenderk{\displaystyle k}hasta el infinito para obtener la suma

12πi|z|=ε1znorte+11(1z)norte+1k1kzk(1z)kdz=12πi|z|=ε1znorte+11(1z)norte+1z/(1z)(1z/(1z))2dz=12πi|z|=ε1znorte1(1z)norte1(12z)2dz.{\displaystyle {\begin{aligned}&{\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {1}{z^{n+1}}}{\frac {1}{(1-z)^{n+1}}}\sum _{k\geq 1}k{\frac {z^{k}}{(1-z)^{k}}}\;dz\\[6pt]={}&{\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {1}{z^{n+1}}}{\frac {1}{(1-z)^{n+1}}}{\frac {z/(1-z)}{(1-z/(1-z))^{2}}}\;dz\\[6pt]={}&{\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {1}{z^{n}}}{\frac {1}{(1-z)^{n}}}{\frac {1}{(1-2z)^{2}}}\;dz.\end{aligned}}}

Ahora ponz(1z)=w{\displaystyle z(1-z)=w}para que (observe que conw=z+{\displaystyle w=z+\cdots }la imagen de|z|=ε{\displaystyle |z|=\varepsilon }conε{\displaystyle \varepsilon }pequeño es otro contorno cerrado parecido a un círculo que hace una vuelta y que ciertamente podemos deformar para obtener otro círculo.|w|=γ{\displaystyle |w|=\gamma })

z=114w2y(12z)2=14w{\displaystyle z={\frac {1-{\sqrt {1-4w}}}{2}}\quad {\text{and}}\quad (1-2z)^{2}=1-4w}

y además

dz=12×12×(4)×(14w)1/2dw=(14w)1/2dw{\displaystyle dz=-{\frac {1}{2}}\times {\frac {1}{2}}\times (-4)\times (1-4w)^{-1/2}\;dw=(1-4w)^{-1/2}\;dw}

para obtener la integral

12πi|w|=γ1wnorte114w(14w)1/2dw=12πi|w|=γ1wnorte1(14w)3/2dw.{\displaystyle {\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {1}{w^{n}}}{\frac {1}{1-4w}}(1-4w)^{-1/2}\;dw={\frac {1}{2\pi i}}\int _{|w|=\gamma }{\frac {1}{w^{n}}}{\frac {1}{(1-4w)^{3/2}}}\;dw.}

Esto se evalúa por inspección como (usar el binomio de Newton )

4norte1(norte1+1/2norte1)=4norte1(norte1/2norte1)=4norte1(norte1)¡q=0norte2(norte1/2q)=2norte1(norte1)¡q=0norte2(2norte2q1)=2norte1(norte1)¡(2norte1)¡2norte1(norte1)¡=norte22norte(2nortenorte)=12norte(2nortenorte).{\displaystyle {\begin{aligned}&4^{n-1}{n-1+1/2 \choose n-1}=4^{n-1}{n-1/2 \choose n-1}={\frac {4^{n-1}}{(n-1)!}}\prod _{q=0}^{n-2}(n-1/2-q)\\={}&{\frac {2^{n-1}}{(n-1)!}}\prod _{q=0}^{n-2}(2n-2q-1)={\frac {2^{n-1}}{(n-1)!}}{\frac {(2n-1)!}{2^{n-1}(n-1)!}}\\[6pt]={}&{\frac {n^{2}}{2n}}{2n \choose n}={\frac {1}{2}}n{2n \choose n}.\end{aligned}}}

Aquí el mapeo de z=0{\displaystyle z=0}aw=0{\displaystyle w=0}determina la elección de la raíz cuadrada . Para las condiciones enϵ{\displaystyle \epsilon } yγ{\displaystyle \gamma } tenemos que para que la serie converja necesitamos |z/(1z)|<1{\displaystyle |z/(1-z)|<1}o ϵ/(1ϵ)<1{\displaystyle \epsilon /(1-\epsilon )<1}o ϵ<1/2.{\displaystyle \epsilon <1/2.} Cuanto más cerca esté el contorno de la imagen de |z|=ϵ{\displaystyle |z|=\epsilon } llega al origen es ϵϵ2{\displaystyle \epsilon -\epsilon ^{2}} por eso elegimos γ<ϵϵ2{\displaystyle \gamma <\epsilon -\epsilon ^{2}} Por ejemplo γ=ϵ2ϵ3.{\displaystyle \gamma =\epsilon ^{2}-\epsilon ^{3}.} Esto también garantiza que γ<1/4{\displaystyle \gamma <1/4}entonces |w|=γ{\displaystyle |w|=\gamma }no intersecta el corte de la rama [1/4,){\displaystyle [1/4,\infty )} (y está contenido en la imagen de |z|=ϵ{\displaystyle |z|=\epsilon }). Por ejemploϵ=1/3{\displaystyle \epsilon =1/3} yγ=2/27{\displaystyle \gamma =2/27}funcionará.

Este ejemplo también se presta a métodos más sencillos, pero se incluyó aquí para demostrar el efecto de sustituir en la variable de integración.

Cálculo mediante series de potencias formales

Podemos usar la regla de cambio de variables 1.8 (5) del texto de Egorychev (página 16) en la integral (recordemos que por el requisito de convergencia los polos enz=1{\displaystyle z=1}yz=1/2{\displaystyle z=1/2}no están dentro del contorno ya queε<1/2{\displaystyle \varepsilon <1/2}):

12πi|z|=ε1znorte1(1z)norte1(12z)2dz=rmisz1znorte1(1z)norte1(12z)2{\displaystyle {\frac {1}{2\pi i}}\int _{|z|=\varepsilon }{\frac {1}{z^{n}}}{\frac {1}{(1-z)^{n}}}{\frac {1}{(1-2z)^{2}}}\;dz={\underset {z}{\mathrm {res} }}{\frac {1}{z^{n}}}{\frac {1}{(1-z)^{n}}}{\frac {1}{(1-2z)^{2}}}}

conA(z)=z(12z)2{\displaystyle A(z)={\frac {z}{(1-2z)^{2}}}}yF(z)=11z.{\displaystyle f(z)={\frac {1}{1-z}}.}Nosotros obtenemosh(z)=z(1z){\displaystyle h(z)=z(1-z)}y encontrar

rmisw1wnorte+1[A(z)F(z)h(z)]|z=gramo(w).{\displaystyle {\underset {w}{\mathrm {res} }}{\frac {1}{w^{n+1}}}\left.\left[{\frac {A(z)}{f(z)h'(z)}}\right]\right|_{z=g(w).}}

congramo{\displaystyle g}lo contrario deh{\displaystyle h}.

Esto se convierte en

rmisw1wnorte+1[z/(12z)2(12z)/(1z)]|z=gramo(w){\displaystyle {\underset {w}{\mathrm {res} }}{\frac {1}{w^{n+1}}}\left.\left[{\frac {z/(1-2z)^{2}}{(1-2z)/(1-z)}}\right]\right|_{z=g(w)}}

o alternativamente

rmisw1wnorte+1[z(1z)(12z)3]|z=gramo(w)=rmisw1wnorte[1(12z)3]|z=gramo(w).{\displaystyle {\underset {w}{\mathrm {res} }}{\frac {1}{w^{n+1}}}\left.\left[{\frac {z(1-z)}{(1-2z)^{3}}}\right]\right|_{z=g(w)}={\underset {w}{\mathrm {res} }}{\frac {1}{w^{n}}}\left.\left[{\frac {1}{(1-2z)^{3}}}\right]\right|_{z=g(w).}}

Observa que(12z)2=14z+4z2=14z(1z)=14w{\displaystyle (1-2z)^{2}=1-4z+4z^{2}=1-4z(1-z)=1-4w} así que esto es

rmisw1wnorte1(14w)3/2{\displaystyle {\underset {w}{\mathrm {res} }}{\frac {1}{w^{n}}}{\frac {1}{(1-4w)^{3/2}}}}

y el resto del cálculo continúa como antes.

  • Hosam Mahmoud, 2022, Historia y ejemplos del método Egorychev
  • Marko Riedel, 2024, Ejemplos computacionales del uso del método de Egorychev para evaluar sumas que involucran tipos de números combinatorios (partes 1 y 2, series de potencias formales y operadores de residuos).
  • Marko Riedel, 2024, Ejemplos computacionales del uso del método de Egorychev para evaluar sumas que involucran tipos de números combinatorios (parte 3, variables complejas)

Referencias

  • Egorychev, GP (1984). Representación integral y el cálculo de sumas combinatorias . American Mathematical Society. ISBN 9780821898093.
  • Riedel, Marko; Mahmoud, Hosam (2023). "Método Egorychev: un tesoro escondido" . La Matemática . 2 (4): 893– 933. doi : 10.1007/s44007-023-00065-y .