Mi pregunta es literalmente el título:
Mostrar $y = \sum_{k=1}^\infty \frac{k^{k-1}}{k!} x^k$ satisface $ye^{-y} = x$ .
He aquí una pequeña motivación para la pregunta. Si no te interesa la motivación, ¡ignora lo que sigue!
Considerando los grafos aleatorios de Erdos-Renyi, en el caso supercrítico $p = c/n$ con $c > 1$ y considerar el tamaño del componente gigante. Podemos demostrar (por ejemplo, Grimmett, Probability on Graphs, Theorem 11.1) que el tamaño es $\theta(p) n$ , donde $\theta(p) \equiv \theta(c)$ es la probabilidad de supervivencia para un Poisson $(c)$ proceso de bifurcación ( $p = c/n$ ), y esto satisface $$ 1 - \theta(c) = \frac{1}{c} \sum_{k=1}^\infty \frac{k^{k-1}}{k!} (ce^{-c})^k.$$ Además, también se puede demostrar (por ejemplo, Frieze-Karonski, Introduction to Random Graphs, Theorem 2.14) que $\theta(c) = 1 - c'/c$ , donde $c'$ es el conjugar de $c$ , es decir, para $c>1$ es la solución $c' < 1$ de $c'e^{-c'} = c e^{-c}$ . Combinando estas dos fórmulas, obtenemos (cancelando el $1/c$ factor) $$ c' = \frac{1}{c} \sum_{k=1}^\infty \frac{k^{k-1}}{k!} (ce^{-c})^k.$$ En mi declaración general de la parte superior, acabo de dejar $y = c'$ y $x = ce^{-c}$ para facilitar la anotación.
Sin embargo, por mucho que lo intente, no he podido encontrar ninguna forma de justificar la igualdad, aparte de, por supuesto, probar los dos resultados sobre el tamaño del gigante por separado y luego equipararlos. Es de esperar que haya una forma de demostrarlo directamente -- después de todo, no es exactamente una expresión fea: después de no mucho tiempo se vuelve básicamente geométrica; lo he intentado en Matlab y con $c = 3$ Obtengo la precisión de la máquina después de sólo $40$ términos en la suma. (Mi supervisora de doctorado también le ha echado un vistazo, pero tampoco ha sido capaz de resolverlo).
0 votos
¿Qué ocurre si se intenta calcular la serie de Taylor para $y$ (utilizando la diferenciación implícita) en torno a $x=0$ ? Usando la aproximación de Stirling, parece que tienes una serie de potencia razonablemente convergente alrededor de $x=0$ .
1 votos
Esto es $y=W(-x)$ donde $W$ es el Función Lambert W .
0 votos
@ThomasAndrews sí, lo vi al buscar preguntas relacionadas. No veo ninguna formulación explícita de la función W de Lambert -- de hecho, una persona dijo "usa la función W de Lambert, ya que una solución explícita es muy complicada" (o algo así). Más arriba he dado una fórmula explícita en forma de serie de potencias con unos coeficientes muy bonitos.