所有离散型随机变量的模拟都是以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一直往后查找,看究竟落在哪一个区间里
优化:将区间按照长度的逆序排序,可以有更大的几率提早结束查找