xorshift128+算法的真正定义是什么

What is the real definition of the xorshift128+ algorithm?

本文关键字：定义 是什么 算法 xorshift128+      更新时间：2023-10-16

我需要一个好的伪随机数生成器（PRNG），目前的技术似乎是xorshift128+算法。不幸的是，我发现了两个不同的版本。在wikipedia:Xorshift上显示为：

uint64_t s[2];
uint64_t xorshift128plus(void) {
    uint64_t x = s[0];
    uint64_t const y = s[1];
    s[0] = y;
    x ^= x << 23; // a
    s[1] = x ^ y ^ (x >> 17) ^ (y >> 26); // b, c
    return s[1] + y;
}

这似乎已经足够直截了当了。此外，编辑日志似乎显示，这个代码片段是由一个名为"Vigna"的用户添加的，该用户可能是"Sebastiano Vigna"，他是关于xorshift128+的论文的作者：Marsaglia的xorshift生成器的进一步加扰。不幸的是，该文件中的实现略有不同：

uint64_t next(void) {
    uint64_t s1 = s[0];
    const uint64_t s0 = s[1];
    s[0] = s0;
    s1 ^= s1 << 23; // a
    s[1] = s1 ^ s0 ^ (s1 >> 18) ^ (s0 >> 5); // b, c
    return s[1] + s0;
}

除了一些不同的名字外，这两个片段除了最后两个移位之外都是相同的。在维基百科版本中，这些变化是17和26，而论文中的变化是18和5。

有人知道哪种算法是"正确的"吗？这有什么不同吗？这显然是一个使用相当广泛的算法，但使用的是哪个版本？

多亏了@Blastfurt，答案似乎是，根据算法作者的说法，最近的一组常数是：23、18和5。显然，这并不重要，但理论上这些数字比他最初使用的数字要好。Sebastiano Vigna在回应V8 Javascript引擎正在转向使用该算法的消息时发表了这些评论。

我使用的实现是：

uint64_t a = s[0];
uint64_t b = s[1];
s[0] = b;
a ^= a << 23;
a ^= a >> 18;
a ^= b;
a ^= b >>  5;
s[1] = a;
return a + b;
相关文章：