%I #25 Sep 16 2020 05:05:42
%S 1,-33,-244,32,-3126,8052,-16808,0,243,103158,-161052,-7808,-371294,
%T 554664,762744,0,-1419858,-8019,-2476100,-100032,4101152,5314716,
%U -6436344,0,3125,12252702,0,-537856,-20511150,-25170552,-28629152,0,39296688,46855314,52541808
%N Dirichlet inverse of A001160, sigma_5.
%H Amiram Eldar, <a href="/A178448/b178448.txt">Table of n, a(n) for n = 1..10000</a>
%F Dirichlet g.f.: 1/(zeta(s)*zeta(s-5)). - _R. J. Mathar_, Mar 10 2011
%F a(n) = Sum_{d|n} mu(n/d)*mu(d)*d^5. - _Ilya Gutkovskiy_, Nov 06 2018
%F Multiplicative with a(p) = -1 - p^5, a(p^2) = p^5, and a(p^e) = 0 for e>=3. - _Amiram Eldar_, Sep 16 2020
%t a[1] = 1; a[n_] := a[n] = -Sum[ DivisorSigma[5, n/d] a[d], {d, Most @ Divisors[n]}]; Table[a[n], {n, 1, 29}] (* _Jean-François Alcover_, Jun 24 2013 *)
%t f[p_, e_] := If[e == 1, -p^5 - 1, If[e == 2, p^5, 0]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Sep 16 2020 *)
%o (PARI) A178448_vec(len)={
%o a063524=vector(len) ; a063524[1] = 1 ;
%o a001160=direuler(p=2,len, 1/(1-p^5*X)/(1-X)) ;
%o dirdiv(a063524,a001160) ;}
%o { A178448_vec(70) } /* _R. J. Mathar_, Mar 10 2011 */
%o (PARI) a(n) = sumdiv(n, d, moebius(n/d)*moebius(d)*d^5); \\ _Michel Marcus_, Nov 06 2018
%o (PARI) for(n=1, 100, print1(direuler(p=2, n, (1 - X)*(1 - p^5*X))[n], ", ")) \\ _Vaclav Kotesovec_, Sep 16 2020
%Y Cf. A001160, A053825, A053826.
%K sign,mult
%O 1,2
%A _R. J. Mathar_, Dec 22 2010
%E More terms from _Amiram Eldar_, Sep 16 2020
|