线性同余生成器 C

Question

• 平台：Linux 3.2.0 x86 (Debian Wheezy)
• 编译器：GCC 4.7.2 (Debian 4.7.2-5)

我正在编写一个线性同余生成器。我读过线性同余生成器生成的高阶位具有更高的周期。我想使用高阶位，但我不知道如何使用。我当前的实现使用了《数值食谱》一书中的“a”和“c”值，因为方程 ax+c 的最大可能值可以表示为 64 位整数。

#define ran32(seed) (uint32_t)((seed) * UINT64_C(1664525) + UINT64_C(1013904223))

我尝试了以下方法来消除低位，但据我所知，它不会返回无符号 32 位整数的整个范围内的值。

#define ran32(seed) (uint32_t)(((seed) * UINT64_C(1664525) + UINT64_C(1013904223)) / UINT64_C(256))

这个想法是用 2^8(256) 积分除以 'foo' >> 8。这个想法的灵感来自我在 rand() 的文档中找到的以下代码段。

static unsigned long next = 1;
int myrand(void)  /* RAND_MAX assumed to be 32767. */
{
    next = next * 1103515245 + 12345;
    return((unsigned)(next/65536) % 32768);
}

void mysrand(unsigned seed)
{
    next = seed;
}

Answer 1

一个问题可能是计算：

seed * 1664525 + 1013904223

生成一个 32 位值 (uint32_t)，然后将其分配给一个 64 位变量。您可能希望将其作为 64 位计算。

但是，除非您经常溢出 64 位数字，否则您的最高有效字节不一定是非常随机的。

我认为您的代码也容易受到字节序问题的影响。

您对表达式最小值较大的分析忽略了 uint32_t 所包含的模数 2³² 算术。