Eu gero um número aleatório, testo se ele é divisível pelos primos abaixo de 5000, se ele não passar, eu adiciono 2 e tento denovo, se ele passar dessa etapa, eu teto-o usando o Teste de Miller-Rabin, com 10 checagens, se ele falhar, ele volta o processo todo, mas a quantidade de checagens pode ser alterada passando o parametro -checks (o numero desejado)
Ex:
pgen -bits 1024 -checks 5

Obrigado por ter me avisado do repo