This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A078988 #38 Feb 25 2023 03:07:46
%S A078988 1,65,4289,283009,18674305,1232221121,81307919681,5365090477825,
%T A078988 354014663616769,23359602708228929,1541379764079492545,
%U A078988 101707704826538279041,6711167138787446924161,442835323455144958715585,29220420180900779828304449,1928104896615996323709378049
%N A078988 Chebyshev sequence with Diophantine property.
%C A078988 Bisection (even part) of A041025.
%C A078988 (4*A078989(n))^2 - 17*a(n)^2 = -1 (Pell -1 equation, see A077232-3).
%C A078988 Starting with a(1), hypotenuses of primitive Pythagorean triples in A195619 and A195620. - _Clark Kimberling_, Sep 22 2011
%H A078988 Colin Barker, <a href="/A078988/b078988.txt">Table of n, a(n) for n = 0..549</a>
%H A078988 A. J. C. Cunningham, <a href="https://archive.org/details/binomialfactoris01cunn/page/n46/mode/1up">Binomial Factorisations</a>, Vols. 1-9, Hodgson, London, 1923-1929. See Vol. 1, page xxxv.
%H A078988 Tanya Khovanova, <a href="http://www.tanyakhovanova.com/RecursiveSequences/RecursiveSequences.html">Recursive Sequences</a>
%H A078988 Giovanni Lucca, <a href="http://forumgeom.fau.edu/FG2019volume19/FG201902index.html">Integer Sequences and Circle Chains Inside a Hyperbola</a>, Forum Geometricorum (2019) Vol. 19, 11-16.
%H A078988 <a href="/index/Ch#Cheby">Index entries for sequences related to Chebyshev polynomials.</a>
%H A078988 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (66,-1).
%F A078988 G.f.: (1-x)/(1-66*x+x^2).
%F A078988 a(n) = T(2*n+1, sqrt(17))/sqrt(17) = ((-1)^n)*S(2*n, 8*i) = S(n, 66) - S(n-1, 66) with i^2=-1 and T(n, x), resp. S(n, x), Chebyshev's polynomials of the first, resp. second, kind. See A053120 and A049310.
%F A078988 a(n) = A041025(2*n).
%F A078988 a(n) = 66*a(n-1) - a(n-2) for n>1 ; a(0)=1, a(1)=65. - _Philippe Deléham_, Nov 18 2008
%e A078988 (x,y) = (4,1), (268,65), (17684,4289), ... give the positive integer solutions to x^2 - 17*y^2 =-1.
%t A078988 CoefficientList[Series[(1-x)/(1-66x+x^2), {x,0,20}], x] (* _Michael De Vlieger_, Apr 15 2019 *)
%t A078988 LinearRecurrence[{66,-1}, {1,65}, 21] (* _G. C. Greubel_, Aug 01 2019 *)
%o A078988 (PARI) Vec((1-x)/(1-66*x+x^2) + O(x^20)) \\ _Colin Barker_, Jun 15 2015
%o A078988 (Magma) I:=[1, 65]; [n le 2 select I[n] else 66*Self(n-1) - Self(n-2): n in [1..20]]; // _G. C. Greubel_, Aug 01 2019
%o A078988 (Sage) ((1-x)/(1-66*x+x^2)).series(x, 20).coefficients(x, sparse=False) # _G. C. Greubel_, Aug 01 2019
%o A078988 (GAP) a:=[1,65];; for n in [3..20] do a[n]:=66*a[n-1]-a[n-2]; od; a; # _G. C. Greubel_, Aug 01 2019
%Y A078988 Row 66 of array A094954.
%Y A078988 Cf. A097316 for S(n, 66).
%Y A078988 Row 4 of array A188647.
%K A078988 nonn,easy
%O A078988 0,2
%A A078988 _Wolfdieter Lang_, Jan 10 2003