import Text.Printf
f::Int->Double
f n=1-(product$map(\x->1-x/365)[1..fromIntegral$n-1])
main=mapM(s(printf"n=%d, p=%.8f\n")f)[2..100]where s f g x=f x$g x