Don’t assume NumPy.vectorize is faster
A look at solving the birthday paradox with different vectorized functions