OFFSET
1,8
COMMENTS
Let b(n, k) be the sequence that counts the divisors of n that are k-th powers of squarefree numbers. Then, b(n, 1) = A034444(n), b(n, 2) = A323308(n), b(n, 3) = A368248(n). b(n, k) is multiplicative with b(p^e, k) = 2 if e >= k, and 1 otherwise. The asymptotic mean of b(n, k) for k >= 2 is lim_{m->oo} (1/m) * Sum_{n=1..m} b(n, k) = zeta(k)/zeta(2*k). Since a(n) = Sum_{k>=2} (b(n, k) - 1), the formula for the asymptotic mean of this sequence follows (see the Formula section).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
a[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, 1 + Total[2^Accumulate[Count[e, #] & /@ Range[Max[e], 1, -1]] - 1] - 2^Length[e]]; a[1] = 0; Array[a, 100]
PROG
(PARI) a(n) = {my(f = factor(n), e, m, h, c); if(n == 1, 0, e = f[, 2]; m = vecmax(e); h = vector(m); for(i = 1, m, c = 0; for(j = 1, #e, if(e[j] == (m+1-i), c++)); h[i] = c); for(i = 2, m, h[i] += h[i-1]); for(i = 1, m, h[i] = 2^h[i]-1); 1 + vecsum(h) - 1<<#e); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Dec 19 2023
STATUS
approved