python - 随机数算法

Question

我需要按组生成随机数:100、500、1000 和 10000 个数字制服和高斯分布。这对于创建一些直方图和其他统计资料是必要的。

条件是不使用python自带的随机函数，所以想着用这个方法(线性同余生成器):Xn+1 ≡ (aXn + c) mod m。这里我需要 4 个变量。有人可以告诉我如何实现这个算法吗？我假设 m 变量第一次是 100

Answer 1

您已经知道这是 Linear congruential generator , 那么读起来有什么难的呢？

它告诉你你已经知道的公式以及如何选择它们的解释:

The period of a general LCG is at most m, and for some choices of factor a much less than that. Provided that the offset c is nonzero, the LCG will have a full period for all seed values if and only if:

1) c and m are relatively prime

2) a - 1 is divisible by all prime factors of m

3) a - 1 is a multiple of 4 if m is a multiple of 4.

他们甚至在下表中为您提供了这些值的一些示例。这足以实现一个简单的功能:

def LCG(seed, n, a=1664525, c=1013904223, m=2**32):
    numbers = []
    for i in xrange(n):
        seed = (a * seed + c) % m
        numbers.append(seed)

    return numbers

print LCG(3, 5)