Random Number系列：离散随机变量的分布模拟(Inverse Transform方法)

所有离散型随机变量的模拟都是以uniform distribution (0,1) U为基础的

如果离散随机变量的分布为：

P{X=xj}=pj, j=0,1,... （∑ pj )=1）

令X= x0, if U<p0

x1, if p0<=U<p1

x2, if p1<=U<p2

xj, if ∑ (i=0 to j-1) pi<=U<∑(i=0 to j) pi

则

P{X=xj} = P{  ∑ (i=0 to j-1) pi<=U<∑(i=0 to j) pi } = pj（由uniform的性质得到）

写成代码其实就是从x0一直往后查找，看究竟落在哪一个区间里

优化：将区间按照长度的逆序排序，可以有更大的几率提早结束查找