OFFSET
1,1
COMMENTS
LINKS
MATHEMATICA
ok[p_]:= Reduce[Mod[x^30 - 2, p] == 0, x, Integers] =!= False; Select[Prime[Range[300]], ok] (* Vincenzo Librandi, Sep 14 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(1300) | exists(t){x : x in ResidueClassRing(p) | x^30 eq 2}]; // Vincenzo Librandi, Sep 14 2012
(PARI) forprime(p=2, 1300, if(0!=#polrootsmod(x^30-2, p), print1(p, ", "))); \\ Joerg Arndt, Apr 19 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved
