cp's OEIS Frontend

A097315 Pell equation solutions (3*b(n))^2 - 10*a(n)^2 = -1 with b(n) = A097314(n), n >= 0.

%I A097315 #61 Apr 24 2025 17:39:48
%S A097315 1,37,1405,53353,2026009,76934989,2921503573,110940200785,
%T A097315 4212806126257,159975692596981,6074863512559021,230684837784645817,
%U A097315 8759948972303982025,332647376109766671133,12631840343198829521029,479677285665445755127969,18215105014943739865341793,691694313282196669127860165
%N A097315 Pell equation solutions (3*b(n))^2 - 10*a(n)^2 = -1 with b(n) = A097314(n), n >= 0.
%C A097315 Hypotenuses of primitive Pythagorean triples in A195616 and A195617. - _Clark Kimberling_, Sep 22 2011
%H A097315 Indranil Ghosh, <a href="/A097315/b097315.txt">Table of n, a(n) for n = 0..631</a>
%H A097315 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 A097315 Tanya Khovanova, <a href="http://www.tanyakhovanova.com/RecursiveSequences/RecursiveSequences.html">Recursive Sequences</a>
%H A097315 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 A097315 <a href="/index/Ch#Cheby">Index entries for sequences related to Chebyshev polynomials.</a>
%H A097315 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (38,-1).
%F A097315 a(n) = S(n, 38) - S(n-1, 38) = T(2*n+1, sqrt(10))/sqrt(10), with Chebyshev polynomials of the second and first kind. See A049310 for the triangle of S(n, x)= U(n, x/2) coefficients. S(-1, x) := 0 =: U(-1, x); and A053120 for the T-triangle.
%F A097315 a(n) = ((-1)^n)*S(2*n, 6*i) with the imaginary unit i and Chebyshev polynomials S(n, x) with coefficients shown in A049310.
%F A097315 G.f.: (1-x)/(1-38*x+x^2).
%F A097315 a(n) = 38*a(n-1) - a(n-2) for n > 1. - _Philippe Deléham_, Nov 18 2008
%F A097315 a(n) = sqrt(2+(19-6*sqrt(10))^(1+2*n) + (19+6*sqrt(10))^(1+2*n))/(2*sqrt(10)). - _Gerry Martens_, Jun 04 2015
%F A097315 a(n) = A078987(n) - A078987(n-1). - _R. J. Mathar_, Dec 05 2015
%F A097315 a(n) = A005668(2*n+1). - _Michael Somos_, Feb 24 2023
%F A097315 E.g.f.: exp(19*x)*(10*cosh(6*sqrt(10)*x) + 3*sqrt(10)*sinh(6*sqrt(10)*x))/10. - _Stefano Spezia_, Apr 24 2025
%e A097315 (x,y) = (3,1), (117,37), (4443,1405), ... give the positive integer solutions to x^2 - 10*y^2 = -1.
%e A097315 G.f. = 1 + 37*x + 1405*x^2 + 53353*x^3 + ... - _Michael Somos_, Feb 24 2023
%t A097315 CoefficientList[Series[(1-x)/(1-38x+x^2), {x,0,20}], x] (* _Michael De Vlieger_, Feb 04 2017 *)
%t A097315 LinearRecurrence[{38,-1}, {1,37}, 21] (* _G. C. Greubel_, Aug 01 2019 *)
%o A097315 (PARI) Vec((1-x)/(1-38*x+x^2) + O(x^20)) \\ _Michel Marcus_, Jun 04 2015
%o A097315 (Magma) I:=[1, 37]; [n le 2 select I[n] else 38*Self(n-1) - Self(n-2): n in [1..30]]; // _G. C. Greubel_, Aug 01 2019
%o A097315 (Sage) ((1-x)/(1-38*x+x^2)).series(x, 20).coefficients(x, sparse=False) # _G. C. Greubel_, Aug 01 2019
%o A097315 (GAP) a:=[1,37];; for n in [3..20] do a[n]:=38*a[n-1]-a[n-2]; od; a; # _G. C. Greubel_, Aug 01 2019
%o A097315 (Python)
%o A097315 from itertools import islice
%o A097315 def A097315_gen(): # generator of terms
%o A097315     x, y = 30, 10
%o A097315     while True:
%o A097315         yield y//10
%o A097315         x, y = x*19+y*60, x*6+y*19
%o A097315 A097315_list = list(islice(A097315_gen(),20)) # _Chai Wah Wu_, Apr 24 2025
%Y A097315 Cf. A078987 (partial sums), A097314, A195616, A195617, A049310, A053120, A005668.
%Y A097315 Row 3 of array A188647.
%Y A097315 Cf. A221874.
%Y A097315 Cf. similar sequences listed in A238379.
%K A097315 nonn,easy
%O A097315 0,2
%A A097315 _Wolfdieter Lang_, Aug 31 2004
%E A097315 Typo in recurrence formula corrected by Laurent Bonaventure (bonave(AT)free.fr), Oct 03 2010
%E A097315 More terms added by _Indranil Ghosh_, Feb 04 2017