OFFSET
1,16
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
MATHEMATICA
Table[Last[Select[Divisors[n], # <= Sqrt[n] && IntegerQ[Sqrt[#]] &]], {n, 100}]
PROG
(PARI) a(n) = {my(m=1); fordiv(sqrtint(n/core(n)), d, if(d^4 <= n, m=max(m, d))); m^2} \\ Andrew Howroyd, Oct 11 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Oct 11 2023
STATUS
approved