A358040 -id:A358040 - cp's OEIS frontend

A366438 The number of divisors of the exponentially odd numbers (A268335).

1, 2, 2, 2, 4, 2, 4, 4, 2, 2, 4, 4, 2, 2, 4, 4, 2, 8, 4, 4, 2, 8, 2, 6, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 2, 4, 2, 8, 4, 8, 4, 4, 2, 2, 4, 4, 8, 2, 4, 8, 2, 2, 4, 4, 8, 2, 4, 2, 4, 4, 4, 8, 2, 4, 4, 4, 4, 12, 2, 2, 8, 2, 8, 8, 4, 2, 2, 8, 4, 2, 8, 4, 4, 4, 16, 4

Offset: 1

Amiram Eldar, Oct 10 2023

1 is the only odd term in this sequence.

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

Cf. A000005, A268335, A366439.

Similar sequences: A048691, A072048, A076400, A358040, A363194, A363195.

f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[AllTrue[e, OddQ], Times @@ (e + 1), Nothing]]; f[1] = 1; Array[f, 150]

lista(max) = for(k = 1, max, my(e = factor(k)[, 2], isexpodd = 1); for(i = 1, #e, if(!(e[i] % 2), isexpodd = 0; break)); if(isexpodd, print1(vecprod(apply(x -> x+1, e)), ", ")));

from math import prod
from itertools import count, islice
from sympy import factorint
def A366438_gen(): # generator of terms
    for n in count(1):
        f = factorint(n).values()
        if all(e&1 for e in f):
            yield prod(e+1 for e in f)
A366438_list = list(islice(A366438_gen(),30)) # Chai Wah Wu, Oct 10 2023

a(n) = A000005(A268335(n)).

A366440 The sum of divisors of the cubefree numbers (A004709).

Original entry on oeis.org

1, 3, 4, 7, 6, 12, 8, 13, 18, 12, 28, 14, 24, 24, 18, 39, 20, 42, 32, 36, 24, 31, 42, 56, 30, 72, 32, 48, 54, 48, 91, 38, 60, 56, 42, 96, 44, 84, 78, 72, 48, 57, 93, 72, 98, 54, 72, 80, 90, 60, 168, 62, 96, 104, 84, 144, 68, 126, 96, 144, 72, 74, 114, 124, 140

Offset: 1

Amiram Eldar, Oct 10 2023

Amiram Eldar, Table of n, a(n) for n = 1..10000
Rafael Jakimczuk and Matilde Lalín, Asymptotics of sums of divisor functions over sequences with restricted factorization structure, Notes on Number Theory and Discrete Mathematics, Vol. 28, No. 4 (2022), pp. 617-634; p. 619, eq. (2).

Cf. A000203, A004709, A358040.
Cf. A002117, A082020.
Similar sequences: A062822, A065764, A180114, A362986, A366439.

f[p_, e_] := (p^(e+1)-1)/(p-1); s[n_] := Module[{fct = FactorInteger[n]}, If[AllTrue[fct[[;;, 2]], # < 3 &], Times @@ f @@@ fct, Nothing]]; s[1] = 1; Array[s, 100]

lista(max) = for(k = 1, max, my(f = factor(k), iscubefree = 1); for(i = 1, #f~, if(f[i, 2] > 2, iscubefree = 0; break)); if(iscubefree, print1(sigma(f), ", ")));

from sympy import mobius, integer_nthroot, divisor_sigma
def A366440(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 divisor_sigma(m) # Chai Wah Wu, Aug 06 2024

a(n) = A000203(A004709(n)).
Sum_{k=1..n} a(k) ~ c * n^2, where c = 15*zeta(3)/(2*Pi^2) = A082020 * A002117 / 2 = 0.913453711751... .
The asymptotic mean of the abundancy index of the cubefree numbers: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A004709(k) = 15/Pi^2 = 1.519817... (A082020).

A366536 The number of unitary divisors of the cubefree numbers (A004709).

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 2, 4, 2, 4, 2, 4, 4, 2, 4, 2, 4, 4, 4, 2, 2, 4, 4, 2, 8, 2, 4, 4, 4, 4, 2, 4, 4, 2, 8, 2, 4, 4, 4, 2, 2, 4, 4, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 8, 2, 4, 4, 8, 2, 2, 4, 4, 4, 4, 8, 2, 4, 2, 8, 4, 4, 4, 2, 8, 4, 4, 4, 4, 4, 2, 4, 4, 4, 2, 8

Offset: 1

Amiram Eldar, Oct 12 2023

The number of unitary divisors of the squarefree numbers (A005117) is the same as the number of divisors of the squarefree numbers (A072048), because all the divisors of a squarefree number are unitary.

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

Cf. A004709, A005117, A034444, A072048, A077610, A358040, A366440, A366537.

Similar sequences: A366534, A366538.

f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[AllTrue[e, # < 3 &], 2^Length[e], Nothing]]; f[1] = 1; Array[f, 150]

lista(max) = for(k = 1, max, my(e = factor(k)[, 2], iscubefree = 1); for(i = 1, #e, if(e[i] > 2, iscubefree = 0; break)); if(iscubefree, print1(2^(#e), ", ")));

from sympy.ntheory.factor_ import udivisor_count
from sympy import mobius, integer_nthroot
def A366536(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_count(m) # Chai Wah Wu, Aug 05 2024

a(n) = A034444(A004709(n)).

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

Original entry on oeis.org

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

Amiram Eldar, Oct 12 2023

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

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

Similar sequences: A034676, A366535, A366539.

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]

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]), ", ")));

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

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... .

A358039 a(n) is the Euler totient function phi applied to the n-th cubefree number.

Original entry on oeis.org

1, 1, 2, 2, 4, 2, 6, 6, 4, 10, 4, 12, 6, 8, 16, 6, 18, 8, 12, 10, 22, 20, 12, 12, 28, 8, 30, 20, 16, 24, 12, 36, 18, 24, 40, 12, 42, 20, 24, 22, 46, 42, 20, 32, 24, 52, 40, 36, 28, 58, 16, 60, 30, 36, 48, 20, 66, 32, 44, 24, 70, 72, 36, 40, 36, 60, 24, 78, 40

Offset: 1

Amiram Eldar, Oct 29 2022

nonn

The analogous sequence with squarefree numbers is A049200.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Zhu Weiyi, On the cube free number sequences, Smarandache Notions J., Vol. 14 (2004), pp. 199-202.

Cf. A000010, A002117, A004709, A049200, A358040.

EulerPhi[Select[Range[100], Max[FactorInteger[#][[;; , 2]]] < 3 &]]

from sympy import mobius, integer_nthroot, totient
def A358039(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 totient(m) # Chai Wah Wu, Aug 06 2024

a(n) = A000010(A004709(n)).
Sum_{k=1..n} a(k) = (c/(2*zeta(3)))*n^2 + O(n^(3/2+eps)), where c = Product_{p prime} (1 - (p+1)/(p^3+p^2+1)) = 0.62583324412633345811... (Weiyi, 2004).
From Amiram Eldar, Oct 09 2023: (Start)
Sum_{n>=1} 1/(A004709(n)*a(n)) = Product_{p prime} (1 + (p^2+1)/((p-1)*p^3)) = 2.14437852780769816048... .
Sum_{n>=1} 1/a(n)^2 = Product_{p prime} (1 + (p^2+1)/((p-1)^2*p^2)) = 3.26032746607943673536... . (End)

A366441 The number of divisors of the 5-rough numbers (A007310).

Original entry on oeis.org

1, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 4, 2, 2, 2, 2, 3, 2, 4, 2, 2, 4, 2, 2, 2, 4, 2, 2, 4, 2, 4, 4, 2, 2, 2, 2, 2, 2, 4, 4, 3, 4, 2, 2, 4, 2, 2, 4, 4, 2, 2, 4, 2, 4, 2, 2, 3, 2, 6, 2, 2, 4, 4, 2, 2, 2, 2, 4, 4, 4, 2, 4, 4, 4, 2, 2, 2, 2, 4, 2, 2, 6, 4, 2, 4, 2, 4

Offset: 1

Amiram Eldar, Oct 10 2023

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

Cf. A000005, A001620, A007310, A366442.

Similar sequences: A048691, A072048, A076400, A099774, A358040, A363194, A363195.

a[n_] := DivisorSigma[0, 2*Floor[3*n/2] - 1]; Array[a, 100]

PARI
```
a(n) = numdiv((3*n)\2 << 1 - 1)
```

from sympy import divisor_count
def A366441(n): return divisor_count((n+(n>>1)<<1)-1) # Chai Wah Wu, Oct 10 2023

a(n) = A000005(A007310(n)).

Sum_{k=1..n} a(k) ~ (log(n) + 2*gamma - 1 + 2*log(6)) * n / 3, where gamma is Euler's constant (A001620).

cp's OEIS Frontend

A366438 The number of divisors of the exponentially odd numbers (A268335).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A366440 The sum of divisors of the cubefree numbers (A004709).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A366536 The number of unitary divisors of the cubefree numbers (A004709).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

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

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A358039 a(n) is the Euler totient function phi applied to the n-th cubefree number.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Python

Formula

A366441 The number of divisors of the 5-rough numbers (A007310).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula