OFFSET
1,1
COMMENTS
Complement of A042966 relative to A000040. Coincides for the first 96 terms with the sequence of primes p such that x^49 = 2 has no solution mod p (first divergence is at 4999, cf. A059667). - Klaus Brockhaus, Feb 04 2001
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
x^7 = 2 has no solution mod 29, so 29 is in the sequence.
8^7 = 2097152 and (2097152 - 2)/31 = 67650, so 31 is not in the sequence.
MATHEMATICA
sevPow2ModPQ[p_] := Reduce[Mod[x^7 - 2, p] == 0, x, Integers] == False; Select[Prime[Range[700]], sevPow2ModPQ] (* Vincenzo Librandi, Sep 19 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(3000) | forall{x: x in ResidueClassRing(p) | x^7 ne 2}]; // Vincenzo Librandi, Aug 21 2012
(Magma) [p: p in PrimesUpTo(2600) | not exists{x : x in ResidueClassRing(p) | x^7 eq 2} ]; // Vincenzo Librandi, Sep 19 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved