OFFSET
0,3
COMMENTS
a(n*(n+1)/2) = a(A000217(n)) = 1.
LINKS
Joerg Arndt and Alois P. Heinz, Table of n, a(n) for n = 0..140
FORMULA
a(n) = A238858(n,Re(n-floor((sqrt(8*n-7)+1)/2))).
MAPLE
b:= proc(n, i, t) option remember; `if`(n=0, 1, expand(add(
`if`(j<i, x, 1) *b(n-1, j, t+`if`(j>i, 1, 0)), j=0..t+1)))
end:
a:= n-> (p-> coeff(p, x, degree(p)))(b(n, -1$2)):
seq(a(n), n=0..40);
MATHEMATICA
b[n_, i_, t_] := b[n, i, t] = If[n == 0, 1, Expand[Sum[If[j<i, x, 1] *b[n-1, j, t + If[j>i, 1, 0]], {j, 0, t+1}]]]; a[n_] := Function[{p}, Coefficient[p, x, Exponent[ p, x ]]][b[n, -1, -1]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Feb 13 2015, after Maple *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Joerg Arndt and Alois P. Heinz, May 01 2014
STATUS
approved