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.

A178702 Coefficients of the sum 1+ x/(1-x) + x^2/(1-x^2) + x^3/ ( (1-x)(1-x^2)) + x^4/ ( (1-x)(1-x^3) ) + x^5/ ( (1-x)(1-x^4) ) + x^5 /((1-x^2)(1-x^3)) + x^6/ ( (1-x)(1-x^2)(1-x^3)) + ...

This page as a plain text file.
%I A178702 #25 Sep 20 2014 02:46:37
%S A178702 1,1,2,2,4,6,7,10,14,20,24,32,40,54,69,86,106,135,165,206,256,311,378,
%T A178702 460,555,670,808,970,1156,1380,1638,1938,2296,2706,3188,3752,4390,
%U A178702 5136,6000,6990,8128,9444,10944,12672,14659,16904,19476,22420,25753,29550,33873,38759,44306
%N A178702 Coefficients of the sum 1+ x/(1-x) + x^2/(1-x^2) + x^3/ ( (1-x)(1-x^2)) + x^4/ ( (1-x)(1-x^3) ) +  x^5/  ( (1-x)(1-x^4) ) + x^5 /((1-x^2)(1-x^3)) + x^6/ ( (1-x)(1-x^2)(1-x^3)) + ...
%C A178702 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.
%C A178702 This sequence is an upper bound for A098859.
%C A178702 The initial values of the sequence were verified by several members of the Sequence Fans mailing list.
%C A178702  |S_n| = A144328(n) (if the offset of the latter is changed to 0). - _Benoit Jubin_, Dec 13 2010.
%t A178702 <<Combinatorica`
%t A178702 f[set_]:=x^Total[set]/Product[1-x^set[[i]],{i,1,Length[set]}]
%t A178702 gf=1;
%t A178702 For[n=1,n<=20,n++,Print["n= ",n];
%t A178702 For[i=1,i<=Length[Partitions[n]],i++,
%t A178702 t=Tally[Partitions[n][[i]]];
%t A178702 sum=0;
%t A178702 For[j=1,j<=Length[t],j++,
%t A178702 If[t[[j]][[2]]==1,sum++]];
%t A178702 If[sum==Length[t]
%t A178702 &&Length[t]==Floor[(-1+Sqrt[1+8 n])/2],
%t A178702 gf=gf+f[Partitions[n][[i]]];Print[CoefficientList[Series[gf,{x,0,n}],x]];
%t A178702 ]]]
%Y A178702 Cf. A098859, A144328.
%K A178702 nonn
%O A178702 0,3
%A A178702 _David S. Newman_, Dec 27 2010
%E A178702 More terms from _Alois P. Heinz_, Jun 25 2011