M <- 20000
N <- 12^8
q <- factorial(N)/(factorial(N-M)*N^M)
1-q