Я работаю над программой, которая запускает симуляцию Монте-Карло; в частности, я использую алгоритм Метрополис. Программа должна генерировать, возможно, миллиарды «случайных» чисел. Я знаю, что твистер Мерсенна очень популярен для симуляции Монте-Карло, но я хотел бы убедиться, что я сею генератор наилучшим образом.

В настоящее время я вычисляю 32-битное начальное число, используя следующий метод:

mt19937_64 prng; //pseudo random number generator
unsigned long seed; //store seed so that every run can follow the same sequence
unsigned char seed_count; //to help keep seeds from repeating because of temporal proximity

unsigned long genSeed() {
    return (  static_cast<unsigned long>(time(NULL))      << 16 )
         | ( (static_cast<unsigned long>(clock()) & 0xFF) << 8  )
         | ( (static_cast<unsigned long>(seed_count++) & 0xFF) );
}

//...

seed = genSeed();
prng.seed(seed);

У меня есть ощущение, что есть гораздо лучшие способы обеспечить неповторяющиеся новые семена, и я вполне уверен, что mt19937_64 можно посеять с более чем 32-битными. У кого-нибудь есть предложения?