OFFSET
0,3
COMMENTS
After the first term, row sums of triangle A111669. - Emanuele Munarini, Dec 05 2017
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..140
FORMULA
G.f.: (1 - G(0) )/(1-x) where G(k) = 1 - 1/(1-Fibonacci(k)*x)/(1-x/(x-1/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 17 2013
G.f.: 1/(x*(1-x)*G(0)) - 1/x where G(k) = 1 - x/(x - 1/(1 + 1/(x*Fibonacci(k)-1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Feb 13 2013
EXAMPLE
A(x) = 1 + x/(1-x) + x^2/((1-x)*(1-x)) + x^3/((1-x)*(1-x)*(1-2*x)) +
x^4/((1-x)*(1-x)*(1-2*x)(1-3*x)) + x^5/((1-x)*(1-x)*(1-2*x)*(1-3*x)*(1-5*x)) + x^6/((1-x)*(1-x)*(1-2*x)*(1-3*x)*(1-5*x)*(1-8*x)) +...
MAPLE
b:= proc(n, m) option remember; `if`(n=0, 1,
(<<0|1>, <1|1>>^m)[1, 2]*b(n-1, m)+b(n-1, m+1))
end:
a:= n-> b(n, 0):
seq(a(n), n=0..25); # Alois P. Heinz, Aug 08 2021
MATHEMATICA
b[n_, m_] := b[n, m] = If[n == 0, 1,
MatrixPower[{{0, 1}, {1, 1}}, m][[1, 2]]*b[n-1, m]+b[n-1, m+1]];
a[n_] := b[n, 0];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Sep 07 2022, after Alois P. Heinz *)
PROG
(PARI) {a(n)=polcoeff(sum(k=0, n, x^k/prod(j=0, k, 1-fibonacci(j)*x+x*O(x^n))), n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 07 2007
STATUS
approved