OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
g:= proc(n) local p;
p:= 1;
do
p:= nextprime(p);
if numtheory:-quadres(n, p) = -1 then return p fi
od
end proc:
map(g, select(isprime, [seq(i, i=5..10000, 8)])); # Robert Israel, Apr 17 2023
MATHEMATICA
f[n_] := Block[{k = 2}, While[ JacobiSymbol[n, Prime[k]] == 1, k++ ]; Prime[k]]; f /@ Select[ Prime[ Range[435]], Mod[ #, 8] == 5 &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jun 24 2004
STATUS
approved