A continuación explico cómo enumerar todas las cubiertas sin mala pares de forma explícita. El número exacto resulta ser $$3668033946384704437729512814619767610579526911188666362431432294400.$$
Dividiendo por $52!$, me sale aproximadamente el $$0.045476282331.$$
Vamos a escribir una relación de recurrencia para $f(H,k)$, que es el número de parciales cubiertas con $t$ tarjetas apareciendo $H_t$ veces, la última carta haciendo su $k$th apariencia. El número total de cubiertas es, a continuación,$f((0,0,0,0,13),4)$.
El caso base es $f((12,1,0,0,0),1) = 52$. Ahora supongamos que se nos ha dado algunos de los $H$ $k$ tal que $H_k \geq 1$. La primera forma de $H'$ mediante el movimiento de un "rango" de$k$$k-1$. Ahora, para cada una de las $t \neq k - 1$, podemos extender cada una de las $f(H_k,t)$ parcial cubiertas de una nueva cubierta en $H_k(k-1) \cdot (4-(k-1))$ maneras. Para $t = k-1$, podemos extender la $f(H_k,k-1)$ parcial de cubiertas en sólo $(H_k(k-1)-1) \cdot (4-(k-1))$ formas (desde un rango).
Sólo hay 9481 pares de $H,k$ a considerar, por lo que incluso mi pésima implementación de python que utiliza una tabla hash para mantener un seguimiento de todas las parejas es lo suficientemente rápido.
Michael Lugo método da $(1-3/51)^{51} \approx 0.0454176$, que está muy cerca de la respuesta correcta. La probabilidad correcta es ligeramente más grande ya que si nos imaginamos a un proceso en el que se revelan las cartas una por una, la probabilidad de que la siguiente carta es mala, comienza a $3/51$ pero eventualmente disminuye (aproximadamente).
Su aproximación a $e^{-3} \approx 0.049787$ es demasiado grande, y eso es debido a que $(1-3/n)^n$ enfoques $e^{-3}$ de los de abajo.
leonbloy la aproximación también es demasiado pequeño, por razones similares. Podemos imaginar un proceso en el que los rangos se reveló uno por uno. Después de que el primer rango se revela, un menor número de posiciones adyacentes, y por lo que la probabilidad disminuye.