login
A370897
Partial alternating sums of the number of abelian groups sequence (A000688).
1
1, 0, 1, -1, 0, -1, 0, -3, -1, -2, -1, -3, -2, -3, -2, -7, -6, -8, -7, -9, -8, -9, -8, -11, -9, -10, -7, -9, -8, -9, -8, -15, -14, -15, -14, -18, -17, -18, -17, -20, -19, -20, -19, -21, -19, -20, -19, -24, -22, -24, -23, -25, -24, -27, -26, -29, -28, -29, -28
OFFSET
1,8
LINKS
László Tóth, Alternating Sums Concerning Multiplicative Arithmetic Functions, Journal of Integer Sequences, Vol. 20 (2017), Article 17.2.1.
FORMULA
a(n) = Sum_{k=1..n} (-1)^(k+1) * A000688(k).
a(n) = k_1 * A021002 * n + k_2 * A084892 * n^(1/2) + k_3 * A084893 * n^(1/3) + O(n^(1/4 + eps)), where eps > 0 is arbitrarily small, k_j = -1 + 2 * Product_{i>=1} (1 - 1/2^(i/j)), k_1 = 2*A048651 - 1 = -0.422423809826..., k_2 = -0.924973966404..., and k_3 = -0.991478298912... (Tóth, 2017).
MATHEMATICA
f[n_] := Times @@ (PartitionsP[Last[#]] & /@ FactorInteger[n]); f[1] = 1; Accumulate[Array[(-1)^(#+1) * f[#] &, 100]]
PROG
(PARI) f(n) = vecprod(apply(numbpart, factor(n)[, 2]));
lista(kmax) = {my(s = 0); for(k = 1, kmax, s += (-1)^(k+1) * f(k); print1(s, ", "))};
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Amiram Eldar, Mar 05 2024
STATUS
approved