OFFSET
1,1
COMMENTS
LINKS
MATHEMATICA
ad[1] = 0; ad[n_] := n * Total@(Last[#]/First[#]& /@ FactorInteger[n]); Select[Range[1, 4000, 2], (d = ad[DivisorSigma[1, #^2]] - #^2) > 0 && Divisible[d, 3] &] (* Amiram Eldar, Sep 18 2021 *)
PROG
CROSSREFS
Cf. also A342923.
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 18 2021
STATUS
approved