Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #33 Jun 14 2022 16:45:21
%S 2,3,3,5,6,7,8,8,10,11,11,13,14,15,14,17,18,19,20,21,22,23,23,24,26,
%T 27,27,29,30,31,32,31,34,35,35,37,38,39,38,41,42,43,44,45,46,47,48,48,
%U 50,51,51,53,54,55,56,56,58,59,59,61,62,63,63,65,66,67
%N Largest quadratic nonresidue modulo n (with n >= 3).
%C The largest nonnegative integer less than n which is not a square modulo n.
%C If p is a prime congruent 3 modulo 4 then a(p) = p-1 since -1 is not a quadratic residue for such primes.
%H Robert Israel, <a href="/A334819/b334819.txt">Table of n, a(n) for n = 3..10000</a>
%e The squares modulo 3 are 0 and 1. Therefore a(3) = 2. The nonsquares modulo 4 are 2 and 3 which makes a(4) = 3. Modulo 5 we have 0, 1 and 4 as squares giving a(5) = 3. 0, 1 and 4 are also the squares modulo 6 resulting in a(6) = 5. Since 10007 is a prime of the form 4*k + 3, a(10007) = 10006.
%p f:= proc(n) local k;
%p for k from n-1 by -1 do if numtheory:-msqrt(k,n)=FAIL then return k fi
%p od
%p end proc:
%p map(f, [$3..100]); # _Robert Israel_, May 14 2020
%t a[n_] := Module[{r}, For[r = n-1, r >= 1, r--, If[!IntegerQ[Sqrt[Mod[r, n]] ], Return[r]]]];
%t a /@ Range[3, 100] (* _Jean-François Alcover_, Aug 15 2020 *)
%o (PARI) a(n) = forstep(r = n - 1, 1, -1, if(!issquare(Mod(r, n)), return(r)))
%Y Cf. A020649, A047210 (the largest square modulo n), A192450 (a(n)=n-1).
%K nonn,easy
%O 3,1
%A _Peter Schorn_, May 12 2020