|
|
A005064
|
|
Sum of cubes of primes dividing n.
|
|
22
|
|
|
0, 8, 27, 8, 125, 35, 343, 8, 27, 133, 1331, 35, 2197, 351, 152, 8, 4913, 35, 6859, 133, 370, 1339, 12167, 35, 125, 2205, 27, 351, 24389, 160, 29791, 8, 1358, 4921, 468, 35, 50653, 6867, 2224, 133, 68921, 378, 79507, 1339, 152, 12175, 103823, 35, 343, 133, 4940, 2205, 148877, 35, 1456, 351, 6886, 24397, 205379, 160
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Additive with a(p^e) = p^3.
G.f.: Sum_{k>=1} prime(k)^3*x^prime(k)/(1 - x^prime(k)). - Ilya Gutkovskiy, Dec 24 2016
(End)
|
|
MATHEMATICA
|
f[p_, e_] := p^3; a[n_] := Plus @@ f @@@ FactorInteger[n]; a[1] = 0; Array[a, 100] (* Amiram Eldar, Jun 20 2022 *)
|
|
PROG
|
(Python)
from sympy import primefactors
def a(n): return sum(p**3 for p in primefactors(n))
(PARI) a(n) = my(f=factor(n)); sum(k=1, #f~, f[k, 1]^3); \\ Michel Marcus, Jul 11 2017
|
|
CROSSREFS
|
Sum of the k-th powers of the primes dividing n for k=0..10 : A001221 (k=0), A008472 (k=1), A005063 (k=2), this sequence (k=3), A005065 (k=4), A351193 (k=5), A351194 (k=6), A351195 (k=7), this sequence (k=8), A351197 (k=9), A351198 (k=10).
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|