|
|
A160896
|
|
a(n) = ((2^b-1)/phi(n))*Sum_{d|n} Moebius(n/d)*d^(b-1) for b = 6.
|
|
1
|
|
|
63, 1953, 7623, 31248, 49203, 236313, 176463, 499968, 617463, 1525293, 1014615, 3781008, 1949283, 5470353, 5953563, 7999488, 5590683, 19141353, 8666343, 24404688, 21352023, 31453065, 18431343, 60496128, 30751875, 60427773, 50014503, 87525648, 46150083, 184560453
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
Sum_{k=1..n} a(k) ~ c * n^5, where c = (63/5) * Product_{p prime} (1 + (p^4-1)/((p-1)*p^5)) = 23.9347523175... .
Sum_{k>=1} 1/a(k) = (zeta(4)*zeta(5)/63) * Product_{p prime} (1 - 2/p^5 + 1/p^9) = 0.01658573169... . (End)
|
|
MATHEMATICA
|
f[p_, e_] := p^(4*e - 4)*(1 + p + p^2 + p^3 + p^4); a[1] = 63; a[n_] := 63 * Times @@ f @@@ FactorInteger[n]; Array[a, 25] (* Amiram Eldar, Nov 08 2022 *)
|
|
PROG
|
(PARI) a(n) = {my(f = factor(n)); 63 * prod(i = 1, #f~, (f[i, 1]^4 + f[i, 1]^3 + f[i, 1]^2 + f[i, 1] + 1)*f[i, 1]^(4*f[i, 2] - 4)); } \\ Amiram Eldar, Nov 08 2022
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|