#2 is double primorial, tested with gmplib mpz_probab_prime_p(k, 10)

To be completed...