OFFSET
0,4
COMMENTS
See A225768 for motivation and references.
PROG
(PARI) (a, b=5)->{#factor(x^b+a)~==1&for(n=1, 9e9, ispseudoprime(n^b+a)&return(n)); a==1 || print1("/*"factor(x^b+a)"*/")} \\ For illustrative purpose only. The polynomial is factored to avoid an infinite search loop when it is composite. But a factored polynomial can yield a prime when all factors but one equal 1. This happens for b=4, n=4, cf. A225766.
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Jul 25 2013
STATUS
approved