pb <- function(n) 1 - factorial(365)/(factorial(365-n)*365^n)