La respuesta es casi seguro que la densidad de la limitación no existe. Sin entrar en los detalles de la prueba que me permiten dar una heurística argumento que se basa en cómo el OP probable que generó su gráfico de la frecuencia relativa de feliz números.
Deje $Y_n$ ser la r.v. distribuidos de manera uniforme entre los números enteros en el intervalo de $[0,10^n -1]$ (que es $Y_n$ elige una al azar $n$-dígito entero). Si $X_i$ denota la r.v. para el dígito de $10^i$ en $Y_n$,, a continuación,$s(Y_n) = \sum\limits_{i=0}^{n-1} s(X_i)$.
Supongo que la forma en que se genera el gráfico era la primera vez calculada la distribución de $s(Y_n)$ (esto se puede hacer de forma recursiva), a continuación, calculadas $\mathbb{P}\big(s(Y_n) \text{ is happy}\big)$. Esto daría a la densidad relativa de feliz números de entre todos los $n$ dígitos enteros.
El estudio de la distribución de $s(Y_n)$ nos puede decir mucho. Su equivalente a rodar $n$ veces un 10 colindado mueren con caras $0,1,4,\dots, 81$ y la búsqueda de la suma. Su distribución es Gaussiana como $n$ se hace grande por el teorema del límite central. Lo que es más importante la mayor parte de la distribución se concentra cerca de la media, que es $28.5n$. Esto implica que la densidad feliz números de entre todos los $n$dígitos enteros depende casi exclusivamente de la distribución de los números happy cerca de $28.5n$.
Por ejemplo, hay un pico en el gráfico a alrededor de $n = 400$ acerca de $.185$ densidad. El cálculo de la densidad de feliz números dentro de una desviación estándar de la media de $s(Y_{400})$ obtenemos una densidad de .1911 (el intervalo miré a se $[10916,11884]$). Si asumimos $s(Y_{400})$ es "exactamente" normalmente distribuida y se estimó la densidad de esta manera se podría obtener una mejor aproximación.
Esto significa recoger $n$ s.t. la media de $s(Y_n)$ tierras en el intervalo de $[10^{400},10^{401}-1]$, entonces la densidad de feliz números entre $n$dígitos enteros debe ser de alrededor de $.185$. Probablemente algunas de las opciones de $n$ le dará densidades estrictamente mayor que $.185$ y algunos estrictamente menor. Esto me ha llevado a sospechar que por recorrer en este proceso, la parte superior de la densidad de feliz números pueden ser $1$, y una baja densidad de $0$.
El artículo de Joe Silverman mencionados son de mi propiedad. En ella intento dar por encima de la heurística de un riguroso fundación. Todavía es un borrador y que solo ha sido revisado por uno de mis estudiantes de posgrado, por lo que no voy a decir que es sin duda correcto, aunque estoy muy seguro de que es. He estado trabajando en él durante las últimas semanas, al ver tu pregunta en MO me decidí a seguir adelante y subir un borrador. En el uso que tengo un promedio de argumento para decir que si usted se encuentra experimentalmente un gran intervalo de $n$dígitos enteros ($n$ suficientemente grande) que contienen números happy con la densidad de $d$, entonces la parte superior de la densidad de los números happy es, al menos,$d(1 - o(1))$. Que es donde la parte superior de la densidad de $\geq .18$ y menor densidad de $ \leq .12$ proviene.