Question

我想设计一个彩票调度程序，我需要一个非常好的（伪）随机数生成器，类似于LCG，但我想知道是否还有其他更好的选择？我特意寻找用C语言编写的随机生成器。

LCG代码：

unsigned long lcg_rand(unsigned long a)
{
  return (a * 279470273UL) % 4294967291UL;
}

此外，我想知道srand()是否可用于此目的或不准确？

Answer 1

如果你需要简单但体面的质量，我会使用64位LCG的高32位（或更少），可能还有一个应用于输出的回火功能。这样做时，我复制了Mersenne Twister中使用的缓和函数。我不建议实际使用Mersenne Twister，因为它比其他PRNG具有更多的复杂性和内部状态而没有明显更好的品质。

以下是一些示例代码：

static uint32_t temper(uint32_t x)
{
    x ^= x>>11;
    x ^= x<<7 & 0x9D2C5680;
    x ^= x<<15 & 0xEFC60000;
    x ^= x>>18;
    return x;
}
uint32_t lcg64_temper(uint64_t *seed)
{
    *seed = 6364136223846793005ULL * *seed + 1;
    return temper(*seed >> 32);
}

Answer 2

Mersenne Twister是一种选择。另一个选项是Subtract with carry

Answer 3

C rand()函数的大多数实现都使用 LGC 的变体。 rand()，就像任何计算机化的随机生成器不是真正随机的，它只是伪随机的。使用srand()可以改善随机性，但不能使其完美。这取决于srand()中使用的种子的变化和随机性。例如，如果您在rand()中使用相同的相同种子调用srand() n次，则结果将相同。但是如果你每次都调用srand(clock())（并且调用之间经过的时间大于clock()中的滴答时间段），那么你将有一个改进的随机生成器。

这是一个简单的代码示例，其中clock()和支持函数 NotRecentlyUsed（）（对于 min 和 max的小样本）：

#include <ansi_c.h>

#define _UNIQUE_

int randomGenerator(int min, int max);
int NotUsedRecently (int number);

int main(void)
{
    int i=0;
    for(i=0;i<1000;i++)
    {
        printf("%d,\t", randomGenerator(0, 20));
            if(i%20 == 0) printf("\n");
    }
    getchar();
    return 0;   
}

//////////////////////////////////////////////////////
//generates pseudo random numbers between min and max
//If it is desired to use this without a guarantee of uniqueness
//for a specified span of values, un-define _UNIQUE_
//
int randomGenerator(int min, int max)
{
    int random, trying;

    trying = 1;         
    while(trying)
    {
        srand(clock());
        random = (rand()/32767.0)*(max+1);
        ((random >= min)
#ifdef _UNIQUE_
            && NotUsedRecently(random) 
#endif
            ) ? (trying = 0) : (trying = 1);
    }

    return random;
}

//This function is used to guarantee that a span of n generated values
//will not be the same. To adjust the span, change the index in 
//static int recent[n];  Note: it is required that n < (max-min)
//otherwise an infinite loop will occur
int NotUsedRecently (int number)
{
    static int recent[20];//Use No greater value for index than max - min
    int i,j;
    int notUsed = 1;

    for(i=0;i<(sizeof(recent)/sizeof(recent[0]));i++)  (number != recent[i]) ? (notUsed==notUsed) : (notUsed=0, i=(sizeof(recent)/sizeof(recent[0])));
    if(notUsed) 
    {
        for(j=(sizeof(recent)/sizeof(recent[0]));j>1;j--)
        {
            recent[j-1] = recent[j-2];
        }
        recent[j-1] = number;
    }
    return notUsed;
}

与彩票调度程序的LCG相比，更好的（伪）随机数发生器是什么？

3 个答案: