OFFSET
1,1
COMMENTS
LINKS
Bruno Berselli, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Prime[Range[PrimePi[500]]], ! MemberQ[PowerMod[Range[#], 18, #], Mod[2, #]] &]
ok[p_]:= Reduce[Mod[x^18 - 2, p] == 0, x, Integers] == False; Select[Prime[Range[80]], ok] (* Vincenzo Librandi, Sep 20 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(500) | forall{x: x in ResidueClassRing(p) | x^18 ne 2}];
(Magma) [p: p in PrimesUpTo(400) | not exists{x : x in ResidueClassRing(p) | x^18 eq 2} ]; // Vincenzo Librandi, Sep 20 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Sep 13 2012
STATUS
approved