giorni <- c(1:365)
prove <- 100000
fass <- rep(NA,prove)
for (i in 1:prove)
{
    out <- rep(0,365)
    n <- 0
    repeat {
       ris <- sample(giorni,size=1)
       out[ris] <- out[ris] + 1
       n <- n + 1
       if (out[ris]>1) break
    }
    fass[i] <- n
}
table(fass)/prove
plot(table(fass)/prove, xlab="n. persone da intervistare", ylab="probabilità")