Estou trabalhando em um programa que executa a simulação de Monte Carlo; especificamente, estou usando um algoritmo Metropolis. O programa precisa gerar bilhões de números "aleatórios". Sei que o twister de Mersenne é muito popular na simulação de Monte Carlo, mas gostaria de ter certeza de que estou semeando o gerador da melhor maneira possível.

Atualmente, estou computando uma semente de 32 bits usando o seguinte método:

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);

Tenho a sensação de que existem maneiras muito melhores de garantir novas sementes não repetidas, e tenho certeza que o mt19937_64 pode ser semeado com mais de 32 bits. Alguém tem alguma sugestão?