Найти где-то, надо посмотреть как в пакетах это делается, либо взять готовый.
Но стало интересно как там делают это.
Оно работает очень быстро, это не считается проблемой, значит там есть тривиальное решение.
Сначала берётся случайное нечётное натуральное, а потом уже проверяется на простоту вероятностными тестами. Если тест провален, число увеличивается на 2 или 4 в зависимости на какую цифру оно оканчивается.