OFFSET
1,1
COMMENTS
Coincides with sequence of primes p such that x^54 = 2 has a solution mod p for the first 167 terms (and then diverges).
LINKS
MATHEMATICA
ok[p_]:= Reduce[Mod[x^18- 2, p] == 0, x, Integers]=!=False; Select[Prime[Range[200]], ok] (* Vincenzo Librandi, Sep 13 2012 *)
PROG
(PARI) forprime(p=2, 2000, if([]~!=polrootsmod(x^18-2, p), print1(p, ", "))); print();
/* Joerg Arndt, Jul 27 2011 */
(Magma) [p: p in PrimesUpTo(1000) | exists(t){x : x in ResidueClassRing(p) | x^18 eq 2}]; // Vincenzo Librandi, Sep 13 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved