Update pn_generator.py

updated with factorial primes