OFFSET
0,3
LINKS
Colin Barker, Table of n, a(n) for n = 0..500
Index entries for linear recurrences with constant coefficients, signature (99,-99,1).
FORMULA
O.g.f.: x/((1 - x)*(1 - 98*x + x^2)).
E.g.f.: ((5-2*sqrt(6))*exp((5-2*sqrt(6))^2*x) + (5+2*sqrt(6))*exp((5+2*sqrt(6))^2*x) - 10*exp(x))/960.
a(n) = 99*a(n-1) - 99*a(n-2) + a(n-3) for n>2.
a(n) = 98*a(n-1) - a(n-2) + 1 for n>1.
a(n) = a(-n-1) = ((5+2*sqrt(6))^(2*n+1) + (5-2*sqrt(6))^(2*n+1))/960 - 1/96.
a(n) = floor((5+2*sqrt(6))^(2*n+1)/960).
a(n)*a(n-2) = a(n-1)*(a(n-1)-1) for n>1.
Lim_{i -> infinity} a(i)/a(i-1) = (5 + 2*sqrt(6))^2.
From the closed form: a(n) + a(-n) = A108741(n).
a(n) = A123479(n)/20 for n>0.
a(n) = A045502(n)/40.
MAPLE
P:=proc(q) local a, b, c, n; a:=0; b:=1; print(a); print(b); for n from 1 to q do
c:=98*b-a+1; a:=b; b:=c; print(b); od; end: P(100); # Paolo P. Lava, Nov 30 2016
MATHEMATICA
CoefficientList[x/(1 - 99 x + 99 x^2 - x^3) + O[x]^20, x]
LinearRecurrence[{99, -99, 1}, {0, 1, 99}, 20] (* Harvey P. Dale, Aug 22 2020 *)
PROG
(PARI) concat(0, Vec(1/(1-99*x+99*x^2-x^3) + O(x^20)))
(Sage) gf = x/((1-x)*(1-98*x+x^2)); print(taylor(gf, x, 0, 20).list())
(Maxima) makelist(coeff(taylor(x/((1-x)*(1-98*x+x^2)), x, 0, n), x, n), n, 0, 20);
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Nov 24 2016
STATUS
approved