0

所以我有一个使用 Mersenne Twister 的自定义随机器类（我使用的代码改编自这个站点）。一切似乎都运行良好，直到我开始测试不同的种子（我通常使用 42 作为种子，以确保每次运行我的程序时，结果都是相同的，因此我可以看到代码更改如何影响事物）。

事实证明，无论我选择什么种子，代码每次都会产生完全相同的数字序列。显然我做错了什么，但我不知道是什么。这是我的种子功能：

void Randomizer::Seed(unsigned long int Seed)
{
    int ii;
    x[0] = Seed & 0xffffffffUL;
    for (ii = 0; ii < N; ii++)
    {
        x[ii] = (1812433253UL * (x[ii - 1] ^ (x[ii - 1] >> 30)) + ii);
        x[ii] &= 0xffffffffUL;
    }
}

这是我的 Rand() 函数

unsigned long int Randomizer::Rand()
{
    unsigned long int Result;
    unsigned long int a;
    int ii;

    // Refill x if exhausted
    if (Next == N)
    {
        Next = 0;

        for (ii = 0; ii < N - 1; ii++)
        {
            Result = (x[ii] & U) | x[ii + 1] & L;
            a = (Result & 0x1UL) ? A : 0x0UL;
            x[ii] = x[( ii + M) % N] ^ (Result >> 1) ^ a;
        }

        Result = (x[N - 1] & U) | x[0] & L;
        a = (Result & 0x1UL) ? A : 0x0UL;
        x[N - 1] = x[M - 1] ^ (Result >> 1) ^ a;
    }
    Result = x[Next++];

    //Improves distribution
    Result ^= (Result >> 11);
    Result ^= (Result << 7) & 0x9d2c5680UL;
    Result ^= (Result << 15) & 0xefc60000UL;
    Result ^= (Result >> 18);

    return Result;
}

各种值是：

#define A 0x9908b0dfUL
#define U 0x80000000UL
#define L 0x7fffffffUL

int Randomizer::N = 624;
int Randomizer::M = 397;
int Randomizer::Next = 0;
unsigned long Randomizer::x[624];

谁能帮我弄清楚为什么不同的种子不会产生不同的数字序列？

4

1 回答 1

5

您的 Seed() 函数分配给x[0]，然后在 开始循环ii=0，它用未定义的值覆盖x[0]（它引用x[-1]）。从 1 开始循环，您可能已经准备就绪。

编写自己的随机化器是危险的。为什么？很难做到正确（见上文），很难知道你是否做对了，如果做错了，依赖正确分布的随机数的事情将无法正常工作。希望那东西不是尾巴很重要的密码学或统计建模……考虑使用std::random，或者如果您还没有使用 C++11，请使用boost::random。

于 2012-06-11T19:37:25.020 回答