He aquí una aproximación por simulación: En un millón de sesiones de 70 tiradas del dado según la pregunta 1, he encontrado que es casi seguro obtener al menos 20 aciertos. [Nótese que he etiquetado las caras del dado con los números de aciertos correspondientes; R trata los tres 0
's como caras diferentes].
set.seed(926); m = 10^6
die = c(0,0,0,1,1,2); rolls=70
hits = replicate(m, sum(sample(die,rolls,rep=T)))
mean(hits >= 20)
[1] 0.999998
Si se pretendía hacer 30 tiradas del dado (como se mencionó en un comentario), entonces es un poco más del 50%.
set.seed(926); m = 10^6
die = c(0,0,0,1,1,2); rolls=30
hits = replicate(m, sum(sample(die,rolls,rep=T)))
mean(hits >= 20)
[1] 0.540681
En cualquier caso, te sugiero que pruebes el método CLT propuesto por @awkward. Setenta tiradas, o incluso 30, deberían ser suficientes para obtener una aproximación razonablemente buena. aproximación. Luego puedes comparar los resultados del CLT con los de la la simulación.
Nota: Aquí hay un histograma (azul) para el número de aciertos en 70 tiradas, junto con el curva normal mejor ajustada (naranja). El ajuste no es perfecto, pero es lo suficientemente bueno suficiente para mostrar que casi no hay probabilidad por debajo de 20.
# for rolls=70
mu = mean(hits); sg = sd(hits); mu; sg
[1] 46.66131
[1] 6.235426
hist(hits, prob=T, col="skyblue2", ylim=c(0,.07))
curve(dnorm(x, mu, sg), add=T, col="orange", lwd=2)