A261144 -id:A261144 - 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

A344296 Numbers with at least as many prime factors (counted with multiplicity) as half their sum of prime indices.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 9, 10, 12, 16, 18, 20, 24, 27, 28, 30, 32, 36, 40, 48, 54, 56, 60, 64, 72, 80, 81, 84, 88, 90, 96, 100, 108, 112, 120, 128, 144, 160, 162, 168, 176, 180, 192, 200, 208, 216, 224, 240, 243, 252, 256, 264, 270, 280, 288, 300, 320, 324, 336, 352

Offset: 1

Gus Wiseman, May 16 2021

nonn

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.

These are the Heinz numbers of certain partitions counted by A025065, but different from palindromic partitions, which have Heinz numbers A265640.

			The sequence of terms together with their prime indices begins:
      1: {}            30: {1,2,3}
      2: {1}           32: {1,1,1,1,1}
      3: {2}           36: {1,1,2,2}
      4: {1,1}         40: {1,1,1,3}
      6: {1,2}         48: {1,1,1,1,2}
      8: {1,1,1}       54: {1,2,2,2}
      9: {2,2}         56: {1,1,1,4}
     10: {1,3}         60: {1,1,2,3}
     12: {1,1,2}       64: {1,1,1,1,1,1}
     16: {1,1,1,1}     72: {1,1,1,2,2}
     18: {1,2,2}       80: {1,1,1,1,3}
     20: {1,1,3}       81: {2,2,2,2}
     24: {1,1,1,2}     84: {1,1,2,4}
     27: {2,2,2}       88: {1,1,1,5}
     28: {1,1,4}       90: {1,2,2,3}

The case with difference at least 1 is A322136.
The case of equality is A340387, counted by A000041 or A035363.
The opposite version is A344291, counted by A110618.
The conjugate version is A344414, with even-weight case A344416.
A025065 counts palindromic partitions, ranked by A265640.
A056239 adds up prime indices, row sums of A112798.
A300061 lists numbers whose sum of prime indices is even.
Cf. A001399, A002865, A025147, A027336, A036036, A067712, A244990, A261144, A325691, A344293, A344295.

Select[Range[100],PrimeOmega[#]>=Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]/2&]

A056239(a(n)) <= 2*A001222(a(n)).

a(n) = A322136(n)/4.

A209862 Permutation of nonnegative integers which maps A209642 into ascending order (A209641).

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 11, 13, 14, 15, 16, 17, 18, 20, 24, 19, 21, 25, 22, 26, 28, 23, 27, 29, 30, 31, 32, 33, 34, 36, 40, 48, 35, 37, 41, 49, 38, 42, 50, 44, 52, 56, 39, 43, 51, 45, 53, 57, 46, 54, 58, 60, 47, 55, 59, 61, 62, 63, 64, 65, 66, 68, 72, 80, 96, 67, 69, 73, 81, 97, 70, 74, 82, 98, 76, 84, 100, 88, 104, 112, 71, 75, 83

Offset: 0

Antti Karttunen, Mar 24 2012

nonn

Conjecture: For all n, a(A054429(n)) = A054429(a(n)), i.e. A054429 acts as a homomorphism (automorphism) of the cyclic group generated by this permutation. This implies also a weaker conjecture given in A209860.
From Gus Wiseman, Aug 24 2021: (Start)
As a triangle with row lengths 2^n, T(n,k) for n > 0 appears (verified up to n = 2^15) to be the unique nonnegative integer whose binary indices are the k-th subset of {1..n} containing n. Here, a binary index of n (row n of A048793) is any position of a 1 in its reversed binary expansion, and sets are sorted first by length, then lexicographically. For example, the triangle begins:
   1
   2  3
   4  5  6  7
   8  9 10 12 11 13 14 15
  16 17 18 20 24 19 21 25 22 26 28 23 27 29 30 31
Mathematica: Table[Total[2^(Append[#,n]-1)]&/@Subsets[Range[n-1]],{n,5}]
Row lengths are A000079 (shifted right). Also Column k = 1.
Row sums are A010036.
Using reverse-lexicographic order gives A059893.
Using lexicographic order gives A059894.
Taking binary indices to prime indices gives A339195 (or A019565).
The ordering of sets is A344084.
A version using Heinz numbers is A344085.
(End)

			From _Gus Wiseman_, Aug 24 2021: (Start)
The terms, their binary expansions, and their binary indices begin:
   0:      ~ {}
   1:    1 ~ {1}
   2:   10 ~ {2}
   3:   11 ~ {1,2}
   4:  100 ~ {3}
   5:  101 ~ {1,3}
   6:  110 ~ {2,3}
   7:  111 ~ {1,2,3}
   8: 1000 ~ {4}
   9: 1001 ~ {1,4}
  10: 1010 ~ {2,4}
  12: 1100 ~ {3,4}
  11: 1011 ~ {1,2,4}
  13: 1101 ~ {1,3,4}
  14: 1110 ~ {2,3,4}
  15: 1111 ~ {1,2,3,4}
(End)

Antti Karttunen, Table of n, a(n) for n = 0..32767
Index entries for sequences that are permutations of the natural numbers

Inverse permutation: A209861. Cf. A209860, A209863, A209864, A209865, A209866, A209867, A209868.

Cf. A010036, A026793, A048793, A111059, A147655, A246688, A246867, A261144, A272020, A294648, A339360.

a(n) = A209859(A036044(A209641(n))) = A209859(A056539(A209641(n))).

A339195 Triangle of squarefree numbers grouped by greatest prime factor, read by rows.

Original entry on oeis.org

1, 2, 3, 6, 5, 10, 15, 30, 7, 14, 21, 35, 42, 70, 105, 210, 11, 22, 33, 55, 66, 77, 110, 154, 165, 231, 330, 385, 462, 770, 1155, 2310, 13, 26, 39, 65, 78, 91, 130, 143, 182, 195, 273, 286, 390, 429, 455, 546, 715, 858, 910, 1001, 1365, 1430, 2002, 2145, 2730, 3003, 4290, 5005, 6006, 10010, 15015, 30030

Offset: 0

Gus Wiseman, Dec 02 2020

Also Heinz numbers of subsets of {1..n} that contain n if n>0, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
A019565 in its triangle form, with each row's terms in increasing order. - Peter Munn, Feb 26 2021
From David James Sycamore, Jan 09 2025: (Start)
Alternative definition, with offset = 1: a(1) = 1. For n>1 if a(n-1) = A002110(k), a(n) = prime(k+1). Otherwise a(n) is the smallest novel squarefree number whose prime factors have already occurred as previous terms.
Permutation of A005117, Squarefree version A379746. (End)

			Triangle begins:
   1
   2
   3   6
   5  10  15  30
   7  14  21  35  42  70  105  210

Alois P. Heinz, Rows n = 0..14, flattened
Michael De Vlieger, Plot p | a(n) at (x,y) = (n,pi(p)), n = 0..2047, 12X vertical exaggeration.
Michael De Vlieger, Fan style binary tree showing a(n), n = 0..2047, with a color function related to the order of a(n) in A019565.
Michael De Vlieger, Fan style binary tree showing a(n), n = 0..2047, with a color function showing 1 in gray, primes in red, primorials in bright green, even squarefree semiprimes in yellow, odd squarefree semiprimes in light green, thereafter, progressively deeper green related to omega(a(n)) = m until m >= 6.

A011782 gives row lengths.
A339360 gives row sums.
A008578 (shifted) is column k = 1.
A100484 is column k = 2.
A001748 is column k = 3.
A002110 is column k = 2^(n-1).
A070826 is column k = 2^(n-1) - 1.
A209862 takes prime indices to binary indices in these terms.
A246867 groups squarefree numbers by Heinz weight, with row sums A147655.
A261144 divides the n-th row by prime(n), with row sums A054640.
A339116 is the restriction to semiprimes, with row sums A339194.
A005117 lists squarefree numbers, ordered lexicographically by prime factors: A019565.
A006881 lists squarefree semiprimes.
A072047 counts prime factors of squarefree numbers.
A319246 is the sum of prime indices of the n-th squarefree number.
A329631 lists prime indices of squarefree numbers, reversed: A319247.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
Cf. A001221, A014342, A014466, A019565, A098350, A112798, A320656, A326882, A338901, A379770.
Cf. A379746.

T:= proc(n) option remember; `if`(n=0, 1, (p-> map(
      x-> x*p, {seq(T(i), i=0..n-1)})[])(ithprime(n)))
    end:
seq(T(n), n=0..6);  # Alois P. Heinz, Jan 08 2025

Table[Prime[n]*Sort[Times@@Prime/@#&/@Subsets[Range[n-1]]],{n,5}]

For n > 1, T(n,k) = prime(n) * A261144(n-1,k).

a(n) = A019565(A379770(n)). - Michael De Vlieger, Jan 08 2025

Row n=0 (term 1) prepended by Alois P. Heinz, Jan 08 2025

A344293 5-smooth numbers n whose sum of prime indices A056239(n) is at least twice the number of prime indices A001222(n).

Original entry on oeis.org

1, 3, 5, 9, 10, 15, 25, 27, 30, 45, 50, 75, 81, 90, 100, 125, 135, 150, 225, 243, 250, 270, 300, 375, 405, 450, 500, 625, 675, 729, 750, 810, 900, 1000, 1125, 1215, 1250, 1350, 1500, 1875, 2025, 2187, 2250, 2430, 2500, 2700, 3000, 3125, 3375, 3645, 3750, 4050

Offset: 1

Gus Wiseman, May 16 2021

nonn

A number is 5-smooth if its prime divisors are all <= 5.

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.

			The sequence of terms together with their prime indices begins:
       1: {}            125: {3,3,3}
       3: {2}           135: {2,2,2,3}
       5: {3}           150: {1,2,3,3}
       9: {2,2}         225: {2,2,3,3}
      10: {1,3}         243: {2,2,2,2,2}
      15: {2,3}         250: {1,3,3,3}
      25: {3,3}         270: {1,2,2,2,3}
      27: {2,2,2}       300: {1,1,2,3,3}
      30: {1,2,3}       375: {2,3,3,3}
      45: {2,2,3}       405: {2,2,2,2,3}
      50: {1,3,3}       450: {1,2,2,3,3}
      75: {2,3,3}       500: {1,1,3,3,3}
      81: {2,2,2,2}     625: {3,3,3,3}
      90: {1,2,2,3}     675: {2,2,2,3,3}
     100: {1,1,3,3}     729: {2,2,2,2,2,2}

Allowing any number of parts and sum gives A051037, counted by A001399.
These are Heinz numbers of the partitions counted by A266755.
Allowing parts > 5 gives A344291, counted by A110618.
The non-3-smooth case is A344294, counted by A325691.
Requiring the sum of prime indices to be even gives A344295.
A000070 counts non-multigraphical partitions, ranked by A344292.
A025065 counts partitions of n with >= n/2 parts, ranked by A344296.
A035363 counts partitions of n with n/2 parts, ranked by A340387.
A056239 adds up prime indices, row sums of A112798.
A300061 ranks partitions of even numbers, with 5-smooth case A344297.
Cf. A000041, A000244, A026811, A080193, A244990, A261144, A279622.

Select[Range[1000],PrimeOmega[#]<=Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]/2&&Max@@First/@FactorInteger[#]<=5&]

Intersection of A051037 and A344291.

A111059 a(n) = Product_{k=1..n} A005117(k), the product of the first n squarefree positive integers.

Original entry on oeis.org

1, 2, 6, 30, 180, 1260, 12600, 138600, 1801800, 25225200, 378378000, 6432426000, 122216094000, 2566537974000, 56463835428000, 1298668214844000, 33765373585944000, 979195833992376000, 29375875019771280000

Offset: 1

Leroy Quet, Oct 07 2005

nonn

Do all terms belong to A242031 (weakly decreasing prime signature)? - Gus Wiseman, May 14 2021

			Since the first 6 squarefree positive integers are 1, 2, 3, 5, 6, 7, the 6th term of the sequence is 1*2*3*5*6*7 = 1260.
From _Gus Wiseman_, May 14 2021: (Start)
The sequence of terms together with their prime signatures begins:
             1: ()
             2: (1)
             6: (1,1)
            30: (1,1,1)
           180: (2,2,1)
          1260: (2,2,1,1)
         12600: (3,2,2,1)
        138600: (3,2,2,1,1)
       1801800: (3,2,2,1,1,1)
      25225200: (4,2,2,2,1,1)
     378378000: (4,3,3,2,1,1)
    6432426000: (4,3,3,2,1,1,1)
  122216094000: (4,3,3,2,1,1,1,1)
(End)

Charles R Greathouse IV, Table of n, a(n) for n = 1..417

A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A072047 applies Omega to each squarefree number.
A246867 groups squarefree numbers by Heinz weight (row sums: A147655).
A261144 groups squarefree numbers by smoothness (row sums: A054640).
A319246 gives the sum of prime indices of each squarefree number.
A329631 lists prime indices of squarefree numbers (reversed: A319247).
Cf. A001221, A002110, A014466, A070826, A338901, A339360.

Rest[FoldList[Times,1,Select[Range[40],SquareFreeQ]]] (* Harvey P. Dale, Jun 14 2011 *)

m=30;k=1;for(n=1,m,if(issquarefree(n),print1(k=k*n,",")))

More terms from Klaus Brockhaus, Oct 08 2005

A344294 5-smooth but not 3-smooth numbers k such that A056239(k) >= 2*A001222(k).

Original entry on oeis.org

5, 10, 15, 25, 30, 45, 50, 75, 90, 100, 125, 135, 150, 225, 250, 270, 300, 375, 405, 450, 500, 625, 675, 750, 810, 900, 1000, 1125, 1215, 1250, 1350, 1500, 1875, 2025, 2250, 2430, 2500, 2700, 3000, 3125, 3375, 3645, 3750, 4050, 4500, 5000, 5625, 6075, 6250

Offset: 1

Gus Wiseman, May 16 2021

nonn

A number is d-smooth iff its prime divisors are all <= d.

A prime index of k is a number m such that prime(m) divides k, and the multiset of prime indices of k is row k of A112798. This row has length A001222(k) and sum A056239(k).

			The sequence of terms together with their prime indices begins:
       5: {3}           270: {1,2,2,2,3}
      10: {1,3}         300: {1,1,2,3,3}
      15: {2,3}         375: {2,3,3,3}
      25: {3,3}         405: {2,2,2,2,3}
      30: {1,2,3}       450: {1,2,2,3,3}
      45: {2,2,3}       500: {1,1,3,3,3}
      50: {1,3,3}       625: {3,3,3,3}
      75: {2,3,3}       675: {2,2,2,3,3}
      90: {1,2,2,3}     750: {1,2,3,3,3}
     100: {1,1,3,3}     810: {1,2,2,2,2,3}
     125: {3,3,3}       900: {1,1,2,2,3,3}
     135: {2,2,2,3}    1000: {1,1,1,3,3,3}
     150: {1,2,3,3}    1125: {2,2,3,3,3}
     225: {2,2,3,3}    1215: {2,2,2,2,2,3}
     250: {1,3,3,3}    1250: {1,3,3,3,3}

Allowing any number of parts and sum gives A080193, counted by A069905.
The partitions with these Heinz numbers are counted by A325691.
Relaxing the smoothness conditions gives A344291, counted by A110618.
Allowing 3-smoothness gives A344293, counted by A266755.
A025065 counts partitions of n with at least n/2 parts, ranked by A344296.
A035363 counts partitions of n whose length is n/2, ranked by A340387.
A051037 lists 5-smooth numbers (complement: A279622).
A056239 adds up prime indices, row sums of A112798.
A257993 gives the least gap of the partition with Heinz number n.
A300061 lists numbers with even sum of prime indices (5-smooth: A344297).
A342050/A342051 list Heinz numbers of partitions with even/odd least gap.
Cf. A000041, A001399, A002865, A027336, A244990, A261144, A344295.

Select[Range[1000],PrimeOmega[#]<=Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]/2&&Max@@First/@FactorInteger[#]==5&]

Intersection of A080193 and A344291.

A344297 Heinz numbers of integer partitions of even numbers with no part greater than 3.

Original entry on oeis.org

1, 3, 4, 9, 10, 12, 16, 25, 27, 30, 36, 40, 48, 64, 75, 81, 90, 100, 108, 120, 144, 160, 192, 225, 243, 250, 256, 270, 300, 324, 360, 400, 432, 480, 576, 625, 640, 675, 729, 750, 768, 810, 900, 972, 1000, 1024, 1080, 1200, 1296, 1440, 1600, 1728, 1875, 1920

Offset: 1

Gus Wiseman, May 16 2021

nonn

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

			The sequence of terms together with their prime indices begins:
       1: {}                 81: {2,2,2,2}
       3: {2}                90: {1,2,2,3}
       4: {1,1}             100: {1,1,3,3}
       9: {2,2}             108: {1,1,2,2,2}
      10: {1,3}             120: {1,1,1,2,3}
      12: {1,1,2}           144: {1,1,1,1,2,2}
      16: {1,1,1,1}         160: {1,1,1,1,1,3}
      25: {3,3}             192: {1,1,1,1,1,1,2}
      27: {2,2,2}           225: {2,2,3,3}
      30: {1,2,3}           243: {2,2,2,2,2}
      36: {1,1,2,2}         250: {1,3,3,3}
      40: {1,1,1,3}         256: {1,1,1,1,1,1,1,1}
      48: {1,1,1,1,2}       270: {1,2,2,2,3}
      64: {1,1,1,1,1,1}     300: {1,1,2,3,3}
      75: {2,3,3}           324: {1,1,2,2,2,2}

These partitions are counted by A007980.
Including partitions of odd numbers gives A051037 (complement: A279622).
Allowing parts > 3 gives A300061.
A001358 lists semiprimes.
A035363 counts partitions whose length is half their sum, ranked by A340387.
A056239 adds up prime indices, row sums of A112798.
Cf. A001399, A002620, A030229, A080193, A244990, A261144, A266755, A344291, A344293, A344294.

Select[Range[1000],EvenQ[Total[Cases[FactorInteger[#],{p_,k_}:>k*PrimePi[p]]]]&&Max@@First/@FactorInteger[#]<=Prime[3]&]

Intersection of A051037 and A300061.

A339360 Sum of all squarefree numbers with greatest prime factor prime(n).

Original entry on oeis.org

1, 2, 9, 60, 504, 6336, 89856, 1645056, 33094656, 801239040, 24246190080, 777550233600, 29697402470400, 1250501433753600, 55083063155097600, 2649111037319577600, 143390180403000115200, 8619643674791667302400, 534710099148093259776000, 36412881178052121329664000

Offset: 0

Gus Wiseman, Dec 04 2020

nonn

			The initial terms are:
   1 = 1,
   2 = 2,
   9 = 3 + 6,
  60 = 5 + 10 + 15 + 30.

Robert Israel, Table of n, a(n) for n = 0..349

A010036 takes prime indices here to binary indices, row sums of A209862.
A048672 takes prime indices to binary indices in squarefree numbers.
A054640 divides the n-th term by prime(n), row sums of A261144.
A072047 counts prime factors of squarefree numbers.
A339194 is the restriction to semiprimes, row sums of A339116.
A339195 has this as row sums.
A002110 lists primorials.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A056239 is the sum of prime indices of n (Heinz weight).
A246867 groups squarefree numbers by weight, with row sums A147655.
A319246 is the sum of prime indices of the n-th squarefree number.
A319247 lists reversed prime indices of squarefree numbers.
A329631 lists prime indices of squarefree numbers.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
Cf. A001221, A014466, A098350, A168472, A282935, A302590, A326878, A338901.

f:= proc(n) local i;
  `if`(n=0, 1, ithprime(n)) *mul(1+ithprime(i),i=1..n-1)
end proc:
map(f, [$0..20]); # Robert Israel, Dec 08 2020

Table[Sum[Times@@Prime/@stn,{stn,Select[Subsets[Range[n]],MemberQ[#,n]&]}],{n,10}]

For n >= 1, a(n) = A054640(n-1) * prime(n).

a(0)=1 prepended by Alois P. Heinz, Jan 08 2025

A344295 Heinz numbers of partitions of 2*n with at most n parts, none greater than 3, for some n.

Original entry on oeis.org

1, 3, 9, 10, 25, 27, 30, 75, 81, 90, 100, 225, 243, 250, 270, 300, 625, 675, 729, 750, 810, 900, 1000, 1875, 2025, 2187, 2250, 2430, 2500, 2700, 3000, 5625, 6075, 6250, 6561, 6750, 7290, 7500, 8100, 9000, 10000, 15625, 16875, 18225, 18750, 19683, 20250, 21870

Offset: 1

Gus Wiseman, May 15 2021

nonn

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

			The sequence of terms together with their prime indices begins:
      1: {}
      3: {2}
      9: {2,2}
     10: {1,3}
     25: {3,3}
     27: {2,2,2}
     30: {1,2,3}
     75: {2,3,3}
     81: {2,2,2,2}
     90: {1,2,2,3}
    100: {1,1,3,3}
    225: {2,2,3,3}
    243: {2,2,2,2,2}
    250: {1,3,3,3}
    270: {1,2,2,2,3}
    300: {1,1,2,3,3}

These partitions are counted by A001399.
Allowing any number of parts and sum gives A051037.
Allowing parts > 3 and any length gives A300061.
Not requiring the sum of prime indices to be even gives A344293.
Allowing any number of parts (but still with even sum) gives A344297.
Allowing parts > 3 gives A344413.
A001358 lists semiprimes.
A025065 counts partitions of n with at least n/2 parts, ranked by A344296.
A035363 counts partitions of n of length n/2, ranked by A340387.
A056239 adds up prime indices, row sums of A112798.
A110618 counts partitions of n with at most n/2 parts, ranked by A344291.
A344414 counts partitions of n with all parts >= n/2, ranked by A344296.
Cf. A030229, A080193, A244990, A261144, A266755, A279622, A334433, A344294.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[1000],EvenQ[Total[primeMS[#]]]&&PrimeOmega[#]<=Total[primeMS[#]]/2&&Max@@primeMS[#]<=3&]

Intersection of A300061 (even Heinz weight), A344291 (Omega > half Heinz weight), and A051037 (5-smooth).

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

A344296 Numbers with at least as many prime factors (counted with multiplicity) as half their sum of prime indices.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A209862 Permutation of nonnegative integers which maps A209642 into ascending order (A209641).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A339195 Triangle of squarefree numbers grouped by greatest prime factor, read by rows.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A344293 5-smooth numbers n whose sum of prime indices A056239(n) is at least twice the number of prime indices A001222(n).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A111059 a(n) = Product_{k=1..n} A005117(k), the product of the first n squarefree positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A344294 5-smooth but not 3-smooth numbers k such that A056239(k) >= 2*A001222(k).

Views

Author