A368251 The number of nonsquarefree divisors of n that are powers of squarefree numbers (A072777).

0, 0, 0, 1, 0, 0, 0, 2, 1, 0, 0, 1, 0, 0, 0, 3, 0, 1, 0, 1, 0, 0, 0, 2, 1, 0, 2, 1, 0, 0, 0, 4, 0, 0, 0, 3, 0, 0, 0, 2, 0, 0, 0, 1, 1, 0, 0, 3, 1, 1, 0, 1, 0, 2, 0, 2, 0, 0, 0, 1, 0, 0, 1, 5, 0, 0, 0, 1, 0, 0, 0, 4, 0, 0, 1, 1, 0, 0, 0, 3, 3, 0, 0, 1, 0, 0, 0

Offset: 1

Author

Amiram Eldar, Dec 19 2023

Keywords

• nonn

Comments

First differs from A046660 and A066301 at n = 36, and from A183094 at n = 72.

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

Crossrefs

Cf. A005117, A048105, A072777, A327527, A368250.

Cf. A034444, A323308, A368248.

Cf. A046660, A066301, A183094.

Programs

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]

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);}

Formula

a(n) = A327527(n) - A034444(n).

a(n) = 0 if and only if n is squarefree (A005117).

Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{k>=2} (zeta(k)/zeta(2*k) - 1) = 0.848633... (A368250).