Even faster implementation with slow_primes lib