mpz_probab_prime_p (10 reps) from gmplib used for the 4 PRPs. It does some trial divisions, then some Miller-Rabin probabilistic primality tests.

To be completed...