OFFSET
1,1
COMMENTS
This is not the same as "Primes congruent to 1 mod 59", A216315. The first missing number is A216315(27) = 11447. - Zak Seidov, Sep 03 2012
LINKS
Zak Seidov, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Prime[Range[PrimePi[21000]]], ! MemberQ[PowerMod[Range[#], 59, #], Mod[2, #]] &] (* Bruno Berselli, Sep 20 2012 *)
ok[p_] := Reduce[Mod[x^59 - 2, p] == 0, x, Integers] == False; Select[Prime[Range[2500]], ok] (* Vincenzo Librandi, Sep 20 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(21000) | forall{x: x in ResidueClassRing(p) | x^59 ne 2}]; // Bruno Berselli, Sep 20 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Klaus Brockhaus, Jan 25 2001
STATUS
approved