The OEIS is supported by the many generous donors to the OEIS Foundation.
%I #18 Sep 08 2022 08:44:59
%S 2,3,11,43,59,83,107,131,179,227,251,257,281,283,307,347,419,443,467,
%T 491,499,563,587,617,643,659,683,691,739,811,827,881,947,971,1019,
%U 1049,1051,1091,1097,1163,1187,1193,1217,1259,1283,1307,1427,1451,1459,1481
%N Primes p such that x^24 = -2 has a solution mod p.
%C Complement of A216743 relative to A000040. - _Vincenzo Librandi_, Sep 17 2012
%H Vincenzo Librandi, <a href="/A051081/b051081.txt">Table of n, a(n) for n = 1..1000</a>
%p isA051081 := proc(p) local x; for x from 0 to p-1 do if (x^24 mod p) = (-2 mod p) then RETURN(true) ; fi; od: RETURN(false) ; end: for i from 1 to 300 do p := ithprime(i) ; if isA051081(p) then printf("%d,",p) ; fi; od: # _R. J. Mathar_, Oct 15 2008
%t ok[p_]:= Reduce[Mod[x^24 + 2, p] == 0, x, Integers] =!= False; Select[Prime[Range[500]], ok] (* _Vincenzo Librandi_, Sep 15 2012 *)
%o (PARI) /* see A051071 */
%o (Magma) [p: p in PrimesUpTo(1500) | exists(t){x : x in ResidueClassRing(p) | x^24 eq - 2}]; // _Vincenzo Librandi_, Sep 15 2012
%K nonn,easy
%O 1,1
%A _N. J. A. Sloane_
%E More terms from _R. J. Mathar_, Oct 15 2008