OFFSET
1,2
COMMENTS
Numbers of the form p^(m*(4m-1)) and p^(m*(4m+1)) are terms of the sequence, where p is prime. p^2*q are terms of the sequence, where p and q are prime and p^2 > q > p.
Complement of A072497 in the positive integers. - Robert Israel, Dec 10 2024
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local s, F, d, p;
s:= n^2;
F:= sort(convert(numtheory:-divisors(s), list));
p:= 1:
for d in F do
p:= p*d;
if p > s then return false
elif p = s then return true
fi
od;
end proc:
select(filter, [$1..5000]); # Robert Israel, Dec 10 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Baltic, Aug 04 2002
EXTENSIONS
More terms from Sean A. Irvine, Sep 23 2024
STATUS
approved
