A112798 -id:A112798 - cp's OEIS frontend

A344616 Alternating sum of the integer partition with Heinz number n.

0, 1, 2, 0, 3, 1, 4, 1, 0, 2, 5, 2, 6, 3, 1, 0, 7, 1, 8, 3, 2, 4, 9, 1, 0, 5, 2, 4, 10, 2, 11, 1, 3, 6, 1, 0, 12, 7, 4, 2, 13, 3, 14, 5, 3, 8, 15, 2, 0, 1, 5, 6, 16, 1, 2, 3, 6, 9, 17, 1, 18, 10, 4, 0, 3, 4, 19, 7, 7, 2, 20, 1, 21, 11, 2, 8, 1, 5, 22, 3, 0, 12

Offset: 1

Gus Wiseman, Jun 03 2021

nonn

The alternating sum of a partition (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i, which is equal to the number of odd parts in the conjugate partition.
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 the reverse-alternating sum of the prime indices of n.

			The partition (6,4,3,2,2) has Heinz number 4095 and conjugate (5,5,3,2,1,1), so a(4095) = 5.

Alois P. Heinz, Table of n, a(n) for n = 1..20000

Positions of nonzeros are A000037.
Positions of 0's are A000290.
The version for prime factors is A071321 (reverse: A071322).
A version for compositions is A124754.
The version for prime multiplicities is A316523.
The reverse version is A316524, with sign A344617.
A000041 counts partitions of 2n with alternating sum 0.
A056239 adds up prime indices, row sums of A112798.
A103919 counts partitions by sum and alternating sum.
A335433 ranks separable partitions.
A335448 ranks inseparable partitions.
A344606 counts wiggly permutations of prime indices with twins.
A344610 counts partitions by sum and positive reverse-alternating sum.
A344612 counts partitions by sum and reverse-alternating sum.
A344618 gives reverse-alternating sums of standard compositions.
Cf. A000070, A001222, A026424, A028260, A116406, A119899, A343938, A344607, A344608, A344609, A344619, A344653, A344739.

a:= n-> (l-> -add(l[i]*(-1)^i, i=1..nops(l)))(sort(map(
    i-> numtheory[pi](i[1])$i[2], ifactors(n)[2]), `>`)):
seq(a(n), n=1..82);  # Alois P. Heinz, Jun 04 2021

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
Table[ats[Reverse[primeMS[n]]],{n,100}]

a(n) = A257991(A122111(n)).

A057427(a(n)) = A049240(n).

A335448 Numbers whose prime indices are inseparable.

Original entry on oeis.org

4, 8, 9, 16, 24, 25, 27, 32, 40, 48, 49, 54, 56, 64, 80, 81, 88, 96, 104, 112, 121, 125, 128, 135, 136, 144, 152, 160, 162, 169, 176, 184, 189, 192, 208, 224, 232, 240, 243, 248, 250, 256, 272, 288, 289, 296, 297, 304, 320, 324, 328, 336, 343, 344, 351, 352

Offset: 1

Gus Wiseman, Jun 21 2020

nonn

First differs from A212164 in lacking 72.
First differs from A293243 in lacking 72.
No terms are squarefree.
Also Heinz numbers of inseparable partitions (A325535). The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
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 also numbers that can be written as a product of prime numbers, each different from the last but not necessarily different from those prior to the last.
A multiset is inseparable iff its maximal multiplicity is greater than one plus the sum of its remaining multiplicities.

			The sequence of terms together with their prime indices begins:
   4: {1,1}
   8: {1,1,1}
   9: {2,2}
  16: {1,1,1,1}
  24: {1,1,1,2}
  25: {3,3}
  27: {2,2,2}
  32: {1,1,1,1,1}
  40: {1,1,1,3}
  48: {1,1,1,1,2}
  49: {4,4}
  54: {1,2,2,2}
  56: {1,1,1,4}
  64: {1,1,1,1,1,1}
  80: {1,1,1,1,3}
  81: {2,2,2,2}
  88: {1,1,1,5}
  96: {1,1,1,1,1,2}

Complement of A335433.
Separations are counted by A003242 and A335452 and ranked by A333489.
Permutations of prime indices are counted by A008480.
Inseparable partitions are counted by A325535.
Strict permutations of prime indices are counted by A335489.
Cf. A000670, A000961, A005117, A056239, A112798, A181796, A261962, A333221, A335451.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],Select[Permutations[primeMS[#]],!MatchQ[#,{_,x_,x_,_}]&]=={}&]

A335433 Numbers whose multiset of prime indices is separable.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 02 2020

nonn

First differs from A212167 in having 72.
Includes all squarefree numbers A005117.
A multiset is separable if it has a permutation that is an anti-run, meaning there are no adjacent equal parts.
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 separable partitions (A325534). The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Also numbers that cannot be written as a product of prime numbers, each different from the last but not necessarily different from those prior to the last.

			The sequence of terms together with their prime indices begins:
      1: {}          20: {1,1,3}       39: {2,6}
      2: {1}         21: {2,4}         41: {13}
      3: {2}         22: {1,5}         42: {1,2,4}
      5: {3}         23: {9}           43: {14}
      6: {1,2}       26: {1,6}         44: {1,1,5}
      7: {4}         28: {1,1,4}       45: {2,2,3}
     10: {1,3}       29: {10}          46: {1,9}
     11: {5}         30: {1,2,3}       47: {15}
     12: {1,1,2}     31: {11}          50: {1,3,3}
     13: {6}         33: {2,5}         51: {2,7}
     14: {1,4}       34: {1,7}         52: {1,1,6}
     15: {2,3}       35: {3,4}         53: {16}
     17: {7}         36: {1,1,2,2}     55: {3,5}
     18: {1,2,2}     37: {12}          57: {2,8}
     19: {8}         38: {1,8}         58: {1,10}

The version for a multiset with prescribed multiplicities is A335127.
Separable factorizations are counted by A335434.
The complement is A335448.
Separations are counted by A003242 and A335452 and ranked by A333489.
Permutations of prime indices are counted by A008480.
Inseparable partitions are counted by A325535.
Strict permutations of prime indices are counted by A335489.
Cf. A000961, A005117, A056239, A112798, A114938, A181796, A292884, A335407, A335451, A335516.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],Select[Permutations[primeMS[#]],!MatchQ[#,{_,x_,x_,_}]&]!={}&]

A181821 a(n) = smallest integer with factorization as Product p(i)^e(i) such that Product p(e(i)) = n.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 30, 36, 24, 32, 60, 64, 48, 72, 210, 128, 180, 256, 120, 144, 96, 512, 420, 216, 192, 900, 240, 1024, 360, 2048, 2310, 288, 384, 432, 1260, 4096, 768, 576, 840, 8192, 720, 16384, 480, 1800, 1536, 32768, 4620, 1296, 1080, 1152, 960, 65536

Offset: 1

Matthew Vandermast, Dec 07 2010

nonn

A permutation of A025487. a(n) is the member m of A025487 such that A181819(m) = n. a(n) is also the member of A025487 whose prime signature is conjugate to the prime signature of A108951(n).

If n = Product_i prime(e(i)) with the e(i) weakly decreasing, then a(n) = Product_i prime(i)^e(i). For example, 90 = prime(3) * prime(2) * prime(2) * prime(1), so a(90) = prime(1)^3 * prime(2)^2 * prime(3)^2 * prime(4)^1 = 12600. - Gus Wiseman, Jan 02 2019

			The canonical factorization of 24 is 2^3*3^1. Therefore, p(e(i)) = prime(3)*prime(1)(i.e., A000040(3)*A000040(1)), which equals 5*2 = 10. Since 24 is the smallest integer for which p(e(i)) = 10, a(10) = 24.

Alois P. Heinz, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Conjugate Partition

Other rearrangements of A025487 include A036035, A059901, A063008, A077569, A085988, A086141, A087443, A108951, A181822.
Cf. A046523, A181819, A181820.
Cf. A000040, A001221, A001222, A002110, A056239, A071625, A112798, A118914, A122111, A124859, A182850, A305936, A361808.

a:= n-> (l-> mul(ithprime(i)^l[i], i=1..nops(l)))(sort(map(i->
             numtheory[pi](i[1])$i[2], ifactors(n)[2]), `>`)):
seq(a(n), n=1..70);  # Alois P. Heinz, Sep 05 2018

With[{s = Array[If[# == 1, 1, Times @@ Map[Prime@ Last@ # &, FactorInteger@ #]] &, 2^16]}, Array[First@ FirstPosition[s, #] &, LengthWhile[Differences@ Union@ s, # == 1 &]]] (* Michael De Vlieger, Dec 17 2018 *)
Table[Times@@MapIndexed[Prime[#2[[1]]]^#1&,Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]],{n,30}] (* Gus Wiseman, Jan 02 2019 *)

A181821(n) = { my(f=factor(n),p=0,m=1); forstep(i=#f~,1,-1,while(f[i,2], f[i,2]--; m *= (p=nextprime(p+1))^primepi(f[i,1]))); (m); }; \\ Antti Karttunen, Dec 10 2018

from math import prod
from sympy import prime, primepi, factorint
def A181821(n): return prod(prime(i)**e for i, e in enumerate(sorted(map(primepi,factorint(n,multiple=True)),reverse=True),1)) # Chai Wah Wu, Sep 15 2023

If A108951(n) = Product p(i)^e(i), then a(n) = Product A002110(e(i)). I.e., a(n) = A108951(A181819(A108951(n))).
a(A181819(n)) = A046523(n). - [See also A124859]. Antti Karttunen, Dec 10 2018
a(n) = A025487(A361808(n)). - Pontus von Brömssen, Mar 25 2023
a(n) = A108951(A122111(n)). - Antti Karttunen, Sep 15 2023

Definition corrected by Gus Wiseman, Jan 02 2019

A106529 Numbers having k prime factors (counted with multiplicity), the largest of which is the k-th prime.

Original entry on oeis.org

2, 6, 9, 20, 30, 45, 50, 56, 75, 84, 125, 126, 140, 176, 189, 196, 210, 264, 294, 315, 350, 396, 416, 440, 441, 490, 525, 594, 616, 624, 660, 686, 735, 875, 891, 924, 936, 968, 990, 1029, 1040, 1088, 1100, 1225, 1386, 1404, 1452, 1456, 1485, 1540, 1560

Offset: 1

Matthew Ryan (mattryan1994(AT)hotmail.com), May 30 2005

nonn

It seems that the ratio between successive terms tends to 1 as n increases, meaning perhaps that most numbers are in this sequence.
The number of terms that have the k-th prime as their largest prime factor is A000984(k), the k-th central binomial coefficient. E.g., 6 and 9 are the A000984(2)=2 terms in {a(n)} that have prime(2)=3 as their largest prime factor.
The sequence contains the positive integers m such that the rank of the partition B(m) = 0. For m >= 2, B(m) is defined as the partition obtained by taking the prime decomposition of m and replacing each prime factor p with its index i (i.e., i-th prime = p); also B(1) = the empty partition. For example, B(350) = B(2*5^2*7) = [1,3,3,4]. B is a bijection between the positive integers and the set of all partitions. The rank of a partition P is the largest part of P minus the number of parts of P. - Emeric Deutsch, May 09 2015
Also Heinz numbers of balanced partitions, counted by A047993. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Feb 08 2021

			a(7)=50 because 50=2*5*5 is, for k=3, the product of k primes, the largest of which is the k-th prime, and 50 is the 7th such number.

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Cf. A000984.
A001222 counts prime factors.
A056239 adds up prime indices.
A061395 selects maximum prime index.
A112798 lists the prime indices of each positive integer.
Other balance-related sequences:
- A010054 counts balanced strict partitions.
- A047993 counts balanced partitions.
- A090858 counts partitions of rank 1.
- A098124 counts balanced compositions.
- A340596 counts co-balanced factorizations.
- A340598 counts balanced set partitions.
- A340599 counts alt-balanced factorizations.
- A340600 counts unlabeled balanced multiset partitions.
- A340653 counts balanced factorizations.
Cf. A006141, A064174, A096401, A117409, A168659, A324522, A340654, A340655, A340656, A340657.

with(numtheory): a := proc (n) options operator, arrow: pi(max(factorset(n)))-bigomega(n) end proc: A := {}: for i from 2 to 1600 do if a(i) = 0 then A := `union`(A, {i}) else  end if end do: A; # Emeric Deutsch, May 09 2015

Select[Range@ 1560, PrimePi@ FactorInteger[#][[-1, 1]] == PrimeOmega@ # &] (* Michael De Vlieger, May 09 2015 *)

For all terms, A001222(a(n)) = A061395(a(n)). - Gus Wiseman, Feb 08 2021

A182850 a(n) = number of iterations that n requires to reach a fixed point under the x -> A181819(x) map.

Original entry on oeis.org

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

Offset: 1

Matthew Vandermast, Jan 04 2011

nonn

The fixed points of the x -> A181819(x) map are 1 and 2. Note that the x -> A000005(x) map has the same fixed points, and that A000005(n) = A181819(n) iff n is cubefree (cf. A004709). Under the x -> A181819(x) map, it seems significantly easier to generalize about which kinds of integers take a given number of iterations to reach a fixed point than under the x -> A000005(x) map.

Also the number of steps in the reduction of the multiset of prime factors of n wherein one repeatedly takes the multiset of multiplicities. For example, the a(90) = 5 steps are {2,3,3,5} -> {1,1,2} -> {1,2} -> {1,1} -> {2} -> {1}. - Gus Wiseman, May 13 2018

			A181819(6) = 4; A181819(4) = 3; A181819(3) = 2; A181819(2) = 2. Therefore, a(6) = 3, a(4) = 2, a(3)= 1, and a(2) = 0.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Fixed Point
Eric Weisstein's World of Mathematics, Map

A182857 gives values of n where a(n) increases to a record.

Cf. A000961, A001222, A003434, A005117, A007916, A036459, A112798, A130091, A181819, A182851-A182858, A238748, A304455, A304464, A304465.

a182850 n = length $ takeWhile (`notElem` [1,2]) $ iterate a181819 n
-- Reinhard Zumkeller, Mar 26 2012

Table[If[n<=2,0,Length[FixedPointList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]]]]-1],{n,100}] (* Gus Wiseman, May 13 2018 *)

;; With memoization-macro definec.
(definec (A182850 n) (if (<= n 2) 0 (+ 1 (A182850 (A181819 n))))) ;; Antti Karttunen, Feb 05 2016

For n > 2, a(n) = a(A181819(n)) + 1.
a(n) = 0 iff n equals 1 or 2.
a(n) = 1 iff n is an odd prime (A000040(n) for n>1).
a(n) = 2 iff n is a composite perfect prime power (A025475(n) for n>1).
a(n) = 3 iff n is a squarefree composite integer or a power of a squarefree composite integer (cf. A182853).
a(n) = 4 iff n's prime signature a) contains at least two distinct numbers, and b) contains no number that occurs less often than any other number (cf. A182854).

A326568 Denominator of the average of the multiset of prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 2, 2, 1, 1, 3, 1, 3, 1, 1, 1, 4, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 3, 1, 3, 3, 1, 1, 5, 1, 3, 2, 3, 1, 4, 1, 4, 1, 2, 1, 4, 1, 1, 3, 1, 2, 3, 1, 1, 2, 3, 1, 5, 1, 2, 3, 3, 2, 1, 1, 5, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 2, 1, 2, 6, 1, 1, 1, 1, 1, 3, 1, 4, 1, 2, 1, 5

Offset: 2

Gus Wiseman, Jul 13 2019

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 prime indices of 12 are {1,1,2}, with average 4/3, so a(12) = 3.

Antti Karttunen, Table of n, a(n) for n = 2..65537
Index entries for sequences related to prime indices in the factorization of n.

a(n) is a divisor of Omega(n) = A001222(n).
Positions of 1's are A316413.
Cf. A001414, A056239, A067629, A112798, A123528/A123529, A289508, A289509, A290103, A316428, A326567, A326620.

Table[Denominator[Sum[q[[2]]*PrimePi[q[[1]]],{q,FactorInteger[n]}]/PrimeOmega[n]],{n,2,100}]

A326568(n) = { my(f=factor(n)); denominator(sum(i=1,#f~,f[i,2]*primepi(f[i,1]))/bigomega(n)); }; \\ Antti Karttunen, Jan 28 2025

Starting offset corrected from 0 to 2 and data section extended to a(108) by Antti Karttunen, Jan 28 2025

A326567 Numerator of the average of the multiset of prime indices of n.

Original entry on oeis.org

1, 2, 1, 3, 3, 4, 1, 2, 2, 5, 4, 6, 5, 5, 1, 7, 5, 8, 5, 3, 3, 9, 5, 3, 7, 2, 2, 10, 2, 11, 1, 7, 4, 7, 3, 12, 9, 4, 3, 13, 7, 14, 7, 7, 5, 15, 6, 4, 7, 9, 8, 16, 7, 4, 7, 5, 11, 17, 7, 18, 6, 8, 1, 9, 8, 19, 3, 11, 8, 20, 7, 21, 13, 8, 10, 9, 3, 22, 7, 2, 7

Offset: 2

Gus Wiseman, Jul 13 2019

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 prime indices of 12 are {1,1,2}, with average 4/3, so a(12) = 4.

Cf. A001222, A001414, A056239, A067629, A112798, A123528/A123529, A289508, A289509, A290103, A316413, A326568.

Table[Numerator[Sum[q[[2]]*PrimePi[q[[1]]],{q,FactorInteger[n]}]/PrimeOmega[n]],{n,2,100}]

A036036 Triangle read by rows in which row n lists all the parts of all reversed partitions of n, sorted first by length and then lexicographically.

Original entry on oeis.org

1, 2, 1, 1, 3, 1, 2, 1, 1, 1, 4, 1, 3, 2, 2, 1, 1, 2, 1, 1, 1, 1, 5, 1, 4, 2, 3, 1, 1, 3, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 6, 1, 5, 2, 4, 3, 3, 1, 1, 4, 1, 2, 3, 2, 2, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 7, 1, 6, 2, 5, 3, 4, 1, 1, 5, 1, 2, 4, 1, 3, 3, 2, 2, 3, 1, 1, 1

Offset: 1

N. J. A. Sloane

First differs from A334442 for reversed partitions of 9. Namely, this sequence has (1,4,4) before (2,2,5), while A334442 has (2,2,5) before (1,4,4). - Gus Wiseman, May 07 2020
This is the "Abramowitz and Stegun" ordering of the partitions, referenced in numerous other sequences. The partitions are in reverse order of the conjugates of the partitions in Mathematica order (A080577). Each partition is the conjugate of the corresponding partition in Maple order (A080576). - Franklin T. Adams-Watters, Oct 18 2006
The "Abramowitz and Stegun" ordering of the partitions is the graded reflected colexicographic ordering of the partitions. - Daniel Forgues, Jan 19 2011
The "Abramowitz and Stegun" ordering of partitions has been traced back to C. F. Hindenburg, 1779, in the Knuth reference, p. 38. See the Hindenburg link, pp. 77-5 with the listing of the partitions for n=10.  This is also mentioned in the P. Luschny link. - Wolfdieter Lang, Apr 04 2011
The "Abramowitz and Stegun" order used here means that the partitions of a given number are listed by increasing number of (nonzero) parts, then by increasing lexicographical order with parts in (weakly) indecreasing order. This differs from n=9 on from A334442 which considers reverse lexicographic order of parts in (weakly) decreasing order. - M. F. Hasler, Jul 12 2015, corrected thanks to Gus Wiseman, May 14 2020
This is the Abramowitz-Stegun ordering of reversed partitions (finite weakly increasing sequences of positive integers). The same ordering of non-reversed partitions is A334301. - Gus Wiseman, May 07 2020

			1
2; 1,1
3; 1,2; 1,1,1
4; 1,3; 2,2; 1,1,2; 1,1,1,1
5; 1,4; 2,3; 1,1,3; 1,2,2; 1,1,1,2; 1,1,1,1,1;
6; 1,5; 2,4; 3,3; 1,1,4; 1,2,3; 2,2,2; 1,1,1,3; 1,1,2,2; 1,1,1,1,2; 1,1,1,1,1,1;
...

Abramowitz and Stegun, Handbook, p. 831, column labeled "pi".
D. Knuth, The Art of Computer Programming, Vol. 4, fascicle 3, 7.2.1.4, Addison-Wesley, 2005.

Franklin T. Adams-Watters, First 20 rows, flattened.
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy]. (uses Flash)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions.
C. F. Hindenburg, Infinitinomii Dignitatum Exponentis Indeterminati, Goettingen 1779.
Peter Luschny, Counting with partitions.
OEIS Wiki, Orderings of partitions (a comparison).
Wikiversity, Lexicographic and colexicographic order

Cf. A036037-A036040.
See A036037 for the graded colexicographic ordering.
See A080576 for the Maple (graded reflected lexicographic) ordering.
See A080577 for the Mathematica (graded reverse lexicographic) ordering.
See A193073 for the graded lexicographic ordering.
See A228100 for the Fenner-Loizou (binary tree) ordering.
The version ignoring length is A026791.
Same as A036037 with partitions reversed.
The lengths of these partitions are A036043.
The number of distinct parts is A103921.
The corresponding ordering of compositions is A124734.
Showing partitions as Heinz numbers gives A185974.
The version for non-reversed partitions is A334301.
Lexicographically ordered reversed partitions are A026791.
Sorting reversed partitions by Heinz number gives A112798.
The version for revlex instead of lex is A334302.
The version for revlex instead of colex is A334442.
Cf. A000041, A211992, A228531, A296774, A334433, A334435, A334436, A334437, A334438, A334439, A334440, A334441.

Join@@Table[Sort[Reverse/@IntegerPartitions[n]],{n,0,8}] (* Gus Wiseman, May 07 2020 *)
- or -
colen[f_,c_]:=OrderedQ[{Reverse[f],Reverse[c]}];
Reverse/@Join@@Table[Sort[IntegerPartitions[n],colen],{n,0,8}] (* Gus Wiseman, May 07 2020 *)

T036036(n,k)=k&&return(T036036(n)[k]);concat(partitions(n))
\\ If 2nd arg "k" is not given, return the n-th row as a vector. Assumes PARI version >= 2.7.1. See A193073 for "hand made" code.
concat(vector(8,n,T036036(n))) \\ to get the "flattened" sequence
\\ M. F. Hasler, Jul 12 2015

Edited by Daniel Forgues, Jan 21 2011
Edited by M. F. Hasler, Jul 12 2015
Name corrected by Gus Wiseman, May 12 2020

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

Original entry on oeis.org

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

cp's OEIS Frontend

A344616 Alternating sum of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A335448 Numbers whose prime indices are inseparable.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A335433 Numbers whose multiset of prime indices is separable.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A181821 a(n) = smallest integer with factorization as Product p(i)^e(i) such that Product p(e(i)) = n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A106529 Numbers having k prime factors (counted with multiplicity), the largest of which is the k-th prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A182850 a(n) = number of iterations that n requires to reach a fixed point under the x -> A181819(x) map.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Scheme

Formula

A326568 Denominator of the average of the multiset of prime indices of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs