OFFSET
1,1
COMMENTS
This sequence is a subsequence of A135844.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(p) isprime(p) and nops([msolve(x^5-x-1, p)])=5 end proc:
select(filter, [seq(seq(10*k+j, j=[3, 7, 9]), k=0..10000)]); # Robert Israel, Jul 03 2018
MATHEMATICA
a = {}; Do[If[Mod[Prime[n], 1, 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
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Dec 01 2007
STATUS
approved