OFFSET
1,5
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..20000
MAPLE
a:= proc(n) local k; for k from (s-> `if`(s^2<n, s+1, s))(isqrt(n))
while not issqr(irem(k^2, n)) do od; isqrt(irem(k^2, n))
end:
seq(a(n), n=1..120);
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Alois P. Heinz, Feb 13 2019
STATUS
approved