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-4 of 4 results.

A117433 Number of planar partitions of n with all part sizes distinct.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 11, 15, 21, 35, 41, 59, 75, 103, 149, 187, 243, 321, 413, 527, 735, 895, 1165, 1467, 1885, 2335, 2997, 3853, 4765, 5977, 7473, 9269, 11531, 14255, 17537, 22201, 26897, 33233, 40613, 50027, 60637, 74459, 89963, 109751, 134407, 162117, 195859
Offset: 0

Views

Author

Franklin T. Adams-Watters, Mar 16 2006, Apr 01 2008

Keywords

Comments

Matches A072706 for n < 10, since a unimodal composition into distinct parts can be placed uniquely as a hook. Starting with n = 10, additional partitions are possible (starting with [4,3|2,1] and [4,2|3,1]).

Examples

			From _Gus Wiseman_, Nov 15 2018: (Start)
The a(10) = 35 strict plane partitions (A = 10):
  A  64  73  82  532  91  541  631  721  4321
.
  9  54  63  72  432  8  53  71  431  7  43  52  61  421  6  42  51
  1  1   1   1   1    2  2   2   2    3  21  3   3   3    4  31  4
.
  7  6  5  43  42  5  41
  2  3  4  2   3   3  3
  1  1  1  1   1   2  2
.
  4
  3
  2
  1
(End)
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) b(n, i):= `if`(n=0, [1], `if`(i<1, [], zip((x, y)
          -> x+y, b(n, i-1), `if`(i>n, [], [0, b(n-i, i-1)[]]), 0)))
        end:
    g:= proc(n) g(n):= `if`(n<2, 1, (n-1)*g(n-2) +g(n-1)) end:
    a:= proc(n) b(n, n); add(%[i]*g(i-1), i=1..nops(%)) end:
    seq(a(n), n=0..60);  # Alois P. Heinz, Nov 18 2012
  • Mathematica
    prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    Table[Length[Select[multsubs[Tuples[Range[n],2],n],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#],UnsameQ@@DeleteCases[Join@@prs2mat[#],0],And@@(OrderedQ[#,Greater]&/@prs2mat[#]),And@@(OrderedQ[#,Greater]&/@Transpose[prs2mat[#]])]&]],{n,5}] (* Gus Wiseman, Nov 15 2018 *)
    zip[f_, x_List, y_List, z_] := With[{m = Max[Length[x], Length[y]]}, f[PadRight[x, m, z], PadRight[y, m, z]]];
    b[n_, i_] := b[n, i] = If[n == 0, {1}, If[i < 1, {}, zip[Plus, b[n, i - 1], If[i > n, {}, Join[{0}, b[n - i, i - 1]]], 0]]];
    g[n_] := g[n] = If[n < 2, 1, (n - 1)*g[n - 2] + g[n - 1]];
    a[n_] := With[{bn = b[n, n]}, Sum[bn[[i]]*g[i - 1], {i, 1, Length[bn]}]];
    Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Dec 05 2023, after Alois P. Heinz *)

Formula

a(n) = Sum_{k=1..floor((sqrt(8*n+1)-1)/2)} A000085(k)*A008289(n,k).

A321662 Number of non-isomorphic multiset partitions of weight n whose incidence matrix has all distinct entries.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 13, 15, 23, 33, 49, 59, 83, 101, 133, 281, 321, 477, 655, 941, 1249, 1795, 2241, 3039, 3867, 5047, 6257, 8063, 11459, 13891, 18165, 23149, 29975, 37885, 49197, 61829, 89877, 109165, 145673, 185671, 246131, 310325, 408799, 514485, 668017, 871383
Offset: 0

Views

Author

Gus Wiseman, Nov 15 2018

Keywords

Comments

The incidence matrix of a multiset partition has entry (i, j) equal to the multiplicity of vertex i in part j.
Also the number of positive integer matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, with all different entries.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(3) = 3 through a(7) = 15 multiset partitions:
  {{111}}    {{1111}}    {{11111}}    {{111111}}      {{1111111}}
  {{122}}    {{1222}}    {{11222}}    {{112222}}      {{1112222}}
  {{1}{11}}  {{1}{111}}  {{12222}}    {{122222}}      {{1122222}}
                         {{1}{1111}}  {{122333}}      {{1222222}}
                         {{11}{111}}  {{1}{11111}}    {{1223333}}
                                      {{11}{1111}}    {{1}{111111}}
                                      {{1}{11222}}    {{11}{11111}}
                                      {{11}{1222}}    {{111}{1111}}
                                      {{112}{222}}    {{1}{112222}}
                                      {{122}{222}}    {{11}{12222}}
                                      {{2}{11222}}    {{112}{2222}}
                                      {{22}{1222}}    {{122}{2222}}
                                      {{1}{11}{111}}  {{2}{112222}}
                                                      {{22}{12222}}
                                                      {{1}{11}{1111}}
		

Crossrefs

Programs

  • Mathematica
    (* b = A121860 *) b[n_] := Sum[n!/(d! (n/d)!), {d, Divisors[n]}];
    (* c = A008289 *) c[n_, k_] := c[n, k] = If[n < k || k < 1, 0, If[n == 1, 1, c[n - k, k] + c[n - k, k - 1]]];
    a[n_] := If[n == 0, 1, Sum[ (b[k] + b[k + 1] - 2) c[n, k], {k, 1, n}]];
    a /@ Range[0, 45] (* Jean-François Alcover, Sep 14 2019 *)
  • PARI
    \\ here b(n) is A121860(n).
    b(n)={sumdiv(n, d, n!/(d!*(n/d)!))}
    seq(n)={my(B=vector((sqrtint(8*(n+1))+1)\2, n, if(n==1, 1, b(n-1)+b(n)-2))); apply(p->sum(i=0, poldegree(p), B[i+1]*polcoef(p, i)), Vec(prod(k=1, n, 1 + x^k*y + O(x*x^n))))} \\ Andrew Howroyd, Nov 16 2018

Formula

a(n) = Sum_{k>=1} (A121860(k) + A121860(k+1) - 2)*A008289(n,k) for n > 0. - Andrew Howroyd, Nov 17 2018

Extensions

Terms a(11) and beyond from Andrew Howroyd, Nov 16 2018

A321660 Number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns, whose entries are all distinct.

Original entry on oeis.org

1, 1, 1, 5, 5, 9, 45, 49, 85, 125, 233, 273, 417, 529, 745, 2573, 2861, 4761, 6837, 10489, 14317, 22637, 28289, 40041, 52041, 70177, 88561, 117605, 234773, 274761, 407469, 553681, 792613, 1052525, 1493033, 1959009, 3135537, 3904129, 5475673, 7173725, 9853325
Offset: 0

Views

Author

Gus Wiseman, Nov 15 2018

Keywords

Examples

			The a(5) = 9 matrices:
  [5] [4 1] [3 2] [2 3] [1 4]
.
  [4] [3] [2] [1]
  [1] [2] [3] [4]
		

Crossrefs

Programs

  • Mathematica
    prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    Table[Length[Select[multsubs[Tuples[Range[n],2],n],And[Union[First/@#]==Range[Max@@First/@#],Union[Last/@#]==Range[Max@@Last/@#],UnsameQ@@Join@@prs2mat[#]]&]],{n,5}]
  • PARI
    seq(n)={my(B=vector((sqrtint(8*(n+1))+1)\2, n, if(n==1, 1, (n-1)!*numdiv(n-1) + n!*(numdiv(n) - 2)))); apply(p->sum(i=0, poldegree(p), B[i+1]*polcoef(p, i)), Vec(prod(k=1, n, 1 + x^k*y + O(x*x^n))))} \\ Andrew Howroyd, Nov 16 2018

Formula

a(n) = Sum_{k>=1} (k!*A000005(k) + (k+1)!*(A000005(k+1) - 2))*A008289(n,k) for n > 0. - Andrew Howroyd, Nov 17 2018

Extensions

Terms a(11) and beyond from Andrew Howroyd, Nov 16 2018

A321661 Number of non-isomorphic multiset partitions of weight n where the nonzero entries of the incidence matrix are all distinct.

Original entry on oeis.org

1, 1, 1, 4, 4, 7, 22, 25, 40, 58, 186, 204, 347, 478, 734, 2033, 2402, 3814, 5464, 8142, 11058, 30142, 34437, 55940, 77794, 116954, 156465, 229462, 533612, 640544, 994922, 1397896, 2048316, 2778750, 3987432, 5292293, 11921070, 14076550, 21802928, 29917842, 44080285
Offset: 0

Views

Author

Gus Wiseman, Nov 15 2018

Keywords

Comments

The incidence matrix of a multiset partition has entry (i, j) equal to the multiplicity of vertex i in part j.
Also the number of positive integer matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, whose nonzero entries are all distinct.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(6) = 22 multiset partitions:
  {{1}}  {{11}}  {{111}}    {{1111}}    {{11111}}    {{111111}}
                 {{122}}    {{1222}}    {{11222}}    {{112222}}
                 {{1}{11}}  {{1}{111}}  {{12222}}    {{122222}}
                 {{1}{22}}  {{1}{222}}  {{1}{1111}}  {{122333}}
                                        {{11}{111}}  {{1}{11111}}
                                        {{11}{222}}  {{11}{1111}}
                                        {{1}{2222}}  {{1}{11222}}
                                                     {{11}{1222}}
                                                     {{11}{2222}}
                                                     {{112}{222}}
                                                     {{11}{2333}}
                                                     {{1}{22222}}
                                                     {{122}{222}}
                                                     {{1}{22333}}
                                                     {{122}{333}}
                                                     {{2}{11222}}
                                                     {{22}{1222}}
                                                     {{1}{11}{111}}
                                                     {{1}{11}{222}}
                                                     {{1}{22}{222}}
                                                     {{1}{22}{333}}
                                                     {{2}{11}{222}}
		

Crossrefs

Programs

  • PARI
    \\ here b(n) is A059849(n).
    b(n)={sum(k=0, n, stirling(n,k,1)*sum(i=0, k, stirling(k,i,2))^2)}
    seq(n)={my(B=vector((sqrtint(8*(n+1))+1)\2, n, b(n-1))); apply(p->sum(i=0, poldegree(p), B[i+1]*polcoef(p,i)), Vec(prod(k=1, n, 1 + x^k*y + O(x*x^n))))} \\ Andrew Howroyd, Nov 16 2018

Formula

a(n) = Sum_{k>=1} A059849(k)*A008289(n,k) for n > 0. - Andrew Howroyd, Nov 17 2018

Extensions

Terms a(11) and beyond from Andrew Howroyd, Nov 16 2018
Showing 1-4 of 4 results.