A385042 -id:A385042 - cp's OEIS frontend

A385043 The sum of the unitary divisors of n whose exponents in their prime factorizations are all powers of 2 (A138302).

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

Offset: 1

Amiram Eldar, Jun 16 2025

The number of these divisors is A385042(n), and the largest of them is A367168(n).

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

The unitary analog of A353900.
Cf. A005117, A034448, A138302, A209229, A385042.
The sum of unitary divisors of n that are: A092261 (squarefree), A192066 (odd), A358346 (exponentially odd), A358347 (square), A360720 (powerful), A371242 (cubefree), A380396 (cube), A383763 (exponentially squarefree), this sequence (exponentially 2^n), A385045 (5-rough), A385046 (3-smooth), A385047 (power of 2), A385048 (cubefull), A385049 (biquadratefree).

f[p_, e_] := If[e == 2^IntegerExponent[e, 2], p^e + 1, 1]; a[ 1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i, 2] == 1<

Multiplicative with a(p^e) = p^(A209229(e)) + 1.
a(n) <= A034448(n), with equality if and only if n is in A138302.
a(n) <= A353900(n), with equality if and only if n is squarefree (A005117).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (1/(p*(p+1)) + Sum_{k>=2} (1/p^(2^k)-1/p^(2^k-1))) = 1.21427559551509410114... .

A385044 The number of unitary divisors of n that are 5-rough numbers (A007310).

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Jun 16 2025

The sum of these divisors is A385045(n), and the largest of them is A065330(n).

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

The unitary analog of A035218.
Cf. A007310, A034444, A065330, A385045.
Cf. A001620, A013661, A306016.
The number of unitary divisors of n that are: A000034 (power of 2), A055076 (exponentially odd), A056624 (square), A056671 (squarefree), A068068 (odd), A323308 (powerful), A365498 (cubefree), A365499 (biquadratefree), A368248 (cubefull), A380395 (cube), A382488 (3-smooth), A385042 (exponentially 2^n), this sequence (5-rough).

f[p_, e_] := If[p <= 3, 1, 2]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

a(n) = vecprod(apply(x -> if(x <= 3, 1, 2), factor(n)[, 1]));

Multiplicative with a(p^e) = 1 if p <= 3, and 2 if p >= 5.
a(n) = A034444(n)/A382488(n).
a(n) <= A034444(n), with equality if and only if n is 5-rough.
a(n) <= A035218(n).
Dirichlet g.f.: (zeta(s)^2/zeta(2*s)) * (1/((1+1/2^s)*(1+1/3^s))).
Sum_{k=1..n} a(k) ~ (n / (2 * zeta(2))) *(log(n) + 2*gamma - 1 + log(2)/3 + log(3)/4 - 2*zeta'(2)/zeta(2)), where gamma is Euler's constant (A001620).

cp's OEIS Frontend

A385043 The sum of the unitary divisors of n whose exponents in their prime factorizations are all powers of 2 (A138302).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula