OFFSET
1,4
COMMENTS
Dirichlet convolution of A010051(n) and n^6. - Wesley Ivan Hurt, Jul 15 2025
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
FORMULA
a(A000040(n)) = 1.
From Wesley Ivan Hurt, Jul 15 2025: (Start)
a(n) = Sum_{d|n} c(d) * (n/d)^6, where c = A010051.
a(p^k) = p^(6*k-6) for p prime and k>=1. (End)
EXAMPLE
a(6) = 793; a(6) = 6^6 * Sum_{p|6, p prime} 1/p^6 = 46656 * (1/2^6 + 1/3^6) = 793.
MATHEMATICA
Array[#^6*DivisorSum[#, 1/#^6 &, PrimeQ] &, 50] (* Wesley Ivan Hurt, Jul 15 2025 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Feb 05 2022
STATUS
approved
