OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1)=3037 since p=3037 is the least prime of the form 4k+1 for which sp-(floor(sqrt(sp)))^2 is not a square for s=1..16, but 17p-(floor(sqrt(17p)))^2 is a square (for p=3037 it is 100).
MAPLE
filter:= proc(p) local s;
if not isprime(p) then return false fi;
for s from 1 to 17 do
if issqr(s*p - floor(sqrt(s*p))^2) then return evalb(s=17) fi
od;
false
end proc:
select(filter, [seq(i, i=1..10000, 4)]); # Robert Israel, Jan 22 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Sep 30 2008
STATUS
approved