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

A255906 Number of collections of nonempty multisets with a total of n objects having color set {1,...,k} for some k<=n.

Original entry on oeis.org

1, 1, 4, 16, 76, 400, 2356, 15200, 106644, 806320, 6526580, 56231024, 513207740, 4941362512, 50013751812, 530481210672, 5880285873060, 67954587978448, 816935340368068, 10196643652651664, 131904973822724540, 1765645473517011568, 24420203895517396180
Offset: 0

Views

Author

Alois P. Heinz, Mar 10 2015

Keywords

Comments

Number of multiset partitions of normal multisets of size n, where a multiset is normal if it spans an initial interval of positive integers. - Gus Wiseman, Jul 30 2018

Examples

			a(0) = 1: {}.
a(1) = 1: {{1}}.
a(2) = 4: {{1},{1}}, {{1,1}}, {{1},{2}}, {{1,2}}.
a(3) = 16: {{1},{1},{1}}, {{1},{1,1}}, {{1,1,1}}, {{1},{1},{2}}, {{1},{2},{2}}, {{1},{1,2}}, {{1},{2,2}}, {{2},{1,1}}, {{2},{1,2}}, {{1,1,2}}, {{1,2,2}}, {{1},{2},{3}}, {{1},{2,3}}, {{2},{1,3}}, {{3},{1,2}}, {{1,2,3}}.

Links

Crossrefs

Row sums of A255903. Also row sums of A317532.
Cf. A007716, A034691, A116539, A317073.

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]]]];
allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
Table[Length[Join@@mps/@allnorm[n]],{n,6}] (* Gus Wiseman, Jul 30 2018 *)
  • PARI
    R(n, k)={Vec(-1 + 1/prod(j=1, n, (1 - x^j + O(x*x^n))^binomial(k+j-1, j) ))}
seq(n) = {concat([1], sum(k=1, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) ))} \\ Andrew Howroyd, Sep 23 2023

Formula

a(n) = Sum_{k=0..n} A255903(n,k).

A317533 Regular triangle read rows: T(n,k) = number of non-isomorphic multiset partitions of size n and length k.

Original entry on oeis.org

1, 2, 2, 3, 4, 3, 5, 14, 9, 5, 7, 28, 33, 16, 7, 11, 69, 104, 74, 29, 11, 15, 134, 294, 263, 142, 47, 15, 22, 285, 801, 948, 599, 263, 77, 22, 30, 536, 2081, 3058, 2425, 1214, 453, 118, 30, 42, 1050, 5212, 9769, 9276, 5552, 2322, 761, 181, 42, 56, 1918, 12645, 29538, 34172, 23770, 11545, 4179, 1223, 267, 56
Offset: 1

Views

Author

Gus Wiseman, Jul 30 2018

Keywords

Examples

			Non-isomorphic representatives of the T(3,2) = 4 multiset partitions:
  {{1},{1,1}}
  {{1},{1,2}}
  {{1},{2,2}}
  {{1},{2,3}}
Triangle begins:
    1
    2    2
    3    4    3
    5   14    9    5
    7   28   33   16    7
   11   69  104   74   29   11
   15  134  294  263  142   47   15

Links

Crossrefs

Row sums are A007716. First and last columns are both A000041.
Cf. A034691, A255397, A255903, A255906, A317532, A334550.

Programs

  • Mathematica
    permcount[v_List] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
c[p_List, q_List, k_] := SeriesCoefficient[1/Product[(1 - x^LCM[p[[i]], q[[j]]])^GCD[p[[i]], q[[j]]], {j, 1, Length[q]}, {i, 1, Length[p]}], {x, 0, k}];
M[m_, n_, k_] := Module[{s = 0}, Do[Do[s += permcount[p]*permcount[q]*c[p, q, k], {q, IntegerPartitions[n]}], {p, IntegerPartitions[m]}]; s/(m!*n!)];
T[n_, k_] := M[k, n, n] - M[k - 1, n, n];
Table[T[n, k], {n, 1, 11}, {k, 1, n}] // Flatten (* Jean-François Alcover, Feb 08 2020, after Andrew Howroyd *)
  • PARI
    \\ See A318795 for definition of M.
T(n,k)={M(k, n, n) - M(k-1, n, n)}
for(n=1, 10, for(k=1, n, print1(T(n,k),", "));print) \\ Andrew Howroyd, Dec 28 2019
  • PARI
    \\ Faster version.
permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
K(q, t, n)={1/prod(j=1, #q, (1-x^lcm(t, q[j]) + O(x*x^n))^gcd(t, q[j]))}
G(m,n)={my(s=0); forpart(q=m, s+=permcount(q)*exp(sum(t=1, n, (K(q, t, n)-1)/t) + O(x*x^n))); s/m!}
A(n,m=n)={my(p=sum(k=0, m, G(k,n)*y^k)*(1-y)); matrix(n, m, n, k, polcoef(polcoef(p, n, x), k, y))}
{ my(T=A(10)); for(n=1, #T, print(T[n,1..n])) } \\ Andrew Howroyd, Aug 30 2020

Extensions

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

A318284 Number of multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 5, 9, 7, 7, 11, 11, 12, 16, 15, 15, 26, 22, 21, 29, 19, 30, 36, 31, 30, 66, 38, 42, 52, 56, 52, 47, 45, 57, 92, 77, 67, 77, 74, 101, 98, 135, 64, 137, 97, 176, 135, 109, 109, 118, 105, 231, 249, 97, 141, 181, 139, 297, 198, 385, 195, 269
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2018

Keywords

Examples

			The a(12) = 11 multiset partitions of {1,1,2,3}:
  {{1,1,2,3}}
  {{1},{1,2,3}}
  {{2},{1,1,3}}
  {{3},{1,1,2}}
  {{1,1},{2,3}}
  {{1,2},{1,3}}
  {{1},{1},{2,3}}
  {{1},{2},{1,3}}
  {{1},{3},{1,2}}
  {{2},{3},{1,1}}
  {{1},{1},{2},{3}}

Links

Crossrefs

Cf. A001055, A001970, A007716, A181821, A219727, A255906, A269134, A305936.
Cf. A317532, A318283, A318285, A318286, A318360, A318361.

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[facs[Times@@Prime/@nrmptn[n]]],{n,60}]
  • PARI
    permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i,2], j, primepi(f[i,1]))))}
count(sig)={my(n=vecsum(sig), A=O(x*x^vecmax(sig)), s=0); forpart(p=n, my(q=1/prod(i=1, #p, 1 - x^p[i] + A)); s+=prod(i=1, #sig, polcoef(q,sig[i]))*permcount(p)); s/n!}
a(n)={if(n==1, 1, my(s=sig(n)); if(#s==1, numbpart(s[1]), count(sig(n))))} \\ Andrew Howroyd, Dec 10 2018

Formula

a(n) = A001055(A181821(n)).
a(prime(n)^k) = A219727(n,k). - Andrew Howroyd, Dec 10 2018

A381718 Number of normal multiset partitions of weight n into sets with distinct sums.

Original entry on oeis.org

1, 1, 2, 6, 23, 106, 549, 3184, 20353, 141615, 1063399, 8554800, 73281988, 665141182, 6369920854, 64133095134, 676690490875, 7462023572238, 85786458777923, 1025956348473929, 12739037494941490
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(3) = 6 multiset partitions:
  {{1}}  {{1,2}}    {{1,2,3}}
         {{1},{2}}  {{1},{1,2}}
                    {{1},{2,3}}
                    {{2},{1,2}}
                    {{2},{1,3}}
                    {{1},{2},{3}}
The a(4) = 23 factorizations:
  2*3*6  5*30    3*30    2*30    210
         10*15   6*15    6*10    2*105
         2*5*15  2*3*15  2*3*10  3*70
         3*5*10                  5*42
                                 7*30
                                 6*35
                                 10*21
                                 2*3*35
                                 2*5*21
                                 2*7*15
                                 3*5*14
                                 2*3*5*7

Crossrefs

For distinct blocks instead of sums we have A116539, see A050326.
Without distinct sums we have A116540 (normal set multipartitions).
Twice-partitions of this type are counted by A279785.
Without strict blocks we have A326519.
Factorizations of this type are counted by A381633.
For constant instead of strict blocks we have A382203.
For distinct sizes instead of sums we have A382428, non-strict blocks A326517.
For equal instead of distinct block-sums we have A382429, non-strict blocks A326518.
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.
Cf. A000110, A007716, A050320, A255903, A275780, A317532, A317583, A321469, A381635, A382204, A382214, A382216.

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

Extensions

a(10)-a(11) from Robert Price, Mar 31 2025
a(12)-a(20) from Christian Sievers, Apr 05 2025

A317752 Number of multiset partitions of normal multisets of size n such that the blocks have empty intersection.

Original entry on oeis.org

0, 1, 8, 49, 305, 1984, 13686, 100124, 776885, 6386677, 55532358, 509549386, 4921352952, 49899820572, 529807799836, 5876162077537, 67928460444139, 816764249684450, 10195486840926032, 131896905499007474, 1765587483656124106, 24419774819813602870
Offset: 1

Views

Author

Gus Wiseman, Aug 06 2018

Keywords

Comments

A multiset is normal if it spans an initial interval of positive integers.

Examples

			The a(3) = 8 multiset partitions with empty intersection:
  {{2},{1,1}}
  {{1},{2,2}}
  {{1},{2,3}}
  {{2},{1,3}}
  {{3},{1,2}}
  {{1},{1},{2}}
  {{1},{2},{2}}
  {{1},{2},{3}}

Links

Crossrefs

Cf. A007716, A255903, A255906, A317073, A281116, A317077, A317532, A317533.
Cf. A317748, A317751, A317755, A317757, A317776.

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]]]];
allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
Table[Length[Join@@Table[Select[mps[m],Intersection@@#=={}&],{m,allnorm[n]}]],{n,6}]
  • PARI
    P(n,k)={1/prod(i=1, n, (1 - x^i*y + O(x*x^n))^binomial(k+i-1, k-1))}
R(n,k)={my(p=P(n,k), q=p/(1-y+O(y*y^n))); Vec(sum(i=2, n, polcoef(p,i,y) + polcoef(q,i,y)*sum(j=1, n\i, (-1)^j*binomial(k,j)*x^(i*j))), -n)}
seq(n)={sum(k=2, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) )} \\ Andrew Howroyd, Feb 05 2021

Extensions

Terms a(9) and beyond from Andrew Howroyd, Feb 05 2021

A317755 Number of multiset partitions of strongly normal multisets of size n such that the blocks have empty intersection.

Original entry on oeis.org

0, 1, 6, 30, 130, 629, 2930, 15019, 78224, 438626, 2548481
Offset: 1

Views

Author

Gus Wiseman, Aug 06 2018

Keywords

Comments

A multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities.

Examples

			The a(3) = 6 strongly normal multiset partitions with empty intersection:
  {{2},{1,1}}
  {{1},{2,3}}
  {{2},{1,3}}
  {{3},{1,2}}
  {{1},{1},{2}}
  {{1},{2},{3}}

Crossrefs

Cf. A007716, A035310, A255906, A317073, A281116, A317077, A317078, A317532, A317533.
Cf. A317748, A317751, A317752, A317757, A317775.

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]]]];
strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
Table[Length[Join@@Table[Select[mps[m],Intersection@@#=={}&],{m,strnorm[n]}]],{n,6}]

Extensions

a(10)-a(11) from Robert Price, May 08 2021

A055884 Euler transform of partition triangle A008284.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 4, 4, 5, 1, 4, 8, 7, 7, 1, 6, 12, 16, 12, 11, 1, 6, 17, 25, 28, 19, 15, 1, 8, 22, 43, 49, 48, 30, 22, 1, 8, 30, 58, 87, 88, 77, 45, 30, 1, 10, 36, 87, 134, 167, 151, 122, 67, 42, 1, 10, 45, 113, 207, 270, 296, 247, 185, 97, 56, 1, 12, 54, 155, 295, 448, 510, 507, 394, 278, 139, 77
Offset: 1

Views

Author

Christian G. Bower, Jun 09 2000

Keywords

Comments

Number of multiset partitions of length-k integer partitions of n. - Gus Wiseman, Nov 09 2018

Examples

			From _Gus Wiseman_, Nov 09 2018: (Start)
Triangle begins:
   1
   1   2
   1   2   3
   1   4   4   5
   1   4   8   7   7
   1   6  12  16  12  11
   1   6  17  25  28  19  15
   1   8  22  43  49  48  30  22
   1   8  30  58  87  88  77  45  30
   ...
The fifth row {1, 4, 8, 7, 7} counts the following multiset partitions:
  {{5}}   {{1,4}}     {{1,1,3}}       {{1,1,1,2}}         {{1,1,1,1,1}}
          {{2,3}}     {{1,2,2}}      {{1},{1,1,2}}       {{1},{1,1,1,1}}
         {{1},{4}}   {{1},{1,3}}     {{1,1},{1,2}}       {{1,1},{1,1,1}}
         {{2},{3}}   {{1},{2,2}}     {{2},{1,1,1}}      {{1},{1},{1,1,1}}
                     {{2},{1,2}}    {{1},{1},{1,2}}     {{1},{1,1},{1,1}}
                     {{3},{1,1}}    {{1},{2},{1,1}}    {{1},{1},{1},{1,1}}
                    {{1},{1},{3}}  {{1},{1},{1},{2}}  {{1},{1},{1},{1},{1}}
                    {{1},{2},{2}}
(End)

Links

Crossrefs

Row sums give A001970.
Main diagonal gives A000041.
Columns k=1-2 give: A057427, A052928.
T(n+2,n+1) gives A000070.
T(2n,n) gives A360468.
Cf. A055885, A055886, A360742.
Cf. A000219, A007716, A008284, A255906, A317449, A317532, A317533, A320796, A320801, A320808.

Programs

  • Maple
    h:= proc(n, i) option remember; expand(`if`(n=0, 1,
      `if`(i<1, 0, h(n, i-1)+x*h(n-i, min(n-i, i)))))
    end:
g:= proc(n, i, j) option remember; expand(`if`(j=0, 1, `if`(i<0, 0, add(
      g(n, i-1, j-k)*x^(i*k)*binomial(coeff(h(n$2), x, i)+k-1, k), k=0..j))))
    end:
b:= proc(n, i) option remember; expand(`if`(n=0, 1,
     `if`(i<1, 0, add(b(n-i*j, i-1)*g(i$2, j), j=0..n/i))))
    end:
T:= (n, k)-> coeff(b(n$2), x, k):
seq(seq(T(n,k), k=1..n), n=1..12);  # Alois P. Heinz, Feb 17 2023
  • 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]]]];
Table[Length[Join@@mps/@IntegerPartitions[n,{k}]],{n,5},{k,n}] (* Gus Wiseman, Nov 09 2018 *)

A382216 Number of normal multisets of size n that can be partitioned into a set of sets with distinct sums.

Original entry on oeis.org

1, 1, 1, 3, 5, 11, 23, 48, 101, 208, 434
Offset: 0

Views

Author

Gus Wiseman, Mar 29 2025

Keywords

Comments

We call a multiset normal iff it covers an initial interval of positive integers. The size of a multiset is the number of elements, counting multiplicity.

Examples

			The multiset {1,2,2,3,3} can be partitioned into a set of sets with distinct sums in 4 ways:
  {{2,3},{1,2,3}}
  {{2},{3},{1,2,3}}
  {{2},{1,3},{2,3}}
  {{1},{2},{3},{2,3}}
so is counted under a(5).
The multisets counted by A382214 but not by A382216 are:
  {1,1,1,1,2,2,3,3,3}
  {1,1,2,2,2,2,3,3,3}
The a(1) = 1 through a(5) = 11 multisets:
  {1}  {1,2}  {1,1,2}  {1,1,2,2}  {1,1,1,2,3}
              {1,2,2}  {1,1,2,3}  {1,1,2,2,3}
              {1,2,3}  {1,2,2,3}  {1,1,2,3,3}
                       {1,2,3,3}  {1,1,2,3,4}
                       {1,2,3,4}  {1,2,2,2,3}
                                  {1,2,2,3,3}
                                  {1,2,2,3,4}
                                  {1,2,3,3,3}
                                  {1,2,3,3,4}
                                  {1,2,3,4,4}
                                  {1,2,3,4,5}

Crossrefs

Twice-partitions of this type are counted by A279785, without distinct sums A358914.
Factorizations of this type are counted by A381633, without distinct sums A050326.
Normal multiset partitions of this type are counted by A381718, A116539.
The complement is counted by A382202.
Without distinct sums we have A382214, complement A292432.
The case of a unique choice is counted by A382459, without distinct sums A382458.
For Heinz numbers: A293243, A381806, A382075, A382200.
For integer partitions: A381990, A381992, A382077, A382078.
Strong version: A382523, A382430, A381996, A292444.
Normal multiset partitions: A034691, A035310, A255906.
Set systems: A050342, A296120, A318361.
Set multipartitions: A089259, A270995, A296119, A318360.
Cf. A000110, A000670, A007716, A050320, A116540, A255903, A275780, A317532, A326518, A326519, A382429, A382460.

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

A321449 Regular triangle read by rows where T(n,k) is the number of twice-partitions of n with a combined total of k parts.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 1, 2, 3, 0, 1, 4, 5, 5, 0, 1, 4, 8, 8, 7, 0, 1, 6, 13, 19, 16, 11, 0, 1, 6, 17, 27, 32, 24, 15, 0, 1, 8, 24, 47, 61, 62, 41, 22, 0, 1, 8, 30, 63, 99, 111, 100, 61, 30, 0, 1, 10, 38, 94, 158, 209, 210, 170, 95, 42, 0, 1, 10, 45, 119, 229, 328, 382, 348, 259, 136, 56
Offset: 0

Views

Author

Gus Wiseman, Nov 10 2018

Keywords

Comments

A twice partition of n (A063834) is a choice of an integer partition of each part in an integer partition of n.

Examples

			Triangle begins:
   1
   0   1
   0   1   2
   0   1   2   3
   0   1   4   5   5
   0   1   4   8   8   7
   0   1   6  13  19  16  11
   0   1   6  17  27  32  24  15
   0   1   8  24  47  61  62  41  22
   0   1   8  30  63  99 111 100  61  30
The sixth row {0, 1, 6, 13, 19, 16, 11} counts the following twice-partitions:
  (6)  (33)    (222)      (2211)        (21111)          (111111)
       (42)    (321)      (3111)        (1111)(2)        (111)(111)
       (51)    (411)      (111)(3)      (111)(21)        (1111)(11)
       (3)(3)  (21)(3)    (211)(2)      (21)(111)        (11111)(1)
       (4)(2)  (22)(2)    (21)(21)      (211)(11)        (11)(11)(11)
       (5)(1)  (31)(2)    (22)(11)      (2111)(1)        (111)(11)(1)
               (3)(21)    (221)(1)      (11)(11)(2)      (1111)(1)(1)
               (32)(1)    (3)(111)      (111)(2)(1)      (11)(11)(1)(1)
               (4)(11)    (31)(11)      (11)(2)(11)      (111)(1)(1)(1)
               (41)(1)    (311)(1)      (2)(11)(11)      (11)(1)(1)(1)(1)
               (2)(2)(2)  (11)(2)(2)    (21)(11)(1)      (1)(1)(1)(1)(1)(1)
               (3)(2)(1)  (2)(11)(2)    (211)(1)(1)
               (4)(1)(1)  (21)(2)(1)    (11)(2)(1)(1)
                          (2)(2)(11)    (2)(11)(1)(1)
                          (22)(1)(1)    (21)(1)(1)(1)
                          (3)(11)(1)    (2)(1)(1)(1)(1)
                          (31)(1)(1)
                          (2)(2)(1)(1)
                          (3)(1)(1)(1)

Links

Crossrefs

Row sums are A063834. Last column is A000041.
Cf. A000219, A001970, A007716, A008284, A055884, A289501, A317449, A317532, A317533, A320801, A320808.

Programs

  • Maple
    g:= proc(n, i) option remember; `if`(n=0 or i=1, x^n,
      g(n, i-1)+ `if`(i>n, 0, expand(g(n-i, i)*x)))
    end:
b:= proc(n, i) option remember; `if`(n=0 or i=1, x^n,
      b(n, i-1)+ `if`(i>n, 0, expand(b(n-i, i)*g(i$2))))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n$2)):
seq(T(n), n=0..12);  # Alois P. Heinz, Nov 11 2018
  • Mathematica
    Table[Length[Join@@Table[Select[Tuples[IntegerPartitions/@ptn],Length[Join@@#]==k&],{ptn,IntegerPartitions[n]}]],{n,0,10},{k,0,n}]
(* Second program: *)
g[n_, i_] := g[n, i] = If[n == 0 || i == 1, x^n,
     g[n, i - 1] + If[i > n, 0, Expand[g[n - i, i]*x]]];
b[n_, i_] := b[n, i] = If[n == 0 || i == 1, x^n,
     b[n, i - 1] + If[i > n, 0, Expand[b[n - i, i]*g[i, i]]]];
T[n_] := CoefficientList[b[n, n], x];
T /@ Range[0, 12] // Flatten (* Jean-François Alcover, May 20 2021, after Alois P. Heinz *)

Formula

O.g.f.: Product_{n >= 0} 1/(1 - x^n * (Sum_{0 <= k <= n} A008284(n,k) * t^k)).

A382214 Number of normal multisets of size n that can be partitioned into a set of sets.

Original entry on oeis.org

1, 1, 1, 3, 5, 11, 23, 48, 101, 210, 436, 894
Offset: 0

Views

Author

Gus Wiseman, Mar 29 2025

Keywords

Comments

First differs from A382216 at a(9) = 210, A382216(9) = 208.
We call a multiset or multiset partition normal iff it covers an initial interval of positive integers. The size of a multiset is the number of elements, counting multiplicity.

Examples

			The normal multiset {1,1,1,1,2,2,3,3,3} has partition {{1},{3},{1,2},{1,3},{1,2,3}}, so is counted under a(9).
The a(1) = 1 through a(5) = 11 multisets:
  {1}  {1,2}  {1,1,2}  {1,1,2,2}  {1,1,1,2,3}
              {1,2,2}  {1,1,2,3}  {1,1,2,2,3}
              {1,2,3}  {1,2,2,3}  {1,1,2,3,3}
                       {1,2,3,3}  {1,1,2,3,4}
                       {1,2,3,4}  {1,2,2,2,3}
                                  {1,2,2,3,3}
                                  {1,2,2,3,4}
                                  {1,2,3,3,3}
                                  {1,2,3,3,4}
                                  {1,2,3,4,4}
                                  {1,2,3,4,5}

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 A292432.
Twice-partitions of this type are counted by A358914, distinct sums A279785.
The strong version is A381996, complement A292444.
For integer partitions we have A382077, ranks A382200, complement A382078, ranks A293243.
For distinct sums we have A382216, complement A382202.
The case of a unique choice is counted by A382458, distinct sums A382459.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count factorizations, strict A045778.
Normal multiset partitions: A034691, A035310, A255906.
Set systems: A050342, A296120, A318361.
Set multipartitions: A089259, A270995, A296119, A318360.
Cf. A000110, A007716, A050320, A116540, A255903, A317532, A326517, A326519, A381990, A381992, A382428, A382460.

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[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]& /@ sps[Range[Length[mset]]]];
Table[Length[Select[allnorm[n],Select[mps[#], UnsameQ@@#&&And@@UnsameQ@@@#&]!={}&]],{n,0,5}]
Showing 1-10 of 22 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.