OFFSET
1,1
COMMENTS
Differs from A051085 first at the 541st entry, at p=15809. - R. J. Mathar, Oct 14 2008
From Christopher J. Smyth, Jul 24 2009: (Start)
Differs from A163183 (primes dividing 2^j+1 for some odd j) at the 827th entry, at p=25601. See comment at A163186 for explanation.
Sequence is union of A163183 and A163186 (primes p such that the equation x^64 = -2 mod p has a solution, and ord_p(-2) is even).
(End)
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
MATHEMATICA
ok[p_]:= Reduce[Mod[x^64 + 2, p] == 0, x, Integers] =!= False; Select[Prime[Range[400]], ok] (* Vincenzo Librandi, Sep 16 2012 *)
PROG
(PARI)
forprime(p=2, 2000, if([]~!=polrootsmod(x^64+2, p), print1(p, ", "))); print();
/* Joerg Arndt, Jun 24 2012 */
(Magma) [p: p in PrimesUpTo(1200) | exists(t){x : x in ResidueClassRing(p) | x^64 eq - 2}]; // Vincenzo Librandi, Sep 16 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
More terms from Joerg Arndt, Jul 27 2011
STATUS
approved