A360539 -id:A360539 - cp's OEIS frontend

A366077 The number of prime factors of the cubefree part of n (A360539), counted with multiplicity.

0, 1, 1, 2, 1, 2, 1, 0, 2, 2, 1, 3, 1, 2, 2, 0, 1, 3, 1, 3, 2, 2, 1, 1, 2, 2, 0, 3, 1, 3, 1, 0, 2, 2, 2, 4, 1, 2, 2, 1, 1, 3, 1, 3, 3, 2, 1, 1, 2, 3, 2, 3, 1, 1, 2, 1, 2, 2, 1, 4, 1, 2, 3, 0, 2, 3, 1, 3, 2, 3, 1, 2, 1, 2, 3, 3, 2, 3, 1, 1, 0, 2, 1, 4, 2, 2, 2

Offset: 1

Amiram Eldar, Sep 28 2023

The sum of exponents smaller than 3 in the prime factorization of n.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Rafael Jakimczuk and Matilde Lalín, Sums of omega(n) and Omega(n) over the k-free parts and k-full parts of some particular sequences, Integers, Vol. 22 (2022), Article #A113.

Cf. A001222, A004709, A005117, A036966, A056169, A077761, A085541, A085548, A360539, A366076, A366078.

f[p_, e_] := If[e < 3, e, 0]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]

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

a(n) = A001222(A360539(n)).
a(n) = A001222(n) - A366076(n).
Additive with a(p^e) = e if e <= 2, and a(p^e) = 0 for e >= 3.
a(n) >= 0, with equality if and only if n is cubefull (A036966).
a(n) <= A001222(n), with equality if and only if n is cubefree (A004709).
a(n) >= A366078(n), with equality if and only if n is squarefree (A005117).
Sum_{k=1..m} a(k) ~ n * (log(log(n)) + B + C), where B is Mertens's constant (A077761) and C = Sum_{p prime} (p-2)/p^3 = A085548 - 2 * A085541 = 0.10272214144217842566... .

A366078 The number of distinct prime factors of the cubefree part of n (A360539).

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 0, 1, 2, 1, 2, 1, 2, 2, 0, 1, 2, 1, 2, 2, 2, 1, 1, 1, 2, 0, 2, 1, 3, 1, 0, 2, 2, 2, 2, 1, 2, 2, 1, 1, 3, 1, 2, 2, 2, 1, 1, 1, 2, 2, 2, 1, 1, 2, 1, 2, 2, 1, 3, 1, 2, 2, 0, 2, 3, 1, 2, 2, 3, 1, 1, 1, 2, 2, 2, 2, 3, 1, 1, 0, 2, 1, 3, 2, 2, 2

Offset: 1

Amiram Eldar, Sep 28 2023

The number of exponents smaller than 3 in the prime factorization of n.

The number of prime factors of the cubefree part of n (A360539), counted with multiplicity is A366077(n).

Amiram Eldar, Table of n, a(n) for n = 1..10000
Rafael Jakimczuk and Matilde Lalín, Sums of omega(n) and Omega(n) over the k-free parts and k-full parts of some particular sequences, Integers, Vol. 22 (2022), Article #A113.

Cf. A001221, A004709, A005117, A036966, A077761, A085541, A295659, A360539, A366077.

f[p_, e_] := If[e < 3, 1, 0]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]

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

a(n) = A001221(A360539(n)).
a(n) = A001221(n) - A295659(n).
Additive with a(p^e) = 1 if e <= 2, and a(p^e) = 0 for e >= 3.
a(n) >= 0, with equality if and only if n is cubefull (A036966).
a(n) <= A001221(n), with equality if and only if n is cubefree (A004709).
a(n) <= A366077(n), with equality if and only if n is squarefree (A005117).
Sum_{k=1..m} a(k) = n * (log(log(n)) + B - C) + O(n/log(n)), where B is Mertens's constant (A077761) and C = Sum_{p prime} 1/p^3 = 0.174762... (A085541).

A366147 The number of divisors of the cubefree part of n (A360539).

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 1, 3, 4, 2, 6, 2, 4, 4, 1, 2, 6, 2, 6, 4, 4, 2, 2, 3, 4, 1, 6, 2, 8, 2, 1, 4, 4, 4, 9, 2, 4, 4, 2, 2, 8, 2, 6, 6, 4, 2, 2, 3, 6, 4, 6, 2, 2, 4, 2, 4, 4, 2, 12, 2, 4, 6, 1, 4, 8, 2, 6, 4, 8, 2, 3, 2, 4, 6, 6, 4, 8, 2, 2, 1, 4, 2, 12, 4, 4, 4

Offset: 1

Amiram Eldar, Oct 01 2023

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

Cf. A000005, A004709, A036966, A056671, A360539, A366077, A366078, A366145, A366148.

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

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

a(n) = A000005(A360539(n)).
a(n) = A000005(n)/A366145(n).
Multiplicative with a(p^e) = e+1 if e <= 2, and 1 otherwise.
a(n) >= 1, with equality if and only if n is cubefull (A036966).
a(n) <= A000005(n), with equality if and only if n is cubefree (A004709).
Dirichlet g.f.: zeta(s) * Product_{p prime} (1 + 1/p^s + 1/p^(2*s) - 2/p^(3*s)).
From Vaclav Kotesovec, Oct 01 2023: (Start)
Dirichlet g.f.: zeta(s)^2 * Product_{p prime} (1 - 3/p^(3*s) + 2/p^(4*s)).
Let f(s) = Product_{p prime} (1 - 3/p^(3*s) + 2/p^(4*s)).
Sum_{k=1..n} a(k) ~ f(1) * n * (log(n) + 2*gamma - 1 + f'(1)/f(1)), where
f(1) = Product_{p prime} (1 - 3/p^3 + 2/p^4) = 0.66219033176371496870504912254207846719824904470940603905284774924086...,
f'(1) = f(1) * Sum_{p prime} (9*p - 8) * log(p) / (p^4 - 3*p + 2) = f(1) * 1.04316863044761953555286128194165251303791613504188623828521117799260...
and gamma is the Euler-Mascheroni constant A001620. (End)

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

A384049 The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is cubefree.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 7, 9, 10, 11, 12, 13, 14, 15, 15, 17, 18, 19, 20, 21, 22, 23, 21, 25, 26, 26, 28, 29, 30, 31, 31, 33, 34, 35, 36, 37, 38, 39, 35, 41, 42, 43, 44, 45, 46, 47, 45, 49, 50, 51, 52, 53, 52, 55, 49, 57, 58, 59, 60, 61, 62, 63, 63, 65, 66, 67, 68

Offset: 1

Amiram Eldar, May 18 2025

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

Unitary analog of A254926.
Cf. A004709, A065468, A360539, A360540, A384046, A384047.
The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is: A047994 (1), A384048 (squarefree), this sequence (cubefree), A384050 (powerful), A384051 (cubefull), A384052 (square), A384053 (cube), A384054 (exponentially odd), A384055 (odd), A384056 (power of 2), A384057 (3-smooth), A384058 (5-rough).

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

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

Multiplicative with a(p^e) = p^e if e <= 2, and p^e - 1 if e >= 3.
a(n) = n * A047994(n) / A384051(n).
a(n) = A047994(A360540(n)) * A360539(n).
Dirichlet g.f.: zeta(s-1) * zeta(s) * Product_{p prime} (1 - 1/p^s - 1/p^(3*s) + 1/p^(4*s-1)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (1 - 1/(p^5*(p+1))) = 0.988504... (A065468).

A360540 a(n) is the cubefull part of n: the largest divisor of n that is a cubefull number (A036966).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 16, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 27, 1, 1, 1, 1, 32, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 16, 1, 1, 1, 1, 1, 27, 1, 8, 1, 1, 1, 1, 1, 1, 1, 64, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 16, 81, 1, 1, 1, 1

Offset: 1

Amiram Eldar, Feb 11 2023

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

Cf. A004709, A036966, A057521, A360539.

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

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

a(n) = 1 if and only if n is a cubefree number (A004709).
a(n) = n if and only if n is a cubefull number (A036966).
a(n) <= A057521(n) with equality if and only if n is in A337050.
a(n) = n/A360539(n).
Multiplicative with a(p^e) = p^e if e >= 3, and 1 otherwise.
Dirichlet g.f.: zeta(s-1) * Product_{p prime} (1 - p^(1-s) + p^(-s) - p^(1-3*s) - p^(1-2*s) + p^(-2*s) + p^(3-3*s)).

A384051 The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is a cubefull number.

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, May 18 2025

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

Unitary analog of A384040.
Cf. A036966, A360539, A360540, A384046, A384047.
The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is: A047994 (1), A384048 (squarefree), A384049 (cubefree), A384050 (powerful), this sequence (cubefull), A384052 (square), A384053 (cube), A384054 (exponentially odd), A384055 (odd), A384056 (power of 2), A384057 (3-smooth), A384058 (5-rough).

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

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

Multiplicative with a(p^e) = p^e-1 if e <= 2, and p^e if e >= 3.
a(n) = n * A047994(n) / A384049(n).
a(n) = A047994(A360539(n)) * A360540(n).
Dirichlet g.f.: zeta(s-1) * Product_{p prime} (1 - 1/p^s - 1/p^(2*s) + 1/p^(2*s-1) + 1/p^(3*s-1)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (1 - 1/p^2 + 1/p^3 - 1/p^4 + 1/p^5) = 0.714093594477970831206... .

A367168 The largest unitary divisor of n that is a term of A138302.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 1, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 3, 25, 26, 1, 28, 29, 30, 31, 1, 33, 34, 35, 36, 37, 38, 39, 5, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 2, 55, 7, 57, 58, 59, 60, 61, 62, 63, 1, 65, 66, 67, 68, 69, 70

Offset: 1

Amiram Eldar, Nov 07 2023

First differs from A056192 at n = 32 and from A270418 at n = 128.

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

Cf. A048298, A065465, A138302, A367169, A367170, A367171.
Cf. A056192, A270418.
Similar sequences: A350388, A350389, A360539.

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

from math import prod
from sympy import factorint
def A367168(n): return prod(p**e for p,e in factorint(n).items() if not(e&-e)^e) # Chai Wah Wu, Nov 10 2023

Multiplicative with a(p^e) = p^A048298(e).
a(n) <= n, with equality if and only if n is in A138302.
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = 0.881513... (A065465).

A382903 The largest cubefree unitary divisor of the n-th biquadratefree number.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 1, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 3, 25, 26, 1, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 5, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 2, 55, 7, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 9, 73

Offset: 1

Amiram Eldar, Apr 08 2025

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

Cf. A046100, A013662, A360539.

Similar sequences: A382902, A382904, A382905, A382906.

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

list(lim) = {my(f); print1(1, ", "); for(k = 2, lim, f = factor(k); if(vecmax(f[, 2]) < 4, print1(prod(i = 1, #f~, f[i, 1]^if(f[i, 2] < 3, f[i, 2], 0)), ", ")));}

a(n) = A360539(A046100(n)).
a(n) = A382902(n)^3 / A046100(n)^2.
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = zeta(4)^2 * Product_{p prime} (1 - 1/p^3 + 1/p^6 - 1/p^7) = 0.98420942697128846925... .

A384040 The number of integers k from 1 to n such that gcd(n,k) is a cubefull number.

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, May 18 2025

The number of integers k from 1 to n such that the cubefree part (A360539) of gcd(n,k) is 1.

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

Cf. A005117, A036966, A078429, A360539, A254926.

The number of integers k from 1 to n such that gcd(n,k) is: A026741 (odd), A062570 (power of 2), A063659 (squarefree), A078429 (cube), A116512 (power of a prime), A117494 (prime), A126246 (1 or 2), A206369 (square), A254926 (cubefree), A372671 (3-smooth), A384039 (powerful), this sequence (cubefull), A384041 (exponentially odd), A384042 (5-rough).

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

Multiplicative with a(p^e) = (p^3-p^2+1)*p^(e-3) if e >= 3, p*(p-1) if e = 2, and p-1 otherwise.
a(n) >= A384039(n), with equality if and only if n is squarefree (A005117).
Dirichlet g.f.: zeta(s-1) * Product_{p prime} (1 - 1/p^s + 1/p^(3*s)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (1 - 1/p^2 + 1/p^6) = 0.62159731307414305346... .

cp's OEIS Frontend

A366077 The number of prime factors of the cubefree part of n (A360539), counted with multiplicity.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A366078 The number of distinct prime factors of the cubefree part of n (A360539).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A366147 The number of divisors of the cubefree part of n (A360539).

Views

Author

Keywords

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

A384049 The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is cubefree.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A360540 a(n) is the cubefull part of n: the largest divisor of n that is a cubefull number (A036966).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A384051 The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is a cubefull number.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A367168 The largest unitary divisor of n that is a term of A138302.

Views

Author

Keywords

Comments

Links

Crossrefs