OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..102
FORMULA
MAPLE
b:= proc(n) option remember; local i; if n=0 then 1 else 0;
for i while i*(i+1)/2<=n do %+b(n-i*(i+1)/2) od; % fi
end:
a:= n-> b(n*(n+1)/2):
seq(a(n), n=0..20); # Alois P. Heinz, Feb 05 2018
MATHEMATICA
b[n_] := b[n] = If[n==0, 1, Sum[If[IntegerQ[Sqrt[8j+1]], b[n-j], 0], {j, 1, n}]];
a[n_] := b[n(n+1)/2];
PROG
(PARI) {a(n)=polcoeff(1/(1-sum(r=1, n+1, x^(r*(r+1)/2)+x*O(x^(n*(n+1)/2)))), n*(n+1)/2)}
for(n=0, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Apr 14 2013
STATUS
approved