在 Python 中重新创建一个 RNG

Question

我正在尝试重新创建一个在 Matlab 中使用的 乘法同余算法（找到它here，第 9.2 节，阅读它以了解一些上下文）

基本上，我有 4 个变量

a = 13
c = 0
m = 31

最后一个，我无法格式化为代码，因为它破坏了下标：

x₀= 1

有了这四个，我想重现这样的东西：

x_k+1= a *x_k+ c mod m

到目前为止，我有：

a = 13
c = 0
m = 31
base = 2 # Does not work with 0 for some reason
x = int(x='1', base=base)

for i in range(8):
    rand = a * x + c % m
    base += 1
    x = int(x='1', base=base)  
    print(rand)

我的输出是：

13, 13, 13, 13, 13, 13, 13, 13...

即使，根据文档，应该是：

1, 13, 14, 27, 10, 6, 16, 22...

所以我不确定我是否完全误解了这个问题，或者我想要完成的事情在 Python 中是不可能的，或者我只是发疯了，但我非常需要一些建议。非常感谢任何见解。

Answer 1

看起来问题是公式应该是

x_k+1 = (a * x_k + c) mod m

这在文档中并不清楚，但是如果您使用 c mod m 就像他们拥有的那样，您每次都会执行相同的操作（0mod31）。这段python代码运行正常：

a = 13
c = 0
m = 31

x = [1]

for i in range(1, 9):
    x.append( (a * x[i-1] + c)%m )
    print x[i],

13 14 27 10 6 16 22 7