OFFSET
1,1
LINKS
Index entries for linear recurrences with constant coefficients, signature (1,0,34,-34,0,-1,1).
FORMULA
a(n) = 34*a(n-3)-a(n-6)-2312 for n > 6; a(1)=169, a(2)=289, a(3)=625, a(4)=2809, a(5)=7225, a(6)=18769.
G.f.: x*(169+120*x+336*x^2-3562*x^3+336*x^4+120*x^5+169*x^6)/((1-x)*(1-34*x^3+x^6)).
Limit_{n -> oo} a(n)/a(n-3) = (17+12*sqrt(2)).
Limit_{n -> oo} a(n)/a(n-1) = ((19+6*sqrt(2))/17)^2 for n mod 3 = {0, 2}.
Limit_{n -> oo} a(n)/a(n-1) = ((387+182*sqrt(2))/17^2)^2 for n mod 3 = 1.
EXAMPLE
625 = 25^2 is of the form k^2+(k+17)^2 with k = 7: 7^2+24^2 = 625. Hence 625 is in the sequence.
MATHEMATICA
LinearRecurrence[{1, 0, 34, -34, 0, -1, 1}, {169, 289, 625, 2809, 7225, 18769, 93025}, 30] (* Harvey P. Dale, Apr 22 2022 *)
PROG
(PARI) {forstep(n=-5, 1600000, [1, 3], if(issquare(a=2*n*(n+17)+289), print1(a, ", ")))}
CROSSREFS
Equals A155923^2. Cf. A156160 (first trisection), A156161 (second trisection), A156162 (third trisection).
KEYWORD
nonn,easy
AUTHOR
Klaus Brockhaus, Feb 09 2009
EXTENSIONS
G.f. corrected, fourth comment and cross-references edited by Klaus Brockhaus, Sep 23 2009
STATUS
approved