login
Number of partitions of n in which any two parts differ by at most 8.
4

%I #13 May 20 2018 11:34:52

%S 1,1,2,3,5,7,11,15,22,30,42,55,75,96,127,161,208,260,330,407,509,621,

%T 765,925,1127,1350,1627,1934,2310,2725,3227,3782,4446,5178,6044,7000,

%U 8122,9355,10791,12370,14195,16196,18494,21012,23887,27029,30596,34492,38894

%N Number of partitions of n in which any two parts differ by at most 8.

%H Alois P. Heinz, <a href="/A218510/b218510.txt">Table of n, a(n) for n = 0..1000</a>

%F G.f.: 1 + Sum_{j>0} x^j / Product_{i=0..8} (1-x^(i+j)).

%p b:= proc(n, i, k) option remember; `if`(n<0 or k<0, 0,

%p `if`(n=0, 1, `if`(i<1, 0, b(n, i-1, k-1) +b(n-i, i, k))))

%p end:

%p a:= n-> `if`(n=0, 1, 0) +add(b(n-i, i, 8), i=1..n):

%p seq(a(n), n=0..80);

%t b[n_, i_, k_] := b[n, i, k] = If[n < 0 || k < 0, 0, If[n == 0, 1, If[i < 1, 0, b[n, i - 1, k - 1] + b[n - i, i, k]]]];

%t a[n_] := If[n == 0, 1, 0] + Sum[b[n - i, i, 8], {i, 1, n}];

%t Table[a[n], {n, 0, 80}] (* _Jean-François Alcover_, May 20 2018, after _Alois P. Heinz_ *)

%Y Column k=8 of A194621.

%K nonn

%O 0,3

%A _Alois P. Heinz_, Oct 31 2012