cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 22 results. Next

A318360 Number of set multipartitions (multisets of sets) of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 5, 3, 2, 1, 6, 1, 2, 3, 15, 1, 9, 1, 6, 3, 2, 1, 21, 4, 2, 16, 6, 1, 10, 1, 52, 3, 2, 4, 35, 1, 2, 3, 22, 1, 10, 1, 6, 19, 2, 1, 83, 5, 13, 3, 6, 1, 66, 4, 22, 3, 2, 1, 41, 1, 2, 20, 203, 4, 10, 1, 6, 3, 14, 1, 153, 1, 2, 26, 6, 5, 10, 1
Offset: 1

Views

Author

Gus Wiseman, Aug 24 2018

Keywords

Examples

			The a(12) = 6 set multipartitions of {1,1,2,3}:
  {{1},{1,2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{1},{1},{2},{3}}
		

Crossrefs

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
    Table[Length[sqfacs[Times@@Prime/@nrmptn[n]]],{n,80}]
  • PARI
    permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
    sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i,2], j, primepi(f[i,1]))))}
    count(sig)={my(n=vecsum(sig), s=0); forpart(p=n, my(q=prod(i=1, #p, 1 + x^p[i] + O(x*x^n))); s+=prod(i=1, #sig, polcoef(q,sig[i]))*permcount(p)); s/n!}
    a(n)={if(n==1, 1, my(s=sig(n)); if(#s<=2, if(#s==1, 1, min(s[1],s[2])+1), count(sig(n))))} \\ Andrew Howroyd, Dec 10 2018

Formula

a(n) = A050320(A181821(n)).
From Andrew Howroyd, Dec 10 2018:(Start)
a(p) = 1 for prime(p).
a(prime(i)*prime(j)) = min(i,j) + 1.
a(prime(n)^k) = A188392(n,k). (End)

A318284 Number of multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 5, 9, 7, 7, 11, 11, 12, 16, 15, 15, 26, 22, 21, 29, 19, 30, 36, 31, 30, 66, 38, 42, 52, 56, 52, 47, 45, 57, 92, 77, 67, 77, 74, 101, 98, 135, 64, 137, 97, 176, 135, 109, 109, 118, 105, 231, 249, 97, 141, 181, 139, 297, 198, 385, 195, 269
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			The a(12) = 11 multiset partitions of {1,1,2,3}:
  {{1,1,2,3}}
  {{1},{1,2,3}}
  {{2},{1,1,3}}
  {{3},{1,1,2}}
  {{1,1},{2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{2},{3},{1,1}}
  {{1},{1},{2},{3}}
		

Crossrefs

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[facs[Times@@Prime/@nrmptn[n]]],{n,60}]
  • PARI
    permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
    sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i,2], j, primepi(f[i,1]))))}
    count(sig)={my(n=vecsum(sig), A=O(x*x^vecmax(sig)), s=0); forpart(p=n, my(q=1/prod(i=1, #p, 1 - x^p[i] + A)); s+=prod(i=1, #sig, polcoef(q,sig[i]))*permcount(p)); s/n!}
    a(n)={if(n==1, 1, my(s=sig(n)); if(#s==1, numbpart(s[1]), count(sig(n))))} \\ Andrew Howroyd, Dec 10 2018

Formula

a(n) = A001055(A181821(n)).
a(prime(n)^k) = A219727(n,k). - Andrew Howroyd, Dec 10 2018

A305936 Irregular triangle whose n-th row is the multiset spanning an initial interval of positive integers with multiplicities equal to the n-th row of A296150 (the prime indices of n in weakly decreasing order).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 3, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 2, 3, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 2, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			Row 90 is {1,1,1,2,2,3,3,4} because 90 = prime(3)*prime(2)*prime(2)*prime(1).
Triangle begins:
   1:
   2:  1
   3:  1  1
   4:  1  2
   5:  1  1  1
   6:  1  1  2
   7:  1  1  1  1
   8:  1  2  3
   9:  1  1  2  2
  10:  1  1  1  2
  11:  1  1  1  1  1
  12:  1  1  2  3
  13:  1  1  1  1  1  1
		

Crossrefs

Row lengths are A056239. Number of distinct elements in row n is A001222(n). Number of distinct multiplicities in row n is A001221(n).

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Array[nrmptn,30]

A318361 Number of strict set multipartitions (sets of sets) of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 0, 2, 0, 1, 0, 5, 1, 0, 0, 4, 0, 0, 0, 15, 0, 5, 0, 1, 0, 0, 0, 16, 0, 0, 8, 0, 0, 2, 0, 52, 0, 0, 0, 23, 0, 0, 0, 7, 0, 0, 0, 0, 5, 0, 0, 68, 0, 1, 0, 0, 0, 40, 0, 1, 0, 0, 0, 14, 0, 0, 1, 203, 0, 0, 0, 0, 0, 0, 0, 111, 0, 0, 4, 0, 0, 0, 0, 41, 80, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Aug 24 2018

Keywords

Examples

			The a(24) = 16 sets of sets with multiset union {1,1,2,3,4}:
  {{1},{1,2,3,4}}
  {{1,2},{1,3,4}}
  {{1,3},{1,2,4}}
  {{1,4},{1,2,3}}
  {{1},{2},{1,3,4}}
  {{1},{3},{1,2,4}}
  {{1},{4},{1,2,3}}
  {{1},{1,2},{3,4}}
  {{1},{1,3},{2,4}}
  {{1},{1,4},{2,3}}
  {{2},{1,3},{1,4}}
  {{3},{1,2},{1,4}}
  {{4},{1,2},{1,3}}
  {{1},{2},{3},{1,4}}
  {{1},{2},{4},{1,3}}
  {{1},{3},{4},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
    Table[Length[sqfacs[Times@@Prime/@nrmptn[n]]],{n,90}]
  • PARI
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i, 2], j, primepi(f[i, 1]))))}
    count(sig)={my(r=0, A=O(x*x^vecmax(sig))); for(n=1, vecsum(sig)+1, my(s=0); forpart(p=n, my(q=prod(i=1, #p, 1 + x^p[i] + A)); s+=prod(i=1, #sig, polcoef(q, sig[i]))*(-1)^#p*permcount(p)); r+=(-1)^n*s/n!); r/2}
    a(n)={if(n==1, 1, my(s=sig(n)); if(#s==1, s[1]==1, count(sig(n))))} \\ Andrew Howroyd, Dec 18 2018

Formula

a(n) = A050326(A181821(n)).
a(prime(n)^k) = A188445(n, k). - Andrew Howroyd, Dec 17 2018

A045779 Number of factorizations of n into distinct factors for some n (image of A045778).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 17, 18, 19, 21, 22, 25, 27, 31, 32, 33, 34, 38, 40, 42, 43, 44, 46, 52, 54, 55, 56, 57, 59, 61, 64, 67, 70, 74, 76, 80, 83, 88, 89, 91, 93, 100, 104, 110, 111, 112, 116, 117, 120, 122, 123, 132, 137, 140, 141, 142, 143, 148
Offset: 1

Views

Author

Keywords

Comments

We may use A045778(k*m) >= A045778(k) for any k, m >= 1 to disprove presence of some positive integer in this sequence. - David A. Corneth, Oct 24 2024

Examples

			From _David A. Corneth_, Oct 24 2024: (Start)
5 is a term as 24 has five factorizations into distinct divisors of 24 namely 24 = 2 * 12 = 3 * 8 = 4 * 6 = 2 * 3 * 4 which is five such factorizations.
11 is not a term. From terms in A025487 only the numbers 2, 4, 6, 8, 12, 16, 24, 30, 32, 36, 48, 60, 64, 72, 96, 128, 256, 512, 1024 have no more than 11 such factorizations. Any multiple of these numbers in A025487 that is not already listed has more than 11 such factorizations which proves 11 is not in this sequence. (End)
		

Crossrefs

Factorizations are A001055, with image A045782, with complement A330976.
Strict factorizations are A045778 with image A045779 and complement A330975.
The least number with A045779(n) strict factorizations is A045780(n).
The least number with n strict factorizations is A330974(n).

Extensions

Name edited by Gus Wiseman, Jan 11 2020

A045780 Least value with A045779(n) factorizations into distinct factors.

Original entry on oeis.org

1, 6, 12, 64, 24, 256, 48, 512, 60, 96, 2048, 144, 210, 120, 216, 180, 384, 288, 16384, 240, 432, 420, 65536, 1536, 360, 480, 900, 864, 3072, 1152, 1296, 2310, 524288, 6144, 960, 720, 840, 2304, 1728, 1080, 1260, 2592, 2097152, 1800, 4608, 24576, 4194304, 1440, 3456
Offset: 1

Views

Author

Keywords

Examples

			From _Gus Wiseman_, Jan 11 2020: (Start)
The strict factorizations of a(n) for n = 1..9:
  ()  (6)    (12)   (64)     (24)     (256)     (48)     (512)     (60)
      (2*3)  (2*6)  (2*32)   (3*8)    (4*64)    (6*8)    (8*64)    (2*30)
             (3*4)  (4*16)   (4*6)    (8*32)    (2*24)   (16*32)   (3*20)
                    (2*4*8)  (2*12)   (2*128)   (3*16)   (2*256)   (4*15)
                             (2*3*4)  (2*4*32)  (4*12)   (4*128)   (5*12)
                                      (2*8*16)  (2*3*8)  (2*4*64)  (6*10)
                                                (2*4*6)  (2*8*32)  (2*5*6)
                                                         (4*8*16)  (3*4*5)
                                                                   (2*3*10)
(End)
30 is not in the sequence even though A045779(30) = 5. As 24 is the smallest k such that A045779(k) = 5 we have a(m) = 24 where m is such that A045779(m) = 5 which turns out to be m = 5 (not every positive integer is in A045779). So a(5) = 24. - _David A. Corneth_, Oct 24 2024
		

Crossrefs

All terms belong to A025487.
The non-strict version is A045783.
The sorted version is A330997.
Factorizations are A001055 with image A045782 and complement A330976.
Strict factorizations are A045778 with image A045779 and complement A330975.
The least number with exactly n strict factorizations is A330974(n).

Extensions

More terms from David A. Corneth, Oct 24 2024

A219585 Number A(n,k) of k-partite partitions of {n}^k into distinct k-tuples; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 5, 5, 2, 1, 1, 15, 40, 17, 2, 1, 1, 52, 457, 364, 46, 3, 1, 1, 203, 6995, 14595, 2897, 123, 4, 1, 1, 877, 136771, 937776, 407287, 21369, 323, 5, 1, 1, 4140, 3299218, 88507276, 107652681, 10200931, 148257, 809, 6, 1
Offset: 0

Views

Author

Alois P. Heinz, Nov 23 2012

Keywords

Comments

A(n,k) is the number of factorizations of m^n into distinct factors where m is a product of k distinct primes. A(2,2) = 5: (2*3)^2 = 36 has 5 factorizations into distinct factors: 36, 3*12, 4*9, 2*18, 2*3*6.

Examples

			A(1,3) = 5: [(1,1,1)], [(1,1,0),(0,0,1)], [(1,0,1),(0,1,0)], [(1,0,0),(0,1,0),(0,0,1)], [(0,1,1),(1,0,0)].
A(3,2) = 17: [(3,3)], [(3,0),(0,3)], [(3,2),(0,1)], [(2,3),(1,0)], [(3,1),(0,2)], [(2,2),(1,1)], [(1,3),(2,0)], [(2,1),(1,2)], [(2,1),(1,1),(0,1)], [(3,0),(0,2),(0,1)], [(2,2),(1,0),(0,1)], [(2,1),(0,2),(1,0)], [(1,2),(2,0),(0,1)], [(1,2),(1,1),(1,0)], [(0,3),(2,0),(1,0)], [(2,0),(1,1),(0,2)], [(2,0),(0,2),(1,0),(0,1)].
Square array A(n,k) begins:
  1,  1,   1,      1,          1,            1,         1, ...
  1,  1,   2,      5,         15,           52,       203, ...
  1,  1,   5,     40,        457,         6995,    136771, ...
  1,  2,  17,    364,      14595,       937776,  88507276, ...
  1,  2,  46,   2897,     407287,    107652681,  ...
  1,  3, 123,  21369,   10200931,  10781201973,  ...
  1,  4, 323, 148257,  233051939,  ...
  1,  5, 809, 970246, 4909342744,  ...
		

Crossrefs

Columns k=0..5 give: A000012, A000009, A219554, A219560, A219561, A219565.
Rows n=0..3 give: A000012, A000110, A094574, A319591.
Cf. A188445, A219727 (partitions of {n}^k into k-tuples), A318286.

Programs

  • Mathematica
    f[n_, k_] := f[n, k] = 1/2 Product[Sum[O[x[j]]^(n+1), {j, 1, k}]+1+ Product[x[j]^i[j], {j, 1, k}], Evaluate[Sequence @@ Table[{i[j], 0, n}, {j, 1, k}]]];
    a[0, ] = a[, 0] = 1; a[n_, k_] := SeriesCoefficient[f[n, k], Sequence @@ Table[{x[j], 0, n}, {j, 1, k}]];
    Table[Print[a[n-k, k]]; a[n-k, k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Dec 11 2013, updated Sep 16 2019 *)
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    D(p, n, k)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); EulerT(v)[n]^k/prod(i=1, #v, i^v[i]*v[i]!)}
    T(n, k)={my(m=n*k+1, q=Vec(exp(intformal(O(x^m) - x^n/(1-x)))/(1+x))); if(n==0, 1, (-1)^m*sum(j=0, m, my(s=0); forpart(p=j, s+=(-1)^#p*D(p, n, k), [1, n]); s*q[#q-j])/2)} \\ Andrew Howroyd, Dec 16 2018

Formula

A(n,k) = [(Product_{j=1..k} x_j)^n] 1/2 * Product_{i_1,...,i_k>=0} (1+Product_{j=1..k} x_j^i_j).

A330974 Least positive integer with n factorizations into distinct factors > 1, and 0 if no such number exists.

Original entry on oeis.org

1, 6, 12, 64, 24, 256, 48, 512, 60, 96, 0, 2048, 0, 144, 210, 120, 216, 180, 384, 0, 288, 16384, 0, 0, 240, 0, 432, 0, 0, 0, 420, 65536, 1536, 360, 0, 0, 0, 480, 0, 900, 0, 864, 3072, 1152, 0, 1296, 0, 0, 0, 0, 0, 2310, 0, 524288, 6144, 960, 720, 0, 840, 0, 2304
Offset: 1

Views

Author

Gus Wiseman, Jan 06 2020

Keywords

Crossrefs

All nonzero terms belong to A025487.
Strict factorizations are A045778, with image A045779.
The version with zeros removed is A045780.
The non-strict version is A330973.
Positions of zeros are A330975.
The sorted version is A330997.

Programs

  • Mathematica
    nn=10;
    fam[n_]:=fam[n]=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[fam[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    nds=Length/@Array[Select[fam[#],UnsameQ@@#&]&,2^nn];
    Table[If[#=={},0,#[[1,1]]]&[Position[nds,i]],{i,nn}]

Extensions

More terms from Jinyuan Wang, Jul 07 2021

A381452 Number of multisets that can be obtained by partitioning the prime indices of n into a set of multisets and taking their sums.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 3, 1, 2, 2, 2, 1, 3, 1, 3, 2, 2, 1, 4, 1, 2, 2, 3, 1, 5, 1, 3, 2, 2, 2, 4, 1, 2, 2, 5, 1, 5, 1, 3, 3, 2, 1, 5, 1, 3, 2, 3, 1, 5, 2, 5, 2, 2, 1, 7, 1, 2, 3, 4, 2, 5, 1, 3, 2, 5, 1, 6, 1, 2, 3, 3, 2, 5, 1, 6, 2, 2, 1, 8, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2025

Keywords

Comments

First differs from A045778 at a(24) = 4, A045778(24) = 5.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a factorization of n into distinct factors > 1.
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.
A multiset partition can be regarded as an arrow in the poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Sets of multisets are generally not transitive. For example, we have arrows: {{1},{2},{1,2}}: {1,1,2,2} -> {1,2,3} and {{1,2},{3}}: {1,2,3} -> {3,3}, but there is no set of multisets {1,1,2,2} -> {3,3}.

Examples

			The prime indices of 24 are {1,1,1,2}, with 5 partitions into a set of multisets:
  {{1,1,1,2}}
  {{1},{1,1,2}}
  {{2},{1,1,1}}
  {{1,1},{1,2}}
  {{1},{2},{1,1}}
with block-sums: {5}, {1,4}, {2,3}, {2,3}, {1,2,2}, of which 4 are distinct, so a(24) = 4.
		

Crossrefs

Before taking sums we had A045778.
If each block is a set we have A381441, before sums A050326.
For distinct block-sums instead of blocks we have A381637, before sums A321469.
Other multiset partitions of prime indices:
- For multisets of constant multisets (A000688) see A381455 (upper), A381453 (lower).
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For set multipartitions (A050320) see A381078 (upper), A381454 (lower).
- For sets of constant multisets (A050361) see A381715.
- For set systems with distinct sums (A381633) see A381634, zeros A293243.
- For sets of constant multisets with distinct sums (A381635) see A381716, A381636.
More on sets of multisets: A261049, A317776, A317775, A296118, A318286.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Table[Length[Union[Sort[Total/@#]&/@Select[mps[prix[n]],UnsameQ@@#&]]],{n,100}]

Formula

a(A002110(n)) = A066723(n).

A330997 Sorted list containing the least number with each possible nonzero number of factorizations into distinct factors > 1.

Original entry on oeis.org

1, 6, 12, 24, 48, 60, 64, 96, 120, 144, 180, 210, 216, 240, 256, 288, 360, 384, 420, 432, 480, 512, 720, 840, 864, 900, 960, 1080, 1152, 1260, 1296, 1440, 1536, 1680, 1728, 1800, 2048, 2160, 2304, 2310, 2520, 2592, 2880, 3072, 3360, 3456, 3600, 3840, 4320
Offset: 1

Views

Author

Gus Wiseman, Jan 06 2020

Keywords

Examples

			The strict factorizations of a(n) for n = 1..9.
  {}  6    12   24     48     60      64     96      120
      2*3  2*6  3*8    6*8    2*30    2*32   2*48    2*60
           3*4  4*6    2*24   3*20    4*16   3*32    3*40
                2*12   3*16   4*15    2*4*8  4*24    4*30
                2*3*4  4*12   5*12           6*16    5*24
                       2*3*8  6*10           8*12    6*20
                       2*4*6  2*5*6          2*6*8   8*15
                              3*4*5          3*4*8   10*12
                              2*3*10         2*3*16  3*5*8
                                             2*4*12  4*5*6
                                                     2*3*20
                                                     2*4*15
                                                     2*5*12
                                                     2*6*10
                                                     3*4*10
                                                     2*3*4*5
		

Crossrefs

All terms belong to A025487.
Strict factorizations are A045778, with image A045779.
The unsorted version is A045780.
The non-strict version is A330972.
The least number with n strict factorizations is A330974.

Programs

  • Mathematica
    nn=1000;
    strfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strfacs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
    nds=Length/@Array[strfacs,nn];
    Table[Position[nds,i][[1,1]],{i,First/@Gather[nds]}]
Showing 1-10 of 22 results. Next