OFFSET
1,6
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
MAPLE
b:= proc(n) option remember; true end: l:= [0, 1]: for k to 100 do b(l[1]):= false; l:= [l[2], l[1]+l[2]] od: aa:= proc(n, i) option remember; local g, h; if n=0 then [1, 0] elif i=0 or n<0 then [0, 0] else g:= aa(n, i-1); h:= aa(n-i, i); [g[1]+h[1], g[2]+h[2] +`if`(b(i), h[1], 0)] fi end: a:= n-> aa(n, n)[2]: seq(a(n), n=1..60); # Alois P. Heinz, Jul 28 2009
MATHEMATICA
Clear[b]; b[_] = True; l = {0, 1}; For[k=1, k <= 100, k++, b[l[[1]]] = False; l = {l[[2]], l[[1]] + l[[2]]}]; aa[n_, i_] := aa[n, i] = Module[{g, h}, If[n==0, {1, 0}, If[i==0 || n<0, {0, 0}, g = aa[n, i-1]; h = aa[n-i, i]; {g[[1]] + h[[1]], g[[2]] + h[[2]] + If[b[i], h[[1]], 0]}] ]]; a[n_] := aa[n, n][[2]]; Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Jul 30 2015, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Omar E. Pol, Sep 11 2008
EXTENSIONS
More terms from Alois P. Heinz, Jul 28 2009
STATUS
approved