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.

Previous Showing 31-40 of 78 results. Next

A382079 Number of integer partitions of n that can be partitioned into a set of sets in exactly one way.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 3, 4, 6, 5, 10, 9, 13, 14, 21, 20, 32, 31, 42, 47, 63, 62, 90, 94, 117, 138, 170, 186, 235, 260, 315, 363, 429, 493, 588, 674, 795, 901, 1060, 1209, 1431, 1608, 1896, 2152, 2515, 2854, 3310, 3734, 4368, 4905, 5686
Offset: 0

Views

Author

Gus Wiseman, Mar 20 2025

Keywords

Examples

			The unique multiset partition for (3222111) is {{1},{2},{1,2},{1,2,3}}.
The a(1) = 1 through a(12) = 13 partitions:
  1  2  3  4    5    6     7    8      9      A      B      C
           211  221  411   322  332    441    433    443    552
                311  2211  331  422    522    442    533    633
                           511  611    711    622    551    822
                                3311   42111  811    722    A11
                                32111         3322   911    4422
                                              4411   42221  5511
                                              32221  53111  33321
                                              43111  62111  52221
                                              52111         54111
                                                            63111
                                                            72111
                                                            3222111
		

Crossrefs

Normal multiset partitions of this type are counted by A116539, see A381718.
These partitions are ranked by A293511.
MM-numbers of these multiset partitions (sets of sets) are A302494, see A302478, A382201.
Twice-partitions of this type (sets of sets) are counted by A358914, see A279785.
For at least one choice we have A382077 (ranks A382200), see A381992 (ranks A382075).
For no choices we have A382078 (ranks A293243), see A381990 (ranks A381806).
For distinct block-sums instead of blocks we have A382460, ranked by A381870.
Set multipartitions: A089259, A116540, A270995, A296119, A318360.
A000041 counts integer partitions, strict A000009.
A050320 counts multiset partitions of prime indices into sets.
A050326 counts multiset partitions of prime indices into distinct sets, see A381633.
A265947 counts refinement-ordered pairs of integer partitions.

Programs

  • Mathematica
    ssfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[ssfacs[n/d],Min@@#>d&],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
    Table[Length[Select[IntegerPartitions[n],Length[ssfacs[Times@@Prime/@#]]==1&]],{n,0,15}]

Extensions

a(21)-a(50) from Bert Dobbelaere, Mar 29 2025

A336139 Number of ways to choose a strict composition of each part of a strict composition of n.

Original entry on oeis.org

1, 1, 1, 5, 9, 17, 45, 81, 181, 397, 965, 1729, 3673, 7313, 15401, 34065, 68617, 135069, 266701, 556969, 1061921, 2434385, 4436157, 9120869, 17811665, 35651301, 68949549, 136796317, 283612973, 537616261, 1039994921, 2081261717, 3980842425, 7723253181, 15027216049
Offset: 0

Views

Author

Gus Wiseman, Jul 16 2020

Keywords

Comments

A strict composition of n is a finite sequence of distinct positive integers summing to n.

Examples

			The a(1) = 1 through a(5) = 17 splittings:
  (1)  (2)  (3)      (4)        (5)
            (1,2)    (1,3)      (1,4)
            (2,1)    (3,1)      (2,3)
            (1),(2)  (1),(3)    (3,2)
            (2),(1)  (3),(1)    (4,1)
                     (1),(1,2)  (1),(4)
                     (1),(2,1)  (2),(3)
                     (1,2),(1)  (3),(2)
                     (2,1),(1)  (4),(1)
                                (1),(1,3)
                                (1,2),(2)
                                (1),(3,1)
                                (1,3),(1)
                                (2),(1,2)
                                (2,1),(2)
                                (2),(2,1)
                                (3,1),(1)
		

Crossrefs

The version for partitions is A063834.
Row sums of A072574.
The version for non-strict compositions is A133494.
The version for strict partitions is A279785.
Multiset partitions of partitions are A001970.
Strict compositions are A032020.
Taking a composition of each part of a partition: A075900.
Taking a composition of each part of a strict partition: A304961.
Taking a strict composition of each part of a composition: A307068.
Splittings of partitions are A323583.
Compositions of parts of strict compositions are A336127.
Set partitions of strict compositions are A336140.

Programs

  • Mathematica
    strs[n_]:=Join@@Permutations/@Select[IntegerPartitions[n],UnsameQ@@#&];
    Table[Length[Join@@Table[Tuples[strs/@ctn],{ctn,strs[n]}]],{n,0,15}]

A358824 Number of twice-partitions of n of odd length.

Original entry on oeis.org

0, 1, 2, 4, 7, 15, 32, 61, 121, 260, 498, 967, 1890, 3603, 6839, 12972, 23883, 44636, 82705, 150904, 275635, 501737, 905498, 1628293, 2922580, 5224991, 9296414, 16482995, 29125140, 51287098, 90171414, 157704275, 275419984, 479683837, 833154673, 1442550486, 2493570655
Offset: 0

Views

Author

Gus Wiseman, Dec 03 2022

Keywords

Comments

A twice-partition of n is a sequence of integer partitions, one of each part of an integer partition of n.

Examples

			The a(1) = 1 through a(5) = 15 twice-partitions:
  (1)  (2)   (3)        (4)         (5)
       (11)  (21)       (22)        (32)
             (111)      (31)        (41)
             (1)(1)(1)  (211)       (221)
                        (1111)      (311)
                        (2)(1)(1)   (2111)
                        (11)(1)(1)  (11111)
                                    (2)(2)(1)
                                    (3)(1)(1)
                                    (11)(2)(1)
                                    (2)(11)(1)
                                    (21)(1)(1)
                                    (11)(11)(1)
                                    (111)(1)(1)
                                    (1)(1)(1)(1)(1)
		

Crossrefs

The version for set partitions is A024429.
For odd lengths (instead of length) we have A358334.
The case of odd parts also is A358823.
The case of odd sums also is A358826.
The case of odd lengths also is A358834.
For multiset partitions of integer partitions: A358837, ranked by A026424.
A000009 counts partitions into odd parts.
A027193 counts partitions of odd length.
A063834 counts twice-partitions, strict A296122, row-sums of A321449.
A078408 counts odd-length partitions into odd parts.
A300301 aerated counts twice-partitions with odd sums and parts.

Programs

  • Mathematica
    twiptn[n_]:=Join@@Table[Tuples[IntegerPartitions/@ptn],{ptn,IntegerPartitions[n]}];
    Table[Length[Select[twiptn[n],OddQ[Length[#]]&]],{n,0,10}]
  • PARI
    R(u,y) = {1/prod(k=1, #u, 1 - u[k]*y*x^k + O(x*x^#u))}
    seq(n) = {my(u=vector(n,k,numbpart(k))); Vec(R(u, 1) - R(u, -1), -(n+1))/2} \\ Andrew Howroyd, Dec 30 2022

Formula

G.f.: ((1/Product_{k>=1} (1-A000041(k)*x^k)) - (1/Product_{k>=1} (1+A000041(k)*x^k)))/2. - Andrew Howroyd, Dec 30 2022

Extensions

Terms a(26) and beyond from Andrew Howroyd, Dec 30 2022

A381996 Number of non-isomorphic multisets of size n that can be partitioned into a set of sets.

Original entry on oeis.org

1, 1, 1, 2, 3, 4, 6, 9, 13, 18, 25, 34, 47
Offset: 0

Views

Author

Gus Wiseman, Mar 31 2025

Keywords

Comments

First differs from A382523 at a(12) = 47, A382523(12) = 45.
We call a multiset non-isomorphic iff it covers an initial interval of positive integers with weakly decreasing multiplicities. The size of a multiset is the number of elements, counting multiplicity.

Examples

			Differs from A382523 in counting the following under a(12):
  {1,1,1,1,1,1,2,2,3,3,4,5} with partition {{1},{1,2},{1,3},{1,4},{1,5},{1,2,3}}
  {1,1,1,1,2,2,2,2,3,3,3,3} with partition {{1},{2},{3},{1,2},{1,3},{2,3},{1,2,3}}
		

Crossrefs

Factorizations of this type are counted by A050326, distinct sums A381633.
Normal multiset partitions of this type are counted by A116539, distinct sums A381718.
The complement is counted by A292444.
Twice-partitions of this type are counted by A358914, distinct sums A279785.
For integer partitions we have A382077, ranks A382200, complement A382078, ranks A293243.
Weak version is A382214, complement A292432, distinct sums A382216, complement A382202.
For distinct sums we have A382523, complement A382430.
Normal multiset partitions: A034691, A035310, A116540, A255906.
Set systems: A050342, A296120, A318361.
Set multipartitions: A089259, A270995, A296119, A318360.

Programs

  • Mathematica
    strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n];
    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[strnorm[n], Select[mps[#], UnsameQ@@#&&And@@UnsameQ@@@#&]!={}&]], {n,0,5}]

A382204 Number of normal multiset partitions of weight n into constant blocks with a common sum.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 7, 5, 8, 8, 10, 8, 15, 9, 14, 15, 17, 13, 22, 14, 25, 21, 23, 19, 34, 24, 29, 28, 37, 27, 45, 29, 44, 38, 43, 43, 59, 40, 51, 48, 69, 48, 71, 52, 73, 69, 72, 61, 93, 72, 91, 77, 99, 78, 105, 95, 119, 95, 113, 96, 146, 107, 126, 123, 151, 130
Offset: 0

Views

Author

Gus Wiseman, Mar 26 2025

Keywords

Comments

We call a multiset or multiset partition normal iff it covers an initial interval of positive integers. The weight of a multiset partition is the sum of sizes of its blocks.

Examples

			The a(1) = 1 through a(6) = 7 multiset partitions:
  {1} {11}   {111}     {1111}       {11111}         {111111}
      {1}{1} {2}{11}   {11}{11}     {2}{11}{11}     {111}{111}
             {1}{1}{1} {2}{2}{11}   {2}{2}{2}{11}   {22}{1111}
                       {1}{1}{1}{1} {1}{1}{1}{1}{1} {11}{11}{11}
                                                    {2}{2}{11}{11}
                                                    {2}{2}{2}{2}{11}
                                                    {1}{1}{1}{1}{1}{1}
The a(1) = 1 through a(7) = 5 factorizations:
  2  4    8      16       32         64           128
     2*2  3*4    4*4      3*4*4      8*8          3*4*4*4
          2*2*2  3*3*4    3*3*3*4    9*16         3*3*3*4*4
                 2*2*2*2  2*2*2*2*2  4*4*4        3*3*3*3*3*4
                                     3*3*4*4      2*2*2*2*2*2*2
                                     3*3*3*3*4
                                     2*2*2*2*2*2
		

Crossrefs

Without a common sum we have A055887.
Twice-partitions of this type are counted by A279789.
Without constant blocks we have A326518.
For distinct block-sums and strict blocks we have A381718.
Factorizations of this type are counted by A381995.
For distinct instead of equal block-sums we have A382203.
For strict instead of constant blocks we have A382429.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count multiset partitions of prime indices, strict A045778.
A089259 counts set multipartitions of integer partitions.
A255906 counts normal multiset partitions, row sums of A317532.
A321469 counts multiset partitions with distinct block-sums, ranks A326535.
Normal multiset partitions: A035310, A304969, A356945.
Set multipartitions: A116540, A270995, A296119, A318360.
Set multipartitions with distinct sums: A279785, A381806, A381870.
Constant blocks with distinct sums: A381635, A381636, A381716.

Programs

  • Mathematica
    allnorm[n_Integer]:=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[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Table[Length[Join@@(Select[mps[#],SameQ@@Total/@#&&And@@SameQ@@@#&]&/@allnorm[n])],{n,0,5}]
  • PARI
    h(s,x)=my(t=0,p=1,k=1);while(s%k==0,p*=1/(1-x^(s/k))-1;t+=p;k+=1);t
    lista(n)=Vec(1+sum(s=1,n,h(s,x+O(x*x^n)))) \\ Christian Sievers, Apr 05 2025

Formula

G.f.: 1 + Sum_{s>=1} Sum_{k=1..A055874(s)} Product_{v=1..k} (1/(1-x^(s/v)) - 1). - Christian Sievers, Apr 05 2025

Extensions

Terms a(16) and beyond from Christian Sievers, Apr 04 2025

A300436 Number of odd p-trees of weight n (all proper terminal subtrees have odd weight).

Original entry on oeis.org

1, 1, 1, 2, 2, 5, 5, 12, 13, 35, 37, 98, 107, 304, 336, 927, 1037, 3010, 3367, 9585, 10924, 32126, 36438, 105589, 121045, 359691, 412789, 1211214, 1398168, 4188930, 4831708, 14315544, 16636297, 50079792, 58084208, 173370663, 202101971, 611487744, 712709423
Offset: 1

Views

Author

Gus Wiseman, Mar 05 2018

Keywords

Comments

An odd p-tree of weight n > 0 is either a single node (if n = 1) or a finite sequence of at least 3 odd p-trees whose weights are weakly decreasing odd numbers summing to n.

Examples

			The a(7) = 5 odd p-trees: ((ooo)(ooo)o), (((ooo)oo)oo), ((ooooo)oo), ((ooo)oooo), (ooooooo).
		

Crossrefs

Programs

  • Mathematica
    b[n_]:=b[n]=If[n>1,0,1]+Sum[Times@@b/@y,{y,Select[IntegerPartitions[n],Length[#]>1&&And@@OddQ/@#&]}];
    Table[b[n],{n,40}]

Formula

O.g.f: x + Product_{n odd} 1/(1 - a(n)*x^n) - Sum_{n odd} a(n)*x^n. - Gus Wiseman, Aug 27 2018

Extensions

Name corrected by Gus Wiseman, Aug 27 2018

A307068 Expansion of 1/(1 - Sum_{k>=1} k!*x^(k*(k+1)/2) / Product_{j=1..k} (1 - x^j)).

Original entry on oeis.org

1, 1, 2, 6, 14, 34, 88, 216, 532, 1322, 3290, 8142, 20192, 50080, 124144, 307878, 763474, 1893038, 4694060, 11639580, 28861736, 71567206, 177460750, 440037738, 1091134276, 2705618900, 6708953156, 16635775698, 41250705518, 102286806130, 253634237896, 628921097352, 1559496588628
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 22 2019

Keywords

Comments

Invert transform of A032020.
Number of ways to choose a strict composition of each part of a composition of n. - Gus Wiseman, Jul 18 2020
The Invert transform T(a) of a sequence a is given by T(a)n = Sum_c Product_i a(c_i), where the sum is over all compositions c of n. - _Gus Wiseman, Aug 01 2020

Examples

			From _Gus Wiseman_, Jul 18 2020: (Start)
The a(1) = 1 through a(4) = 14 ways to choose a strict composition of each part of a composition:
    (1)  (2)      (3)          (4)
         (1),(1)  (1,2)        (1,3)
                  (2,1)        (3,1)
                  (1),(2)      (1),(3)
                  (2),(1)      (2),(2)
                  (1),(1),(1)  (3),(1)
                               (1),(1,2)
                               (1),(2,1)
                               (1,2),(1)
                               (2,1),(1)
                               (1),(1),(2)
                               (1),(2),(1)
                               (2),(1),(1)
                               (1),(1),(1),(1)
(End)
		

Crossrefs

The version for partitions is A270995.
Starting with a strict composition gives A336139.
Strict compositions are counted by A032020.
Partitions of each part of a partition are A063834.
Compositions of each part of a partition are A075900.
Compositions of each part of a composition are A133494.
Strict partitions of each part of a strict partition are A279785.
Compositions of each part of a strict partition are A304961.
Strict partitions of each part of a composition are A304969.
Compositions of each part of a strict composition are A336127.
Set partitions of strict compositions are A336140.
Strict compositions of each part of a partition are A336141.

Programs

  • Magma
    m:=80;
    R:=PowerSeriesRing(Integers(), m);
    Coefficients(R!( 1/(1 - (&+[Factorial(k)*x^Binomial(k+1,2)/(&*[ 1-x^j: j in [1..k]]): k in [1..m+2]]) ) )); // G. C. Greubel, Jan 25 2024
    
  • Maple
    T:= proc(n, k) option remember; `if`(k<0 or n<0, 0,
          `if`(k=0, `if`(n=0, 1, 0), T(n-k, k) +k*T(n-k, k-1)))
        end:
    g:= proc(n) option remember; add(T(n, k), k=0..floor((sqrt(8*n+1)-1)/2)) end:
    a:= proc(n) option remember; `if`(n<1, 1,
          add(a(n-i)*g(i), i=1..n))
        end:
    seq(a(n), n=0..32);  # Alois P. Heinz, Dec 16 2022
  • Mathematica
    nmax = 32; CoefficientList[Series[1/(1 - Sum[k!*x^(k*(k+1)/2)/Product[ (1-x^j), {j,k}], {k,nmax}]), {x, 0, nmax}], x]
  • SageMath
    m=80;
    def p(x, j): return product(1-x^k for k in range(1,j+1))
    def f(x): return 1/(1 - sum(factorial(j)*x^binomial(j+1,2)/p(x,j) for j in range(1, m+3)) )
    def A307068_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( f(x) ).list()
    A307068_list(m) # G. C. Greubel, Jan 25 2024

Formula

a(0) = 1; a(n) = Sum_{k=1..n} A032020(k)*a(n-k).

A330459 Number of set partitions of set-systems with total sum n.

Original entry on oeis.org

1, 1, 1, 4, 6, 11, 26, 42, 78, 148, 280, 481, 867, 1569, 2742, 4933, 8493, 14857, 25925, 44877, 77022, 132511, 226449, 385396, 657314, 1111115, 1875708, 3157379, 5309439, 8885889, 14861478, 24760339, 41162971, 68328959, 113099231, 186926116, 308230044
Offset: 0

Views

Author

Gus Wiseman, Dec 17 2019

Keywords

Comments

Number of sets of disjoint nonempty sets of nonempty sets of positive integers with total sum n.

Examples

			The a(6) = 26 partitions:
  ((6))  ((15))      ((123))          ((1)(2)(12))
         ((24))      ((1)(14))        ((1))((2)(12))
         ((1)(5))    ((1)(23))        ((12))((1)(2))
         ((2)(4))    ((2)(13))        ((2))((1)(12))
         ((1))((5))  ((3)(12))        ((1))((2))((12))
         ((2))((4))  ((1))((14))
                     ((1))((23))
                     ((1)(2)(3))
                     ((2))((13))
                     ((3))((12))
                     ((1))((2)(3))
                     ((2))((1)(3))
                     ((3))((1)(2))
                     ((1))((2))((3))
		

Crossrefs

Programs

  • Mathematica
    ppl[n_,k_]:=Switch[k,0,{n},1,IntegerPartitions[n],_,Join@@Table[Union[Sort/@Tuples[ppl[#,k-1]&/@ptn]],{ptn,IntegerPartitions[n]}]];
    Table[Length[Select[ppl[n,3],And[UnsameQ@@Join@@#,And@@UnsameQ@@@Join@@#]&]],{n,0,10}]
  • PARI
    \\ here L is A000009 and BellP is A000110 as series.
    L(n)={eta(x^2 + O(x*x^n))/eta(x + O(x*x^n))}
    BellP(n)={serlaplace(exp( exp(x + O(x*x^n)) - 1))}
    seq(n)={my(c=L(n), b=BellP(n), v=Vec(prod(k=1, n, (1 + x^k*y + O(x*x^n))^polcoef(c, k)))); vector(#v, n, my(r=v[n]); sum(k=0, n-1, polcoeff(b,k)*polcoef(r,k)))} \\ Andrew Howroyd, Dec 29 2019

Formula

a(n) = Sum_k A330462(n,k) * A000110(k).

Extensions

Terms a(18) and beyond from Andrew Howroyd, Dec 29 2019

A382429 Number of normal multiset partitions of weight n into sets with a common sum.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 13, 26, 57, 113, 283, 854, 2401, 6998, 24072, 85061, 308956, 1190518, 4770078, 19949106, 87059592
Offset: 0

Views

Author

Gus Wiseman, Mar 26 2025

Keywords

Comments

We call a multiset or multiset partition normal iff it covers an initial interval of positive integers. The weight of a multiset partition is the sum of sizes of its blocks.

Examples

			The a(1) = 1 through a(6) = 13 partitions:
  {1} {12}   {123}     {1234}       {12345}         {123456}
      {1}{1} {3}{12}   {12}{12}     {24}{123}       {123}{123}
             {1}{1}{1} {14}{23}     {34}{124}       {125}{134}
                       {3}{3}{12}   {3}{12}{12}     {135}{234}
                       {1}{1}{1}{1} {5}{14}{23}     {145}{235}
                                    {3}{3}{3}{12}   {12}{12}{12}
                                    {1}{1}{1}{1}{1} {14}{14}{23}
                                                    {14}{23}{23}
                                                    {16}{25}{34}
                                                    {3}{3}{12}{12}
                                                    {5}{5}{14}{23}
                                                    {3}{3}{3}{3}{12}
                                                    {1}{1}{1}{1}{1}{1}
The corresponding factorizations:
  2  6    30     210      2310       30030
     2*2  5*6    6*6      21*30      30*30
          2*2*2  14*15    35*42      6*6*6
                 5*5*6    5*6*6      66*70
                 2*2*2*2  5*5*5*6    110*105
                          11*14*15   154*165
                          2*2*2*2*2  5*5*6*6
                                     14*14*15
                                     14*15*15
                                     26*33*35
                                     5*5*5*5*6
                                     11*11*14*15
                                     2*2*2*2*2*2
		

Crossrefs

Without the common sum we have A116540 (normal set multipartitions).
Twice-partitions of this type are counted by A279788.
For common sizes instead of sums we have A317583.
Without strict blocks we have A326518, non-strict blocks A326517.
For a common length instead of sum we have A331638.
For distinct instead of equal block-sums we have A381718.
Factorizations of this type are counted by A382080.
For distinct block-sums and constant blocks we have A382203.
For constant instead of strict blocks we have A382204.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count multiset partitions of prime indices, strict A045778.
A321469 counts multiset partitions with distinct block-sums, ranks A326535.
Normal multiset partitions: A035310, A255906, A304969, A317532.
Set multipartitions: A089259, A116539, A270995, A296119, A318360.
Set multipartitions with distinct sums: A279785, A381806, A381870.
Constant blocks with distinct sums: A381635, A381636, A381716.

Programs

  • Mathematica
    allnorm[n_Integer]:=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[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Table[Length[Join@@(Select[mps[#],SameQ@@Total/@#&&And@@UnsameQ@@@#&]&/@allnorm[n])],{n,0,5}]

Extensions

a(11) from Robert Price, Mar 30 2025
a(12)-a(20) from Christian Sievers, Apr 06 2025

A330462 Triangle read by rows where T(n,k) is the number of k-element sets of nonempty sets of positive integers with total sum n.

Original entry on oeis.org

1, 0, 1, 0, 1, 0, 0, 2, 1, 0, 0, 2, 2, 0, 0, 0, 3, 4, 0, 0, 0, 0, 4, 6, 2, 0, 0, 0, 0, 5, 11, 3, 0, 0, 0, 0, 0, 6, 16, 8, 0, 0, 0, 0, 0, 0, 8, 25, 15, 1, 0, 0, 0, 0, 0, 0, 10, 35, 28, 4, 0, 0, 0, 0, 0, 0, 0, 12, 52, 46, 9, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Dec 18 2019

Keywords

Examples

			Triangle begins:
  1
  0  1
  0  1  0
  0  2  1  0
  0  2  2  0  0
  0  3  4  0  0  0
  0  4  6  2  0  0  0
  0  5 11  3  0  0  0  0
  0  6 16  8  0  0  0  0  0
  0  8 25 15  1  0  0  0  0  0
  0 10 35 28  4  0  0  0  0  0  0
  ...
Row n = 7 counts the following set-systems:
  {{7}}      {{1},{6}}      {{1},{2},{4}}
  {{1,6}}    {{2},{5}}      {{1},{2},{1,3}}
  {{2,5}}    {{3},{4}}      {{1},{3},{1,2}}
  {{3,4}}    {{1},{1,5}}
  {{1,2,4}}  {{1},{2,4}}
             {{2},{1,4}}
             {{2},{2,3}}
             {{3},{1,3}}
             {{4},{1,2}}
             {{1},{1,2,3}}
             {{1,2},{1,3}}
		

Crossrefs

Programs

  • Mathematica
    ppl[n_,k_]:=Switch[k,0,{n},1,IntegerPartitions[n],_,Join@@Table[Union[Sort/@Tuples[ppl[#,k-1]&/@ptn]],{ptn,IntegerPartitions[n]}]];
    Table[Length[Select[ppl[n,2],And[UnsameQ@@#,And@@UnsameQ@@@#,Length[#]==k]&]],{n,0,10},{k,0,n}]
  • PARI
    L(n)={eta(x^2 + O(x*x^n))/eta(x + O(x*x^n))}
    A(n)={my(c=L(n), v=Vec(prod(k=1, n, (1 + x^k*y + O(x*x^n))^polcoef(c,k)))); vector(#v, n, Vecrev(v[n],n))}
    {my(T=A(12)); for(n=1, #T, print(T[n]))} \\ Andrew Howroyd, Dec 29 2019

Formula

G.f.: Product_{j>=1} (1 + y*x^j)^A000009(j). - Andrew Howroyd, Dec 29 2019
Previous Showing 31-40 of 78 results. Next