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 26 results. Next

A006472 a(n) = n!*(n-1)!/2^(n-1).

Original entry on oeis.org

1, 1, 3, 18, 180, 2700, 56700, 1587600, 57153600, 2571912000, 141455160000, 9336040560000, 728211163680000, 66267215894880000, 6958057668962400000, 834966920275488000000, 113555501157466368000000, 17373991677092354304000000, 2970952576782792585984000000
Offset: 1

Views

Author

Keywords

Comments

Product of first (n-1) positive triangular numbers. - Amarnath Murthy, May 19 2002, corrected by Alex Ratushnyak, Dec 03 2013
Number of ways of transforming n distinguishable objects into n singletons via a sequence of n-1 refinements. Example: a(3)=3 because we have XYZ->X|YZ->X|Y|Z, XYZ->Y|XZ->X|Y|Z and XYZ->Z|XY->X|Y|Z. - Emeric Deutsch, Jan 23 2005
In other words, a(n) is the number of maximal chains in the lattice of set partitions of {1, ..., n} ordered by refinement. - Gus Wiseman, Jul 22 2018
From David Callan, Aug 27 2009: (Start)
With offset 0, a(n) = number of unordered increasing full binary trees of 2n edges with leaf set {n,n+1,...,2n}, where full binary means each nonleaf vertex has two children, increasing means the vertices are labeled 0,1,2,...,2n and each child is greater than its parent, unordered might as well mean ordered and each pair of sibling vertices is increasing left to right. For example, a(2)=3 counts the trees with edge lists {01,02,13,14}, {01,03,12,14}, {01,04,12,13}.
PROOF. Given such a tree of size n, to produce a tree of size n+1, two new leaves must be added to the leaf n. Choose any two of the leaf set {n+1,...,2n,2n+1,2n+2} for the new leaves and use the rest to replace the old leaves n+1,...,2n, maintaining relative order. Thus each tree of size n yields (n+2)-choose-2 trees of the next size up. Since the ratio a(n+1)/a(n)=(n+2)-choose-2, the result follows by induction.
Without the condition on the leaves, these trees are counted by the reduced tangent numbers A002105. (End)
a(n) = Sum(M(t)N(t)), where summation is over all rooted trees t with n vertices, M(t) is the number of ways to take apart t by sequentially removing terminal edges (see A206494) and N(t) is the number of ways to build up t from the one-vertex tree by adding successively edges to the existing vertices (the Connes-Moscovici weight; see A206496). See Remark on p. 3801 of the Hoffman reference. Example: a(3) = 3; indeed, there are two rooted trees with 3 vertices: t' = the path r-a-b and t" = V; we have M(t')=N(t')=1 and M(t") =1, N(t")=2, leading to M(t')N(t') + M(t")N(t")=3. - Emeric Deutsch, Jul 20 2012
Number of coalescence sequences or labeled histories for n lineages: the number of sequences by which n distinguishable leaves can coalesce to a single sequence. The coalescence process merges pairs of lineages into new lineages, labeling each newly formed lineage l by a subset of the n initial lineages corresponding to the union of all initial lineages that feed into lineage l. - Noah A Rosenberg, Jan 28 2019
Conjecture: For n > 1, n divides 2*a(n-1) + 4 if and only if n is prime. - Werner Schulte, Oct 04 2020
For a proof of the above conjecture see Himane. The list of primes p such that p^2 divides (2*a(p-1) + 4) (analog of A007540 - Wilson primes) begins [239, 24049, ...]. - Peter Bala, Nov 06 2024
a(n) is the number of maximal chains in the poset of set of permutations of {1, ..., n} ordered by containment. - Rajesh Kumar Mohapatra, Sep 03 2025

Examples

			From _Gus Wiseman_, Jul 22 2018: (Start)
The a(3) = 3 maximal chains in the lattice of set partitions of {1,2,3}:
  {{1},{2},{3}} < {{1},{2,3}} < {{1,2,3}}
  {{1},{2},{3}} < {{2},{1,3}} < {{1,2,3}}
  {{1},{2},{3}} < {{3},{1,2}} < {{1,2,3}} (End)
From _Rajesh Kumar Mohapatra_, Sep 03 2025: (Start)
The a(3) = 3 maximal chains in the poset of the set of permutations of {1,2,3}:
  {(1)(2)(3)} < (12)(3) < (123)}
  {(1)(2)(3)} < (1)(23) < (123)}
  {(1)(2)(3)} < (13)(2) < (132)} (End)
		

References

  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 148.
  • László Lovász, Combinatorial Problems and Exercises, North-Holland, 1979, p. 165.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Mike Steel, Phylogeny: Discrete and Random Processes in Evolution, SIAM, 2016, p. 47.

Crossrefs

Cf. A000110, A000258, A002846, A005121, A213427, A317145, A363679 (sum of reciprocals).
For the type B and D analogs, see A001044 and A123385.

Programs

  • Magma
    [Factorial(n)*Factorial(n-1)/2^(n-1): n in [1..20]]; // Vincenzo Librandi, Aug 23 2018
    
  • Maple
    A006472 := n -> n!*(n-1)!/2^(n-1):
  • Mathematica
    FoldList[Times,1,Accumulate[Range[20]]] (* Harvey P. Dale, Jan 10 2013 *)
  • PARI
    a(n) = n*(n-1)!^2/2^(n-1) \\ Charles R Greathouse IV, May 18 2015
    
  • Python
    from math import factorial
    def A006472(n): return n*factorial(n-1)**2 >> n-1 # Chai Wah Wu, Jun 22 2022

Formula

a(n) = a(n-1)*A000217(n-1).
a(n) = A010790(n-1)/2^(n-1).
a(n) = polygorial(n, 3) = (A000142(n)/A000079(n))*A000142(n+1) = (n!/2^n)*Product_{i=0..n-1} (i+2) = (n!/2^n)*Pochhammer(2, n) = (n!^2/2^n)*(n+1) = polygorial(n, 4)/2^n*(n+1). - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
a(n-1) = (-1)^(n+1)/(n^2*det(M_n)) where M_n is the matrix M_(i, j) = abs(1/i - 1/j). - Benoit Cloitre, Aug 21 2003
From Ilya Gutkovskiy, Dec 15 2016: (Start)
a(n) ~ 4*Pi*n^(2*n)/(2^n*exp(2*n)).
Sum_{n>=1} 1/a(n) = BesselI(1,2*sqrt(2))/sqrt(2) = 2.3948330992734... (End)
D-finite with recurrence 2*a(n) -n*(n-1)*a(n-1)=0. - R. J. Mathar, May 02 2022
Sum_{n>=1} (-1)^(n+1)/a(n) = BesselJ(1,2*sqrt(2))/sqrt(2). - Amiram Eldar, Jun 25 2022
From Rajesh Kumar Mohapatra, Sep 03 2025: (Start)
a(n) = A331955(n,n)
a(n) = A331956(n,n)
a(n) = A375835(n,n)
a(n) = A375837(n,n) (End)

A005121 Number of ultradissimilarity relations on an n-set.

Original entry on oeis.org

1, 1, 4, 32, 436, 9012, 262760, 10270696, 518277560, 32795928016, 2542945605432, 237106822506952, 26173354092593696, 3375693096567983232, 502995942483693043200, 85750135569136650473360, 16583651916595710735271248, 3611157196483089769387182064, 879518067472225603327860638128
Offset: 1

Views

Author

Keywords

Comments

First column in A154960. - Mats Granvik, Jan 18 2009
Number of chains from minimum to maximum in the lattice of set partitions of {1, ..., n} ordered by refinement. - Gus Wiseman, Jul 22 2018

Examples

			From _Gus Wiseman_, Jul 22 2018: (Start)
The (3) = 4 chains from minimum to maximum in the lattice of set partitions of {1,2,3}:
  {{1},{2},{3}} < {{1,2,3}}
  {{1},{2},{3}} < {{1},{2,3}} < {{1,2,3}}
  {{1},{2},{3}} < {{2},{1,3}} < {{1,2,3}}
  {{1},{2},{3}} < {{3},{1,2}} < {{1,2,3}}
(End)
		

References

  • L. Babai and T. Lengyel, A convergence criterion for recurrent sequences with application to the partition lattice, Analysis 12 (1992), 109-119.
  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 316-321.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    a[1] = 1; a[n_] := a[n] = Sum[StirlingS2[n, k]*a[k], {k, 1, n-1}]; Array[a, 19]
    (* Jean-François Alcover, Jun 24 2011, after Vladeta Jovovic *)
  • PARI
    {a(n) = local(A); if( n<1, 0, for(k=1, n, A = truncate(A) + x*O(x^k); A = x - A + subst(A, x, exp(x + x*O(x^k)) - 1)); n! * polcoeff(A, n))} /* Michael Somos, Sep 22 2007 */

Formula

a(n) = Sum_{i=1..n-1} N_i(n), where N_k(m) = Sum_{j=k..m-1} Stirling2(m, j)*N_{k-1}(j), m=3..n, k=2..m-1; N_1(2)=N_1(3)=...=N_1(n)=1.
a(n) = Sum_{k=1..n-1} Stirling2(n, k)*a(k) [Lengyel]. - Vladeta Jovovic, Apr 16 2003
E.g.f. satisfies Z(z) = 1/2 * (Z(exp(z)-1) - z). [Lengyel]
Asymptotic growth: a(n) ~ C_L*(n!)^2*(2log(2))^(-n)*n^(-1-1/3*log(2)) (Babai and Lengyel), with C_L = 1.0986858055... = A086053 [Flajolet and Salvy].
Sum_{k>=1} a(k-1)/fallfac(n,k) = -1/n^2 + 2*Sum_{k>=1} a(k-1)/n^k, with the falling factorials fallfac(n,k) = Product_{j=0..k-1}(n-j). - Vaclav Kotesovec, Aug 04 2015

Extensions

More terms from Vladeta Jovovic, Apr 16 2003

A318812 Number of total multiset partitions of the multiset of prime indices of n. Number of total factorizations of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 6, 1, 3, 1, 3, 1, 1, 1, 11, 1, 1, 2, 3, 1, 4, 1, 20, 1, 1, 1, 15, 1, 1, 1, 11, 1, 4, 1, 3, 3, 1, 1, 51, 1, 3, 1, 3, 1, 11, 1, 11, 1, 1, 1, 21, 1, 1, 3, 90, 1, 4, 1, 3, 1, 4, 1, 80, 1, 1, 3, 3, 1, 4, 1, 51, 6, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Sep 04 2018

Keywords

Comments

A total multiset partition of m is either m itself or a total multiset partition of a multiset partition of m that is neither minimal nor maximal.
a(n) depends only on the prime signature of n. - Andrew Howroyd, Dec 30 2019

Examples

			The a(24) = 11 total multiset partitions:
  {1,1,1,2}
  {{1},{1,1,2}}
  {{2},{1,1,1}}
  {{1,1},{1,2}}
  {{1},{1},{1,2}}
  {{1},{2},{1,1}}
  {{{1}},{{1},{1,2}}}
  {{{1}},{{2},{1,1}}}
  {{{2}},{{1},{1,1}}}
  {{{1,2}},{{1},{1}}}
  {{{1,1}},{{1},{2}}}
The a(24) = 11 total factorizations:
  24,
  (2*12), (3*8), (4*6),
  (2*2*6), (2*3*4),
  ((2)*(2*6)), ((6)*(2*2)), ((2)*(3*4)), ((3)*(2*4)), ((4)*(2*3)).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    totfac[n_]:=1+Sum[totfac[Times@@Prime/@f],{f,Select[facs[n],1
    				
  • PARI
    MultEulerT(u)={my(v=vector(#u)); v[1]=1; for(k=2, #u, forstep(j=#v\k*k, k, -k, my(i=j, e=0); while(i%k==0, i/=k; e++; v[j]+=binomial(e+u[k]-1, e)*v[i]))); v}
    seq(n)={my(v=vector(n, i, isprime(i)), u=vector(n), m=logint(n,2)+1); for(r=1, m, u += v*sum(j=r, m, (-1)^(j-r)*binomial(j-1, r-1)); v=MultEulerT(v)); u[1]=1; u} \\ Andrew Howroyd, Dec 30 2019

Formula

a(product of n distinct primes) = A005121(n).
a(prime^n) = A318813(n).

A317144 Number of refinement-ordered pairs of factorizations of n into factors > 1.

Original entry on oeis.org

1, 1, 1, 3, 1, 3, 1, 6, 3, 3, 1, 9, 1, 3, 3, 14, 1, 9, 1, 9, 3, 3, 1, 23, 3, 3, 6, 9, 1, 12, 1, 26, 3, 3, 3, 31, 1, 3, 3, 23, 1, 12, 1, 9, 9, 3, 1, 56, 3, 9, 3, 9, 1, 23, 3, 23, 3, 3, 1, 41, 1, 3, 9, 55, 3, 12, 1, 9, 3, 12, 1, 82, 1, 3, 9, 9, 3, 12, 1, 56, 14
Offset: 1

Views

Author

Gus Wiseman, Jul 22 2018

Keywords

Comments

If x and y are factorizations of the same integer and it is possible to produce x by further factoring the factors of y, flattening, and sorting, then x <= y.
As this is a sequence computed from exponents in factorization of n, distinct values of a(n) in this sequence can be found by computing a(A025487(k)) for k >= 0. - David A. Corneth, Jul 30 2018

Examples

			The a(12) = 9 ordered pairs:
  (2*2*3) <= (12)
  (2*2*3) <= (2*6)
  (2*2*3) <= (3*4)
  (2*2*3) <= (2*2*3)
    (2*6) <= (12)
    (2*6) <= (2*6)
    (3*4) <= (12)
    (3*4) <= (3*4)
     (12) <= (12)
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    faccaps[fac_]:=Union[Sort/@Apply[Times,mps[fac],{2}]];
    Table[Sum[Length[faccaps[fac]],{fac,facs[n]}],{n,100}]

Formula

a(n) >= A001055(n) + floor(A000005(n) / 2) - 1. - David A. Corneth, Jul 30 2018

A330976 Numbers that are not the number of factorizations into factors > 1 of any positive integer.

Original entry on oeis.org

6, 8, 10, 13, 14, 17, 18, 20, 23, 24, 25, 27, 28, 32, 33, 34, 35, 37, 39, 40, 41, 43, 44, 46, 48, 49, 50, 51, 53, 54, 55, 58, 59, 60, 61, 62, 63, 65, 68, 69, 70, 71, 72, 73, 75, 76, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 93, 94, 95, 96, 99
Offset: 1

Views

Author

Gus Wiseman, Jan 07 2020

Keywords

Comments

Warning: I have only confirmed the first eight terms. The rest are derived from A045782. - Gus Wiseman, Jan 07 2020

Crossrefs

Complement of A045782.
The strict version is A330975.
Factorizations are A001055, with image A045782.
Strict factorizations are A045778, with image A045779.
The least number with n factorizations is A330973(n).

Programs

  • Mathematica
    nn=15;
    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[fam[#]&,2^nn];
    Complement[Range[nn],nds]

A318813 Number of balanced reduced multisystems with n atoms all equal to 1.

Original entry on oeis.org

1, 1, 2, 6, 20, 90, 468, 2910, 20644, 165874, 1484344, 14653890, 158136988, 1852077284, 23394406084, 317018563806, 4587391330992, 70598570456104, 1151382852200680, 19835976878704628, 359963038816096924, 6863033015330999110, 137156667020252478684, 2867083618970831936826
Offset: 1

Views

Author

Gus Wiseman, Sep 04 2018

Keywords

Comments

For n > 1, also the number of balanced reduced multisystems whose atoms are an integer partition of n with at least one part > 1. A balanced reduced multisystem is either a finite multiset, or a multiset partition with at least two parts, not all of which are singletons, of a balanced reduced multisystem. - Gus Wiseman, Dec 31 2019

Examples

			The a(5) = 20 balanced reduced multisystems (with n written in place of 1^n):
  5  (14)  (23)  (113)      (122)      (1112)
                 ((1)(13))  ((1)(22))  ((1)(112))
                 ((3)(11))  ((2)(12))  ((2)(111))
                                       ((11)(12))
                                       ((1)(1)(12))
                                       ((1)(2)(11))
                                       (((1))((1)(12)))
                                       (((1))((2)(11)))
                                       (((2))((1)(11)))
                                       (((12))((1)(1)))
                                       (((11))((1)(2)))
		

Crossrefs

Programs

  • Mathematica
    normize[m_]:=m/.Rule@@@Table[{Union[m][[i]],i},{i,Length[Union[m]]}];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    totfact[n_]:=totfact[n]=1+Sum[totfact[Times@@Prime/@normize[f]],{f,Select[facs[n],1
    				
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    seq(n)={my(v=vector(n, i, i==1), u=vector(n)); for(r=1, #v, u += v*sum(j=r, #v, (-1)^(j-r)*binomial(j-1, r-1)); v=EulerT(v)); u} \\ Andrew Howroyd, Dec 30 2019

Formula

a(n > 1) = A330679(n)/2. - Gus Wiseman, Dec 31 2019

Extensions

Terms a(14) and beyond from Andrew Howroyd, Dec 30 2019
Terminology corrected by Gus Wiseman, Dec 31 2019

A330977 Numbers whose number of factorizations into factors > 1 (A001055) is a power of 2.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 25, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 65, 67, 68, 69, 71, 72, 73, 74, 75, 76, 77, 79, 82, 83, 85, 86, 87
Offset: 1

Views

Author

Gus Wiseman, Jan 07 2020

Keywords

Comments

The complement starts: 8, 16, 24, 27, 30, 32, 36, 40.

Examples

			Factorizations of n = 1, 4, 12, 72:
  ()  (4)    (12)     (72)
      (2*2)  (2*6)    (8*9)
             (3*4)    (2*36)
             (2*2*3)  (3*24)
                      (4*18)
                      (6*12)
                      (2*4*9)
                      (2*6*6)
                      (3*3*8)
                      (3*4*6)
                      (2*2*18)
                      (2*3*12)
                      (2*2*2*9)
                      (2*2*3*6)
                      (2*3*3*4)
                      (2*2*2*3*3)
		

Crossrefs

The same for strict integer partitions is A331022.
Factorizations are A001055, with image A045782.
The least number with exactly n factorizations is A330973(n).
The least number with exactly 2^n factorizations is A330989(n).
Numbers whose inverse prime shadow belongs to this sequence are A330990.
Numbers with a prime number of factorizations are A330991.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Select[Range[100],IntegerQ[Log[2,Length[facs[#]]]]&]

A080688 Resort the index of A064553 using A080444 and maintaining ascending order within each grouping: seen as a triangle read by rows, the n-th row contains the A001055(n) numbers m with A064553(m)=n.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 6, 11, 13, 8, 10, 17, 9, 19, 14, 23, 29, 12, 15, 22, 31, 37, 26, 41, 21, 43, 16, 20, 25, 34, 47, 53, 18, 33, 38, 59, 61, 28, 35, 46, 67, 39, 71, 58, 73, 79, 24, 30, 44, 51, 55, 62, 83, 49, 89, 74, 97, 27, 57, 101, 52, 65, 82
Offset: 1

Views

Author

Alford Arnold, Mar 23 2003

Keywords

Comments

The number 12 can be written as 3*2*2, 4*3, 6*2 and 12 corresponding to each of the four values (12,15,22,31) in the example. Note that A001055(12) = 4. Since A001055(n) depends only on the least prime signature, the values 1,2,4,6,8,12,16,24,30,32,36,... A025487 are of special interest when counting multisets. (see for example, A035310 and A035341).
A064553(T(n,k)) = A080444(n,k) = n for k=1..A001055(n); T(n,1) = A064554(n); T(n,A001055(n)) = A064554(n). - Reinhard Zumkeller, Oct 01 2012
Row n is the sorted list of shifted Heinz numbers of factorizations of n into factors > 1, where the shifted Heinz number of a factorization (y_1, ..., y_k) is prime(y_1 - 1) * ... * prime(y_k - 1). - Gus Wiseman, Sep 05 2018

Examples

			a(18),a(19),a(20) and a(21) are 12,15,22 and 31 because A064553(12,15,22,31) = (12,12,12,12) similarly, A064553(36,45,66,76,93,95,118,121,149) = (36,36,36,36,36,36,36,36,36)
From _Gus Wiseman_, Sep 05 2018: (Start)
Triangle begins:
   1
   2
   3
   4  5
   7
   6 11
  13
   8 10 17
   9 19
  14 23
  29
  12 15 22 31
  37
  26 41
  21 43
  16 20 25 34 47
Corresponding triangle of factorizations begins:
  (),
  (2),
  (3),
  (2*2), (4),
  (5),
  (2*3), (6),
  (7),
  (2*2*2), (2*4), (8),
  (3*3), (9),
  (2*5), (10),
  (11),
  (2*2*3), (3*4), (2*6), (12).
(End)
		

Crossrefs

Programs

  • Haskell
    a080688 n k = a080688_row n !! (k-1)
    a080688_row n = map (+ 1) $ take (a001055 n) $
                    elemIndices n $ map fromInteger a064553_list
    a080688_tabl = map a080688_row [1..]
    a080688_list = concat a080688_tabl
    -- Reinhard Zumkeller, Oct 01 2012
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[facs[n/d],Min@@#1>=d&],{d,Rest[Divisors[n]]}]];
    Table[Sort[Table[Times@@Prime/@(f-1),{f,facs[n]}]],{n,20}] (* Gus Wiseman, Sep 05 2018 *)

Extensions

More terms from Sean A. Irvine, Oct 05 2011
Keyword tabf added and definition complemented accordingly by Reinhard Zumkeller, Oct 01 2012

A330935 Irregular triangle read by rows where T(n,k) is the number of length-k chains from minimum to maximum in the poset of factorizations of n into factors > 1, ordered by refinement.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 04 2020

Keywords

Comments

This poset is equivalent to the poset of multiset partitions of the prime indices of n, ordered by refinement.

Examples

			Triangle begins:
   1:          16: 0 1 3 2    31: 1            46: 0 1
   2: 1        17: 1          32: 0 1 5 8 4    47: 1
   3: 1        18: 0 1 2      33: 0 1          48: 0 1 10 23 15
   4: 0 1      19: 1          34: 0 1          49: 0 1
   5: 1        20: 0 1 2      35: 0 1          50: 0 1 2
   6: 0 1      21: 0 1        36: 0 1 7 7      51: 0 1
   7: 1        22: 0 1        37: 1            52: 0 1 2
   8: 0 1 1    23: 1          38: 0 1          53: 1
   9: 0 1      24: 0 1 5 5    39: 0 1          54: 0 1 5 5
  10: 0 1      25: 0 1        40: 0 1 5 5      55: 0 1
  11: 1        26: 0 1        41: 1            56: 0 1 5 5
  12: 0 1 2    27: 0 1 1      42: 0 1 3        57: 0 1
  13: 1        28: 0 1 2      43: 1            58: 0 1
  14: 0 1      29: 1          44: 0 1 2        59: 1
  15: 0 1      30: 0 1 3      45: 0 1 2        60: 0 1 9 11
Row n = 48 counts the following chains (minimum and maximum not shown):
  ()  (6*8)      (2*3*8)->(6*8)       (2*2*2*6)->(2*4*6)->(6*8)
      (2*24)     (2*4*6)->(6*8)       (2*2*3*4)->(2*3*8)->(6*8)
      (3*16)     (2*3*8)->(2*24)      (2*2*3*4)->(2*4*6)->(6*8)
      (4*12)     (2*3*8)->(3*16)      (2*2*2*6)->(2*4*6)->(2*24)
      (2*3*8)    (2*4*6)->(2*24)      (2*2*2*6)->(2*4*6)->(4*12)
      (2*4*6)    (2*4*6)->(4*12)      (2*2*3*4)->(2*3*8)->(2*24)
      (3*4*4)    (3*4*4)->(3*16)      (2*2*3*4)->(2*3*8)->(3*16)
      (2*2*12)   (3*4*4)->(4*12)      (2*2*3*4)->(2*4*6)->(2*24)
      (2*2*2*6)  (2*2*12)->(2*24)     (2*2*3*4)->(2*4*6)->(4*12)
      (2*2*3*4)  (2*2*12)->(4*12)     (2*2*3*4)->(3*4*4)->(3*16)
                 (2*2*2*6)->(6*8)     (2*2*3*4)->(3*4*4)->(4*12)
                 (2*2*3*4)->(6*8)     (2*2*2*6)->(2*2*12)->(2*24)
                 (2*2*2*6)->(2*24)    (2*2*2*6)->(2*2*12)->(4*12)
                 (2*2*2*6)->(4*12)    (2*2*3*4)->(2*2*12)->(2*24)
                 (2*2*3*4)->(2*24)    (2*2*3*4)->(2*2*12)->(4*12)
                 (2*2*3*4)->(3*16)
                 (2*2*3*4)->(4*12)
                 (2*2*2*6)->(2*4*6)
                 (2*2*3*4)->(2*3*8)
                 (2*2*3*4)->(2*4*6)
                 (2*2*3*4)->(3*4*4)
                 (2*2*2*6)->(2*2*12)
                 (2*2*3*4)->(2*2*12)
		

Crossrefs

Row lengths are A001222.
Row sums are A317176.
Column k = 1 is A010051.
Column k = 2 is A066247.
Column k = 3 is A330936.
Final terms of each row are A317145.
The version for set partitions is A008826, with row sums A005121.
The version for integer partitions is A330785, with row sums A213427.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    upfacs[q_]:=Union[Sort/@Join@@@Tuples[facs/@q]];
    paths[eds_,start_,end_]:=If[start==end,Prepend[#,{}],#]&[Join@@Table[Prepend[#,e]&/@paths[eds,Last[e],end],{e,Select[eds,First[#]==start&]}]];
    Table[Length[Select[paths[Join@@Table[{y,#}&/@DeleteCases[upfacs[y],y],{y,facs[n]}],{n},First[facs[n]]],Length[#]==k-1&]],{n,100},{k,PrimeOmega[n]}]

Formula

T(2^n,k) = A330785(n,k).
T(n,1) + T(n,2) = 1.

A330665 Number of balanced reduced multisystems of maximal depth whose atoms are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 3, 1, 5, 1, 1, 1, 7, 1, 1, 1, 5, 1, 3, 1, 2, 2, 1, 1, 16, 1, 2, 1, 2, 1, 5, 1, 5, 1, 1, 1, 11, 1, 1, 2, 16, 1, 3, 1, 2, 1, 3, 1, 27, 1, 1, 2, 2, 1, 3, 1, 16, 2, 1, 1, 11, 1
Offset: 1

Views

Author

Gus Wiseman, Dec 27 2019

Keywords

Comments

First differs from A317145 at a(32) = 5, A317145(32) = 4.
A balanced reduced multisystem is either a finite multiset, or a multiset partition with at least two parts, not all of which are singletons, of a balanced reduced multisystem.
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 series/singleton-reduced factorizations of n with Omega(n) levels of parentheses. See A001055, A050336, A050338, A050340, etc.

Examples

			The a(n) multisystems for n = 2, 6, 12, 24, 48:
  {1}  {1,2}  {{1},{1,2}}  {{{1}},{{1},{1,2}}}  {{{{1}}},{{{1}},{{1},{1,2}}}}
              {{2},{1,1}}  {{{1,1}},{{1},{2}}}  {{{{1}}},{{{1,1}},{{1},{2}}}}
                           {{{1}},{{2},{1,1}}}  {{{{1},{1}}},{{{1}},{{1,2}}}}
                           {{{1,2}},{{1},{1}}}  {{{{1},{1,1}}},{{{1}},{{2}}}}
                           {{{2}},{{1},{1,1}}}  {{{{1,1}}},{{{1}},{{1},{2}}}}
                                                {{{{1}}},{{{1}},{{2},{1,1}}}}
                                                {{{{1}}},{{{1,2}},{{1},{1}}}}
                                                {{{{1},{1}}},{{{2}},{{1,1}}}}
                                                {{{{1},{1,2}}},{{{1}},{{1}}}}
                                                {{{{1,1}}},{{{2}},{{1},{1}}}}
                                                {{{{1}}},{{{2}},{{1},{1,1}}}}
                                                {{{{1},{2}}},{{{1}},{{1,1}}}}
                                                {{{{1,2}}},{{{1}},{{1},{1}}}}
                                                {{{{2}}},{{{1}},{{1},{1,1}}}}
                                                {{{{2}}},{{{1,1}},{{1},{1}}}}
                                                {{{{2},{1,1}}},{{{1}},{{1}}}}
		

Crossrefs

The last nonzero term in row n of A330667 is a(n).
The chain version is A317145.
The non-maximal version is A318812.
Unlabeled versions are A330664 and A330663.
Other labeled versions are A330675 (strongly normal) and A330676 (normal).

Programs

  • Mathematica
    primeMS[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[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    totm[m_]:=Prepend[Join@@Table[totm[p],{p,Select[mps[m],1
    				

Formula

a(2^n) = A000111(n - 1).
a(product of n distinct primes) = A006472(n).
Showing 1-10 of 26 results. Next