A053149 -id:A053149 - cp's OEIS frontend

A356191 a(n) is the smallest exponentially odd number that is divisible by n.

1, 2, 3, 8, 5, 6, 7, 8, 27, 10, 11, 24, 13, 14, 15, 32, 17, 54, 19, 40, 21, 22, 23, 24, 125, 26, 27, 56, 29, 30, 31, 32, 33, 34, 35, 216, 37, 38, 39, 40, 41, 42, 43, 88, 135, 46, 47, 96, 343, 250, 51, 104, 53, 54, 55, 56, 57, 58, 59, 120, 61, 62, 189, 128, 65

Offset: 1

Amiram Eldar, Jul 29 2022

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

Cf. A064549, A007913, A268335, A336643, A350390.

Similar sequences: A053149, A053143, A053167, A066638, A087320, A087321, A197863, A356192, A356193, A356194.

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

for(n=1, 100, print1(direuler(p=2, n, 1/(1 - p^2*X^2) * (1 + p*X + p^3*X^2 - p^2*X^2))[n], ", ")) \\ Vaclav Kotesovec, Sep 09 2023

Multiplicative with a(p^e) = p^e if e is odd and p^(e+1) otherwise.
a(n) = n iff n is in A268335.
a(n) = A064549(n)/A007913(n).
a(n) = n*A336643(n).
a(n) = n^2/A350390(n).
From Vaclav Kotesovec, Sep 09 2023: (Start)
Let f(s) = Product_{p prime} (1 - p^(6-5*s) + p^(7-5*s) + 2*p^(5-4*s) - p^(6-4*s) + p^(3-3*s) - p^(4-3*s) - 2*p^(2-2*s)).
Sum_{k=1..n} a(k) ~ Pi^2 * f(2) * n^2 / 24 * (log(n) + 3*gamma - 1/2 + 12*zeta'(2)/Pi^2 + f'(2)/f(2)), where
f(2) = Product_{p prime} (1 - 4/p^2 + 4/p^3 - 1/p^4) = A256392 = 0.2177787166195363783230075141194468131307977550013559376482764035236264911...,
f'(2) = f(2) * Sum_{p prime} (11*p - 5) * log(p) / (p^3 + p^2 - 3*p + 1) = f(1) * 4.7165968208567630786609552448708126340725121316268495170070986645608062483...
and gamma is the Euler-Mascheroni constant A001620. (End)

A356192 a(n) is the smallest cubefull exponentially odd number (A335988) that is divisible by n.

Original entry on oeis.org

1, 8, 27, 8, 125, 216, 343, 8, 27, 1000, 1331, 216, 2197, 2744, 3375, 32, 4913, 216, 6859, 1000, 9261, 10648, 12167, 216, 125, 17576, 27, 2744, 24389, 27000, 29791, 32, 35937, 39304, 42875, 216, 50653, 54872, 59319, 1000, 68921, 74088, 79507, 10648, 3375, 97336

Offset: 1

Amiram Eldar, Jul 29 2022

First differs from A053149 and A356193 at n=16.

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

Cf. A001694, A013664, A268335, A335988.

Similar sequences: A053149, A053143, A053167, A066638, A087320, A087321, A197863, A356191, A356193, A356194.

f[p_, e_] := If[OddQ[e], p^Max[e, 3], p^(e + 1)]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50]

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

Multiplicative with a(p^e) = p^max(e,3) if e is odd and p^(e+1) otherwise.
a(n) = n iff n is in A335988.
a(n) = A356191(n) iff n is a powerful number (A001694).
Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + (3*p^2-1)/(p^3*(p^2-1))) = 1.69824776889117043774... .
Sum_{k=1..n} a(k) ~ c * n^4, where c = (zeta(6)/4) * Product_{p prime} (1 - 1/p^2 + 1/p^5 - 2/p^6 + 1/p^8 + 1/p^9 - 1/p^10) = 0.1559368144... . - Amiram Eldar, Nov 13 2022

A048798 Smallest k > 0 such that n*k is a perfect cube.

Original entry on oeis.org

1, 4, 9, 2, 25, 36, 49, 1, 3, 100, 121, 18, 169, 196, 225, 4, 289, 12, 361, 50, 441, 484, 529, 9, 5, 676, 1, 98, 841, 900, 961, 2, 1089, 1156, 1225, 6, 1369, 1444, 1521, 25, 1681, 1764, 1849, 242, 75, 2116, 2209, 36, 7, 20, 2601, 338, 2809, 4, 3025, 49, 3249

Offset: 1

Charles T. Le (charlestle(AT)yahoo.com)

Note that in general the smallest number k(>0) such that n*k is a perfect m-th power (rather obviously) = (the smallest m-th power divisible by n)/n and also (slightly less obviously) =n^(m-1)/(the number of solutions of x^m==0 mod n)^m. - Henry Bottomley, Mar 03 2000

			a(12) = a(2*2*3) = 2*3*3 = 18 since 12*18 = 6^3.
a(28) = a(2*2*7) = 2*7*7 = 98 since 28*98 = 14^3.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..5000 from Peter Kagey)
Krassimir T. Atanassov, On the 22nd, the 23rd and 24th Smarandache Problems, Notes on Number Theory and Discrete Mathematics, Sophia, Bulgaria, Vol. 5, No. 2 (1998), pp. 80-82.
Krassimir T. Atanassov, On Some of Smarandache's Problems, American Research Press, 1999, 16-21.
Henry Bottomley, Some Smarandache-type multiplicative sequences.
Marcela Popescu and Mariana Nicolescu, About the Smarandache Complementary Cubic Function, Smarandache Notions Journal, Vol. 7, No. 1-2-3, 1996, pp. 54-62.
Florentin Smarandache, Only Problems, Not Solutions! (see Unsolved Problem: 28, p. 26).

Cf. A000189, A007913, A053149.
Cf. A254767 (analogous sequence with the restriction that k > n).
Cf. A002117, A013667.

a[n_] := For[k = 1, True, k++, If[ Divisible[c = k^3, n], Return[c/n]]]; Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Sep 03 2012 *)
f[p_, e_] := p^(Mod[-e, 3]); a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Sep 10 2020 *)
With[{cbs=Range[3300]^3},Table[SelectFirst[cbs,Mod[#,n]==0&]/n,{n,60}]] (* Harvey P. Dale, May 10 2024 *)

a(n)=my(f=factor(n));prod(i=1,#f[,1],f[i,1]^(-f[i,2]%3)) \\ Charles R Greathouse IV, Feb 27 2013

a(n)=for(k=1,n^2,if(ispower(k*n,3),return(k)))
vector(100,n,a(n)) \\ Derek Orr, Feb 07 2015

from math import prod
from sympy import factorint
def A048798(n): return prod(p**(-e%3) for p, e in factorint(n).items()) # Chai Wah Wu, Aug 05 2024

a(n) = A053149(n)/n = n^2/A000189(n)^3.
Multiplicative with a(p^e) = p^((-e) mod 3). - Mitch Harris, May 17 2005
Sum_{k=1..n} a(k) ~ c * n^3, where c = (zeta(9)/(3*zeta(3))) * Product_{p prime} (1 - 1/p^2 + 1/p^3) = 0.2079875504... . - Amiram Eldar, Oct 28 2022

More terms from Patrick De Geest, Feb 15 2000

A062378 n divided by largest cubefree factor of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 16, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 9, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 2

Offset: 1

Henry Bottomley, Jun 18 2001

Numerator of n/rad(n)^2, where rad is the squarefree kernel of n (A007947), denominator: A055231. - Reinhard Zumkeller, Dec 10 2002

Antti Karttunen, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.

Cf. A000189, A000578, A007948, A008834, A019555, A048798, A050985, A053149, A053150, A056551, A056552. See A003557 for squares and A062379 for 4th powers.
Differs from A073753 for the first time at n=90, where a(90) = 1, while A073753(90) = 3.
Cf. A007947, A055231.

f[p_, e_] := p^Max[e-2, 0]; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Sep 07 2020 *)

a(n)=my(f=factor(n));prod(i=1,#f~,f[i,1]^max(f[i,2]-2,0)) \\ Charles R Greathouse IV, Aug 08 2013

(define (A062378 n) (/ n (A007948 n)))
(definec (A007948 n) (if (= 1 n) n (* (expt (A020639 n) (min 2 (A067029 n))) (A007948 (A028234 n)))))
;; Antti Karttunen, Nov 28 2017

a(n) = n / A007948(n).
a(n) = A003557(A003557(n)). - Antti Karttunen, Nov 28 2017
Multiplicative with a(p^e) = p^max(e-2, 0). - Amiram Eldar, Sep 07 2020
Dirichlet g.f.: zeta(s-1) * Product_{p prime} (1 - 1/p^(s-1) + 1/p^s - 1/p^(2*s-1) + 1/p^(2*s)). - Amiram Eldar, Dec 07 2023

A356193 a(n) is the smallest cubefull number (A036966) that is a multiple of n.

Original entry on oeis.org

1, 8, 27, 8, 125, 216, 343, 8, 27, 1000, 1331, 216, 2197, 2744, 3375, 16, 4913, 216, 6859, 1000, 9261, 10648, 12167, 216, 125, 17576, 27, 2744, 24389, 27000, 29791, 32, 35937, 39304, 42875, 216, 50653, 54872, 59319, 1000, 68921, 74088, 79507, 10648, 3375, 97336

Offset: 1

Amiram Eldar, Jul 29 2022

First differs from A053149 and A356192 at n=16.

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

Cf. A002117, A036966, A360541.

Similar sequences: A053143, A053149, A053167, A066638, A087320, A087321, A197863, A356191, A356192, A356194.

f[p_, e_] := p^Max[e, 3]; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50]

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

Multiplicative with a(p^e) = p^max(e,3).
a(n) = n iff n is in A036966.
Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + (3*p-2)/(p^3*(p-1))) = 1.76434793373691907811... . - Amiram Eldar, Jul 29 2022
Sum_{k=1..n} a(k) ~ c * n^4, where c = (zeta(3)/4) * Product_{p prime} (1 - 1/p^2 - 1/p^3 + 2/p^5 - 1/p^6 - 1/p^8 + 2/p^9 - 1/p^10) = 0.1559111567... . - Amiram Eldar, Nov 13 2022
a(n) = n * A360541(n). - Amiram Eldar, Sep 01 2023

A356194 a(n) is the smallest multiple of n whose prime factorization exponents are all powers of 2.

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Jul 29 2022

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

Cf. A138302, A353897.

Similar sequences: A053143, A053149, A053167, A066638, A087320, A087321, A197863, A356191, A356192, A356193.

f[p_, e_] := p^(2^Ceiling[Log2[e]]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

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

Multiplicative with a(p^e) = p^(2^ceiling(log_2(e))).

a(n) = n iff n is in A138302.

A056552 Powerfree kernel of cubefree part of n.

Original entry on oeis.org

1, 2, 3, 2, 5, 6, 7, 1, 3, 10, 11, 6, 13, 14, 15, 2, 17, 6, 19, 10, 21, 22, 23, 3, 5, 26, 1, 14, 29, 30, 31, 2, 33, 34, 35, 6, 37, 38, 39, 5, 41, 42, 43, 22, 15, 46, 47, 6, 7, 10, 51, 26, 53, 2, 55, 7, 57, 58, 59, 30, 61, 62, 21, 1, 65, 66, 67, 34, 69, 70, 71, 3, 73, 74, 15, 38, 77

Offset: 1

Henry Bottomley, Jun 25 2000

			a(32) = 2 because cubefree part of 32 is 4 and powerfree kernel of 4 is 2.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.

Cf. A000189, A000578, A007947, A008834, A013664, A019555, A048798, A050985, A053149, A053150, A056551.

f[p_, e_] :=  p^If[Divisible[e, 3], 0, 1]; a[n_] := Times @@ (f @@@ FactorInteger[ n]); Array[a, 100] (* Amiram Eldar, Aug 29 2019 *)

a(n) = my(f=factor(n)); for (k=1, #f~, if (frac(f[k,2]/3), f[k,2] = 1, f[k,2] = 0)); factorback(f); \\ Michel Marcus, Feb 28 2019

a(n) = A007947(A050985(n)) = A019555(A050985(n)) = n/(A053150(n)*A000189(n)) = A019555(n)/A053150(n) = A056551(n)^(1/3).
If n = Product_{j} Pj^Ej then a(n) = Product_{j} Pj^Fj, where Fj = 0 if Ej is 0 or a multiple of 3 and Fj = 1 otherwise.
Multiplicative with a(p^e) = p^(if 3|e, then 0, else 1). - Mitch Harris, Apr 19 2005
Sum_{k=1..n} a(k) ~ c * n^2, where c = (zeta(6)/2) * Product_{p prime} (1 - 1/p^2 + 1/p^3 - 1/p^4) = 0.3480772773... . - Amiram Eldar, Oct 28 2022
Dirichlet g.f.: zeta(3*s) * Product_{p prime} (1 + 1/p^(s-1) + 1/p^(2*s-1)). - Amiram Eldar, Sep 16 2023

A365488 The number of divisors of the smallest number whose cube is divisible by n.

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Sep 05 2023

First differs from A365171 at n = 32.

The number of divisors of the smallest cube divisible by n, A053149(n), is A365489(n).

Amiram Eldar, Table of n, a(n) for n = 1..10000
Vaclav Kotesovec, Graph - the asymptotic ratio (100000 terms)

Cf. A000005, A005117, A019555, A053149, A365171, A365489, A365498.

f[p_, e_] := Ceiling[e/3] + 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
With[{c=Range[200]^3},Table[DivisorSigma[0,Surd[SelectFirst[c,Mod[#,n]==0&],3]],{n,90}]] (* Harvey P. Dale, Sep 15 2024 *)

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

a(n) = A000005(A019555(n)).
Multiplicative with a(p^e) = ceiling(e/3) + 1.
a(n) <= A000005(n) with equality if and only if n is squarefree  (A005117).
Dirichlet g.f.: zeta(s) * zeta(3*s) * Product_{p prime} (1 + 1/p^s - 1/p^(3*s)).
From Vaclav Kotesovec, Sep 06 2023: (Start)
Dirichlet g.f.: zeta(s)^2 * zeta(3*s) * Product_{p prime} (1 - 1/p^(2*s) - 1/p^(3*s) + 1/p^(4*s)).
Let f(s) = Product_{p prime} (1 - 1/p^(2*s) - 1/p^(3*s) + 1/p^(4*s)).
Sum_{k=1..n} a(k) ~ zeta(3) * f(1) * n * (log(n) + 2*gamma - 1 + 3*zeta'(3)/zeta(3) + f'(1)/f(1)), where
f(1) = Product_{p prime} (1 - 1/p^2 - 1/p^3 + 1/p^4) = 0.5358961538283379998085026313185459506482223745141452711510108346133288...,
f'(1) = f(1) * Sum_{p prime} (-4 + 3*p + 2*p^2) * log(p) / (1 - p - p^2 + p^4) = f(1) * 1.4525924794451595590371439593828547341482465114411929136723476679...
and gamma is the Euler-Mascheroni constant A001620. (End)

A056551 Smallest cube divisible by n divided by largest cube which divides n.

Original entry on oeis.org

1, 8, 27, 8, 125, 216, 343, 1, 27, 1000, 1331, 216, 2197, 2744, 3375, 8, 4913, 216, 6859, 1000, 9261, 10648, 12167, 27, 125, 17576, 1, 2744, 24389, 27000, 29791, 8, 35937, 39304, 42875, 216, 50653, 54872, 59319, 125, 68921, 74088, 79507, 10648

Offset: 1

Henry Bottomley, Jun 25 2000

			a(16) = 8 since smallest cube divisible by 16 is 64 and smallest cube which divides 16 is 8 and 64/8 = 8.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.

Cf. A000189, A000578, A008834, A019555, A048798, A050985, A053149, A053150, A056552.

Cf. A002117, A013670, A330596.

f[p_, e_] := p^If[Divisible[e, 3], 0, 1]; a[n_] := (Times @@ (f @@@ FactorInteger[ n]))^3; Array[a, 100] (* Amiram Eldar, Aug 29 2019*)

a(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i,2]%3, f[i,1], 1))^3; } \\ Amiram Eldar, Oct 28 2022

a(n) = A053149(n)/A008834(n) = A048798(n)*A050985(n) = A056552(n)^3.
From Amiram Eldar, Oct 28 2022: (Start)
Multiplicative with a(p^e) = 1 if e is divisible by 3, and a(p^e) = p^3 otherwise.
Sum_{k=1..n} a(k) ~ c * n^4, where c = (zeta(12)/(4*zeta(3))) * Product_{p prime} (1 - 1/p^2 + 1/p^3) = A013670 * A330596 / (4*A002117) = 0.1557163105... . (End)
Dirichlet g.f.: zeta(3*s) * Product_{p prime} (1 + 1/p^(s-3) + 1/p^(2*s-3)). - Amiram Eldar, Sep 16 2023

A365489 The number of divisors of the smallest cube divisible by n.

Original entry on oeis.org

1, 4, 4, 4, 4, 16, 4, 4, 4, 16, 4, 16, 4, 16, 16, 7, 4, 16, 4, 16, 16, 16, 4, 16, 4, 16, 4, 16, 4, 64, 4, 7, 16, 16, 16, 16, 4, 16, 16, 16, 4, 64, 4, 16, 16, 16, 4, 28, 4, 16, 16, 16, 4, 16, 16, 16, 16, 16, 4, 64, 4, 16, 16, 7, 16, 64, 4, 16, 16, 64, 4, 16, 4

Offset: 1

Amiram Eldar, Sep 05 2023

The number of divisors of the cube root of the smallest cube divisible by n, A019555(n), is A365488(n).

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

Cf. A000005, A019555, A053149, A365345, A365488.

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

a(n) = vecprod(apply(x -> 3*((x-1)\3) + 4, factor(n)[, 2]));

a(n) = A000005(A053149(n)).
Multiplicative with a(p^e) = 3*ceiling(e/3) + 1.
Dirichlet g.f.: zeta(s) * zeta(3*s) * Product_{p prime} (1 + 3/p^s - 1/p^(3*s)).

cp's OEIS Frontend

A356191 a(n) is the smallest exponentially odd number that is divisible by n.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A356192 a(n) is the smallest cubefull exponentially odd number (A335988) that is divisible by n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A048798 Smallest k > 0 such that n*k is a perfect cube.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Formula

Extensions

A062378 n divided by largest cubefree factor of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula

A356193 a(n) is the smallest cubefull number (A036966) that is a multiple of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A356194 a(n) is the smallest multiple of n whose prime factorization exponents are all powers of 2.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A056552 Powerfree kernel of cubefree part of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica