%I #13 Sep 08 2022 08:45:02
%S 11,31,41,61,71,101,131,151,181,191,211,251,271,281,311,331,401,421,
%T 461,491,521,541,601,631,661,691,701,751,761,811,821,881,941,991,1021,
%U 1031,1051,1061,1091,1151,1171,1201,1231,1291,1301,1321,1361,1381,1451
%N Primes p such that x^25 = 2 has no solution mod p.
%C Complement of A049557 relative to A000040.
%H Vincenzo Librandi, <a href="/A059313/b059313.txt">Table of n, a(n) for n = 1..1000</a>
%t ok[p_] := Reduce[Mod[x^25 - 2, p] == 0, x, Integers] == False; Select[Prime[Range[250]], ok ] (* _Vincenzo Librandi_, Sep 20 2012 *)
%o (Magma) [p: p in PrimesUpTo(1500) | forall{x: x in ResidueClassRing(p) | x^25 ne 2}]; // _Vincenzo Librandi_, Sep 20 2012
%Y Cf. A000040, A049557.
%K nonn,easy
%O 1,1
%A _Klaus Brockhaus_, Jan 25 2001
|