OFFSET
0,1
COMMENTS
Continued fraction for sqrt(a(n)) has period n.
LINKS
Index entries for linear recurrences with constant coefficients, signature (22, -71, -36, 73, 14, -1).
FORMULA
G.f.: (-x^5 - 91*x^4 + 534*x^3 + 370*x^2 - 181*x + 9)/((x^2-18*x+1)*(x^2-1)*(x^2+4*x-1)). - Michael Somos
a(n) = 22*a(n-1) - 71*a(n-2) - 36*a(n-3) + 73*a(n-4) + 14*a(n-5) - a(n-6).
MATHEMATICA
#[[1]]^2-2#[[1]]+4#[[2]]+5&/@Table[{Fibonacci[3n], Fibonacci[3n+1]}, {n, 0, 20}] (* or *) LinearRecurrence[{22, -71, -36, 73, 14, -1}, {9, 17, 105, 1313, 21385, 374833}, 20] (* Harvey P. Dale, Jul 07 2017 *)
PROG
(PARI) a(n) = fibonacci(3*n)^2 - 2*fibonacci(3*n) + 4*fibonacci(3*n+1) + 5 \\ Andrew Howroyd, Aug 10 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved