Denotar por $q_r(j)$ la probabilidad de que después de $r\geq0$ lanza exactamente $j\geq0$ urnas están ocupadas. A continuación, $p_0(j)=\delta_{0j}$ y
$$q_r(j)={j\over n} q_{r-1}(j)+{n-(j-1)\over n}q_{r-1}(j-1)\qquad(r\geq1)\ .$$
Deje $f_r(j):=n^r q_r(j)$; a continuación,
$$f_r(j)=j f_{r-1}(j)+\bigl(n-(j-1)\bigr)f_{r-1}(j-1)\ .$$
Tenga en cuenta que $r$ no aparece en los coeficientes de esta recursividad. Entonces es fácil ver que, de hecho,
$$f_r(j)=c_{r,j}\prod_{i=0}^{j-1} (n-i)$$
para las constantes $c_{r,j}$ que satisfacen la recursividad
$$c_{r,j}=j c_{r-1,j} + c_{r-1,j-1}\ .$$
Ahora esta es la recursividad de la fórmula de los números de Stirling del segundo tipo, que se denota por a $S(r,j)$ en Stanley combinatoria Enumerativa, vol. I. Uno fácilmente se comprueba que, de hecho,$c_{r,j}=S(r,j)$. Por lo tanto, obtenemos
$$q_r(j)={1\over n^r} S(r,j)\prod_{i=0}^{j-1}(n-i)\ .$$
Al final, el número de bolas en ${\rm urn}_0$ es igual al número de desocupados urnas de entre ${\rm urn}_1,\ldots,{\rm urn}_n$. La probabilidad de $p(k)$ que este número es $k$ es el dado por
$$p(k)=q_n(n-k)={n!\ S(n,n-k)\over k!\ n^n}\ .$$