%I #7 Sep 08 2022 08:46:03
%S 5,7,11,13,23,29,31,37,41,47,53,61,71,79,101,103,109,127,131,149,151,
%T 157,167,173,181,191,197,199,211,223,229,239,251,263,269,271,277,281,
%U 293,311,317,331,349,359,367,373,383,389,397,401,421,431,439,461,463,479
%N Primes p such that x^50 = -2 has no solution mod p.
%C Complement of A051094 relative to A000040.
%H Vincenzo Librandi, <a href="/A216770/b216770.txt">Table of n, a(n) for n = 1..1000</a>
%t ok[p_] := Reduce[Mod[x^50 + 2, p] == 0, x, Integers] == False;Select[Prime[Range[200]], ok]
%o (Magma) [p: p in PrimesUpTo(600) | not exists{x : x in ResidueClassRing(p) | x^50 eq -2}];
%Y Cf. A000040, A051094.
%K nonn,easy
%O 1,1
%A _Vincenzo Librandi_, Sep 16 2012
|