OFFSET
2,1
COMMENTS
The terms are either primes or squares of primes.
LINKS
Robert Israel, Table of n, a(n) for n = 2..10000
MAPLE
f:= proc(n) local F, t;
F:= sort(ifactors(n)[2], (a, b) -> a[1]<b[1]);
if nops(F) >= 2 and F[2, 1] < F[1, 1]^2 then F[2, 1] else F[1, 1]^2 fi
end proc:
map(f, [$2..100]); # Robert Israel, Oct 09 2024
MATHEMATICA
a[n_] := Divisors[n^2][[3]]; Array[a, 100, 2] (* Amiram Eldar, Jan 20 2021 *)
PROG
(PARI) a(n) = divisors(n^2)[3]; \\ Michel Marcus, Jan 20 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Charles Kusniec, Jan 20 2021
STATUS
approved