A059485 -id:A059485 - cp's OEIS frontend

A051037 5-smooth numbers, i.e., numbers whose prime divisors are all <= 5.

1, 2, 3, 4, 5, 6, 8, 9, 10, 12, 15, 16, 18, 20, 24, 25, 27, 30, 32, 36, 40, 45, 48, 50, 54, 60, 64, 72, 75, 80, 81, 90, 96, 100, 108, 120, 125, 128, 135, 144, 150, 160, 162, 180, 192, 200, 216, 225, 240, 243, 250, 256, 270, 288, 300, 320, 324, 360, 375, 384, 400, 405

Offset: 1

Eric W. Weisstein

Sometimes called the Hamming sequence, since Hamming asked for an efficient algorithm to generate the list, in ascending order, of all numbers of the form 2^i*3^j*5^k for i,j,k >= 0. The problem was popularized by Edsger Dijkstra.
Numbers k such that 8*k = EulerPhi(30*k). - Artur Jasinski, Nov 05 2008
Where record values greater than 1 occur in A165704: A165705(n) = A165704(a(n)). - Reinhard Zumkeller, Sep 26 2009
Also called "harmonic whole numbers", see Howard and Longair, 1982, Table I, page 121. - Hugo Pfoertner, Jul 16 2020
Also called ugly numbers, although it is not clear why. - Gus Wiseman, May 21 2021
Some woody bamboo species have extraordinarily long and stable flowering intervals that belong to this sequence. The model by Veller, Nowak & Davis justifies this observation from the evolutionary point of view. - Andrey Zabolotskiy, Jun 27 2021
Also those integers k for which, for every prime p > 5, p^(4*k) - 1 == 0 (mod 240*k). - Federico Provvedi, May 23 2022
As noted in the comments to A085152, Størmer's theorem implies that the only pairs of consecutive integers that appear as consecutive terms of this sequence are (1,2), (2,3), (3,4), (4,5), (5,6), (8,9), (9,10), (15,16), (24,25), and (80,81). These all represent significant musical intervals. - Hal M. Switkay, Dec 05 2022

			From _Gus Wiseman_, May 21 2021: (Start)
The sequence of terms together with their prime indices begins:
      1: {}            25: {3,3}
      2: {1}           27: {2,2,2}
      3: {2}           30: {1,2,3}
      4: {1,1}         32: {1,1,1,1,1}
      5: {3}           36: {1,1,2,2}
      6: {1,2}         40: {1,1,1,3}
      8: {1,1,1}       45: {2,2,3}
      9: {2,2}         48: {1,1,1,1,2}
     10: {1,3}         50: {1,3,3}
     12: {1,1,2}       54: {1,2,2,2}
     15: {2,3}         60: {1,1,2,3}
     16: {1,1,1,1}     64: {1,1,1,1,1,1}
     18: {1,2,2}       72: {1,1,1,2,2}
     20: {1,1,3}       75: {2,3,3}
     24: {1,1,1,2}     80: {1,1,1,1,3}
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Benoit Cloitre, Plot of abs(f(n)-s(n)) vs its mean values (blue) and vs loglog(n) (red).
M. J. Dominus, Infinite Lists in Perl.
Deborah Howard and Malcolm Longair, Harmonic Proportion and Palladio's "Quattro Libri", Journal of the Society of Architectural Historians (1982) 41 (2): 116-143.
Vaclav Kotesovec, Plot of a(n) / (exp((6*log(2)*log(3)*log(5)*n)^(1/3))/sqrt(30)) for n = 1..1200000
Rosetta Code, A collection of computer codes to compute 5-smooth numbers.
Raphael Schumacher, The Formula for the Distribution of the 3-Smooth Numbers, 5-Smooth, 7-Smooth and all other Smooth Numbers, arXiv:1608.06928 [math.NT], 2016.
Sci.math, Ugly numbers.
Carl Veller, Martin A. Nowak and Charles C. Davis, Extended flowering intervals of bamboos evolved by discrete multiplication, Ecology Letters, 18 (2015), 653-659.
Eric Weisstein's World of Mathematics, Smooth Number.
Wikipedia, Regular number.
Wikipedia, Talk:Regular number. Includes a discussion of the name.
Wikipedia, Størmer's theorem.

Subsequences: A003592, A003593, A051916 , A257997.
For p-smooth numbers with other values of p, see A003586, A002473, A051038, A080197, A080681, A080682, A080683.
Cf. A006530, A112757, A112758, A112759, A112763, A112764, A291719.
The partitions with these Heinz numbers are counted by A001399.
The conjugate opposite is A033942, counted by A004250.
The opposite is A059485, counted by A004250.
The non-3-smooth case is A080193, counted by A069905.
The conjugate is A037144, counted by A001399.
The complement is A279622, counted by A035300.
Requiring the sum of prime indices to be even gives A344297.
Cf. A000244, A002182, A002183, A035301, A056239, A112798, A261144, A344293.

import Data.Set (singleton, deleteFindMin, insert)
a051037 n = a051037_list !! (n-1)
a051037_list = f $ singleton 1 where
   f s = y : f (insert (5 * y) $ insert (3 * y) $ insert (2 * y) s')
               where (y, s') = deleteFindMin s
-- Reinhard Zumkeller, May 16 2015

[n: n in [1..500] | PrimeDivisors(n) subset [2,3,5]]; // Bruno Berselli, Sep 24 2012

A051037 := proc(n)
    option remember;
    local a;
    if n = 1 then
        1;
    else
        for a from procname(n-1)+1 do
            numtheory[factorset](a) minus {2, 3,5 } ;
            if % = {} then
                return a;
            end if;
        end do:
    end if;
end proc:
seq(A051037(n),n=1..100) ; # R. J. Mathar, Nov 05 2017

mx = 405; Sort@ Flatten@ Table[ 2^a*3^b*5^c, {a, 0, Log[2, mx]}, {b, 0, Log[3, mx/2^a]}, {c, 0, Log[5, mx/(2^a*3^b)]}] (* Or *)
Select[ Range@ 405, Last@ Map[First, FactorInteger@ #] < 7 &] (* Robert G. Wilson v *)
With[{nn=10},Select[Union[Times@@@Flatten[Table[Tuples[{2,3,5},n],{n,0,nn}],1]],#<=2^nn&]] (* Harvey P. Dale, Feb 28 2022 *)

test(n)= {m=n; forprime(p=2,5, while(m%p==0,m=m/p)); return(m==1)}
for(n=1,500,if(test(n),print1(n",")))

a(n)=local(m); if(n<1,0,n=a(n-1); until(if(m=n, forprime(p=2,5, while(m%p==0,m/=p)); m==1),n++); n)

list(lim)=my(v=List(),s,t); for(i=0,logint(lim\=1,5), t=5^i; for(j=0,logint(lim\t,3), s=t*3^j; while(s<=lim, listput(v,s); s<<=1))); Set(v) \\ Charles R Greathouse IV, Sep 21 2011; updated Sep 19 2016

smooth(P:vec,lim)={ my(v=List([1]),nxt=vector(#P,i,1),indx,t);
while(1, t=vecmin(vector(#P,i,v[nxt[i]]*P[i]),&indx);
if(t>lim,break); if(t>v[#v],listput(v,t)); nxt[indx]++);
Vec(v)
};
smooth([2,3,5], 1e4) \\ Charles R Greathouse IV, Dec 03 2013

is_A051037(n)=n<7||vecmax(factor(n,6)[, 1])<7 \\ M. F. Hasler, Jan 16 2015

def isok(n):
  while n & 1 == 0: n >>= 1
  while n % 3 == 0: n //= 3
  while n % 5 == 0: n //= 5
  return n == 1 #  Darío Clavijo, Dec 30 2022

from sympy import integer_log
def A051037(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x):
        c = n+x
        for i in range(integer_log(x,5)[0]+1):
            for j in range(integer_log(y:=x//5**i,3)[0]+1):
                c -= (y//3**j).bit_length()
        return c
    return bisection(f,n,n) # Chai Wah Wu, Sep 16 2024

# faster for initial segment of sequence
import heapq
from itertools import islice
def A051037gen(): # generator of terms
    v, oldv, h, psmooth_primes, = 1, 0, [1], [2, 3, 5]
    while True:
        v = heapq.heappop(h)
        if v != oldv:
            yield v
            oldv = v
            for p in psmooth_primes:
                    heapq.heappush(h, v*p)
print(list(islice(A051037gen(), 65))) # Michael S. Branicky, Sep 17 2024

Let s(n) = Card(k | a(k)Benoit Cloitre, Dec 30 2001
The characteristic function of this sequence is given by:
  Sum_{n>=1} x^a(n) = Sum_{n>=1} -Möbius(30*n)*x^n/(1-x^n). - Paul D. Hanna, Sep 18 2011
a(n) = A143207(n) / 30. - Reinhard Zumkeller, Sep 13 2011
A204455(15*a(n)) = 15, and only for these numbers. - Wolfdieter Lang, Feb 04 2012
A006530(a(n)) <= 5. - Reinhard Zumkeller, May 16 2015
Sum_{n>=1} 1/a(n) = Product_{primes p <= 5} p/(p-1) = (2*3*5)/(1*2*4) = 15/4. - Amiram Eldar, Sep 22 2020

A357982 Replace prime(k) with A000009(k) in the prime factorization of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 3, 1, 4, 2, 2, 1, 5, 1, 6, 2, 2, 3, 8, 1, 4, 4, 1, 2, 10, 2, 12, 1, 3, 5, 4, 1, 15, 6, 4, 2, 18, 2, 22, 3, 2, 8, 27, 1, 4, 4, 5, 4, 32, 1, 6, 2, 6, 10, 38, 2, 46, 12, 2, 1, 8, 3, 54, 5, 8, 4, 64, 1, 76, 15, 4, 6, 6, 4, 89, 2, 1

Offset: 1

Gus Wiseman, Oct 25 2022

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798. This sequence gives the number of ways to choose a strict partition of each prime index of n.

The indices i, where a(i) = 1, form A003586, and the indices j, where a(j) > 1, form A059485. - Ivan N. Ianakiev, Oct 27 2022

			The a(121) = 9 twice-partitions are: (5)(5), (5)(41), (5)(32), (41)(5), (41)(41), (41)(32), (32)(5), (32)(41), (32)(32).

Other multiplicative sequences: A003961, A357852, A064988, A064989, A357980.
The non-strict version is A299200.
A horizontal version is A357978, non-strict A357977.
A000040 lists the primes.
A056239 adds up prime indices, row-sums of A112798.
Cf. A000041, A000720, A003964, A063834, A076610, A215366, A273873, A296150, A299201-A299203, A357975, A357979, A357983.
Cf. A003586, A059485.

Table[Times@@Cases[FactorInteger[n],{p_,k_}:>PartitionsQ[PrimePi[p]]^k],{n,100}]

f9(n) = polcoeff( prod( k=1, n, 1 + x^k, 1 + x * O(x^n)), n); \\ A000009
a(n) = my(f=factor(n)); for (k=1, #f~, f[k,1] = f9(primepi(f[k,1]))); factorback(f); \\ Michel Marcus, Oct 26 2022

A279622 Numbers with a prime factor greater than 5.

Original entry on oeis.org

7, 11, 13, 14, 17, 19, 21, 22, 23, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 46, 47, 49, 51, 52, 53, 55, 56, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 76, 77, 78, 79, 82, 83, 84, 85, 86, 87, 88, 89, 91, 92, 93, 94

Offset: 1

Vincenzo Librandi, Dec 21 2016

Vincenzo Librandi, Table of n, a(n) for n = 1..9825

Complement of A051037.

Cf. A059485.

[n: n in [1..100] | not PrimeDivisors(n) subset [2, 3, 5]]; // Vincenzo Librandi, Jan 29 2017

fQ[n_]:=!PowerMod[30, n, n] == 0; Select[Range[100], fQ]
Select[Range[100],Max[FactorInteger[#][[;;,1]]]>5&] (* Harvey P. Dale, Feb 28 2023 *)

isok(n) = vecmax(factor(n)[,1]) > 5; \\ Michel Marcus, Dec 21 2016

is(n)=if(n<7, return(0)); n>>=valuation(n,2); n/=3^valuation(n,2) * 5^valuation(n,5); n>1 \\ Charles R Greathouse IV, Dec 22 2016

from sympy import integer_log
def A279622(n):
    def f(x):
        c = n
        for i in range(integer_log(x,5)[0]+1):
            for j in range(integer_log(y:=x//5**i,3)[0]+1):
                c += (y//3**j).bit_length()
        return c
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Sep 16 2024

a(n) = n + O(log^3 n). - Charles R Greathouse IV, Dec 22 2016

A354181 Numbers whose number of divisors is not a 3-smooth number.

Original entry on oeis.org

16, 48, 64, 80, 81, 112, 144, 162, 176, 192, 208, 240, 272, 304, 320, 324, 336, 368, 400, 405, 432, 448, 464, 496, 512, 528, 560, 567, 576, 592, 624, 625, 648, 656, 688, 704, 720, 729, 752, 784, 810, 816, 832, 848, 880, 891, 912, 944, 960, 976, 1008, 1024, 1040

Offset: 1

Amiram Eldar, May 18 2022

nonn

Number whose prime factorization includes an exponent e such that e+1 is in A059485.

The asymptotic density of this sequence is 1 - Product_{p prime} ((1 - 1/p) * Sum_{k in A003586} 1/p^(k-1)) = 0.0512963858... (Hilberdink, 2022).

			16 is a term since A000005(16) = 5 is not a 3-smooth number.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Titus Hilberdink, How often is d(n) a power of a given integer?, Journal of Number Theory, Vol. 236 (2022), pp. 261-279.
Index entries for sequences computed from exponents in factorization of n.

Cf. A000005, A003586, A059485, A036537, A162643.

smoothQ[n_] := n == 2^IntegerExponent[n, 2] * 3^IntegerExponent[n, 3]; Select[Range[1000], !smoothQ[DivisorSigma[0, #]] &]

is(n) = n>>=valuation(n, 2); n/=3^valuation(n, 3); n>1; \\ A059485
isok(m) = is(numdiv(m)); \\ Michel Marcus, May 19 2022

cp's OEIS Frontend

A051037 5-smooth numbers, i.e., numbers whose prime divisors are all <= 5.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Python

Python

Python

Formula

A357982 Replace prime(k) with A000009(k) in the prime factorization of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

A279622 Numbers with a prime factor greater than 5.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Python

Formula

A354181 Numbers whose number of divisors is not a 3-smooth number.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI