Choice distribution now works with earlier Numpy

Fixes #536
I was reading that for these to be equivalent, it had to be `np.searchsorted(np.cumsum(self.p), rng.rand(n), side='right'))`. Is this not true in our case?