A339102 - OEIS

%I #6 Nov 23 2020 21:01:33

%S 1,0,0,0,1,1,1,1,1,3,3,5,5,7,7,15,15,23,29,37,43,57,87,101,137,175,

%T 235,279,363,431,665,757,1015,1257,1683,2069,2645,3199,4063,5607,6759,

%U 8591,10877,13837,17251,22185,26871,33773,41273,56047,66499,85647,104811

%N Number of compositions (ordered partitions) of n into distinct parts >= 4.

%H <a href="/index/Com#comp">Index entries for sequences related to compositions</a>

%F G.f.: Sum_{k>=0} k! * x^(k*(k + 7)/2) / Product_{j=1..k} (1 - x^j).

%e a(9) = 3 because we have [9], [5, 4] and [4, 5].

%p b:= proc(n, i, p) option remember;

%p       `if`(n=0, p!, `if`((i-3)*(i+4)/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 = 52; CoefficientList[Series[Sum[k! x^(k (k + 7)/2)/Product[1 - x^j, {j, 1, k}], {k, 0, nmax}], {x, 0, nmax}], x]

%Y Cf. A008484, A017898, A025149, A032020, A032022, A339101, A339103, A339104.

%K nonn

%O 0,10

%A _Ilya Gutkovskiy_, Nov 23 2020