OFFSET
1,1
REFERENCES
M. Kneser, Quadratische Formen, Springer, 2002; see Hilfssatz 18.3.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
n=3, p = 73, a(3) = q = 5: Legendre(73,5) = -1.
MAPLE
f:= proc(p) local q;
q:= 3:
do
if numtheory:-quadres(p, q) = -1 then return q fi;
q:= nextprime(q);
od;
end proc:
map(f, select(isprime, [seq(p, p=1..10000, 8)])); # Robert Israel, May 06 2019
MATHEMATICA
f[n_] := Prime[ Position[ JacobiSymbol[n, Select[Range[3, n - 1], PrimeQ[ # ] &]], -1][[1, 1]] + 1]; f /@ Select[ Prime[ Range[435]], Mod[ #, 8] == 1 &] (* Robert G. Wilson v, Jun 23 2004 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Jun 19 2004
EXTENSIONS
More terms from Robert G. Wilson v, Jun 23 2004
STATUS
approved