cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-7 of 7 results.

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

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Jun 16 2025

Keywords

Comments

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

Links

Crossrefs

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

Programs

  • Mathematica
    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]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i, 2] == 1<

Formula

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

Views

Author

Amiram Eldar, Jun 16 2025

Keywords

Comments

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

Links

Crossrefs

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

Programs

  • Mathematica
    f[p_, e_] := If[p <= 3, 1, p^e + 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i,1] < 5, 1, f[i, 1]^f[i, 2] + 1));}

Formula

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

Views

Author

Amiram Eldar, Jun 16 2025

Keywords

Comments

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

Links

Crossrefs

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

Programs

  • Mathematica
    f[n_, p_] := If[Divisible[n, p], p^IntegerExponent[n, p] + 1, 1]; a[n_] := f[n, 2]*f[n, 3]; Array[a, 100]
  • PARI
    a(n) = if(n%2, 1, 2^valuation(n, 2)+1) * if(!(n%3), 3^valuation(n, 3)+1, 1);

Formula

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

Views

Author

Amiram Eldar, Jun 16 2025

Keywords

Links

Crossrefs

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

Programs

  • Mathematica
    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);

Formula

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

Views

Author

Amiram Eldar, Jun 16 2025

Keywords

Comments

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

Links

Crossrefs

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

Programs

  • Mathematica
    f[p_, e_] := If[e <= 2, 1, p^e + 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i, 2] < 3, 1, f[i, 1]^f[i, 2] + 1));}

Formula

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

Views

Author

Amiram Eldar, Jun 16 2025

Keywords

Comments

First differs from A383763 at n = 32.
The number of these divisors is A365499(n), and the largest of them is A385007(n).

References

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

Links

Crossrefs

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

Programs

  • Mathematica
    f[p_, e_] := If[e < 4, p^e + 1, 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i, 2] < 4, f[i, 1]^f[i, 2] + 1, 1)); }

Formula

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

A380090 The sum of the unitary divisors of n that are terms in A207481.

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Jan 12 2025

Keywords

Comments

First differs from A371242 at n = 27.

Links

Crossrefs

Cf. A005117, A054743, A077610, A185359, A207481, A371242, A377520, A380087, A380088, A380089.

Programs

  • Mathematica
    f[p_, e_] := If[e <= p, p^e, 0] + 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i,2] <= f[i,1], f[i,1]^f[i,2], 0) + 1);}

Formula

a(n) = A034448(A380088(n)).
Multiplicative with a(p^e) = p^e + 1 if e <= p, and 1 otherwise.
a(n) = 1 if and only if n is in A054743.
a(n) < A034448(n) if and only if n is in A185359.
a(n) = A034448(n) if and only if n is in A207481.
a(n) = A377520(n) if and only if n is squarefree (A005117).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (p^(p+2) + p^(p+1) + p^p - p - 1)/(p^(p+1) * (p+1)) = 1.2078161... .
Showing 1-7 of 7 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.