A041047 -id:A041047 - cp's OEIS frontend

A052918 a(0) = 1, a(1) = 5, a(n+1) = 5*a(n) + a(n-1).

1, 5, 26, 135, 701, 3640, 18901, 98145, 509626, 2646275, 13741001, 71351280, 370497401, 1923838285, 9989688826, 51872282415, 269351100901, 1398627786920, 7262490035501, 37711077964425, 195817879857626

Offset: 0

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

A087130(n)^2 - 29*a(n-1)^2 = 4*(-1)^n, n >= 1. - Gary W. Adamson, Jul 01 2003, corrected Oct 07 2008, corrected by Jianing Song, Feb 01 2019
a(p-1) == 29^((p-1)/2) (mod p), for odd primes p. - Gary W. Adamson, Feb 22 2009 [See A087475 for more info about this congruence. - Jason Yuen, Apr 05 2025]
For more information about this type of recurrence, follow the Khovanova link and see A054413, A086902 and A178765. - Johannes W. Meijer, Jun 12 2010
Binomial transform of A015523. - Johannes W. Meijer, Aug 01 2010
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 5's along the main diagonal and 1's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 08 2011
a(n) equals the number of words of length n on alphabet {0,1,...,5} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
From Michael A. Allen, Feb 15 2023: (Start)
Also called the 5-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence.
a(n) is the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are 5 kinds of squares available. (End)

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Michael A. Allen and Kenneth Edwards, Fence tiling derived identities involving the metallonacci numbers squared or cubed, Fib. Q. 60:5 (2022) 5-17.
D. Birmajer, J. B. Gil, and M. D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, Example 8.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 901
Milan Janjić, Hessenberg Matrices and Integer Sequences, J. Int. Seq. 13 (2010) # 10.7.8, section 3.
Milan Janjić, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
Kai Wang, On k-Fibonacci Sequences And Infinite Series List of Results and Examples, 2020.
Index entries for sequences related to Chebyshev polynomials.
Index entries for linear recurrences with constant coefficients, signature (5,1).

Row 5 of A073133, A172236, and A352361.

Cf. A087130, A099365 (squares), A100237, A175184 (Pisano periods), A201005 (prime subsequence).

a:=[1,5];; for n in [3..30] do a[n]:=5*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Oct 16 2019

I:=[1, 5]; [n le 2 select I[n] else 5*Self(n-1)+Self(n-2): n in [1..30]]; // Vincenzo Librandi, Feb 23 2013

R:=PowerSeriesRing(Integers(), 22); Coefficients(R!( 1/(1 - 5*x - x^2) )); // Marius A. Burtea, Oct 16 2019

spec := [S,{S=Sequence(Union(Z,Z,Z,Z,Z,Prod(Z,Z)))},unlabeled]: seq(combstruct[count](spec,size=n), n=0..30);
a[0]:=1: a[1]:=5: for n from 2 to 26 do a[n]:=5*a[n-1]+a[n-2] od: seq(a[n], n=0..30); # Zerinvary Lajos, Jul 26 2006
with(combinat):a:=n->fibonacci(n,5):seq(a(n),n=1..30); # Zerinvary Lajos, Dec 07 2008

LinearRecurrence[{5, 1}, {1, 5}, 30] (* Vincenzo Librandi, Feb 23 2013 *)
Table[Fibonacci[n+1, 5], {n,0,30}] (* Vladimir Reshetnikov, May 08 2016 *)

Vec(1/(1-5*x-x^2)+O(x^30)) \\ Charles R Greathouse IV, Nov 20 2011

[lucas_number1(n,5,-1) for n in range(1, 22)] # Zerinvary Lajos, Apr 24 2009

G.f.: 1/(1 - 5*x - x^2).
a(3n) = A041047(5n), a(3n+1) = A041047(5n+3), a(3n+2) = 2*A041047(5n+4). - Henry Bottomley, May 10 2000
a(n) = Sum_{alpha=RootOf(-1+5*z+z^2)} (1/29)*(5+2*alpha)*alpha^(-1-n).
a(n-1) = (((5 + sqrt(29))/2)^n - ((5 - sqrt(29))/2)^n)/sqrt(29). - Gary W. Adamson, Jul 01 2003
a(n) = U(n, 5*i/2)*(-i)^n with i^2 = -1 and Chebyshev's U(n, x/2) = S(n, x) polynomials. See triangle A049310.
Let M = {{0, 1}, {1, 5}}, then a(n) is the lower-right term of M^n. - Roger L. Bagula, May 29 2005
a(n) = F(n, 5), the n-th Fibonacci polynomial evaluated at x = 5. - T. D. Noe, Jan 19 2006
a(n) = denominator of n-th convergent to [1, 4, 5, 5, 5, ...], for n > 0. Continued fraction [1, 4, 5, 5, 5, ...] = 0.807417596..., the inradius of a right triangle with legs 2 and 5. n-th convergent = A100237(n)/A052918(n), the first few being: 1/1, 4/5, 21/26, 109/135, 566/701, ... - Gary W. Adamson, Dec 21 2007
From Johannes W. Meijer, Jun 12 2010: (Start)
a(2n+1) = 5*A097781(n), a(2n) = A097835(n).
Limit_{k->oo} a(n+k)/a(k) = (A087130(n) + a(n-1)*sqrt(29))/2.
Limit_{n->oo} A087130(n)/a(n-1) = sqrt(29). (End)
From L. Edson Jeffery, Jan 07 2012: (Start)
Define the 2 X 2 matrix A = {{1, 1}, {5, 4}}. Then:
a(n) is the upper-left term of (1/5)*(A^(n+2) - A^(n+1));
a(n) is the upper-right term of A^(n+1);
a(n) is the lower-left term of (1/5)*A^(n+1);
a(n) is the lower-right term of (Sum_{k=0..n} A^k). (End)
Sum_{n>=0} (-1)^n/(a(n)*a(n+1)) = (sqrt(29) - 5)/2. - Vladimir Shevelev, Feb 23 2013
G.f.: x/(1 - 5*x - x^2) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (m*k + 5 - m + x)/(1 + m*k*x) ) for arbitrary m (a telescoping series). - Peter Bala, May 08 2024

A041019 Denominators of continued fraction convergents to sqrt(13).

Original entry on oeis.org

1, 1, 2, 3, 5, 33, 38, 71, 109, 180, 1189, 1369, 2558, 3927, 6485, 42837, 49322, 92159, 141481, 233640, 1543321, 1776961, 3320282, 5097243, 8417525, 55602393, 64019918, 119622311, 183642229, 303264540, 2003229469, 2306494009, 4309723478, 6616217487, 10925940965

Offset: 0

N. J. A. Sloane

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,36,0,0,0,0,1).

Cf. A010122 (continued fraction for sqrt(13)), A041018 (numerators).

Cf. A041047, A041091, A041151, A041227, A041319, A041427 and A041551. - Johannes W. Meijer, Jun 12 2010

I:=[1, 1, 2, 3, 5, 33, 38, 71, 109, 180]; [n le 10 select I[n] else 36*Self(n-5)+Self(n-10): n in [1..50]]; // Vincenzo Librandi, Dec 10 2013

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[13], n]]], {n, 1, 50}] (* Vladimir Joseph Stephan Orlovsky, Mar 16 2011 *)
CoefficientList[Series[((1 - 2 x + 4 x^2 - 3 x^3 + x^4) (1 + 3 x + 4 x^2 + 2 x^3 + x^4))/(1 - 36 x^5 - x^10), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 10 2013 *)
LinearRecurrence[{0,0,0,0,36,0,0,0,0,1},{1,1,2,3,5,33,38,71,109,180},40] (* Harvey P. Dale, Sep 30 2016 *)

From Johannes W. Meijer, Jun 12 2010: (Start)
a(5*n) = A006190(3*n+1),
a(5*n+1) = (A006190(3*n+2) - A006190(3*n+1))/2,
a(5*n+2) = (A006190(3*n+2) + A006190(3*n+1))/2,
a(5*n+3) = A006190(3*n+2) and a(5*n+4) = A006190(3*n+3)/2. (End)
G.f.: ((1 - 2*x + 4*x^2 - 3*x^3 + x^4)*(1 + 3*x + 4*x^2 + 2*x^3 + x^4))/(1 - 36*x^5 - x^10). - Peter J. C. Moses, Jul 29 2013
a(n) = A010122(n)*a(n-1) + a(n-2), a(0)=1, a(-1)=0. - Paul Weisenhorn, Aug 17 2018

More terms from Vincenzo Librandi, Dec 10 2013

A041046 Numerators of continued fraction convergents to sqrt(29).

Original entry on oeis.org

5, 11, 16, 27, 70, 727, 1524, 2251, 3775, 9801, 101785, 213371, 315156, 528527, 1372210, 14250627, 29873464, 44124091, 73997555, 192119201, 1995189565, 4182498331, 6177687896, 10360186227, 26898060350

Offset: 0

N. J. A. Sloane

From Johannes W. Meijer, Jun 12 2010: (Start)
The terms of this sequence can be constructed with the terms of sequence A087130.
For the terms of the periodical sequence of the continued fraction for sqrt(29) see A010128. We observe that its period is five. The decimal expansion of sqrt(29) is A010484.  (End)

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,140,0,0,0,0,1).

Cf. A041047, A041018, A041046, A041090, A041150, A041226, A041318, A041426, A041550, A010484.

Table[Numerator[FromContinuedFraction[ContinuedFraction[Sqrt[29],n]]],{n,1,50}] (* Vladimir Joseph Stephan Orlovsky, Mar 18 2011 *)
Numerator[Convergents[Sqrt[29], 30]] (* Vincenzo Librandi, Oct 28 2013 *)
LinearRecurrence[ {0,0,0,0,140,0,0,0,0,1},{5,11,16,27,70,727,1524,2251,3775,9801},30] (* Harvey P. Dale, Jun 10 2021 *)

a(5*n) = A087130(3*n+1), a(5*n+1) = (A087130(3*n+2) - A087130(3*n+1))/2, a(5*n+2) = ( A087130(3*n+2) + A087130(3*n+1))/2, a(5*n+3) = A087130(3*n+2) and a(5*n+4) = A087130(3*n+3)/2. - Johannes W. Meijer, Jun 12 2010

G.f.: (5 + 11*x + 16*x^2 + 27*x^3 + 70*x^4 + 27*x^5 - 16*x^6 + 11*x^7 - 5*x^8 + x^9)/(1 - 140*x^5 - x^10) - Peter J. C. Moses, Jul 29 2013

A041091 Denominators of continued fraction convergents to sqrt(53).

Original entry on oeis.org

1, 3, 4, 7, 25, 357, 1096, 1453, 2549, 9100, 129949, 398947, 528896, 927843, 3312425, 47301793, 145217804, 192519597, 337737401, 1205731800, 17217982601, 52859679603, 70077662204, 122937341807, 438889687625, 6267392968557, 19241068593296, 25508461561853

Offset: 0

N. J. A. Sloane

The terms of this sequence can be constructed with the terms of sequence A054413. For the terms of the periodic sequence of the continued fraction for sqrt(53) see A010139. We observe that its period is five. The decimal expansion of sqrt(53) is A010506. - Johannes W. Meijer, Jun 12 2010

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,364,0,0,0,0,1).

Cf. A010506, A041090.

Cf. A041019, A041047, A041151, A041227, A041319, A041427 and A041551. - Johannes W. Meijer, Jun 12 2010

convert(sqrt(53), confrac, 30, cvgts): denom(cvgts); # Wesley Ivan Hurt, Dec 17 2013

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[53], n]]], {n, 1, 50}] (* Vladimir Joseph Stephan Orlovsky, Jun 23 2011 *)
Denominator[Convergents[Sqrt[53], 30]] (* Vincenzo Librandi, Oct 24 2013 *)
LinearRecurrence[{0,0,0,0,364,0,0,0,0,1},{1,3,4,7,25,357,1096,1453,2549,9100},30] (* Harvey P. Dale, Nov 13 2019 *)

a(5*n) = A054413(3*n), a(5*n+1) = (A054413(3*n+1) - A054413(3*n))/2, a(5*n+2)= (A054413(3*n+1) + A054413(3*n))/2, a(5*n+3) = A054413(3*n+1) and a(5*n+4) = A054413(3*n+2)/2. - Johannes W. Meijer, Jun 12 2010

G.f.: -(x^8-3*x^7+4*x^6-7*x^5+25*x^4+7*x^3+4*x^2+3*x+1) / (x^10+364*x^5-1). - Colin Barker, Sep 26 2013

A041151 Denominators of continued fraction convergents to sqrt(85).

Original entry on oeis.org

1, 4, 5, 9, 41, 747, 3029, 3776, 6805, 30996, 564733, 2289928, 2854661, 5144589, 23433017, 426938895, 1731188597, 2158127492, 3889316089, 17715391848, 322766369353, 1308780869260, 1631547238613, 2940328107873, 13392859670105, 244011802169763, 989440068349157

Offset: 0

N. J. A. Sloane

From Johannes W. Meijer, Jun 12 2010: (Start)
The a(n) terms of this sequence can be constructed with the terms of sequence A099371.
For the terms of the periodic sequence of the continued fraction for sqrt(85) see A010158. We observe that its period is five. The decimal expansion of sqrt(85) is A010536. (End)

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,756,0,0,0,0,1).

Cf. A041150, A041019, A041047, A041091, A041151, A041227, A041319, A041427, A041551.

I:=[1, 4, 5, 9, 41, 747, 3029, 3776, 6805, 30996]; [n le 10 select I[n] else 756*Self(n-5)+Self(n-10): n in [1..30]]; // Vincenzo Librandi, Dec 12 2013

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[85], n]]], {n, 1, 50}] (* Vladimir Joseph Stephan Orlovsky, Jun 23 2011 *)
Denominator[Convergents[Sqrt[85], 30]] (* Vincenzo Librandi, Dec 12 2013 *)

From Johannes W. Meijer, Jun 12 2010: (Start)
a(5*n) = A099371(3*n+1), a(5*n+1) = (A099371(3*n+2)-A099371(3*n+1))/2, a(5*n+2) = (A099371(3*n+2)+A099371(3*n+1))/2, a(5*n+3):= A099371(3*n+2) and a(5*n+4) = A099371(3*n+3)/2. (End)
G.f.: -(x^8-4*x^7+5*x^6-9*x^5+41*x^4+9*x^3+5*x^2+4*x+1) / (x^10+756*x^5-1). - Colin Barker, Nov 11 2013
a(n) = 756*a(n-5) + a(n-10). - Vincenzo Librandi, Dec 12 2013

A041227 Denominators of continued fraction convergents to sqrt(125).

Original entry on oeis.org

1, 5, 6, 11, 61, 1353, 6826, 8179, 15005, 83204, 1845493, 9310669, 11156162, 20466831, 113490317, 2517253805, 12699759342, 15217013147, 27916772489, 154800875592, 3433536035513, 17322481053157, 20756017088670, 38078498141827, 211148507797805

Offset: 0

N. J. A. Sloane

The a(n) terms of this sequence can be constructed with the terms of sequence A049666. For the terms of the periodical sequence of the continued fraction for sqrt(125) see A010186. We observe that its period is five. - Johannes W. Meijer, Jun 12 2010

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,1364,0,0,0,0,1).

Cf. A041226, A010186, A041019, A041047, A041091, A041151, A041227, A041319, A041427, A041551.

I:=[1, 5, 6, 11, 61, 1353, 6826, 8179, 15005, 83204]; [n le 10 select I[n] else 1364*Self(n-5)+Self(n-10): n in [1..40]]; // Vincenzo Librandi, Dec 13 2013

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[125], n]]], {n, 1, 50}] (* Vladimir Joseph Stephan Orlovsky, Jun 23 2011 *)
Denominator[Convergents[Sqrt[125], 30]]  (* Vincenzo Librandi, Dec 13 2013 *)
LinearRecurrence[{0,0,0,0,1364,0,0,0,0,1},{1,5,6,11,61,1353,6826,8179,15005,83204},30] (* Harvey P. Dale, Apr 29 2022 *)

a(5*n) = A049666(3*n+1), a(5*n+1) = (A049666(3*n+2) - A049666(3*n+1))/2, a(5*n+2) = (A049666(3*n+2)+A049666(3*n+1))/2, a(5*n+3):= A049666(3*n+2) and a(5*n+4) = A049666(3*n+3)/2. - Johannes W. Meijer, Jun 12 2010
G.f.: -(x^8 -5*x^7 +6*x^6 -11*x^5 +61*x^4 +11*x^3 +6*x^2 +5*x +1) / ((x^2 +4*x -1)*(x^4 -7*x^3 +19*x^2 -3*x +1)*(x^4 +3*x^3 +19*x^2 +7*x +1)). - Colin Barker, Nov 12 2013
a(n) = 1364*a(n-5) + a(n-10). - Vincenzo Librandi, Dec 13 2013

A041319 Denominators of continued fraction convergents to sqrt(173).

Original entry on oeis.org

1, 6, 7, 13, 85, 2223, 13423, 15646, 29069, 190060, 4970629, 30013834, 34984463, 64998297, 424974245, 11114328667, 67110946247, 78225274914, 145336221161, 950242601880, 24851643870041, 150060105822126, 174911749692167, 324971855514293, 2124742882777925

Offset: 0

N. J. A. Sloane

The a(n) terms of this sequence can be constructed with the terms of sequence A140455. For the terms of the periodical sequence of the continued fraction for sqrt(173) see A010217. We observe that its period is five. - Johannes W. Meijer, Jun 12 2010

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0, 0, 0, 0, 2236, 0, 0, 0, 0, 1).

Cf. A041318, A010217, A041019, A041047, A041091, A041151, A041227, A041319, A041427 and A041551.

I:=[1,6,7,13,85,2223,13423,15646,29069,190060]; [n le 10 select I[n] else 2236*Self(n-5)+Self(n-10): n in [1..40]]; // Vincenzo Librandi, Dec 15 2013

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[173], n]]], {n, 1, 50}] (* Vladimir Joseph Stephan Orlovsky, Jun 23 2011 *)
Denominator[Convergents[Sqrt[173], 30]] (* Vincenzo Librandi, Dec 15 2013 *)
LinearRecurrence[{0,0,0,0,2236,0,0,0,0,1},{1,6,7,13,85,2223,13423,15646,29069,190060},30] (* Harvey P. Dale, Sep 19 2020 *)

a(5*n) = A140455(3*n+1), a(5*n+1) = (A140455(3*n+2) - A140455(3*n+1))/2, a(5*n+2) = (A140455(3*n+2)+A140455(3*n+1))/2, a(5*n+3) = A140455(3*n+2) and a(5*n+4) = A140455(3*n+3)/2. - Johannes W. Meijer, Jun 12 2010
G.f.: -(x^8-6*x^7+7*x^6-13*x^5+85*x^4+13*x^3+7*x^2+6*x+1) / (x^10+2236*x^5-1). - Colin Barker, Nov 12 2013
a(n) = 2236*a(n-5) + a(n-10). - Vincenzo Librandi, Dec 15 2013

A041427 Denominators of continued fraction convergents to sqrt(229).

Original entry on oeis.org

1, 7, 8, 15, 113, 3405, 23948, 27353, 51301, 386460, 11645101, 81902167, 93547268, 175449435, 1321693313, 39826248825, 280105435088, 319931683913, 600037119001, 4520191516920, 136205782626601, 957960669903127, 1094166452529728, 2052127122432855

Offset: 0

N. J. A. Sloane

The a(n) terms of this sequence can be constructed with the terms of sequence A154597. For the terms of the periodical sequence of the continued fraction for sqrt(229) see A040213. We observe that its period is five. - Johannes W. Meijer, Jun 12 2010

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0, 0, 0, 0, 3420, 0, 0, 0, 0, 1).

Cf. A041426, A166125, A040213, A041019, A041047, A041091, A041151, A041227, A041319, A041427 and A041551.

I:=[1,7,8,15,113,3405,23948,27353,51301,386460]; [n le 10 select I[n] else 3420*Self(n-5)+Self(n-10): n in [1..40]]; // Vincenzo Librandi, Dec 17 2013

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[229], n]]], {n, 1, 50}] (* Vladimir Joseph Stephan Orlovsky, Jun 23 2011 *)
Denominator[Convergents[Sqrt[229], 30]] (* Vincenzo Librandi, Dec 17 2013 *)
LinearRecurrence[{0,0,0,0,3420,0,0,0,0,1},{1,7,8,15,113,3405,23948,27353,51301,386460},30] (* Harvey P. Dale, Oct 14 2020 *)

a(5*n) = A154597(3*n+1), a(5*n+1) = (A154597(3*n+2) - A154597(3*n+1))/2, a(5*n+2) = (A154597(3*n+2) + A154597(3*n+1))/2, a(5*n+3) = A154597(3*n+2) and a(5*n+4) = A154597(3*n+3)/2. - Johannes W. Meijer, Jun 12 2010
G.f.: -(x^8 -7*x^7 +8*x^6 -15*x^5 +113*x^4 +15*x^3 +8*x^2 +7*x +1) / (x^10 +3420*x^5 -1). - Colin Barker, Nov 12 2013
a(n) = 3420*a(n-5) + a(n-10) for n>9. - Vincenzo Librandi, Dec 17 2013

A041551 Denominators of continued fraction convergents to sqrt(293).

Original entry on oeis.org

1, 8, 9, 17, 145, 4947, 39721, 44668, 84389, 719780, 24556909, 197175052, 221731961, 418907013, 3572988065, 121900501223, 978776997849, 1100677499072, 2079454496921, 17736313474440, 605114112627881, 4858649214497488, 5463763327125369, 10322412541622857

Offset: 0

N. J. A. Sloane

The a(n) terms of this sequence can be constructed with the terms of sequence A178765. For the terms of the periodical sequence of the continued fraction for sqrt(293) see A040275. We observe that its period is five. - Johannes W. Meijer, Jun 12 2010

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0, 0, 0, 0, 4964, 0, 0, 0, 0, 1).

Cf. A041550, A040275, A041019, A041047, A041091, A041151, A041227, A041319, A041427 and A041551.

I:=[1,8,9,17,145,4947,39721,44668,84389,719780]; [n le 10 select I[n] else 4964*Self(n-5)+Self(n-10): n in [1..40]]; // Vincenzo Librandi, Dec 20 2013

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[293], n]]], {n, 1, 50}] (* Vladimir Joseph Stephan Orlovsky, Jun 23 2011 *)
Denominator[Convergents[Sqrt[293 ], 30]] (* Vincenzo Librandi, Dec 20 2013 *)

a(5n) = A178765(3n), a(5n+1) = (A178765(3n+1) - A178765(3n))/2, a(5n+2) = (A178765(3n+1) + A178765(3n))/2, a(5n+3) = A178765(3n+1) and a(5n+4) = A178765(3n+2)/2. - Johannes W. Meijer, Jun 12 2010
G.f.: -(x^8-8*x^7+9*x^6-17*x^5+145*x^4+17*x^3+9*x^2+8*x+1) / (x^10+4964*x^5-1). - Colin Barker, Nov 12 2013
a(n) = 4964*a(n-5) + a(n-10) for n>9. - Vincenzo Librandi, Dec 20 2013

A319750 a(n) is the denominator of the Heron sequence with h(0) = 3.

Original entry on oeis.org

1, 3, 33, 3927, 55602393, 11147016454528647, 448011292165037607943004375755833, 723685043824607606355691108666081531638582859833105061571146291527

Offset: 0

Paul Weisenhorn, Sep 27 2018

The numerators of the Heron sequence are in A319749.
There is the following relationship between the denominator of the Heron sequence and the denominator of the continued fraction A041018(n)/ A041019(n) convergent to sqrt(13).
n even: a(n) = A041019((5*2^n-5)/3).
n  odd: a(n) = A041019((5*2^n-1)/3).
General: all numbers c(n) = A078370(n) = (2*n+1)^2 + 4 have the same relationship between the denominator of the Heron sequence and the denominator of the continued fraction convergent to 2*n+1.
sqrt(c(n)) has the continued fraction [2*n+1; n, 1, 1, n, 4*n+2].
hn(n)^2 - c(n)*hd(n)^2 = 4 for n > 1.

			A078370(2) = 29.
hd(0) = A041047(0) = 1, hd(1) = A041047(3) = 5,
hd(2) = A041047(5) = 135, hd(3) = A041047(13) = 38145.

Cf. A041018, A041019, A078370, A010122, A041047, A319749.

hn[0]:=3: hd[0]:=1:
for n from 1 to 6 do
  hn[n]:=(hn[n-1]^2+13*hd[n-1]^2)/2:
  hd[n]:=hn[n-1]*hd[n-1]:
  printf("%5d%40d%40d\n", n, hn[n], hd[n]):
end do:

def aupton(nn):
    hn, hd, alst = 3, 1, [1]
    for n in range(nn):
        hn, hd = (hn**2 + 13*hd**2)//2, hn*hd
        alst.append(hd)
    return alst
print(aupton(7)) # Michael S. Branicky, Mar 15 2022

h(n) = hn(n)/hd(n), hn(0) = 3, hd(0) = 1.
hn(n+1) = (hn(n)^2 + 13*hd(n)^2)/2.
hd(n+1) = hn(n)*hd(n).
A041018(n) = A010122(n)*A041018(n-1) + A041018(n-2).
A041019(n) = A010122(n)*A041019(n-1) + A041019(n-2).
a(0) = 1, a(1) = 3 and a(n) = 2*T(2^(n-2), 11/2)*a(n-1) for n >= 2, where T(n,x) denotes the n-th Chebyshev polynomial of the first kind. - Peter Bala, Mar 16 2022

a(5) corrected and terms a(6) and a(7) added by Peter Bala, Mar 15 2022

cp's OEIS Frontend

A052918 a(0) = 1, a(1) = 5, a(n+1) = 5*a(n) + a(n-1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Magma

Maple

Mathematica

PARI

Sage

Formula

A041019 Denominators of continued fraction convergents to sqrt(13).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

Formula

Extensions

A041046 Numerators of continued fraction convergents to sqrt(29).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A041091 Denominators of continued fraction convergents to sqrt(53).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A041151 Denominators of continued fraction convergents to sqrt(85).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

Formula

A041227 Denominators of continued fraction convergents to sqrt(125).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

Formula

A041319 Denominators of continued fraction convergents to sqrt(173).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma