|
%I #13 May 21 2018 11:27:58
%S 1,1,2,3,5,7,11,15,22,30,42,56,77,100,133,171,223,282,362,453,572,709,
%T 884,1084,1337,1626,1984,2394,2896,3468,4162,4951,5897,6972,8249,9696,
%U 11402,13330,15586,18131,21090,24417,28264,32580,37541,43097,49449,56544
%N Number of partitions of n in which any two parts differ by at most 10.
%H Alois P. Heinz, <a href="/A218512/b218512.txt">Table of n, a(n) for n = 0..1000</a>
%F G.f.: 1 + Sum_{j>0} x^j / Product_{i=0..10} (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, 10), 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, 10], {i, 1, n}];
%t Table[a[n], {n, 0, 80}] (* _Jean-François Alcover_, May 21 2018, translated from Maple *)
%Y Column k=10 of A194621.
%K nonn
%O 0,3
%A _Alois P. Heinz_, Oct 31 2012
|
