OFFSET
1,1
COMMENTS
LINKS
Colin Barker, Table of n, a(n) for n = 1..500
Index entries for linear recurrences with constant coefficients, signature (99,-99,1).
FORMULA
O.g.f.: x*(7 + 106*x - x^2)/((1 - x)*(1 - 98*x + x^2)).
E.g.f.: (exp((5-2*sqrt(6))^2*x) + exp((5+2*sqrt(6))^2*x) - 14*exp(x))/12 + 1.
a(n) = 99*a(n-1) - 99*a(n-2) + a(n-3) for n>3.
a(n) = 98*a(n-1) - a(n-2) + 112 for n>2.
a(n) = a(-n) = ((5 + 2*sqrt(6))^(2*n) + (5 - 2*sqrt(6))^(2*n))/12 - 7/6.
a(n) = A001079(2*n)/6 - 7/6.
Lim_{n -> infinity} a(n)/a(n-1) = (5 + 2*sqrt(6))^2.
MAPLE
P:=proc(q) local n; for n from 1 to q do if type(sqrt((3*n^2+7*n+4)/2), integer) then print(n); fi; od; end: P(10^9); # Paolo P. Lava, Nov 25 2016
MATHEMATICA
Table[((5 + 2 Sqrt[6])^(2 n) + (5 - 2 Sqrt[6])^(2 n))/12 - 7/6, {n, 1, 20}]
RecurrenceTable[{a[1] == 7, a[2] == 799, a[n] == 98 a[n - 1] - a[n - 2] + 112}, a, {n, 1, 20}]
LinearRecurrence[{99, -99, 1}, {7, 799, 78407}, 20] (* Harvey P. Dale, Oct 18 2024 *)
PROG
(Sage)
def A278438():
a, b = 7, 799
yield a
while True:
yield b
a, b = b, 98*b - a + 112
a = A278438(); print([next(a) for _ in range(15)]) # Peter Luschny, Nov 24 2016
(Magma) Iv:=[7, 799]; [n le 2 select Iv[n] else 98*Self(n-1)-Self(n-2)+112: n in [1..20]];
(PARI) Vec(x*(7 + 106*x - x^2)/((1 - x)*(1 - 98*x + x^2)) + O(x^20)) \\ Colin Barker, Nov 27 2016
CROSSREFS
KEYWORD
nonn,easy,changed
AUTHOR
Bruno Berselli, Nov 23 2016
STATUS
approved