OFFSET
1,2
COMMENTS
We need d and sigma(d) odd which happens precisely when d is an odd square.
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000
EXAMPLE
As 9 * sigma(9) = 9 * (1 + 3 + 9) = 9 * 13 = 117 is odd, 117 is in the sequence.
MATHEMATICA
Select[2Range[0, 9999] + 1, MemberQ[(DivisorSigma[1, #] * # &)/@Divisors[#], #] &] (* Alonso del Arte, Sep 18 2019 *)
PROG
(PARI) upto(n) = {my(res = List()); forstep(i = 1, sqrtnint(n, 4), 2, c = i^2*sigma(i^2); if(c <= n, listput(res, c))); listsort(res, 1); res}
CROSSREFS
KEYWORD
nonn
AUTHOR
David A. Corneth, Sep 18 2019
STATUS
approved