|
|
A321810
|
|
Sum of 6th powers of odd divisors of n.
|
|
11
|
|
|
1, 1, 730, 1, 15626, 730, 117650, 1, 532171, 15626, 1771562, 730, 4826810, 117650, 11406980, 1, 24137570, 532171, 47045882, 15626, 85884500, 1771562, 148035890, 730, 244156251, 4826810, 387952660, 117650, 594823322, 11406980, 887503682
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
LINKS
|
|
|
FORMULA
|
G.f.: Sum_{k>=1} (2*k - 1)^6*x^(2*k-1)/(1 - x^(2*k-1)). - Ilya Gutkovskiy, Dec 22 2018
Multiplicative with a(2^e) = 1 and a(p^e) = (p^(6*e+6)-1)/(p^6-1) for p > 2.
Sum_{k=1..n} a(k) ~ c * n^7, where c = zeta(7)/14 = 0.0720249... . (End)
a(n) + a(n/2)*2^6 = A013954(n) where a(.)=0 for non-integer arguments. - R. J. Mathar, Aug 15 2023
|
|
MATHEMATICA
|
f[2, e_] := 1; f[p_, e_] := (p^(6*e + 6) - 1)/(p^6 - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 02 2022 *)
|
|
PROG
|
(Python)
from sympy import divisor_sigma
def A321810(n): return int(divisor_sigma(n>>(~n&n-1).bit_length(), 6)) # Chai Wah Wu, Jul 16 2022
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,mult
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|