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 51-60 of 99 results. Next

A317653 Number of free pure symmetric multifunctions whose leaves are a normal multiset of size n.

Original entry on oeis.org

1, 3, 34, 602, 14872, 472138, 18323359, 840503724, 44489123726, 2668985463839, 178960530393633, 13263068003965046, 1076580864432281157, 94987639225399100006, 9051397653144246683937, 926407121115738135640677, 101357200280211387377806719, 11804887470887800839909147484
Offset: 1

Views

Author

Gus Wiseman, Aug 03 2018

Keywords

Comments

A multiset is normal if it spans an initial interval of positive integers. A free pure symmetric multifunction f in EPSM is either (case 1) a positive integer, or (case 2) an expression of the form h[g_1, ..., g_k] where k > 0, h is in EPSM, each of the g_i for i = 1, ..., k is in EPSM, and for i < j we have g_i <= g_j under a canonical total ordering of EPSM, such as the Mathematica ordering of expressions.

Examples

			The a(3) = 34 free pure symmetric multifunctions:
1[1[1]], 1[1,1], 1[1][1],
1[2[2]], 1[2,2], 2[1[2]], 2[2[1]], 2[1,2], 1[2][2], 2[1][2], 2[2][1],
1[1[2]], 1[2[1]], 1[1,2], 2[1[1]], 2[1,1], 1[1][2], 1[2][1], 2[1][1],
1[2[3]], 1[3[2]], 1[2,3], 2[1[3]], 2[3[1]], 2[1,3], 3[1[2]], 3[2[1]], 3[1,2], 1[2][3], 2[1][3], 1[3][2], 3[1][2], 2[3][1], 3[2][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]]]];
    exprUsing[m_]:=exprUsing[m]=If[Length[m]==0,{},If[Length[m]==1,{First[m]},Join@@Cases[Union[Table[PR[m[[s]],m[[Complement[Range[Length[m]],s]]]],{s,Take[Subsets[Range[Length[m]]],{2,-2}]}]],PR[h_,g_]:>Join@@Table[Apply@@@Tuples[{exprUsing[h],Union[Sort/@Tuples[exprUsing/@p]]}],{p,mps[g]}]]]];
    got[y_]:=Join@@Table[Table[i,{y[[i]]}],{i,Range[Length[y]]}];
    Table[Sum[Length[exprUsing[got[y]]],{y,Join@@Permutations/@IntegerPartitions[n]}],{n,6}]
  • PARI
    \\ here R(n,1) is A052893.
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    R(n,k)={my(v=[k]); for(n=2, n, my(t=EulerT(v)); v=concat(v, sum(k=1, n-1, v[k]*t[n-k]))); v}
    seq(n)={sum(k=1, n, R(n,k)*sum(r=k, n, binomial(r,k)*(-1)^(r-k)) )} \\ Andrew Howroyd, Sep 14 2018

Extensions

Terms a(8) and beyond from Andrew Howroyd, Sep 14 2018

A318559 Number of combinatory separations of the multiset of prime factors 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, 3, 1, 7, 2, 2, 2, 8, 1, 2, 2, 7, 1, 3, 1, 4, 4, 2, 1, 12, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 8, 1, 2, 4, 11, 2, 3, 1, 4, 2, 3, 1, 15, 1, 2, 4, 4, 2, 3, 1, 12, 5, 2, 1, 8, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Aug 28 2018

Keywords

Comments

A multiset is normal if it spans an initial interval of positive integers. The type of a multiset is the unique normal multiset that has the same sequence of multiplicities when its entries are taken in increasing order. For example the type of 335556 is 112223. A (headless) combinatory separation of a multiset m is a multiset of normal multisets {t_1,...,t_k} such that there exist multisets {s_1,...,s_k} with multiset union m and such that s_i has type t_i for each i = 1...k.

Examples

			The a(60) = 8 combinatory separations of {2,2,3,5}:
  {1123},
  {1,112}, {1,123}, {11,12}, {12,12},
  {1,1,11}, {1,1,12},
  {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]]]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    normize[m_]:=m/.Rule@@@Table[{Union[m][[i]],i},{i,Length[Union[m]]}];
    Table[Length[Union[Sort/@Map[normize,mps[primeMS[n]],{2}]]],{n,100}]

A323792 Number of non-isomorphic weight-n multisets of sets of sets.

Original entry on oeis.org

1, 1, 4, 11, 43, 145, 614, 2549
Offset: 0

Views

Author

Gus Wiseman, Jan 27 2019

Keywords

Comments

All sets and multisets must be finite, and only the outermost may be empty.
The weight of an atom is 1, and the weight of a multiset is the sum of weights of its elements, counting multiplicity.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(3) = 11 multiset partitions:
  {{1}}  {{12}}      {{123}}
         {{1}{2}}    {{1}{12}}
         {{1}}{{1}}  {{1}{23}}
         {{1}}{{2}}  {{1}}{{12}}
                     {{1}}{{23}}
                     {{1}{2}{3}}
                     {{1}}{{1}{2}}
                     {{1}}{{2}{3}}
                     {{1}}{{1}}{{1}}
                     {{1}}{{1}}{{2}}
                     {{1}}{{2}}{{3}}
		

Crossrefs

A330676 Number of balanced reduced multisystems of weight n and maximum depth whose atoms cover an initial interval of positive integers.

Original entry on oeis.org

1, 1, 2, 8, 70, 1012, 21944, 665708, 26917492, 1399033348, 90878863352, 7214384973908, 687197223963640, 77354805301801012, 10158257981179981304, 1539156284259756811748, 266517060496258245459352, 52301515332984084095078308, 11546416513975694879642736152
Offset: 0

Views

Author

Gus Wiseman, Dec 30 2019

Keywords

Comments

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. The weight of an atom is 1, while the weight of a multiset is the sum of weights of its elements.
A finite multiset is normal if it covers an initial interval of positive integers.

Examples

			The a(0) = 1 through a(3) = 8 multisystems:
  {}  {1}  {1,1}  {{1},{1,1}}
           {1,2}  {{1},{1,2}}
                  {{1},{2,2}}
                  {{1},{2,3}}
                  {{2},{1,1}}
                  {{2},{1,2}}
                  {{2},{1,3}}
                  {{3},{1,2}}
		

Crossrefs

Row sums of A330778.
The case with all atoms equal is A000111.
The case with all atoms different is A006472.
The version allowing all depths is A330655.
The unlabeled version is A330663.
The version where the atoms are the prime indices of n is A330665.
The strongly normal version is A330675.
The version where the degrees are the prime indices of n is A330728.
Multiset partitions of normal multisets are A255906.
Series-reduced rooted trees with normal leaves are A316651.

Programs

  • Mathematica
    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]]]];
    totm[m_]:=Prepend[Join@@Table[totm[p],{p,Select[mps[m],1
    				
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    R(n, k)={my(v=vector(n), u=vector(n)); v[1]=k; for(n=1, #v, for(i=n, #v, u[i] += v[i]*(-1)^(i-n)*binomial(i-1, n-1)); v=EulerT(v)); u}
    seq(n)={concat([1], sum(k=1, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k))))} \\ Andrew Howroyd, Dec 30 2020

Extensions

Terms a(8) and beyond from Andrew Howroyd, Dec 30 2019

A382203 Number of normal multiset partitions of weight n into constant multisets with distinct sums.

Original entry on oeis.org

1, 1, 2, 4, 9, 19, 37, 76, 159, 326, 671, 1376, 2815, 5759, 11774, 24083, 49249, 100632, 205490, 419420, 855799, 1745889, 3561867, 7268240, 14836127, 30295633, 61888616
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(4) = 9 multiset partitions:
  {{1}}  {{1,1}}    {{1,1,1}}      {{1,1,1,1}}
         {{1},{2}}  {{1},{1,1}}    {{1},{1,1,1}}
                    {{1},{2,2}}    {{1,1},{2,2}}
                    {{1},{2},{3}}  {{1},{2,2,2}}
                                   {{2},{1,1,1}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{3},{2,2}}
                                   {{1},{2},{3},{4}}
The a(5) = 19 factorizations:
  32  2*16  2*3*27   2*3*5*25  2*3*5*7*11
      4*8   2*4*9    2*3*5*9
      2*81  2*3*8    2*3*5*49
      4*27  2*3*125  2*3*7*25
      9*8   2*9*25
      3*16  2*5*27
            5*4*9
		

Crossrefs

Without distinct sums we have A055887.
Twice-partitions of this type are counted by A279786.
For distinct blocks instead of sums we have A304969.
Without constant blocks we have A326519.
Factorizations of this type are counted by A381635.
For strict instead of constant blocks we have A381718.
For equal instead of distinct block-sums we have A382204.
For equal block-sums and strict 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.
A321469 counts multiset partitions with distinct block-sums, ranks A326535.
Normal multiset partitions: A035310, A116540, A255906, A317532.
Set multipartitions with distinct sums: A279785, A381806, A381870.

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[#],UnsameQ@@Total/@#&&And@@SameQ@@@#&]&/@allnorm[n])],{n,0,5}]

Extensions

a(14)-a(26) from Christian Sievers, Apr 04 2025

A382428 Number of normal multiset partitions of weight n into sets with distinct sizes.

Original entry on oeis.org

1, 1, 1, 6, 8, 35, 292, 673, 2818, 16956, 219772, 636748, 3768505, 20309534, 183403268, 3227600747, 12272598308, 81353466578, 561187259734, 4416808925866, 50303004612136, 1238783066956740, 5566249468690291, 44970939483601100, 330144217684933896, 3131452652308459402
Offset: 0

Views

Author

Gus Wiseman, Mar 29 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(4) = 8 multiset partitions:
  {{1}}  {{1,2}}  {{1,2,3}}    {{1,2,3,4}}
                  {{1},{1,2}}  {{1},{1,2,3}}
                  {{1},{2,3}}  {{1},{2,3,4}}
                  {{2},{1,2}}  {{2},{1,2,3}}
                  {{2},{1,3}}  {{2},{1,3,4}}
                  {{3},{1,2}}  {{3},{1,2,3}}
                               {{3},{1,2,4}}
                               {{4},{1,2,3}}
		

Crossrefs

For distinct sums instead of sizes we have A116539, see A050326.
Without distinct lengths we have A116540 (normal set multipartitions).
Without strict blocks we have A326517, for sum instead of size A326519.
For equal instead of distinct sizes we have A331638.
Twice-partitions of this type are counted by A358830.
For distinct sums instead of sizes we have A381718.
For equal instead of distinct sizes we have A382429.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count factorizations, strict A045778.
Normal multiset partitions: A034691, A035310, A255906.
Set multipartitions: A089259, A270995, A296119, A318360.

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[#],UnsameQ@@Length/@#&&And@@UnsameQ@@@#&]&/@allnorm[n])],{n,0,5}]
  • PARI
    R(n, k)={Vec(prod(j=1, n, 1 + binomial(k, j)*x^j + O(x*x^n)))}
    seq(n)={sum(k=0, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)))} \\ Andrew Howroyd, Mar 31 2025

Extensions

a(10) onwards from Andrew Howroyd, Mar 31 2025

A320459 MM-numbers of labeled multigraphs spanning an initial interval of positive integers.

Original entry on oeis.org

1, 13, 169, 377, 611, 1363, 1937, 2021, 2117, 2197, 4901, 7943, 10933, 16211, 17719, 25181, 26273, 27521, 28561, 28717, 39527, 44603, 56173, 58609, 61393, 63713, 64061, 83291, 86903, 91031, 91039, 94987, 99499, 103259, 141401, 142129, 143663, 146653, 147533
Offset: 1

Views

Author

Gus Wiseman, Oct 13 2018

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 formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of terms together with their multiset multisystems begins:
      1: {}
     13: {{1,2}}
    169: {{1,2},{1,2}}
    377: {{1,2},{1,3}}
    611: {{1,2},{2,3}}
   1363: {{1,3},{2,3}}
   1937: {{1,2},{3,4}}
   2021: {{1,4},{2,3}}
   2117: {{1,3},{2,4}}
   2197: {{1,2},{1,2},{1,2}}
   4901: {{1,2},{1,2},{1,3}}
   7943: {{1,2},{1,2},{2,3}}
  10933: {{1,2},{1,3},{1,3}}
  16211: {{1,2},{1,3},{1,4}}
  17719: {{1,2},{1,3},{2,3}}
  25181: {{1,2},{1,2},{3,4}}
  26273: {{1,2},{1,4},{2,3}}
  27521: {{1,2},{1,3},{2,4}}
  28561: {{1,2},{1,2},{1,2},{1,2}}
  28717: {{1,2},{2,3},{2,3}}
  39527: {{1,3},{1,3},{2,3}}
  44603: {{1,2},{2,3},{2,4}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    normQ[sys_]:=Or[Length[sys]==0,Union@@sys==Range[Max@@Max@@sys]];
    Select[Range[100000],And[normQ[primeMS/@primeMS[#]],And@@(And[SquareFreeQ[#],Length[primeMS[#]]==2]&/@primeMS[#])]&]

A320532 MM-numbers of labeled hypergraphs with multiset edges and no singletons spanning an initial interval of positive integers.

Original entry on oeis.org

1, 7, 13, 19, 37, 53, 61, 89, 91, 113, 131, 133, 151, 161, 223, 247, 251, 259, 281, 299, 311, 329, 359, 371, 377, 427, 437, 463, 481, 503, 593, 611, 623, 659, 667, 689, 703, 719, 721, 791, 793, 827, 851, 863, 893, 917, 923, 953, 1007, 1057, 1069, 1073, 1157
Offset: 1

Views

Author

Gus Wiseman, Oct 14 2018

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 formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of terms together with their multiset multisystems begins:
    1: {}
    7: {{1,1}}
   13: {{1,2}}
   19: {{1,1,1}}
   37: {{1,1,2}}
   53: {{1,1,1,1}}
   61: {{1,2,2}}
   89: {{1,1,1,2}}
   91: {{1,1},{1,2}}
  113: {{1,2,3}}
  131: {{1,1,1,1,1}}
  133: {{1,1},{1,1,1}}
  151: {{1,1,2,2}}
  161: {{1,1},{2,2}}
  223: {{1,1,1,1,2}}
  247: {{1,2},{1,1,1}}
  251: {{1,2,2,2}}
  259: {{1,1},{1,1,2}}
  281: {{1,1,2,3}}
  299: {{1,2},{2,2}}
  311: {{1,1,1,1,1,1}}
  329: {{1,1},{2,3}}
  359: {{1,1,1,2,2}}
  371: {{1,1},{1,1,1,1}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    normQ[sys_]:=Or[Length[sys]==0,Union@@sys==Range[Max@@Max@@sys]];
    Select[Range[1000],And[SquareFreeQ[#],normQ[primeMS/@primeMS[#]],And@@(And[PrimeOmega[#]>1]&/@primeMS[#])]&]

A320533 MM-numbers of labeled multi-hypergraphs with multiset edges and no singletons spanning an initial interval of positive integers.

Original entry on oeis.org

1, 7, 13, 19, 37, 49, 53, 61, 89, 91, 113, 131, 133, 151, 161, 169, 223, 247, 251, 259, 281, 299, 311, 329, 343, 359, 361, 371, 377, 427, 437, 463, 481, 503, 593, 611, 623, 637, 659, 667, 689, 703, 719, 721, 791, 793, 827, 851, 863, 893, 917, 923, 931, 953
Offset: 1

Views

Author

Gus Wiseman, Oct 14 2018

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 formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of terms together with their multiset multisystems begins:
    1: {}
    7: {{1,1}}
   13: {{1,2}}
   19: {{1,1,1}}
   37: {{1,1,2}}
   49: {{1,1},{1,1}}
   53: {{1,1,1,1}}
   61: {{1,2,2}}
   89: {{1,1,1,2}}
   91: {{1,1},{1,2}}
  113: {{1,2,3}}
  131: {{1,1,1,1,1}}
  133: {{1,1},{1,1,1}}
  151: {{1,1,2,2}}
  161: {{1,1},{2,2}}
  169: {{1,2},{1,2}}
  223: {{1,1,1,1,2}}
  247: {{1,2},{1,1,1}}
  251: {{1,2,2,2}}
  259: {{1,1},{1,1,2}}
  281: {{1,1,2,3}}
  299: {{1,2},{2,2}}
  311: {{1,1,1,1,1,1}}
  329: {{1,1},{2,3}}
  343: {{1,1},{1,1},{1,1}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    normQ[sys_]:=Or[Length[sys]==0,Union@@sys==Range[Max@@Max@@sys]];
    Select[Range[1000],And[normQ[primeMS/@primeMS[#]],And@@(And[PrimeOmega[#]>1]&/@primeMS[#])]&]

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

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 0, 3, 2, 0, 0, 5, 4, 0, 0, 0, 7, 11, 1, 0, 0, 0, 11, 20, 6, 0, 0, 0, 0, 15, 40, 16, 0, 0, 0, 0, 0, 22, 68, 40, 3, 0, 0, 0, 0, 0, 30, 120, 91, 11, 0, 0, 0, 0, 0, 0, 42, 195, 186, 41, 0, 0, 0, 0, 0, 0, 0, 56, 320, 367, 105, 3, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Dec 19 2019

Keywords

Examples

			Triangle begins:
  1
  0  1
  0  2  0
  0  3  2  0
  0  5  4  0  0
  0  7 11  1  0  0
  0 11 20  6  0  0  0
  0 15 40 16  0  0  0  0
  0 22 68 40  3  0  0  0  0
  ...
Row n = 5 counts the following sets of multisets:
  {{5}}          {{1},{4}}        {{1},{2},{1,1}}
  {{1,4}}        {{2},{3}}
  {{2,3}}        {{1},{1,3}}
  {{1,1,3}}      {{1},{2,2}}
  {{1,2,2}}      {{2},{1,2}}
  {{1,1,1,2}}    {{3},{1,1}}
  {{1,1,1,1,1}}  {{1},{1,1,2}}
                 {{1,1},{1,2}}
                 {{2},{1,1,1}}
                 {{1},{1,1,1,1}}
                 {{1,1},{1,1,1}}
		

Crossrefs

Row sums are A261049.
Column k = 1 is A000041.
Multisets of multisets are A061260, with row sums A001970.
Sets of sets are A330462, with row sums A050342.
Multisets of sets are A285229, with row sums A089259.
Sets of disjoint sets are A330460, with row sums A294617.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, add(binomial(
           combinat[numbpart](i), j)*expand(b(n-i*j, i-1)*x^j), j=0..n/i)))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n$2)):
    seq(T(n), n=0..14);  # Alois P. Heinz, Dec 30 2019
  • 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@@#,Length[#]==k]&]],{n,0,10},{k,0,n}]
    (* Second program: *)
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i<1, 0, Sum[Binomial[
         PartitionsP[i], j]*Expand[b[n - i*j, i - 1]*x^j], {j, 0, n/i}]]];
    T[n_] := Function[p, Table[Coefficient[p, x, i], {i, 0, n}]][b[n, n]];
    T /@ Range[0, 14] // Flatten (* Jean-François Alcover, May 18 2021, after Alois P. Heinz *)
  • PARI
    A(n)={my(v=Vec(prod(k=1, n, (1 + x^k*y + O(x*x^n))^numbpart(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)^A000041(j). - Andrew Howroyd, Dec 29 2019
Previous Showing 51-60 of 99 results. Next