OFFSET
1,1
COMMENTS
This is not a subsequence of A228058. The terms that do not occur there: 441, 1521, 2401, 3249, 8649, 16641, 28561, 35721, etc., seem all to be squares. Terms of A228058 missing from this sequence are: 117, 245, 333, 425, 549, 605, 637, 657, 725, etc. (See A351574.) Terms that are also in A228058 are given in A386429.
LINKS
MATHEMATICA
ad[1] = 0; ad[n_] := n * Total@(Last[#]/First[#]& /@ FactorInteger[n]); f[n_] := ad[DivisorSigma[1, n]] - n; Select[Range[1, 5725, 2], CompositeQ[#] && EvenQ[f[#]] && Divisible[f[2*#], 3] &] (* Amiram Eldar, Sep 18 2021 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 18 2021
STATUS
approved
