OFFSET
0,3
COMMENTS
For a given positive integer, n, let S_n be the set of partitions of n into distinct parts where the number of parts is maximal for that n. For example, for n=6, the set S_6 consists of just one such partition S_6={1,2,3}. Similarly, for n=7, S_7={1,2,4}, But for n=8, S_8 will contain two partitions S_8= { {1,2,5}, {1,3,4} }. Now form the sum whose general term is x^n divided by the product (1-x^(p_1))...(1-x^(p_i)) where the p's come from the partitions in S_n. The sequence is the sequence of coefficients of this sum.
This sequence is an upper bound for A098859.
The initial values of the sequence were verified by several members of the Sequence Fans mailing list.
|S_n| = A144328(n) (if the offset of the latter is changed to 0). - Benoit Jubin, Dec 13 2010.
MATHEMATICA
<<Combinatorica`
f[set_]:=x^Total[set]/Product[1-x^set[[i]], {i, 1, Length[set]}]
gf=1;
For[n=1, n<=20, n++, Print["n= ", n];
For[i=1, i<=Length[Partitions[n]], i++,
t=Tally[Partitions[n][[i]]];
sum=0;
For[j=1, j<=Length[t], j++,
If[t[[j]][[2]]==1, sum++]];
If[sum==Length[t]
&&Length[t]==Floor[(-1+Sqrt[1+8 n])/2],
gf=gf+f[Partitions[n][[i]]]; Print[CoefficientList[Series[gf, {x, 0, n}], x]];
]]]
CROSSREFS
KEYWORD
nonn
AUTHOR
David S. Newman, Dec 27 2010
EXTENSIONS
More terms from Alois P. Heinz, Jun 25 2011
STATUS
approved