OFFSET
1,2
FORMULA
Multiplicative with a(p^e) = p^(6*(e-1))*(p^2+1)*(p^4+1), e>0.
Dirichlet g.f.: zeta(s-6)*Product_{primes p} (1+p^(4-s)+p^(2-s)+p^(-s)).
Dirichlet convolution of A001014 with the multiplicative sequence 1, 21, 91, 0, 651, 1911, 2451, 0, 0, 13671, 14763, 0, 28731, 51471...
Sum_{k=1..n} a(k) ~ c * n^7 / 7, where c = Product_{primes p} (1 + 1/p^3 + 1/p^5 + 1/p^7) = 1.22847463998021088097249049512949441921891884186337179613337753... - Vaclav Kotesovec, Dec 18 2019
MATHEMATICA
f[p_, e_] := p^(6*(e - 1))*(p^2 + 1)*(p^4 + 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 01 2022 *)
CROSSREFS
KEYWORD
nonn,mult,easy
AUTHOR
R. J. Mathar, Aug 28 2011
STATUS
approved