See python implementation: https://github.com/numpy/numpy/blob/3a3e90973bba4f0865503dc6bffd9f3080080b38/numpy/random/mtrand/mtrand.pyx#L1169
