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

A326211 Number of unsortable normal multiset partitions of weight n.

Original entry on oeis.org

0, 0, 0, 1, 17, 170, 1455, 11678, 92871, 752473
Offset: 0

Views

Author

Gus Wiseman, Jun 19 2019

Keywords

Comments

A multiset partition is normal if it covers an initial interval of positive integers. It is unsortable if no permutation has an ordered concatenation, or equivalently if the concatenation of its lexicographically-ordered parts is not weakly increasing. For example, the multiset partition {{1,2},{1,1,1},{2,2,2}} is sortable because the permutation ((1,1,1),(1,2),(2,2,2)) has concatenation (1,1,1,1,2,2,2,2), which is weakly increasing.

Examples

			The a(3) = 1 and a(4) = 17 multiset partitions:
  {{1,3},{2}}  {{1,1,3},{2}}
               {{1,2},{1,2}}
               {{1,2},{1,3}}
               {{1,2,3},{2}}
               {{1,2,4},{3}}
               {{1,3},{2,2}}
               {{1,3},{2,3}}
               {{1,3},{2,4}}
               {{1,3,3},{2}}
               {{1,3,4},{2}}
               {{1,4},{2,3}}
               {{1},{1,3},{2}}
               {{1},{2,4},{3}}
               {{1,3},{2},{2}}
               {{1,3},{2},{3}}
               {{1,3},{2},{4}}
               {{1,4},{2},{3}}
		

Crossrefs

Unsortable set partitions are A058681.
Sortable normal multiset partitions are A326212.
Non-crossing normal multiset partitions are A324171.
MM-numbers of unsortable multiset partitions are A326258.

Programs

  • Mathematica
    lexsort[f_,c_]:=OrderedQ[PadRight[{f,c}]];
    allnorm[n_]:=If[n<=0,{{}},Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1]];
    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]]]];
    Table[Length[Select[Sort[#,lexsort]&/@Join@@mps/@allnorm[n],!OrderedQ[Join@@#]&]],{n,0,5}]

Formula

A255906(n) = a(n) + A326212(n).

A054391 Number of permutations with certain forbidden subsequences.

Original entry on oeis.org

1, 1, 2, 5, 14, 41, 123, 374, 1147, 3538, 10958, 34042, 105997, 330632, 1032781, 3229714, 10109310, 31667245, 99260192, 311294876, 976709394, 3065676758, 9625674442, 30231524869, 94972205349, 298419158008, 937861780439, 2947969125284, 9267666915326
Offset: 0

Views

Author

N. J. A. Sloane, Elisa Pergola (elisa(AT)dsi.unifi.it), May 21 2000

Keywords

Comments

Hankel transform is [1,1,1,...] = A000012. - Paul Barry, Jan 19 2009
The inverse Motzkin transform apparently yields 1 followed by A000930, which implies a generating function g(x)=1+z/(1-z-z^3) where z=x*A001006(x). - R. J. Mathar, Jul 07 2009
It appears that the infinite set of interpolated sequences between the Motzkin and the Catalan can be generated with a succession of INVERT transforms, given each sequence has two leading 1's. Also, the N-th sequence in the set starting with (N=1, A001006) can be generated from a production matrix of the form "M" in the formula section, such that the main diagonal is (N leading 1's, 0, 0, 0, ...). M with a diagonal of (1, 0, 0, 0, ...) generates A001006, while M with a main diagonal of all 1's is the production matrix for A000108. - Gary W. Adamson, Jul 29 2011
From Gus Wiseman, Jun 22 2019: (Start)
Conjecture: Also the number of non-capturing set partitions of {1..n}. A set partition is capturing if it has two blocks of the form {...x...y...} and {...z...t...} where x < z and y > t or x > z and y < t. This is a weaker condition than nesting, so for example {{1,3,5},{2,4}} is capturing but not nesting. The a(0) = 1 through a(4) = 14 non-capturing set partitions are:
{} {{1}} {{1,2}} {{1,2,3}} {{1,2,3,4}}
{{1},{2}} {{1},{2,3}} {{1},{2,3,4}}
{{1,2},{3}} {{1,2},{3,4}}
{{1,3},{2}} {{1,2,3},{4}}
{{1},{2},{3}} {{1,2,4},{3}}
{{1,3},{2,4}}
{{1,3,4},{2}}
{{1},{2},{3,4}}
{{1},{2,3},{4}}
{{1,2},{3},{4}}
{{1},{2,4},{3}}
{{1,3},{2},{4}}
{{1,4},{2},{3}}
{{1},{2},{3},{4}}
(End)
The above conjecture is true: A partition is non-capturing iff its representation in canonical sequential form avoids the patterns 1221 and 2112. In the context of these partition representations, the pattern 2112 is equivalent to the pattern 12112. Partitions whose canonical sequence form avoid 1221 and 12112 are one of the classes that are handled in the Mansour/Shattuck "Pattern Avoiding Partitions,..." paper. It shows that they are counted by this sequence. - Christian Sievers, Oct 29 2024

Examples

			a(4) = 14, a(5) = 41 since the top row of M^4 = (14, 14, 9, 3, 1), with 41 = (14 + 14 + 9 + 3 + 1).
		

Crossrefs

Interpolates between Motzkin numbers (A001006) and Catalan numbers (A000108). Cf. A005773, A054392, ...
Binomial transform of A224747.

Programs

  • Maple
    c := x->(1-sqrt(1-4*x))/(2*x); a := (x,j)->(x)/((1-4*x)*(c(x))^2*(1-c(x))^(j))*(-x^2*(c(x))^2*(1-c(x))*(x^2*(c(x))^4)^(j)-(1-3*x-2*x^2)*(c(x))^2*(x*(c(x))^2)^(j)+x);
    b := (x,j)->1+(1)/((1-4*x)*c(x)*(1-c(x))^(j))*(-2*x^3*(c(x))^2*(x^2*(c(x))^4)^(j)+(1-3*x-2*x^2)*c(x)*(x*(c(x))^2)^(j)-2*x^2);
    co := (x,j)->(1)/((1-4*x)*(1-c(x))^(j))*(x^2*(x^2*(c(x))^4)^(j)-(1-3*x-2*x^2)*(x*(c(x))^2)^(j)+x^2);
    s := (x,j)->(1-b(x,j)+(-1)^j*sqrt((1-b(x,j))^2-4*a(x,j)*co(x,j)))/(2*a(x,j)); j := 3; series(s(x,j),x=0..60); od; # j=1,2,3,... inf gives A001006, A005773, A054391, A054392, ..., A000108
  • Mathematica
    CoefficientList[Series[1 - 2*x^2/(2*x^2 - 3*x + 1 - Sqrt[1 - 2*x - 3*x^2]), {x, 0, 50}], x] (* G. C. Greubel, Apr 27 2017 *)
  • Maxima
    a(n):=sum((sum((binomial(k,l)*l*sum(binomial(j,1-n-2*l+k+2*j)*binomial(n-1+l-k,j),j,0,n+l-k-1))/(n+l-k-1),l,1,k)),k,1,n-1)+1; /* Vladimir Kruchinin, Oct 31 2011 */
    
  • PARI
    x='x+O('x^66); gf=1-2*x^2/(2*x^2-3*x+1-sqrt(1-2*x-3*x^2)); Vec(gf) \\ Joerg Arndt, Jun 29 2013

Formula

G.f.: 1 - 2*x^2 / (2*x^2 - 3*x + 1 - sqrt(1 - 2*x - 3*x^2)). - Mansour and Shattuck
G.f.: 1/(1-x-x^2/(1-2x-x^2/(1-x-x^2/(1-x-x^2/(1-x-x^2/(1-x-x^2/(1-... (continued fraction) (conjecture). - Paul Barry, Jan 19 2009
From Gary W. Adamson, Jul 29 2011: (Start)
a(n) = upper left term of M^n, a(n+1) = sum of top row terms of M^n; M = an infinite square production matrix as follows with a main diagonal of (1, 1, 1, 0, 0, 0, ...):
1, 1, 0, 0, 0, 0, ...
1, 1, 1, 0, 0, 0, ...
1, 1, 1, 1, 0, 0, ...
1, 1, 1, 0, 1, 0, ...
1, 1, 1, 1, 0, 1, ...
1, 1, 1, 1, 1, 0, ...
... (End)
a(n) = Sum_{k=1..n-1} (sum(l=1..k, (binomial(k,l)*l*sum(j=0..n+l-k-1, binomial(j,1-n-2*l+k+2*j)*binomial(n-1+l-k,j)))/(n+l-k-1))) + 1. - Vladimir Kruchinin, Oct 31 2011
D-finite with recurrence (-n+1)*a(n) + 3*(2*n-3)*a(n-1) + (-8*n+11)*a(n-2) + (-5*n+32)*a(n-3) + (7*n-31)*a(n-4) + 3*(-n+4)*a(n-5)= 0. - R. J. Mathar, Nov 26 2012
G.f.: 1 - x*(2*x^2-3*x+1 + 1/G(0))/(2*(x^3-3*x^2+4*x-1)), where G(k)= 1 + x*(2+3*x)*(4*k+1)/( 4*k+2 - x*(2+3*x)*(4*k+2)*(4*k+3)/(x*(2+3*x)*(4*k+3) + 4*(k+1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jun 29 2013

A326258 MM-numbers of unsortable multiset partitions (with empty parts allowed).

Original entry on oeis.org

145, 169, 215, 290, 338, 355, 377, 395, 430, 435, 473, 481, 505, 507, 535, 559, 565, 580, 645, 667, 676, 695, 710, 725, 754, 790, 793, 803, 815, 841, 845, 860, 865, 869, 870, 905, 923, 946, 962, 965, 989, 995, 1010, 1014, 1015, 1027, 1065, 1070, 1073, 1075
Offset: 1

Views

Author

Gus Wiseman, Jun 22 2019

Keywords

Comments

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 multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n.
A multiset partition is unsortable if no permutation has an ordered concatenation. For example, the multiset partition ((1,2),(1,1,1),(2,2,2)) is sortable because the permutation ((1,1,1),(1,2),(2,2,2)) has concatenation (1,1,1,1,2,2,2,2), which is weakly increasing.

Examples

			The sequence of terms together with their multiset multisystems begins:
  145: {{2},{1,3}}
  169: {{1,2},{1,2}}
  215: {{2},{1,4}}
  290: {{},{2},{1,3}}
  338: {{},{1,2},{1,2}}
  355: {{2},{1,1,3}}
  377: {{1,2},{1,3}}
  395: {{2},{1,5}}
  430: {{},{2},{1,4}}
  435: {{1},{2},{1,3}}
  473: {{3},{1,4}}
  481: {{1,2},{1,1,2}}
  505: {{2},{1,6}}
  507: {{1},{1,2},{1,2}}
  535: {{2},{1,1,4}}
  559: {{1,2},{1,4}}
  565: {{2},{1,2,3}}
  580: {{},{},{2},{1,3}}
  645: {{1},{2},{1,4}}
  667: {{2,2},{1,3}}
		

Crossrefs

Unsortable set partitions are A058681.
Normal unsortable multiset partitions are A326211.
Unsortable digraphs are A326209.
MM-numbers of crossing multiset partitions are A324170.
MM-numbers of nesting multiset partitions are A326256.
MM-numbers of capturing multiset partitions are A326255.

Programs

  • Mathematica
    lexsort[f_,c_]:=OrderedQ[PadRight[{f,c}]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1000],!OrderedQ[Join@@Sort[primeMS/@primeMS[#],lexsort]]&]

A326249 Number of capturing set partitions of {1..n} that are not nesting.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 9, 55, 283, 1324, 5838, 24744
Offset: 0

Views

Author

Gus Wiseman, Jun 20 2019

Keywords

Comments

Capturing is a weaker condition than nesting. A set partition is capturing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < t < y or z < x < y < t, and nesting if it has two blocks of the form {...x,y...}, {...z,t...} where x < z < t < y or z < x < y < t. For example, {{1,3,5},{2,4}} is capturing but not nesting, so is counted under a(5).

Examples

			The a(6) = 9 set partitions:
  {{1},{2,4,6},{3,5}}
  {{1,3,5},{2,4},{6}}
  {{1,3,6},{2,4},{5}}
  {{1,3,6},{2,5},{4}}
  {{1,4,6},{2},{3,5}}
  {{1,4,6},{2,5},{3}}
  {{1,3,5},{2,4,6}}
  {{1,2,4,6},{3,5}}
  {{1,3,5,6},{2,4}}
		

Crossrefs

MM-numbers of capturing, non-nesting multiset partitions are A326260.
Nesting set partitions are A016098.
Capturing set partitions are A326243.
Non-crossing, nesting set partitions are A122880 (conjectured).

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    capXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,y_,_},_,{_,z_,t_,_},_}/;x
    				

A326245 Number of crossing, non-capturing set partitions of {1..n}.

Original entry on oeis.org

0, 0, 0, 0, 1, 7, 34, 141, 537, 1941, 6777, 23096, 77340
Offset: 0

Views

Author

Gus Wiseman, Jun 20 2019

Keywords

Comments

A set partition is crossing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < y < t or z < x < t < y, and capturing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < t < y or z < x < y < t. Capturing is a weaker condition than nesting, so for example {{1,3,5},{2,4}} is capturing but not nesting.

Examples

			The a(4) = 1 and a(5) = 7 set partitions:
  {{1,3},{2,4}}  {{1,2,4},{3,5}}
                 {{1,3},{2,4,5}}
                 {{1},{2,4},{3,5}}
                 {{1,3},{2,4},{5}}
                 {{1,3},{2,5},{4}}
                 {{1,4},{2},{3,5}}
                 {{1,4},{2,5},{3}}
		

Crossrefs

Crossing set partitions are A016098.
Non-capturing set partitions are A054391.
Crossing, capturing set partitions are A326246.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				

A326260 MM-numbers of capturing, non-nesting multiset partitions (with empty parts allowed).

Original entry on oeis.org

2599, 4163, 5198, 6463, 6893, 7291, 7797, 8326, 8507, 9131, 9959, 10396, 10649, 11041, 11639, 12489, 12811, 12926, 12995, 13786, 14237, 14582, 14899, 15157, 15594, 16123, 16403, 16652, 17014, 17063, 17089, 17141, 18101, 18193, 18262, 18643, 18659, 19337, 19389
Offset: 1

Views

Author

Gus Wiseman, Jun 22 2019

Keywords

Comments

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 multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n.
A set partition is capturing if it has two blocks of the form {...x...y...} and {...z...t...} where x < z and y > t or x > z and y < t. It is nesting if it has two blocks of the form {...x,y...} and {...z,t...} where x < z and y > t or x > z and y < t. Capturing is a weaker condition than nesting, so for example {{1,3,5},{2,4}} is capturing but not nesting.

Examples

			The sequence of terms together with their multiset multisystems begins:
   2599: {{2,2},{1,2,3}}
   4163: {{2,2},{1,2,4}}
   5198: {{},{2,2},{1,2,3}}
   6463: {{2,2},{1,1,2,3}}
   6893: {{1,2,2},{1,2,3}}
   7291: {{2,2},{1,2,5}}
   7797: {{1},{2,2},{1,2,3}}
   8326: {{},{2,2},{1,2,4}}
   8507: {{2,3},{1,2,4}}
   9131: {{2,2},{1,2,6}}
   9959: {{2,2},{1,1,2,4}}
  10396: {{},{},{2,2},{1,2,3}}
  10649: {{2,2},{1,2,2,3}}
  11041: {{1,2,2},{1,2,4}}
  11639: {{2,2,2},{1,2,3}}
  12489: {{1},{2,2},{1,2,4}}
  12811: {{2,2},{1,2,7}}
  12926: {{},{2,2},{1,1,2,3}}
  12995: {{2},{2,2},{1,2,3}}
  13786: {{},{1,2,2},{1,2,3}}
		

Crossrefs

Non-nesting set partitions are A000108.
Capturing set partitions are A326243.
Capturing, non-nesting set partitions are A326249.
MM-numbers of nesting multiset partitions are A326256.
MM-numbers of capturing multiset partitions are A326255.

Programs

  • Mathematica
    capXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,y_,_},_,{_,z_,t_,_},_}/;xTable[PrimePi[p],{k}]]]];
    Select[Range[10000],!nesXQ[primeMS/@primeMS[#]]&&capXQ[primeMS/@primeMS[#]]&]

A326332 Number of integer partitions of n with unsortable prime factors.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 2, 3, 5, 9, 14, 22, 33, 50, 71, 100, 140, 196, 265, 360, 480, 641, 842, 1104, 1432, 1855, 2378, 3040, 3858, 4888, 6146, 7708, 9616, 11969, 14818, 18305, 22511, 27629, 33773, 41191, 50069, 60744, 73453, 88645, 106681
Offset: 0

Views

Author

Gus Wiseman, Jun 27 2019

Keywords

Comments

An integer partition has unsortable prime factors if there is no permutation (c_1,...,c_k) of the parts such that the maximum prime factor of c_i is at most the minimum prime factor of c_{i+1}. For example, the partition (27,8,6) is sortable because the permutation (8,6,27) satisfies the condition.

Examples

			The a(12) = 1 through a(17) = 14 partitions:
  (6,6)  (10,3)   (6,6,2)    (6,6,3)      (10,6)         (14,3)
         (6,6,1)  (10,3,1)   (10,3,2)     (6,6,4)        (6,6,5)
                  (6,6,1,1)  (6,6,2,1)    (10,3,3)       (10,4,3)
                             (10,3,1,1)   (6,6,2,2)      (10,6,1)
                             (6,6,1,1,1)  (6,6,3,1)      (6,6,3,2)
                                          (10,3,2,1)     (6,6,4,1)
                                          (6,6,2,1,1)    (10,3,2,2)
                                          (10,3,1,1,1)   (10,3,3,1)
                                          (6,6,1,1,1,1)  (6,6,2,2,1)
                                                         (6,6,3,1,1)
                                                         (10,3,2,1,1)
                                                         (6,6,2,1,1,1)
                                                         (10,3,1,1,1,1)
                                                         (6,6,1,1,1,1,1)
		

Crossrefs

Sortable integer partitions are A326333.
Unsortable set partitions are A058681.
Unsortable normal multiset partitions are A326211.
MM-numbers of unsortable multiset partitions are A326258.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!OrderedQ[Join@@Sort[First/@FactorInteger[#]&/@#,OrderedQ[PadRight[{#1,#2}]]&]]&]],{n,0,20}]

Formula

A000041(n) = a(n) + A326333(n).

A326291 Number of unsortable factorizations of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 3, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Jun 24 2019

Keywords

Comments

A factorization into factors > 1 is unsortable if there is no permutation (c_1,...,c_k) of the factors such that the maximum prime factor of c_i is at most the minimum prime factor of c_{i+1}. For example, the factorization (6*8*27) is sortable because the permutation (8,6,27) satisfies the condition.

Examples

			The a(180) = 10 unsortable factorizations:
  (2*3*3*10)  (5*6*6)   (3*60)
              (2*3*30)  (6*30)
              (2*9*10)  (9*20)
              (3*3*20)  (10*18)
              (3*6*10)
Missing from this list are:
  (2*2*3*3*5)  (2*2*5*9)   (4*5*9)   (2*90)   (180)
               (2*3*5*6)   (2*2*45)  (4*45)
               (3*3*4*5)   (2*5*18)  (5*36)
               (2*2*3*15)  (2*6*15)  (12*15)
                           (3*4*15)
                           (3*5*12)
		

Crossrefs

Unsortable set partitions are A058681.
Unsortable normal multiset partitions are A326211.
MM-numbers of unsortable multiset partitions are A326258.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    lexsort[f_,c_]:=OrderedQ[PadRight[{f,c}]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[facs[n],!OrderedQ[Join@@Sort[primeMS/@#,lexsort]]&]],{n,100}]

A326333 Number of integer partitions of n with sortable prime factors.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 11, 15, 22, 30, 42, 56, 76, 99, 132, 171, 222, 283, 363, 457, 577, 721, 902, 1115, 1379, 1693, 2076, 2530, 3077, 3723, 4500, 5410, 6494, 7765, 9270, 11025, 13089, 15491, 18307, 21569, 25369, 29765, 34869, 40750, 47546, 55361, 64367, 74685, 86529
Offset: 0

Views

Author

Gus Wiseman, Jun 27 2019

Keywords

Comments

An integer partition has sortable prime factors if there is a permutation (c_1,...,c_k) of the parts such that the maximum prime factor of c_i is at most the minimum prime factor of c_{i+1}. For example, the partition (27,8,6) is sortable because the permutation (8,6,27) satisfies the condition.

Crossrefs

Unsortable integer partitions are A326332.
Sortable normal multiset partitions are A326212.
Sortable factorizations are A326334.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],OrderedQ[Join@@Sort[First/@FactorInteger[#]&/@#,OrderedQ[PadRight[{#1,#2}]]&]]&]],{n,0,20}]

Formula

A000041(n) = a(n) + A326332(n).

A326334 Number of sortable factorizations of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 5, 1, 4, 1, 4, 2, 2, 1, 7, 2, 2, 3, 4, 1, 4, 1, 7, 2, 2, 2, 8, 1, 2, 2, 7, 1, 4, 1, 4, 4, 2, 1, 12, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 8, 1, 2, 4, 11, 2, 4, 1, 4, 2, 4, 1, 14, 1, 2, 4, 4, 2, 4, 1, 12, 5, 2, 1, 8, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Jun 27 2019

Keywords

Comments

A factorization into factors > 1 is sortable if there is a permutation (c_1,...,c_k) of the factors such that the maximum prime factor (in the standard factorization of an integer into prime numbers) of c_i is at most the minimum prime factor of c_{i+1}. For example, the factorization (6*8*27) is sortable because the permutation (8,6,27) satisfies the condition.

Examples

			The a(180) = 16 sortable factorizations:
  (2*2*3*3*5)  (2*2*5*9)   (4*5*9)   (2*90)   (180)
               (2*3*5*6)   (2*2*45)  (4*45)
               (3*3*4*5)   (2*5*18)  (5*36)
               (2*2*3*15)  (2*6*15)  (12*15)
                           (3*4*15)
                           (3*5*12)
Missing from this list are the following unsortable factorizations:
  (2*3*3*10)  (5*6*6)   (3*60)
              (2*3*30)  (6*30)
              (2*9*10)  (9*20)
              (3*3*20)  (10*18)
              (3*6*10)
		

Crossrefs

Factorizations are A001055.
Unsortable factorizations are A326291.
Sortable integer partitions are A326333.

Programs

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

Formula

A001055(n) = a(n) + A326291(n).
Showing 1-10 of 11 results. Next