OFFSET
1,1
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
MATHEMATICA
a = {}; Do[poly = PolynomialMod[x^5 - x - 1, Prime[n]]; c = FactorList[poly, Modulus -> Prime[n]]; If[Sum[c[[m]][[2]], {m, 1, Length[c]}] == 6, AppendTo[a, Prime[n]]], {n, 1, 10000}]; a
PROG
(PARI) isok(n)=#factormod(x^5-x-1, n)[, 2]==5;
forprime(n=2, 10^6, if(isok(n), print1(n, ", "))); \\ Joerg Arndt, Dec 07 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Dec 01 2007
EXTENSIONS
Terms a(35) and beyond from G. C. Greubel, Dec 06 2016
STATUS
approved