OFFSET
1,4
COMMENTS
a(n) = 1 when n = 1 or n is prime.
a(p^k) = (p^(3*k+3) - 1)/(p^3 - 1) - p^3 for p is prime.
LINKS
EXAMPLE
a(4) = 65 because 4 has 2 nonprime divisors {1,4} and 1^3 + 4^3 = 65.
MATHEMATICA
Table[DivisorSum[n, #1^3 & , !PrimeQ[#1] & ], {n, 55}]
Table[DivisorSigma[3, n] - DivisorSum[n, #1^3 & , PrimeQ[#1] & ], {n, 55}]
Table[Total[Select[Divisors[n], !PrimeQ[#]&]^3], {n, 60}] (* Harvey P. Dale, Aug 02 2024 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ilya Gutkovskiy, Dec 12 2016
STATUS
approved