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

A050342 Expansion of Product_{m>=1} (1+x^m)^A000009(m).

Original entry on oeis.org

1, 1, 1, 3, 4, 7, 12, 19, 30, 49, 77, 119, 186, 286, 438, 670, 1014, 1528, 2300, 3437, 5119, 7603, 11241, 16564, 24343, 35650, 52058, 75820, 110115, 159510, 230522, 332324, 477994, 686044, 982519, 1404243, 2003063, 2851720, 4052429, 5748440, 8140007, 11507125
Offset: 0

Views

Author

Christian G. Bower, Oct 15 1999

Keywords

Comments

Number of partitions of n into distinct parts with one level of parentheses. Each "part" in parentheses is distinct from all others at the same level. Thus (2+1)+(1) is allowed but (2)+(1+1) and (2+1+1) are not.

Examples

			4=(4)=(3)+(1)=(3+1)=(2+1)+(1).
From _Gus Wiseman_, Oct 11 2018: (Start)
a(n) is the number of set systems (sets of sets) whose multiset union is an integer partition of n. For example, the a(1) = 1 through a(6) = 12 set systems are:
  {{1}}  {{2}}  {{3}}      {{4}}        {{5}}        {{6}}
                {{1,2}}    {{1,3}}      {{1,4}}      {{1,5}}
                {{1},{2}}  {{1},{3}}    {{2,3}}      {{2,4}}
                           {{1},{1,2}}  {{1},{4}}    {{1,2,3}}
                                        {{2},{3}}    {{1},{5}}
                                        {{1},{1,3}}  {{2},{4}}
                                        {{2},{1,2}}  {{1},{1,4}}
                                                     {{1},{2,3}}
                                                     {{2},{1,3}}
                                                     {{3},{1,2}}
                                                     {{1},{2},{3}}
                                                     {{1},{2},{1,2}}
(End)
		

Crossrefs

Programs

  • Maple
    g:= proc(n, i) option remember; `if`(n=0, 1,
          `if`(i<1, 0, g(n, i-1)+`if`(i>n, 0, g(n-i, i-1))))
        end:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(binomial(g(i, i), j)*b(n-i*j, i-1), j=0..n/i)))
        end:
    a:= n-> b(n, n):
    seq(a(n), n=0..50);  # Alois P. Heinz, May 19 2013
  • Mathematica
    g[n_, i_] := g[n, i] = If[n==0, 1, If[i<1, 0, g[n, i-1] + If[i>n, 0, g[n-i, i-1]]]]; b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, Sum[Binomial[g[i, i], j]*b[n-i*j, i-1], {j, 0, n/i}]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Dec 19 2015, after Alois P. Heinz *)
    nn=10;Table[SeriesCoefficient[Product[(1+x^k)^PartitionsQ[k],{k,nn}],{x,0,n}],{n,0,nn}] (* Gus Wiseman, Oct 11 2018 *)

Formula

Weigh transform of A000009.

A279375 Number of set partitions of strict integer partitions of n that have distinct block-sums.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 12, 16, 24, 39, 49, 70, 94, 127, 202, 247, 340, 450, 606, 772, 1169, 1407, 1920, 2454, 3267, 4089, 5469, 7293, 9222, 11884, 15291, 19417, 24890, 31469, 39662, 52619, 64764, 82502, 103576, 131169, 162726, 206015, 254233, 318464, 406262, 499210, 620593, 773673, 957073, 1181593
Offset: 0

Views

Author

Gus Wiseman, Dec 11 2016

Keywords

Comments

Also twice partitioned numbers where all partitions are strict. Also triangles of weight n in the multisystem of strict partitions. Strict partitions are an example of a multisystem that is neither transitive nor partitive nor contractible but is decomposable; see link for details.

Examples

			The a(6)=9 set partitions of strict integer partitions of 6 are: ((6)), ((51)), ((5)(1)), ((42)), ((4)(2)), ((321)), ((32)(1)), ((31)(2)), ((3)(2)(1)). The set partition ((3)(21)) is not counted because its blocks do not have distinct sums.
		

Crossrefs

Programs

  • Mathematica
    nn=20;sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Total[Length[Select[sps[Sort[#]],UnsameQ@@Total/@#&]]&/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,nn}]

A259936 Number of ways to express the integer n as a product of its unitary divisors (A034444).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 2, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 5, 2, 2, 2, 2, 1, 5, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 5
Offset: 1

Views

Author

Geoffrey Critzer, Jul 09 2015

Keywords

Comments

Equivalently, a(n) is the number of ways to express the cyclic group Z_n as a direct sum of its Hall subgroups. A Hall subgroup of a finite group G is a subgroup whose order is coprime to its index.
a(n) is the number of ways to partition the set of distinct prime factors of n.
Also the number of singleton or pairwise coprime factorizations of n. - Gus Wiseman, Sep 24 2019

Examples

			a(60) = 5 because we have: 60 = 4*3*5 = 4*15 = 3*20 = 5*12.
For n = 36, its unitary divisors are 1, 4, 9, 36. From these we obtain 36 either as 1*36 or 4*9, thus a(36) = 2. - _Antti Karttunen_, Oct 21 2017
		

Crossrefs

Differs from A050320 for the first time at n=36.
Differs from A354870 for the first time at n=210, where a(210) = 15, while A354870(210) = 12.
Related classes of factorizations:
- No conditions: A001055
- Strict: A045778
- Constant: A089723
- Distinct multiplicities: A255231
- Singleton or coprime: A259936
- Relatively prime: A281116
- Aperiodic: A303386
- Stable (indivisible): A305149
- Connected: A305193
- Strict relatively prime: A318721
- Uniform: A319269
- Intersecting: A319786
- Constant or distinct factors coprime: A327399
- Constant or relatively prime: A327400
- Coprime: A327517
- Not relatively prime: A327658
- Distinct factors coprime: A327695

Programs

  • Maple
    map(combinat:-bell @ nops @ numtheory:-factorset, [$1..100]); # Robert Israel, Jul 09 2015
  • Mathematica
    Table[BellB[PrimeNu[n]], {n, 1, 75}]
    (* second program *)
    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],Length[#]==1||CoprimeQ@@#&]],{n,100}] (* Gus Wiseman, Sep 24 2019 *)
  • PARI
    a(n) = my(t=omega(n), x='x, m=contfracpnqn(matrix(2, t\2, y, z, if( y==1, -z*x^2, 1 - (z+1)*x)))); polcoeff(1/(1 - x + m[2, 1]/m[1, 1]) + O(x^(t+1)), t) \\ Charles R Greathouse IV, Jun 30 2017

Formula

a(n) = A000110(A001221(n)).
a(n > 1) = A327517(n) + 1. - Gus Wiseman, Sep 24 2019

Extensions

Incorrect comment removed by Antti Karttunen, Jun 11 2022

A294617 Number of ways to choose a set partition of a strict integer partition of n.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 10, 12, 17, 24, 44, 51, 76, 98, 138, 217, 272, 366, 493, 654, 848, 1284, 1560, 2115, 2718, 3610, 4550, 6024, 8230, 10296, 13354, 17144, 21926, 27903, 35556, 44644, 59959, 73456, 94109, 117735, 150078, 185800, 235719, 290818, 365334, 467923
Offset: 0

Views

Author

Gus Wiseman, Nov 05 2017

Keywords

Comments

From Gus Wiseman, Sep 17 2024: (Start)
Also the number of strict integer compositions of n whose leaders, obtained by splitting into maximal increasing subsequences and taking the first term of each, are decreasing. For example, the strict composition (3,6,2,1,4) has maximal increasing subsequences ((3,6),(2),(1,4)), with leaders (3,2,1), so is counted under a(16). The a(0) = 1 through a(7) = 12 compositions are:
() (1) (2) (3) (4) (5) (6) (7)
(1,2) (1,3) (1,4) (1,5) (1,6)
(2,1) (3,1) (2,3) (2,4) (2,5)
(3,2) (4,2) (3,4)
(4,1) (5,1) (4,3)
(1,2,3) (5,2)
(2,1,3) (6,1)
(2,3,1) (1,2,4)
(3,1,2) (2,1,4)
(3,2,1) (2,4,1)
(4,1,2)
(4,2,1)
(End)

Examples

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

Crossrefs

Row sums of A330460 and of A330759.
This is a strict case of A374689, weak version A189076.
A011782 counts compositions, strict A032020.
A238130, A238279, A333755 count compositions by number of runs.

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n>i*(i+1)/2, 0,
          `if`(n=0, combinat[bell](t), b(n, i-1, t)+
          `if`(i>n, 0, b(n-i, min(n-i, i-1), t+1))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..50);  # Alois P. Heinz, Nov 07 2017
  • Mathematica
    Table[Total[BellB[Length[#]]&/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,25}]
    (* Second program: *)
    b[n_, i_, t_] := b[n, i, t] = If[n > i (i + 1)/2, 0, If[n == 0, BellB[t], b[n, i - 1, t] + If[i > n, 0, b[n - i, Min[n - i, i - 1], t + 1]]]];
    a[n_] := b[n, n, 0];
    a /@ Range[0, 50] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)

Formula

A279375(n) <= a(n) <= A279790(n).
G.f.: Sum_{k>=0} Bell(k) * x^(k*(k + 1)/2) / Product_{j=1..k} (1 - x^j). - Ilya Gutkovskiy, Jan 28 2020

A007713 Number of 4-level rooted trees with n leaves.

Original entry on oeis.org

1, 1, 4, 10, 30, 75, 206, 518, 1344, 3357, 8429, 20759, 51044, 123973, 299848, 719197, 1716563, 4070800, 9607797, 22555988, 52718749, 122655485, 284207304, 655894527, 1508046031, 3454808143, 7887768997, 17949709753, 40719611684, 92096461012, 207697731344
Offset: 0

Views

Author

Keywords

Examples

			From _Gus Wiseman_, Oct 11 2018: (Start)
Also the number of multiset partitions of multiset partitions of integer partitions of n. For example, the a(1) = 1 through a(4) = 30 multiset partitions are:
  ((1))  ((2))       ((3))            ((4))
         ((11))      ((12))           ((13))
         ((1)(1))    ((111))          ((22))
         ((1))((1))  ((1)(2))         ((112))
                     ((1)(11))        ((1111))
                     ((1))((2))       ((1)(3))
                     ((1))((11))      ((2)(2))
                     ((1)(1)(1))      ((1)(12))
                     ((1))((1)(1))    ((2)(11))
                     ((1))((1))((1))  ((1)(111))
                                      ((11)(11))
                                      ((1))((3))
                                      ((2))((2))
                                      ((1))((12))
                                      ((1)(1)(2))
                                      ((2))((11))
                                      ((1))((111))
                                      ((1)(1)(11))
                                      ((11))((11))
                                      ((1))((1)(2))
                                      ((2))((1)(1))
                                      ((1))((1)(11))
                                      ((1)(1)(1)(1))
                                      ((11))((1)(1))
                                      ((1))((1))((2))
                                      ((1))((1))((11))
                                      ((1))((1)(1)(1))
                                      ((1)(1))((1)(1))
                                      ((1))((1))((1)(1))
                                      ((1))((1))((1))((1))
(End)
		

Crossrefs

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d,j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: b0:= etr(1): b1:= etr(b0): a:= etr(b1): seq(a(n), n=0..30); # Alois P. Heinz, Sep 08 2008
  • Mathematica
    i[ n_, m_ ] := 1 /; m==1 || n==0; i[ n_, m_ ] := (i[ n, m ]=1/n Sum[ i[ k, m ] Plus @@ ((# i[ #, m-1 ])& /@ Divisors[ n-k ]), {k, 0, n-1} ]) /; n>0 && m>1
    etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[ j]}]*b[n-j], {j, 1, n}]/n]; b]; b0 = etr[Function[1]]; b1 = etr[b0]; a = etr[b1]; Table[a[n], {n, 1, 30}] (* Jean-François Alcover, Mar 05 2015, after Alois P. Heinz *)

Formula

Euler transform applied thrice to all-1's sequence.

A256130 Number T(n,k) of partitions of n into parts of exactly k sorts which are introduced in ascending order; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 3, 4, 1, 0, 5, 12, 7, 1, 0, 7, 30, 33, 11, 1, 0, 11, 72, 130, 77, 16, 1, 0, 15, 160, 463, 438, 157, 22, 1, 0, 22, 351, 1557, 2216, 1223, 289, 29, 1, 0, 30, 743, 5031, 10422, 8331, 2957, 492, 37, 1, 0, 42, 1561, 15877, 46731, 52078, 26073, 6401, 788, 46, 1
Offset: 0

Views

Author

Alois P. Heinz, Mar 15 2015

Keywords

Comments

In general, column k>1 is asymptotic to c*k^n, where c = 1/(k!*Product_{n>=1} (1-1/k^n)) = 1/(k!*QPochhammer[1/k, 1/k]). - Vaclav Kotesovec, Jun 01 2015

Examples

			T(3,1) = 3: 1a1a1a, 2a1a, 3a.
T(3,2) = 4: 1a1a1b, 1a1b1a, 1a1b1b, 2a1b.
T(3,3) = 1: 1a1b1c.
Triangle T(n,k) begins:
  1;
  0,  1;
  0,  2,    1;
  0,  3,    4,     1;
  0,  5,   12,     7,     1;
  0,  7,   30,    33,    11,     1;
  0, 11,   72,   130,    77,    16,     1;
  0, 15,  160,   463,   438,   157,    22,    1;
  0, 22,  351,  1557,  2216,  1223,   289,   29,   1;
  0, 30,  743,  5031, 10422,  8331,  2957,  492,  37,  1;
  0, 42, 1561, 15877, 46731, 52078, 26073, 6401, 788, 46,  1;
  ...
		

Crossrefs

Columns k=0-10 give: A000007, A000041 (for n>0), A258457, A258458, A258459, A258460, A258461, A258462, A258463, A258464, A258465.
Row sums give A258466.
T(2n,n) give A258467.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1, k) +`if`(i>n, 0, k*b(n-i, i, k))))
        end:
    T:= (n, k)-> add(b(n$2, k-i)*(-1)^i/(i!*(k-i)!), i=0..k):
    seq(seq(T(n, k), k=0..n), n=0..10);
  • Mathematica
    b[n_, i_, k_] := b[n, i, k] = If[n==0, 1, If[i<1, 0, b[n, i-1, k] + If[i>n, 0, k*b[n-i, i, k]]]]; T[n_, k_] := Sum[b[n, n, k-i]*(-1)^i/(i!*(k-i)!), {i, 0, k}]; Table[Table[T[n, k], {k, 0, n}], {n, 0, 10}] // Flatten (* Jean-François Alcover, Feb 21 2016, after Alois P. Heinz *)

Formula

T(n,k) = A255970(n,k)/k! = (Sum_{i=0..k} (-1)^i * C(k,i) * A246935(n,k-i)) / A000142(k).

A320449 Number of antichains of sets whose multiset union is an integer partition of n.

Original entry on oeis.org

1, 1, 2, 4, 6, 9, 18, 24, 39, 58, 92, 131, 206
Offset: 0

Views

Author

Gus Wiseman, Oct 12 2018

Keywords

Examples

			The a(1) = 1 through a(7) = 24 antichains:
  {{1}}  {{2}}      {{3}}          {{4}}              {{5}}
         {{1},{1}}  {{1,2}}        {{1,3}}            {{1,4}}
                    {{1},{2}}      {{1},{3}}          {{2,3}}
                    {{1},{1},{1}}  {{2},{2}}          {{1},{4}}
                                   {{1},{1},{2}}      {{2},{3}}
                                   {{1},{1},{1},{1}}  {{1},{1},{3}}
                                                      {{1},{2},{2}}
                                                      {{1},{1},{1},{2}}
                                                      {{1},{1},{1},{1},{1}}
.
  {{6}}                      {{7}}
  {{1,5}}                    {{1,6}}
  {{2,4}}                    {{2,5}}
  {{1,2,3}}                  {{3,4}}
  {{1},{5}}                  {{1,2,4}}
  {{2},{4}}                  {{1},{6}}
  {{3},{3}}                  {{2},{5}}
  {{1},{2,3}}                {{3},{4}}
  {{2},{1,3}}                {{1},{2,4}}
  {{3},{1,2}}                {{2},{1,4}}
  {{1},{1},{4}}              {{4},{1,2}}
  {{1,2},{1,2}}              {{1},{1},{5}}
  {{1},{2},{3}}              {{1,2},{1,3}}
  {{2},{2},{2}}              {{1},{2},{4}}
  {{1},{1},{1},{3}}          {{1},{3},{3}}
  {{1},{1},{2},{2}}          {{2},{2},{3}}
  {{1},{1},{1},{1},{2}}      {{1},{1},{2,3}}
  {{1},{1},{1},{1},{1},{1}}  {{1},{1},{1},{4}}
                             {{1},{1},{2},{3}}
                             {{1},{2},{2},{2}}
                             {{1},{1},{1},{1},{3}}
                             {{1},{1},{1},{2},{2}}
                             {{1},{1},{1},{1},{1},{2}}
                             {{1},{1},{1},{1},{1},{1},{1}}
		

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]]]];
    submultisetQ[M_,N_]:=Or[Length[M]==0,MatchQ[{Sort[List@@M],Sort[List@@N]},{{x_,Z___},{_,x_,W___}}/;submultisetQ[{Z},{W}]]];
    antiQ[s_]:=Select[Tuples[s,2],And[UnsameQ@@#,submultisetQ@@#]&]=={};
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],And[And@@UnsameQ@@@#,antiQ[#]]&]],{n,10}]

A320353 Number of antichains of multisets whose multiset union is an integer partition of n.

Original entry on oeis.org

1, 1, 3, 5, 11, 17, 36, 56, 107, 175, 311, 505, 887
Offset: 0

Views

Author

Gus Wiseman, Oct 11 2018

Keywords

Examples

			The a(1) = 1 through a(5) = 17 antichains:
  {{1}}  {{2}}      {{3}}          {{4}}              {{5}}
         {{1,1}}    {{1,2}}        {{1,3}}            {{1,4}}
         {{1},{1}}  {{1,1,1}}      {{2,2}}            {{2,3}}
                    {{1},{2}}      {{1,1,2}}          {{1,1,3}}
                    {{1},{1},{1}}  {{1},{3}}          {{1,2,2}}
                                   {{2},{2}}          {{1},{4}}
                                   {{1,1,1,1}}        {{2},{3}}
                                   {{2},{1,1}}        {{1,1,1,2}}
                                   {{1,1},{1,1}}      {{1},{2,2}}
                                   {{1},{1},{2}}      {{3},{1,1}}
                                   {{1},{1},{1},{1}}  {{1,1,1,1,1}}
                                                      {{1,1},{1,2}}
                                                      {{1},{1},{3}}
                                                      {{1},{2},{2}}
                                                      {{2},{1,1,1}}
                                                      {{1},{1},{1},{2}}
                                                      {{1},{1},{1},{1},{1}}
		

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]]]];
    submultisetQ[M_,N_]:=Or[Length[M]==0,MatchQ[{Sort[List@@M],Sort[List@@N]},{{x_,Z___},{_,x_,W___}}/;submultisetQ[{Z},{W}]]];
    antiQ[s_]:=Select[Tuples[s,2],And[UnsameQ@@#,submultisetQ@@#]&]=={};
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],antiQ]],{n,8}]

A320355 Number of connected antichains of multisets whose multiset union is an integer partition of n.

Original entry on oeis.org

1, 1, 3, 4, 8, 9, 19, 24, 45, 71, 118, 194, 335
Offset: 0

Views

Author

Gus Wiseman, Oct 11 2018

Keywords

Examples

			The a(1) = 1 through a(5) = 9 clutters:
  {{1}}  {{2}}      {{3}}          {{4}}              {{5}}
         {{1,1}}    {{1,2}}        {{1,3}}            {{1,4}}
         {{1},{1}}  {{1,1,1}}      {{2,2}}            {{2,3}}
                    {{1},{1},{1}}  {{1,1,2}}          {{1,1,3}}
                                   {{2},{2}}          {{1,2,2}}
                                   {{1,1,1,1}}        {{1,1,1,2}}
                                   {{1,1},{1,1}}      {{1,1,1,1,1}}
                                   {{1},{1},{1},{1}}  {{1,1},{1,2}}
                                                      {{1},{1},{1},{1},{1}}
		

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]]]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    submultisetQ[M_,N_]:=Or[Length[M]==0,MatchQ[{Sort[List@@M],Sort[List@@N]},{{x_,Z___},{_,x_,W___}}/;submultisetQ[{Z},{W}]]];
    antiQ[s_]:=Select[Tuples[s,2],And[UnsameQ@@#,submultisetQ@@#]&]=={};
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],And[Length[csm[#]]==1,antiQ[#]]&]],{n,8}]

A320356 Number of strict connected antichains of multisets whose multiset union is an integer partition of n.

Original entry on oeis.org

1, 1, 2, 3, 5, 8, 13, 22, 35, 62, 98, 171, 277
Offset: 0

Views

Author

Gus Wiseman, Oct 11 2018

Keywords

Examples

			The a(1) = 1 through a(6) = 13 clutters:
  {{1}}  {{2}}    {{3}}      {{4}}        {{5}}          {{6}}
         {{1,1}}  {{1,2}}    {{1,3}}      {{1,4}}        {{1,5}}
                  {{1,1,1}}  {{2,2}}      {{2,3}}        {{2,4}}
                             {{1,1,2}}    {{1,1,3}}      {{3,3}}
                             {{1,1,1,1}}  {{1,2,2}}      {{1,1,4}}
                                          {{1,1,1,2}}    {{1,2,3}}
                                          {{1,1,1,1,1}}  {{2,2,2}}
                                          {{1,1},{1,2}}  {{1,1,1,3}}
                                                         {{1,1,2,2}}
                                                         {{1,1,1,1,2}}
                                                         {{1,1},{1,3}}
                                                         {{1,1,1,1,1,1}}
                                                         {{1,2},{1,1,1}}
		

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]]]];
    csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    submultisetQ[M_,N_]:=Or[Length[M]==0,MatchQ[{Sort[List@@M],Sort[List@@N]},{{x_,Z___},{_,x_,W___}}/;submultisetQ[{Z},{W}]]];
    antiQ[s_]:=Select[Tuples[s,2],And[UnsameQ@@#,submultisetQ@@#]&]=={};
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],And[UnsameQ@@#,Length[csm[#]]==1,antiQ[#]]&]],{n,8}]
Showing 1-10 of 18 results. Next