A366537 The sum of unitary divisors of the cubefree numbers (A004709).

1, 3, 4, 5, 6, 12, 8, 10, 18, 12, 20, 14, 24, 24, 18, 30, 20, 30, 32, 36, 24, 26, 42, 40, 30, 72, 32, 48, 54, 48, 50, 38, 60, 56, 42, 96, 44, 60, 60, 72, 48, 50, 78, 72, 70, 54, 72, 80, 90, 60, 120, 62, 96, 80, 84, 144, 68, 90, 96, 144, 72, 74, 114, 104, 100, 96

Offset

1,2

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Formula

a(n) = A034448(A004709(n)).

Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = zeta(3)^2 * Product_{p prime} (1 + 1/p^2 - 2/p^3 + 1/p^4 - 1/p^5) = 1.665430860774244601005... .

The asymptotic mean of the unitary abundancy index of the cubefree numbers: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A004709(k) = c / zeta(3) = 1.38548421160152785073... .

Mathematica

s[n_] := Module[{f = FactorInteger[n], e}, e = f[[;; , 2]]; If[AllTrue[e, # < 3 &], Times @@ (1 + Power @@@ f), Nothing]]; s[1] = 1; Array[s, 100]

Prog

(PARI) lista(max) = for(k = 1, max, my(f = factor(k), e = f[, 2], iscubefree = 1); for(i = 1, #e, if(e[i] > 2, iscubefree = 0; break)); if(iscubefree, print1(prod(i = 1, #e, 1 + f[i, 1]^e[i]), ", ")));
(Python)
from sympy.ntheory.factor_ import udivisor_sigma
from sympy import mobius, integer_nthroot
def A366537(n):
    def f(x): return n+x-sum(mobius(k)*(x//k**3) for k in range(1, integer_nthroot(x, 3)[0]+1))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return udivisor_sigma(m) # Chai Wah Wu, Aug 05 2024

Crossrefs

Cf. A002117, A004709, A005117, A034448, A062822, A077610, A358040, A366440, A366536.

Similar sequences: A034676, A366535, A366539.

Sequence in context: A069184 A181549 A366539 * A241405 A322485 A327668

Adjacent sequences: A366534 A366535 A366536 * A366538 A366539 A366540

Keyword

nonn,easy

Author

Amiram Eldar, Oct 12 2023

Status

approved