如何用c++生成poisson分布的样本？

假定平均数为lambda, rand()为产生（0，1）的随机数函数。

给出算法也可以。

答案或建议：

algorithm poisson random number (Knuth):
init:
Let L ← e-λ, k ← 0 and p ← 1.
do:
k ← k + 1.
Generate uniform random number u in [0,1] and let p ← p × u.
while p ≥ L.
return k - 1.

文章来自：http://www.baikewenda.com/h/0903/a35528.html