import Text.Printf main=mapM(\n->printf"n=%.f, p=%.8f\n"n$1-product[366-n..365::Double]/365**n)[2..100]