任何可能的正数的随机函数，但可能接近于零

Question

我正在寻找一个随机数函数，它可能会返回任何（正）数，但在统计上不太可能远离零。我有一个算法可以满足我的需求，但似乎效率低下：

def my_random():
    if random.randrange(2):
        return 1 + my_random()
    return 0

我知道可以在不递归的情况下重写它，但这并不能真正否定低效率。任何人都可以提出具有更好性能的替代算法吗？

我的实现语言是python，但算法才是最重要的。

Answer 1

如果您有一个函数返回[0, 1) 范围内的随机浮点数，例如Python 的random.random，那么您可以从exponential distribution 中采样

def random_exponential(scale=1.):
    x = random.random()
    return -math.log(1. - x) * scale

（受我刚刚阅读的 NumPy 源代码的启发；如果可以选择使用 NumPy，您也可以致电 np.random.exponential。）

如果您需要随机整数，请考虑截断这些。