He aquí una útil en la fórmula (que parece muy difícil de poner en forma cerrada). Si ya has dibujado $n$ tarjetas, a continuación, $P(\mathrm{success})$ es
$$
\sum_{\spadesuit+\heartsuit+\diamondsuit+\clubsuit=n;\ 0\leq \spadesuit,\heartsuit,\diamondsuit,\clubsuit\leq13}\frac{13-\min(\spadesuit,\heartsuit,\diamondsuit,\clubsuit)}{52-n}\cdot\frac{\binom{n}{\spadesuit, \heartsuit, \diamondsuit, \clubsuit}\binom{52-n}{13-\spadesuit, 13-\heartsuit, 13-\diamondsuit, 13-\clubsuit}(13!)^4}{52!}
$$
(Esta fórmula se utiliza el coeficiente multinomial, que en este caso cuenta el número de formas de elegir los $\spadesuit$ posiciones de picas, $\heartsuit$ de las restantes posiciones de corazones, etc. de la $n$ primeras cartas a dibujar.)
El primer factor es la probabilidad de éxito dada una partición $\spadesuit+\heartsuit+\diamondsuit+\clubsuit=n$ (la partición representa cómo el primer $n$ se distribuyeron tarjetas, sin importar el orden en que entró). El resto es la probabilidad de que la primera $n$ tarjetas de proporcionar esa partición (contando todos los posibles ordenamientos, dividido por el total de posibles ordenamientos de una baraja).
Añadido posterior:
La suma por encima de las condiciones sobre la partición que representa las cartas que ya han sido elaborados. También podríamos suma sobre las particiones que representan las cartas aún no se dibuja. Que pone la informática más eficiente para la segunda mitad de la cubierta. Si $r$ tarjetas, a continuación, $P(\mathrm{success})$ es
$$
\sum_{\spadesuit+\heartsuit+\diamondsuit+\clubsuit=r;\ 0\leq \spadesuit,\heartsuit,\diamondsuit,\clubsuit\leq13}\frac{\max(\spadesuit,\heartsuit,\diamondsuit,\clubsuit)}{r}\cdot\frac{\binom{r}{\spadesuit, \heartsuit, \diamondsuit, \clubsuit}\binom{52-r}{13-\spadesuit, 13-\heartsuit, 13-\diamondsuit, 13-\clubsuit}(13!)^4}{52!}
$$
La evaluación de estas sumas requiere la enumeración de las particiones. Esto se podría hacer con un poco de programación, pero lo hice por $r,n\leq10$ en Excel y encontrado el siguiente $P_n(\mathrm{success})$:
$$
\begin{align}
P_0&=0.25 &
P_1&\approx0.2549 &
P_2&=0.26 &
P_3&\approx0.2653 \\
P_4&\approx0.2686 &
P_5&\approx0.2710 &
P_6&\approx0.2733 &
P_7&\approx0.2762 \\
P_8&\approx0.2788 &
P_9&\approx0.2809 &
P_{10}&\approx0.2830 &
&\cdots\\
&\cdots &
&\cdots &
P_{42}&\approx0.4005 &
P_{43}&\approx0.4122 \\
P_{44}&\approx0.4258 &
P_{45}&\approx0.4392 &
P_{46}&\approx0.4548 &
P_{47}&\approx0.4836 \\
P_{48}&\approx0.5201 &
P_{49}&\approx0.5514 &
P_{50}&\approx0.6176 &
P_{51}&=1
\end{align}
$$