|
%I #6 Nov 23 2020 21:14:01
%S 1,0,0,0,0,1,1,1,1,1,1,3,3,5,5,7,7,9,15,17,23,31,37,45,57,65,101,115,
%T 151,189,255,293,383,451,565,777,921,1157,1469,1855,2311,2865,3495,
%U 4313,5231,7063,8269,10509,12849,16217,19829,25171,30031,37485,45183
%N Number of compositions (ordered partitions) of n into distinct parts >= 5.
%H <a href="/index/Com#comp">Index entries for sequences related to compositions</a>
%F G.f.: Sum_{k>=0} k! * x^(k*(k + 9)/2) / Product_{j=1..k} (1 - x^j).
%e a(11) = 3 because we have [11], [6, 5] and [5, 6].
%p b:= proc(n, i, p) option remember;
%p `if`(n=0, p!, `if`((i-4)*(i+5)/2<n, 0,
%p add(b(n-i*j, i-1, p+j), j=0..min(1, n/i))))
%p end:
%p a:= n-> b(n$2, 0):
%p seq(a(n), n=0..60); # _Alois P. Heinz_, Nov 23 2020
%t nmax = 54; CoefficientList[Series[Sum[k! x^(k (k + 9)/2)/Product[1 - x^j, {j, 1, k}], {k, 0, nmax}], {x, 0, nmax}], x]
%Y Cf. A017899, A025150, A032020, A032022, A185325, A339101, A339102, A339104.
%K nonn
%O 0,12
%A _Ilya Gutkovskiy_, Nov 23 2020
|
