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

A266755 Expansion of 1/((1-x^2)(1-x^3)(1-x^4)).

Original entry on oeis.org

1, 0, 1, 1, 2, 1, 3, 2, 4, 3, 5, 4, 7, 5, 8, 7, 10, 8, 12, 10, 14, 12, 16, 14, 19, 16, 21, 19, 24, 21, 27, 24, 30, 27, 33, 30, 37, 33, 40, 37, 44, 40, 48, 44, 52, 48, 56, 52, 61, 56, 65, 61, 70, 65, 75, 70, 80, 75, 85, 80, 91, 85, 96, 91, 102, 96, 108, 102, 114, 108, 120, 114, 127, 120, 133, 127, 140, 133, 147, 140, 154, 147, 161, 154, 169

Offset: 0

N. J. A. Sloane, Jan 10 2016

This is the same as A005044 but without the three leading zeros. There are so many situations where one wants this sequence rather than A005044 that it seems appropriate for it to have its own entry.
But see A005044 (still the main entry) for numerous applications and references.
Also, Molien series for invariants of finite Coxeter group D_3.
The Molien series for the finite Coxeter group of type D_k (k >= 3) has g.f. = 1/Product_i (1-x^(1+m_i)) where the m_i are [1,3,5,...,2k-3,k-1]. If k is even only even powers of x appear, and we bisect the sequence.
Also, Molien series for invariants of finite Coxeter group A_3. The Molien series for the finite Coxeter group of type A_k (k >= 1) has g.f. = 1/Product_{i=2..k+1} (1-x^i). Note that this is the root system A_k not the alternating group Alt_k.
a(n) is the number of partitions of n into parts 2, 3, and 4. - Joerg Arndt, Apr 16 2017
From Gus Wiseman, May 23 2021: (Start)
Also the number of integer partitions of n into at most n/2 parts, none greater than 3. The case of any maximum is A110618. The case of any length is A001399. The Heinz numbers of these partitions are given by A344293.
For example, the a(2) = 1 through a(13) = 5 partitions are:
    2  3  22  32  33   322  332   333   3322   3332   3333    33322
          31      222  331  2222  3222  3331   32222  33222   33331
                  321       3221  3321  22222  33221  33321   322222
                            3311        32221  33311  222222  332221
                                        33211         322221  333211
                                                      332211
                                                      333111
(End)

			G.f. = 1 + x^2 + x^3 + 2*x^4 + x^5 + 3*x^6 + 2*x^7 + 4*x^8 + ... - _Michael Somos_, Jan 29 2022

J. E. Humphreys, Reflection Groups and Coxeter Groups, Cambridge, 1990. See Table 3.1, page 59.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Sara C. Billey, Matjaž Konvalinka, and Joshua P. Swanson, Asymptotic normality of the major index on standard tableaux, arXiv:1905.00975 [math.CO], 2019.
Index entries for Molien series
Index entries for linear recurrences with constant coefficients, signature (0,1,1,1,-1,-1,-1,0,1).

Molien series for finite Coxeter groups A_1 through A_12 are A059841, A103221, A266755, A008667, A037145, A001996, and A266776-A266781.
Molien series for finite Coxeter groups D_3 through D_12 are A266755, A266769, A266768, A003402, and A266770-A266775.
A variant of A005044.
Cf. A001400 (partial sums).
Cf. A308065.
Number of partitions of n whose Heinz number is in A344293.
A001399 counts partitions with all parts <= 3, ranked by A051037.
A025065 counts partitions of n with >= n/2 parts, ranked by A344296.
A035363 counts partitions of n with n/2 parts, ranked by A340387.
A110618 counts partitions of n into at most n/2 parts, ranked by A344291.
Cf. A000041, A008642, A279622, A325691, A344294, A344297.

I:=[1,0,1,1,2,1,3,2,4]; [n le 9 select I[n] else Self(n-2)+ Self(n-3)+Self(n-4)-Self(n-5)-Self(n-6)-Self(n-7)+Self(n-9): n in [1..100]]; // Vincenzo Librandi, Jan 11 2016

CoefficientList[Series[1/((1-x^2)(1-x^3)(1-x^4)), {x, 0, 100}], x] (* JungHwan Min, Jan 10 2016 *)
LinearRecurrence[{0,1,1,1,-1,-1,-1,0,1}, {1,0,1,1,2,1,3,2,4}, 100] (* Vincenzo Librandi, Jan 11 2016 *)
Table[Length[Select[IntegerPartitions[n],Length[#]<=n/2&&Max@@#<=3&]],{n,0,30}] (* Gus Wiseman, May 23 2021 *)
a[ n_] := Round[(n + 3*(2 - Mod[n,2]))^2/48]; (* Michael Somos, Jan 29 2022 *)

Vec(1/((1-x^2)*(1-x^3)*(1-x^4)) + O(x^100)) \\ Michel Marcus, Jan 11 2016

{a(n) = round((n + 3*(2-n%2))^2/48)}; /* Michael Somos, Jan 29 2022 */

(1/((1-x^2)*(1-x^3)*(1-x^4))).series(x, 100).coefficients(x, sparse=False) # G. C. Greubel, Jun 13 2019

a(n) = a(n-2) + a(n-3) + a(n-4) - a(n-5) - a(n-6) - a(n-7) + a(n-9) for n>8. - Vincenzo Librandi, Jan 11 2016
a(n) = a(-9-n) for all n in Z. a(n) = a(n+3) for all n in 2Z. - Michael Somos, Jan 29 2022
E.g.f.: exp(-x)*(81 - 18*x + exp(2*x)*(107 + 60*x + 6*x^2) + 64*exp(x/2)*cos(sqrt(3)*x/2) + 36*exp(x)*(cos(x) - sin(x)))/288. - Stefano Spezia, Mar 05 2023
For n >= 3, if n is even, a(n) = a(n-3) + floor(n/4) + 1, otherwise a(n) = a(n-3). - Robert FERREOL, Feb 05 2024
a(n) = floor((n^2+9*n+(3*n+9)*(-1)^n+39)/48). - Hoang Xuan Thanh, Jun 03 2025

A344415 Numbers whose greatest prime index is half their sum of prime indices.

Original entry on oeis.org

4, 9, 12, 25, 30, 40, 49, 63, 70, 84, 112, 121, 154, 165, 169, 198, 220, 264, 273, 286, 289, 325, 351, 352, 361, 364, 390, 442, 468, 520, 529, 561, 595, 624, 646, 714, 741, 748, 765, 832, 841, 850, 874, 918, 931, 952, 961, 988, 1020, 1045, 1173, 1197, 1224

Offset: 1

Gus Wiseman, May 19 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.

			The sequence of terms together with their prime indices begins:
       4: {1,1}           198: {1,2,2,5}
       9: {2,2}           220: {1,1,3,5}
      12: {1,1,2}         264: {1,1,1,2,5}
      25: {3,3}           273: {2,4,6}
      30: {1,2,3}         286: {1,5,6}
      40: {1,1,1,3}       289: {7,7}
      49: {4,4}           325: {3,3,6}
      63: {2,2,4}         351: {2,2,2,6}
      70: {1,3,4}         352: {1,1,1,1,1,5}
      84: {1,1,2,4}       361: {8,8}
     112: {1,1,1,1,4}     364: {1,1,4,6}
     121: {5,5}           390: {1,2,3,6}
     154: {1,4,5}         442: {1,6,7}
     165: {2,3,5}         468: {1,1,2,2,6}
     169: {6,6}           520: {1,1,1,3,6}

The partitions with these Heinz numbers are counted by A035363.
The conjugate version is A340387.
This sequence is the case of equality in A344414 and A344416.
A001222 counts prime factors with multiplicity.
A025065 counts palindromic partitions, ranked by A265640.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A058696 counts partitions of even numbers, ranked by A300061.
A301987 lists numbers whose sum of prime indices equals their product.
A322109 ranks partitions of n with no part > n/2, counted by A110618.
A334201 adds up all prime indices except the greatest.
A344291 lists numbers m with A001222(m) <= A056239(m)/2, counted by A110618.
A344296 lists numbers m with A001222(m) >= A056239(m)/2, counted by A025065.
Cf. A000070, A001414, A209816, A301988, A316413, A316428, A320924, A325037, A325038, A325044, A330950, A344293, A344294, A344297.

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

A061395(a(n)) = A056239(a(n))/2.

A344414 Heinz numbers of integer partitions whose sum is at most twice their greatest part.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 17, 19, 20, 21, 22, 23, 25, 26, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 40, 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

Offset: 1

Gus Wiseman, May 19 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:
     2: {1}        20: {1,1,3}    39: {2,6}
     3: {2}        21: {2,4}      40: {1,1,1,3}
     4: {1,1}      22: {1,5}      41: {13}
     5: {3}        23: {9}        42: {1,2,4}
     6: {1,2}      25: {3,3}      43: {14}
     7: {4}        26: {1,6}      44: {1,1,5}
     9: {2,2}      28: {1,1,4}    46: {1,9}
    10: {1,3}      29: {10}       47: {15}
    11: {5}        30: {1,2,3}    49: {4,4}
    12: {1,1,2}    31: {11}       51: {2,7}
    13: {6}        33: {2,5}      52: {1,1,6}
    14: {1,4}      34: {1,7}      53: {16}
    15: {2,3}      35: {3,4}      55: {3,5}
    17: {7}        37: {12}       56: {1,1,1,4}
    19: {8}        38: {1,8}      57: {2,8}
For example, 56 has prime indices {1,1,1,4} and 7 <= 2*4, so 56 is in the sequence. On the other hand, 224 has prime indices {1,1,1,1,1,4} and 9 > 2*4, so 224 is not in the sequence.

These partitions are counted by A025065 but are different from palindromic partitions, which have Heinz numbers A265640.
The opposite even-weight version appears to be A320924, counted by A209816.
The opposite version appears to be A322109, counted by A110618.
The case of equality in the conjugate version is A340387.
The conjugate opposite version is A344291, counted by A110618.
The conjugate opposite 5-smooth case is A344293, counted by A266755.
The conjugate version is A344296, also counted by A025065.
The case of equality is A344415.
The even-weight case is A344416.
A001222 counts prime factors with multiplicity.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A058696 counts partitions of even numbers, ranked by A300061.
A301987 lists numbers whose sum of prime indices equals their product.
A334201 adds up all prime indices except the greatest.
Cf. A001414, A067538, A301988, A316413, A316428, A325037, A325038, A325044, A330950, A344294, A344297.

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

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

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.

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.

A344416 Heinz numbers of integer partitions whose sum is even and is at most twice the greatest part.

Original entry on oeis.org

3, 4, 7, 9, 10, 12, 13, 19, 21, 22, 25, 28, 29, 30, 34, 37, 39, 40, 43, 46, 49, 52, 53, 55, 57, 61, 62, 63, 66, 70, 71, 76, 79, 82, 84, 85, 87, 88, 89, 91, 94, 101, 102, 107, 111, 112, 113, 115, 116, 117, 118, 121, 129, 130, 131, 133, 134, 136, 138, 139, 146

Offset: 1

Gus Wiseman, May 20 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.

Also numbers m whose sum of prime indices A056239(m) is even and is at most twice the greatest prime index A061395(m).

			The sequence of terms together with their prime indices begins:
      3: {2}         37: {12}          71: {20}
      4: {1,1}       39: {2,6}         76: {1,1,8}
      7: {4}         40: {1,1,1,3}     79: {22}
      9: {2,2}       43: {14}          82: {1,13}
     10: {1,3}       46: {1,9}         84: {1,1,2,4}
     12: {1,1,2}     49: {4,4}         85: {3,7}
     13: {6}         52: {1,1,6}       87: {2,10}
     19: {8}         53: {16}          88: {1,1,1,5}
     21: {2,4}       55: {3,5}         89: {24}
     22: {1,5}       57: {2,8}         91: {4,6}
     25: {3,3}       61: {18}          94: {1,15}
     28: {1,1,4}     62: {1,11}       101: {26}
     29: {10}        63: {2,2,4}      102: {1,2,7}
     30: {1,2,3}     66: {1,2,5}      107: {28}
     34: {1,7}       70: {1,3,4}      111: {2,12}

These partitions are counted by A000070 = even-indexed terms of A025065.
The opposite version appears to be A320924, counted by A209816.
The opposite version with odd weights allowed appears to be A322109.
The conjugate opposite version allowing odds is A344291, counted by A110618.
The conjugate version is A344296, also counted by A025065.
The conjugate opposite version is A344413, counted by A209816.
Allowing odd weight gives A344414.
The case of equality is A344415, counted by A035363.
A001222 counts prime factors with multiplicity.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A058696 counts partitions of even numbers, ranked by A300061.
A265640 lists Heinz numbers of palindromic partitions.
A301987 lists numbers whose sum of prime indices equals their product.
A334201 adds up all prime indices except the greatest.
A340387 lists Heinz numbers of partitions whose sum is twice their length.
Cf. A001414, A074761, A316413, A316428, A325037, A325038, A325044, A330950, A344293, A344294, A344297.

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

Intersection of A300061 and A344414.

A344413 Numbers n whose sum of prime indices A056239(n) is even and is at least twice the number of prime factors A001222(n).

Original entry on oeis.org

1, 3, 7, 9, 10, 13, 19, 21, 22, 25, 27, 28, 29, 30, 34, 37, 39, 43, 46, 49, 52, 53, 55, 57, 61, 62, 63, 66, 70, 71, 75, 76, 79, 81, 82, 84, 85, 87, 88, 89, 90, 91, 94, 100, 101, 102, 107, 111, 113, 115, 116, 117, 118, 121, 129, 130, 131, 133, 134, 136, 138

Offset: 1

Gus Wiseman, May 19 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.

Also Heinz numbers of integer partitions of even numbers m with at most m/2 parts, counted by A209816 riffled with zeros, or A110618 with odd positions zeroed out.

			The sequence of terms together with their prime indices begins:
      1: {}          37: {12}        75: {2,3,3}
      3: {2}         39: {2,6}       76: {1,1,8}
      7: {4}         43: {14}        79: {22}
      9: {2,2}       46: {1,9}       81: {2,2,2,2}
     10: {1,3}       49: {4,4}       82: {1,13}
     13: {6}         52: {1,1,6}     84: {1,1,2,4}
     19: {8}         53: {16}        85: {3,7}
     21: {2,4}       55: {3,5}       87: {2,10}
     22: {1,5}       57: {2,8}       88: {1,1,1,5}
     25: {3,3}       61: {18}        89: {24}
     27: {2,2,2}     62: {1,11}      90: {1,2,2,3}
     28: {1,1,4}     63: {2,2,4}     91: {4,6}
     29: {10}        66: {1,2,5}     94: {1,15}
     30: {1,2,3}     70: {1,3,4}    100: {1,1,3,3}
     34: {1,7}       71: {20}       101: {26}
For example, 75 has 3 prime indices {2,3,3} with sum 8 >= 2*3, so 75 is in the sequence.

Robert Israel, Table of n, a(n) for n = 1..10000

These are the Heinz numbers of partitions counted by A209816 and A110618.
A subset of A300061 (sum of prime indices is even).
The conjugate version appears to be A320924 (allowing odd weights: A322109).
The case of equality is A340387.
Allowing odd weights gives A344291.
The 5-smooth case is A344295, or A344293 allowing odd weights.
The opposite version allowing odd weights is A344296.
The conjugate opposite version allowing odd weights is A344414.
The case of equality in the conjugate case is A344415.
The conjugate opposite version is A344416, counted by A000070.
A001222 counts prime factors with multiplicity.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A058696 counts partitions of even numbers, ranked by A300061.
A301987 lists numbers whose sum of prime indices equals their product.
A330950 counts partitions of n with Heinz number divisible by n.
A334201 adds up all prime indices except the greatest.
Cf. A001414, A025065, A067538, A301988, A316413, A316428, A325037, A325038, A325044, A338914, A344294, A344297.

filter:= proc(n) local F,a,t;
  F:= ifactors(n)[2];
  a:= add((numtheory:-pi(t[1])-2)*t[2],t=F);
  a::even and a >= 0
end proc:
select(filter, [$1..300]); # Robert Israel, Oct 10 2024

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

Members m of A300061 such that A056239(m) >= 2*A001222(m).

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

A344292 Numbers m whose sum of prime indices A056239(m) is even and is at most twice the number of prime factors counted with multiplicity A001222(m).

Original entry on oeis.org

1, 3, 4, 9, 10, 12, 16, 27, 28, 30, 36, 40, 48, 64, 81, 84, 88, 90, 100, 108, 112, 120, 144, 160, 192, 208, 243, 252, 256, 264, 270, 280, 300, 324, 336, 352, 360, 400, 432, 448, 480, 544, 576, 624, 640, 729, 756, 768, 784, 792, 810, 832, 840, 880, 900, 972

Offset: 1

Gus Wiseman, May 22 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.

Also Heinz numbers of integer partitions of even numbers m with at least m/2 parts, counted by A000070 riffled with 0's, or A025065 with odd positions zeroed out.

			The sequence of terms together with their prime indices begins:
       1: {}                 84: {1,1,2,4}
       3: {2}                88: {1,1,1,5}
       4: {1,1}              90: {1,2,2,3}
       9: {2,2}             100: {1,1,3,3}
      10: {1,3}             108: {1,1,2,2,2}
      12: {1,1,2}           112: {1,1,1,1,4}
      16: {1,1,1,1}         120: {1,1,1,2,3}
      27: {2,2,2}           144: {1,1,1,1,2,2}
      28: {1,1,4}           160: {1,1,1,1,1,3}
      30: {1,2,3}           192: {1,1,1,1,1,1,2}
      36: {1,1,2,2}         208: {1,1,1,1,6}
      40: {1,1,1,3}         243: {2,2,2,2,2}
      48: {1,1,1,1,2}       252: {1,1,2,2,4}
      64: {1,1,1,1,1,1}     256: {1,1,1,1,1,1,1,1}
      81: {2,2,2,2}         264: {1,1,1,2,5}

These are the Heinz numbers of partitions counted by A000070 and A025065.
A subset of A300061 (sum of prime indices is even).
The conjugate opposite version is A320924, counted by A209816.
The conjugate opposite version allowing odds is A322109, counted by A110618.
The case of equality is A340387, counted by A000041.
The opposite version allowing odd weights is A344291, counted by A110618.
Allowing odd weights gives A344296, counted by A025065.
The opposite version is A344413, counted by A209816.
The conjugate version allowing odd weights is A344414, counted by A025065.
The case of equality in the conjugate case is A344415, counted by A035363.
The conjugate version is A344416, counted by A000070.
A001222 counts prime factors with multiplicity.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A058696 counts partitions of even numbers, ranked by A300061.
A301987 lists numbers whose sum of prime indices equals their product.
A330950 counts partitions of n with Heinz number divisible by n.
A334201 adds up all prime indices except the greatest.
Cf. A001414, A067538, A316413, A316428, A325037, A325038, A325044, A338914, A344294, A344297.

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

Members m of A300061 such that A056239(m) <= 2*A001222(m).

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

A266755 Expansion of 1/((1-x^2)*(1-x^3)*(1-x^4)).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Sage

Formula

A344415 Numbers whose greatest prime index is half their sum of prime indices.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A344414 Heinz numbers of integer partitions whose sum is at most twice their greatest part.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

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

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A344416 Heinz numbers of integer partitions whose sum is even and is at most twice the greatest part.

A266755 Expansion of 1/((1-x^2)(1-x^3)(1-x^4)).