|
%I #19 Oct 31 2020 04:35:56
%S 7,14,17,31,34,41,49,62,79,82,89,98,103,113,119,127,137,151,158,161,
%T 178,199,206,217,223,226,233,238,254,257,271,274,281,287,289,302,322,
%U 329,343,353,367,391,398,401,434,439,446,449,463,466,487,497,511,514,521,527,542
%N Numbers whose least quadratic nonresidue (A020649) is 3.
%C n such that n is not divisible by 4, all primes dividing n are in A038873, and at least one prime dividing n is in A003630. - _Robert Israel_, Jul 19 2017
%H Robert Israel, <a href="/A025021/b025021.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/QuadraticNonresidue.html">Quadratic Nonresidue</a>.
%p select(t -> numtheory:-quadres(2,t) = 1 and numtheory:-quadres(3,t)=-1, [$1..1000]); # _Robert Israel_, Jul 19 2017
%t Select[Range[500], Min @ Complement[Range[# - 1], Mod[Range[#/2]^2, #]] == 3 &] (* _Amiram Eldar_, Oct 31 2020 *)
%o (PARI) residue(n,m)={local(r);r=0;for(i=1,floor(m/2),if(i^2%m==n,r=1));r}
%o isA025021(n)=residue(2,n) && !residue(3,n) \\ _Michael B. Porter_, Apr 18 2010
%Y Cf. A020649, A025020, A025022, A025023, A025024, A025025, A025026, A025027, A025028, A025029.
%Y Cf. A003630, A038873.
%K nonn
%O 1,1
%A _David W. Wilson_
|
