%I #37 Jan 13 2022 01:31:58
%S 1,1,2,1,2,3,1,4,3,5,1,4,7,5,7,1,6,10,13,7,11,1,6,14,20,19,11,15,1,8,
%T 18,33,32,31,15,22,1,8,25,43,56,54,43,22,30,1,10,29,66,81,99,78,64,30,
%U 42,1,10,37,83,126,150,148,118,88,42,56,1,12,44,114,174,246,235,230,166,124,56,77
%N Triangle read by rows: T(n,k) is the number of plane partitions of n containing exactly k parts.
%C First column is 1, representing the single-part {{n}}, last column is P(n), since the all-ones plane partitions form the Ferrers-Young plots of the (linear) partitions of n.
%C A plane partition of n is a two-dimensional table (or matrix) with nonnegative elements summing up to n, and nonincreasing rows and columns. (Zero rows and columns are ignored.) - _M. F. Hasler_, Sep 22 2018
%H Alois P. Heinz, <a href="/A091298/b091298.txt">Rows n = 1..50</a>
%H A. Rovenchak, <a href="http://arxiv.org/abs/1401.4367">Enumeration of plane partitions with a restricted number of parts</a>, arXiv preprint arXiv:1401.4367 [math-ph], 2014.
%H E. W. Weisstein, <a href="http://mathworld.wolfram.com/PlanePartition.html">Plane partition</a>.
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Plane_partition">Plane partition</a>.
%e This plane partition of n=7: {{3,1,1},{2}} contains 4 parts: 3,1,1,2.
%e Triangle T(n,k) begins:
%e 1;
%e 1, 2;
%e 1, 2, 3;
%e 1, 4, 3, 5;
%e 1, 4, 7, 5, 7;
%e 1, 6, 10, 13, 7, 11;
%e 1, 6, 14, 20, 19, 11, 15;
%e 1, 8, 18, 33, 32, 31, 15, 22;
%e 1, 8, 25, 43, 56, 54, 43, 22, 30;
%e 1, 10, 29, 66, 81, 99, 78, 64, 30, 42;
%e ...
%t (* see A089924 for "planepartition" *) Table[Length /@ Split[Sort[Length /@ Flatten /@ planepartitions[n]]], {n, 16}]
%o (PARI) A091298(n,k)=sum(i=1,#n=PlanePartitions(n),sum(j=1,#n[i],#n[i][j])==k)
%o PlanePartitions(n,L=0,PP=List())={ n<2&&return([if(n,[[1]],[])]); for(N=1,n, my(P=apply(Vecrev, if(L, select(p->vecmin(L-Vecrev(p,#L))>=0, partitions(N,L[1],#L)), partitions(N)))); if(N<n, for(i=1, #P, my(pp = PlanePartitions(n-N,P[i])); for(j=1, #pp, listput(PP, concat([P[i]], pp[j])))), for(i=1,#P, listput(PP, [P[i]])))); Set(PP)} \\ _M. F. Hasler_, Sep 24 2018
%Y Row sums give A000219.
%Y Cf. A090539, A092288, A089924.
%Y Column 1 is A000012. Column 2 is A052928. Diagonal and subdiagonal are A000041.
%K nonn,tabl
%O 1,3
%A _Wouter Meeussen_, Feb 24 2004