M <- 20000
N <- 12^8
lq <- lfactorial(N) - lfactorial(N-M) - M*log(N)
1-exp(lq)