OFFSET
1,3
COMMENTS
Positive terms are of the form (m^2-9)/16 where m runs over the elements of A077443 that are congruent to 5 modulo 8. Correspondingly, for n>1, sqrt(16*a(n)+9) form a subsequence of A077443, while sqrt(8*a(n)+1) form a subsequence of A077442 with indices congruent to 2,3 modulo 4. [Max Alekseyev]
LINKS
Index entries for linear recurrences with constant coefficients, signature (1, 1154, -1154, -1, 1).
FORMULA
For n>5, a(n) = 1154*a(n-2) - a(n-4) + 396.
For n>6, a(n) = a(n-1) + 1154*a(n-2) - 1154*a(n-3) - a(n-4) + a(n-5).
For n>1, a(n) = 1/64 * ( (9 + 4* sqrt(2)*(-1)^n)*(1+sqrt(2))^(4*n-6) + (9 - 4* sqrt(2)*(-1)^n)*(1-sqrt(2))^(4*n-6) - 22).
a(n) = floor ( 1/64 * (9 + 4*sqrt(2)*(-1)^n) * (1+sqrt(2))^(4*n-6) ).
G.f.: (x^5 + 9*x^4 + 376*x^3 + 9*x^2 + x)/((1 - x)*(x^2 - 34*x + 1)*(x^2 + 34*x + 1)). [corrected by Peter Luschny, Apr 04 2019]
Lim (n -> Infinity, a(2n+1)/a(2n)) = (1/49)*(3649+2580*sqrt(2)).
Lim (n -> Infinity, a(2n)/a(2n-1)) = (1/49)*(193+132*sqrt(2)).
EXAMPLE
MATHEMATICA
LinearRecurrence[{1, 1154, -1154, -1, 1} , {0, 1, 10, 1540, 11935, 1777555}, 17] (* first term 0 corrected by Georg Fischer, Apr 02 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Richard Choulet, Oct 11 2007; Ant King, Nov 04 2011
EXTENSIONS
Entry revised by N. J. A. Sloane, Nov 06 2011
Term 0 prepended and entry revised accordingly by Max Alekseyev, Nov 06 2011
STATUS
approved