OFFSET
0,3
COMMENTS
Multiplicative because A001160 is. - Andrew Howroyd, Jul 23 2018
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
FORMULA
a(n) = A145095(n)/504 for n > 0.
G.f.: phi_{6, 1}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}.
If p is a prime, a(p) = p^6 + p = A131472(p). - Seiichi Manyama, Feb 10 2017
a(n) = n*A001160(n) for n > 0. - Seiichi Manyama, Feb 18 2017
Sum_{k=1..n} a(k) ~ zeta(6) * n^7 / 7. - Amiram Eldar, Sep 06 2023
From Amiram Eldar, Oct 30 2023: (Start)
Multiplicative with a(p^e) = p^e * (p^(5*e+5)-1)/(p^5-1).
Dirichlet g.f.: zeta(s-1)*zeta(s-6). (End)
EXAMPLE
a(6) = 1^6*6^1 + 2^6*3^1 + 3^6*2^1 + 6^6*1^1 = 48312.
MATHEMATICA
terms = 30;
E2[x_] = 1 - 24*Sum[k*x^k/(1 - x^k), {k, 1, terms}];
E4[x_] = 1 + 240*Sum[k^3*x^k/(1 - x^k), {k, 1, terms}];
E6[x_] = 1 - 504*Sum[k^5*x^k/(1 - x^k), {k, 1, terms}];
(E4[x]^2 - E2[x]*E6[x])/1008 + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 26 2018 *)
Table[n * DivisorSigma[5, n], {n, 0, 30}] (* Amiram Eldar, Sep 06 2023 *)
PROG
(PARI) a(n) = if(n < 1, 0, n * sigma(n, 5)); \\ Andrew Howroyd, Jul 23 2018
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Seiichi Manyama, Feb 05 2017
STATUS
approved