Permítanme comentar la conexión con Combinatoria analítica. Es se sabe que todas las endofunciones en $[n]$ son conjuntos de ciclos de árboles etiquetados etiquetados. Esto se deduce del hecho de que cuando iteramos $f$ obtenemos un camino que termina en un ciclo y se documenta en Cartografía aleatoria Estadísticas de P. Flajolet, donde se derivan las asintóticas. La construcción también apareció en este MSE enlace . La clase combinatoria clase combinatoria viene dada entonces por
$$\def\textsc#1{\dosc#1\csod} \def\dosc#1#2\csod{{\rm #1{\small #2}}} \mathcal{F} = \textsc{SET}(\textsc{CYC}(\mathcal{T})) \quad\text{where}\quad \mathcal{T} = \mathcal{Z} \textsc{SET}(\mathcal{T}).$$
Ahora bien, en el presente caso no puede haber ciclos de dos o más elementos porque los valores de esos ciclos no cumplirían la condición de que $f(f(x)) = f(f(f(x))).$ Esto deja conjuntos de árboles, siendo la raíz del del árbol es un punto fijo. Obsérvese que cualquier hoja con un camino a la raíz incluyendo la raíz que contiene más de tres nodos también rompe el requisito de que $f(f(x)) = f(f(f(x))).$ Esto restringe la clase de clase de árboles a aquellos en los que el camino de cada hoja a la raíz incluyendo la raíz, contiene como máximo tres nodos. Ahora, para caracterizar estos árboles son, primero, un nodo único (punto fijo) o, segundo un nodo único (también un punto fijo) con un conjunto de subárboles adjuntos a él, que a su vez son hojas o tienen un conjunto de hojas que tienen un conjunto de hojas. Con la clase combinatoria $\mathcal{T}$ de estos árboles obtenemos así
$$\mathcal{T} = \mathcal{Z} + \mathcal{Z}\textsc{SET}_{\ge 1} (\mathcal{Z} + \mathcal{Z}\textsc{SET}_{\ge 1}(\mathcal{Z})).$$
La clase combinatoria deseada $\mathcal{F}$ es un bosque de estos árboles dado por
$$\mathcal{F} = \textsc{SET}(\mathcal{T}) = \textsc{SET}(\mathcal{Z} + \mathcal{Z}\textsc{SET}_{\ge 1} (\mathcal{Z} + \mathcal{Z}\textsc{SET}_{\ge 1}(\mathcal{Z}))).$$
Obsérvese que la distinción entre los nodos que no tienen subárboles y los nodos con un conjunto de subárboles como característica que que aparece durante el estudio de los problemas dados no es esencial aquí y podemos utilizar el hecho conveniente de que
$$\mathcal{E} + \textsc{SET}_{\ge 1}(\mathcal{Q}) = \textsc{SET}(\mathcal{Q}).$$
Así pues, tenemos para la clase en cuestión que es
$$\mathcal{F} = \textsc{SET}(\mathcal{Z}\textsc{SET}(\mathcal{Z} \textsc{SET}(\mathcal{Z}))).$$
Lo que ocurre aquí es que para los árboles enraizados de altura máxima $h$ tenemos $\mathcal{T}_{\le 0} = \mathcal{Z}$ y para $h\ge 1$
$$\mathcal{T}_{\le h} = \mathcal{Z} \textsc{SET}(\mathcal{T}_{\le h-1}).$$
Obtenemos instantáneamente el EGF
$$F(z) = \exp(z\exp(z\exp(z))).$$
Extrayendo los coeficientes aquí encontramos
$$n! [z^n] F(z) = n! [z^n] \sum_{q=0}^n \frac{1}{q!} z^q \exp(qz\exp(z)) \\ = n! \sum_{q=0}^n \frac{1}{q!} [z^{n-q}] \exp(qz\exp(z)) \\ = n! \sum_{q=0}^n \frac{1}{q!} [z^{n-q}] \sum_{p=0}^{n-q} \frac{1}{p!} q^p z^p \exp(pz) \\ = n! \sum_{q=0}^n \frac{1}{q!} \sum_{p=0}^{n-q} \frac{1}{p!} q^p [z^{n-q-p}] \exp(pz) \\ = n! \sum_{q=0}^n \frac{1}{q!} \sum_{p=0}^{n-q} \frac{1}{p!} q^p \frac{p^{n-q-p}}{(n-q-p)!}.$$
Esta es la fórmula que se presenta en OEIS A000949 . Al parecer, han optado por simplificar omitiendo el término de $q=0$ ( $q^p=1$ sólo cuando $p=0$ pero entonces $p^{n-q-p} = 0$ ) y extrayendo el término para $q=n$ (que se simplifica a $1$ ) para obtener
$$1 + n! \sum_{q=1}^{n-1} \frac{1}{q!} \sum_{p=0}^{n-q} \frac{1}{p!} q^p \frac{p^{n-q-p}}{(n-q-p)!}.$$