A092261 -id:A092261 - 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... .

A385045 The sum of the unitary divisors of n that are 5-rough numbers (A007310).

Original entry on oeis.org

1, 1, 1, 1, 6, 1, 8, 1, 1, 6, 12, 1, 14, 8, 6, 1, 18, 1, 20, 6, 8, 12, 24, 1, 26, 14, 1, 8, 30, 6, 32, 1, 12, 18, 48, 1, 38, 20, 14, 6, 42, 8, 44, 12, 6, 24, 48, 1, 50, 26, 18, 14, 54, 1, 72, 8, 20, 30, 60, 6, 62, 32, 8, 1, 84, 12, 68, 18, 24, 48, 72, 1, 74, 38

Offset: 1

Amiram Eldar, Jun 16 2025

First differs from A186099 at n = 25; a(25) = 26, while A186099(25) = 31.

The number of these divisors is A385044(n), and the largest of them is A065330(n).

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

The unitary analog of A186099.
Cf. A002117, A007310, A034448, A065330, A385044.
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), A385043 (exponentially 2^n), this sequence (5-rough), A385046 (3-smooth), A385047 (power of 2), A385048 (cubefull), A385049 (biquadratefree).

f[p_, e_] := If[p <= 3, 1, p^e + 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,1] < 5, 1, f[i, 1]^f[i, 2] + 1));}

Multiplicative with a(p^e) = 1 if p <= 3, and p^e + 1 if p >= 5.
a(n) = A034448(n)/A385046(n).
a(n) <= A034448(n), with equality if and only if n is 5-rough.
a(n) <= A186099(n).
Dirichlet g.f.: (zeta(s)*zeta(s-1)/zeta(2*s-1)) * ((1-1/2^(s-1))/(1-1/2^(2*s-1))) * ((1-1/3^(s-1))/(1-1/3^(2*s-1))).
Sum_{k=1..n} a(k) ~ c * n^2, where c = 3*Pi^2/(91*zeta(3)) = 0.270679... .

A385046 The sum of the unitary divisors of n that are 3-smooth numbers (A003586).

Original entry on oeis.org

1, 3, 4, 5, 1, 12, 1, 9, 10, 3, 1, 20, 1, 3, 4, 17, 1, 30, 1, 5, 4, 3, 1, 36, 1, 3, 28, 5, 1, 12, 1, 33, 4, 3, 1, 50, 1, 3, 4, 9, 1, 12, 1, 5, 10, 3, 1, 68, 1, 3, 4, 5, 1, 84, 1, 9, 4, 3, 1, 20, 1, 3, 10, 65, 1, 12, 1, 5, 4, 3, 1, 90, 1, 3, 4, 5, 1, 12, 1, 17

Offset: 1

Amiram Eldar, Jun 16 2025

The number of these divisors is A382488(n), and the largest of them is A065331(n).

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

The unitary analog of A072079.
Cf. A001620, A003586, A034448, A065331, A082633, A382488.
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), A385043 (exponentially 2^n), A385045 (5-rough), this sequence (3-smooth), A385047 (power of 2), A385048 (cubefull), A385049 (biquadratefree).

f[n_, p_] := If[Divisible[n, p], p^IntegerExponent[n, p] + 1, 1]; a[n_] := f[n, 2]*f[n, 3]; Array[a, 100]

a(n) = if(n%2, 1, 2^valuation(n, 2)+1) * if(!(n%3), 3^valuation(n, 3)+1, 1);

Multiplicative with a(p^e) = p^e + 1 if p <= 3, and 1 if p >= 5.
a(n) = A034448(n)/A385045(n).
a(n) <= A034448(n), with equality if and only if n 3-smooth.
a(n) <= A072079(n).
Dirichlet g.f.: zeta(s) * ((1-1/2^(2*s-1))/(1-1/2^(s-1))) * ((1-1/3^(2*s-1))/(1-1/3^(s-1))).
Sum_{k=1..n} a(k) ~ (n/(6*log(2)*log(3))) * (log(n)^2 + c1*log(n) + c2), where c1 = 2*gamma - 2 + 7*log(2) + 5*log(3) - 2*log(6) = 5.916004..., c2 = 2 - 5*log(2) - 11*log(2)^2/6 - 3*log(3) - 5*log(3)^2/6 + 15*log(2)*log(3)/2 + (5*log(2) + 3*log(3) - 2)*gamma - 2*gamma_1 = 1.957142..., gamma is Euler's constant (A001620), and gamma_1 is the 1st Stieltjes constant (A082633).

A385047 The sum of the unitary divisors of n that are powers of 2.

Original entry on oeis.org

1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 17, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 33, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 17, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 65, 1, 3, 1, 5, 1, 3, 1, 9, 1, 3, 1, 5, 1, 3, 1, 17, 1, 3, 1, 5, 1, 3

Offset: 1

Amiram Eldar, Jun 16 2025

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

The unitary analog of A038712.
Cf. A001620, A006519, A034448.
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), A385043 (exponentially 2^n), A385045 (5-rough), A385046 (3-smooth), this sequence (power of 2), A385048 (cubefull), A385049 (biquadratefree).

a[n_] := If[OddQ[n], 1, 2^IntegerExponent[n, 2] + 1]; Array[a, 100]

PARI
```
a(n) = if(n%2, 1, 2^valuation(n, 2)+1);
```

Multiplicative with a(2^e) = 2^e + 1, and a(p^e) = 1 for an odd prime p.
a(n) = A034448(n) / A192066(n).
a(n) = A059841(n) + A006519(n), i.e., a(n) = A006519(n) + 1 if n is even, and 1 is n is odd.
Dirichlet g.f.: zeta(s) * ((1-1/2^(2*s-1))/(1-1/2^(s-1))).
Sum_{k=1..n} a(k) ~ (n/(2*log(2))) * (log(n) + gamma - 1 + 5*log(2)/2), where gamma is Euler's constant (A001620).

A385048 The sum of the unitary divisors of n that are cubefull numbers (A036966).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 17, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 28, 1, 1, 1, 1, 33, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 17, 1, 1, 1, 1, 1, 28, 1, 9, 1, 1, 1, 1, 1, 1, 1, 65, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 17, 82, 1, 1, 1, 1

Offset: 1

Amiram Eldar, Jun 16 2025

The number of these divisors is A368248(n), and the largest of them is A360540(n).

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

The unitary analog of A385005.
Cf. A034448, A036966, A046100, A360540, A368248.
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), A385043 (exponentially 2^n), A385045 (5-rough), A385046 (3-smooth), A385047 (power of 2), this sequence (cubefull), A385049 (biquadratefree).

f[p_, e_] := If[e <= 2, 1, p^e + 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] < 3, 1, f[i, 1]^f[i, 2] + 1));}

Multiplicative with a(p^e) = 1 if e <= 2, and a(p^e) = p^e + 1 if e >= 3.
a(n) = A034448(n) / A371242(n).
a(n) <= A034448(n), with equality if and only if n is cubefull (A036966).
a(n) <= A385005(n), with equality if and only if n is biquadratefree (A046100).
Dirichlet g.f.: zeta(s)*zeta(s-1)*Product_{p prime} (1 - 1/p^(s-1) + 1/p^(3*s-3) - 1/p^(4*s-3)).

A385049 The sum of the unitary divisors of n that are biquadratefree numbers (A046100).

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Jun 16 2025

First differs from A383763 at n = 32.

The number of these divisors is A365499(n), and the largest of them is A385007(n).

D. Suryanarayana, The number and sum of k-free integers <= x which are prime to n, Indian J. Math., Vol. 11 (1969), pp. 131-139.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Francesco Pappalardi, A survey on k-freeness, Number Theory, Ramanujan Math. Soc. Lect. Notes Ser., Vol. 1 (2003), pp. 71-88.

The unitary analog of A385006.
Cf. A036967, A046100, A365499, A385007.
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), A385043 (exponentially 2^n), A385045 (5-rough), A385046 (3-smooth), A385047 (power of 2), A385048 (cubefull), this sequence (biquadratefree).

f[p_, e_] := If[e < 4, 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] < 4, f[i, 1]^f[i, 2] + 1, 1)); }

Multiplicative with a(p^e) = p^e + 1 for e <= 3, and a(p^e) = 1 for e >= 4.
a(n) = 1 if and only if n is 4-full (A036967).
a(n) <= A034448(n), with equality if and only if n is biquadratefree.
Dirichlet g.f.: zeta(s) * Product_{p prime} (1 + 1/p^(s-1) + 1/p^(2*s-2) - 1/p^(2*s-1) + 1/p^(3*s-3) - 1/p^(3*s-2) - 1/p^(4*s-3)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (1 + 1/(p^2 + p) - 1/p^4) = 1.27769267395905900191... .

A366148 The sum of divisors of the cubefree part of n (A360539).

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Oct 01 2023

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

Cf. A000203, A004709, A036966, A092261, A360539, A366077, A366078, A366146, A366147.

f[p_, e_] := If[e < 3, (p^(e+1)-1)/(p-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] < 3, (f[i,1]^(f[i,2]+1)-1)/(f[i,1]-1), 1));}

a(n) = A000203(A360539(n)).
a(n) = A000203(n)/A366146(n).
Multiplicative with a(p^e) = (p^(e+1)-1)/(p-1) if e <= 2, and 1 otherwise.
a(n) >= 1, with equality if and only if n is cubefull (A036966).
a(n) <= A000203(n), with equality if and only if n is cubefree (A004709).
Dirichlet g.f.: zeta(s) * Product_{p prime} (1 + 1/p^(s-1) + 1/p^(2*s-2) - 1/p^(3*s-2) - 1/p^(3*s-1)). CHECK
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} (p/(p+1) + 1/p - 1/p^4) = 0.63884633697952950095... .

A286875 If n = Product (p_j^k_j) then a(n) = Sum (k_j >= 2, p_j^k_j).

Original entry on oeis.org

0, 0, 0, 4, 0, 0, 0, 8, 9, 0, 0, 4, 0, 0, 0, 16, 0, 9, 0, 4, 0, 0, 0, 8, 25, 0, 27, 4, 0, 0, 0, 32, 0, 0, 0, 13, 0, 0, 0, 8, 0, 0, 0, 4, 9, 0, 0, 16, 49, 25, 0, 4, 0, 27, 0, 8, 0, 0, 0, 4, 0, 0, 9, 64, 0, 0, 0, 4, 0, 0, 0, 17, 0, 0, 25, 4, 0, 0, 0, 16, 81, 0, 0, 4, 0, 0, 0, 8, 0, 9, 0, 4, 0, 0, 0, 32, 0, 49, 9, 29, 0, 0, 0, 8, 0, 0, 0, 31

Offset: 1

Ilya Gutkovskiy, Aug 02 2017

Sum of unitary, proper prime power divisors of n.

			a(360) = a(2^3*3^2*5) = 2^3 + 3^2 = 17.

Antti Karttunen, Table of n, a(n) for n = 1..16384
Index entries for sequences related to sums of divisors.

Cf. A005117, A008475, A222416, A023888, A023889, A034448, A063956, A077610, A092261, A246547, A284117.

Table[DivisorSum[n, # &, CoprimeQ[#, n/#] && PrimePowerQ[#] && !PrimeQ[#] &], {n, 108}]
f[p_, e_] := If[e == 1, 0, p^e]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jul 24 2024 *)

A286875(n) = { my(f=factor(n)); for (i=1, #f~, if(f[i, 2] < 2, f[i, 1] = 0)); vecsum(vector(#f~,i,f[i,1]^f[i,2])); }; \\ Antti Karttunen, Oct 07 2017

from sympy import primefactors, isprime, gcd, divisors
def a(n): return sum(d for d in divisors(n) if gcd(d, n//d)==1 and len(primefactors(d))==1 and not isprime(d))
print([a(n) for n in range(1, 109)]) # Indranil Ghosh, Aug 02 2017

a(n) = Sum_{d|n, d = p^k, p prime, k >= 2, gcd(d, n/d) = 1} d.
a(A246547(k)) = A246547(k).
a(A005117(k)) = 0.
Additive with a(p^e) = p^e if e >= 2, and 0 otherwise. - Amiram Eldar, Jul 24 2024

A343442 If n = Product (p_j^k_j) then a(n) = Product (p_j + 2), with a(1) = 1.

Original entry on oeis.org

1, 4, 5, 4, 7, 20, 9, 4, 5, 28, 13, 20, 15, 36, 35, 4, 19, 20, 21, 28, 45, 52, 25, 20, 7, 60, 5, 36, 31, 140, 33, 4, 65, 76, 63, 20, 39, 84, 75, 28, 43, 180, 45, 52, 35, 100, 49, 20, 9, 28, 95, 60, 55, 20, 91, 36, 105, 124, 61, 140, 63, 132, 45, 4, 105, 260, 69, 76, 125, 252

Offset: 1

Ilya Gutkovskiy, Apr 15 2021

Cf. A000203, A001615, A007429, A007947, A008683, A048250, A062953, A063441, A074816, A092261, A107758, A107759, A319132, A335005.

a[1] = 1; a[n_] := Times @@ ((#[[1]] + 2) & /@ FactorInteger[n]); Table[a[n], {n, 70}]
nmax = 70; CoefficientList[Series[Sum[MoebiusMu[k]^2 DivisorSigma[1, k] x^k/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x] // Rest

a(n) = sumdiv(n, d, moebius(d)^2 * sigma(d)) \\ Andrew Howroyd, Apr 15 2021

G.f.: Sum_{k>=1} mu(k)^2 * sigma(k) * x^k / (1 - x^k), where mu = A008683 and sigma = A000203.
a(n) = Sum_{d|n} mu(d)^2 * sigma(d).
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/(12*zeta(3)) = 0.684216... (A335005). - Amiram Eldar, Nov 13 2022
a(n) = Sum_{d|n} mu(d)^2*psi(d), where psi is A001615. - Ridouane Oudra, Jul 24 2025

A380161 a(n) is the value of the Euler totient function when applied to the powerfree part of n.

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Jan 13 2025

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

Cf. A000010, A005117, A013661, A055231, A056671, A092261, A335851, A380160.

f[p_, e_] := If[e == 1, p-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, f[i, 1]-1, 1));}

a(n) = A000010(A055231(n)).
a(n) >= 1, with equality if and only if n is in A335851.
a(n) <= A000010(n), with equality if and only if n is squarefree (A005117).
Multiplicative with a(p) = p-1, and a(p^e) = 1 if e >= 2.
Dirichlet g.f.: zeta(s) * Product_{p prime} (1 - 2/p^s + 1/p^(s-1) - 1/p^(2*s-1) + 2/p^(2*s)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = zeta(2) * Product_{p prime} (1 - 2/p^5 + 3/p^4 + 1/p^3 - 3/p^2) = 0.46288631864329056778... .

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

A385045 The sum of the unitary divisors of n that are 5-rough numbers (A007310).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A385046 The sum of the unitary divisors of n that are 3-smooth numbers (A003586).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A385047 The sum of the unitary divisors of n that are powers of 2.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A385048 The sum of the unitary divisors of n that are cubefull numbers (A036966).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A385049 The sum of the unitary divisors of n that are biquadratefree numbers (A046100).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A366148 The sum of divisors of the cubefree part of n (A360539).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A286875 If n = Product (p_j^k_j) then a(n) = Sum (k_j >= 2, p_j^k_j).

Views

Author