$p = 1 - ( factorial($n) * c(2**32, $n) / 365**$n )