|
%I #22 Dec 18 2023 01:42:35
%S 1,2,4,2,4,8,4,0,10,8,4,8,4,8,16,-5,4,20,4,8,16,8,4,0,10,8,20,8,4,32,
%T 4,-14,16,8,16,20,4,8,16,0,4,32,4,8,40,8,4,-20,10,20,16,8,4,40,16,0,
%U 16,8,4,32,4,8,40,-28,16,32,4,8,16,32,4,0,4,8,40,8,16,32,4
%N a(n) = Sum_{d|n} (-1)^(n/d+1) * Sum_{j|d} tau(j), where tau = number of divisors (A000005).
%H Antti Karttunen, <a href="/A318768/b318768.txt">Table of n, a(n) for n = 1..65537</a>
%F G.f.: Sum_{k>=1} tau_3(k)*x^k/(1 + x^k), where tau_3() = A007425.
%F L.g.f.: log(Product_{k>=1} (1 + x^k)^(tau_3(k)/k)) = Sum_{n>=1} a(n)*x^n/n.
%F Multiplicative with a(2^e) = 1 + (7-e^2)*e/6, and a(p^e) = binomial(e+3,3) for an odd prime p. - _Amiram Eldar_, Oct 25 2020
%F From _Amiram Eldar_, Dec 18 2023: (Start)
%F Dirichlet g.f.: zeta(s)^4 * (1 - 1/2^(s-1)).
%F Sum_{k=1..n} a(k) ~ (log(2)/2) * n * (log(n)^2 + (8 * gamma - log(2) - 2) * log(n) + 12 * gamma^2 - 8 * gamma + log(2) + 2 - 4 * gamma * log(2) + log(2)^2/3 - 8 * gamma_1), where gamma is Euler's constant (A001620) and gamma_1 is the first Stieltjes constant (A082633). (End)
%t Table[Sum[(-1)^(n/d + 1) Sum[DivisorSigma[0, j], {j, Divisors[d]}], {d, Divisors[n]}], {n, 79}]
%t nmax = 79; Rest[CoefficientList[Series[Sum[DivisorSum[k, DivisorSigma[0, #] &] x^k/(1 + x^k), {k, 1, nmax}], {x, 0, nmax}], x]]
%t nmax = 79; Rest[CoefficientList[Series[Log[Product[(1 + x^k)^(DivisorSum[k, DivisorSigma[0, #] &]/k), {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]]
%t f[p_, e_] := If[p == 2, 1 + (7 - e^2)*e/6, Binomial[e + 3, 3]]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* _Amiram Eldar_, Oct 25 2020 *)
%o (PARI) a(n) = sumdiv(n, d, (-1)^(n/d+1) * sumdiv(d, j, numdiv(j))); \\ _Michel Marcus_, Sep 04 2018
%Y Cf. A000005, A007425, A007426, A051062 (positions of 0's), A288571.
%Y Cf. A001620, A082633.
%K sign,mult,easy
%O 1,2
%A _Ilya Gutkovskiy_, Sep 03 2018
|
