cp's OEIS Frontend

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.

Showing 1-10 of 10 results.

A054413 a(n) = 7*a(n-1) + a(n-2), with a(0)=1 and a(1)=7.

Original entry on oeis.org

1, 7, 50, 357, 2549, 18200, 129949, 927843, 6624850, 47301793, 337737401, 2411463600, 17217982601, 122937341807, 877779375250, 6267392968557, 44749530155149, 319514104054600, 2281348258537349, 16288951913816043, 116304011655249650, 830417033500563593
Offset: 0

Views

Author

Henry Bottomley, May 10 2000

Keywords

Comments

In general, sequences with recurrence a(n) = k*a(n-1) + a(n-2) and a(0)=1 (and a(-1)=0) have the generating function 1/(1-k*x-x^2). If k is odd (k>=3) they satisfy a(3n) = b(5n), a(3n+1) = b(5n+3), a(3n+2) = 2*b(5n+4) where b(n) is the sequence of denominators of continued fraction convergents to sqrt(k^2+4). [If k is even then a(n) is the sequence of denominators of continued fraction convergents to sqrt(k^2/4+1).]
a(p-1) == 53^((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]
From Johannes W. Meijer, Jun 12 2010: (Start)
For the sequence given above k=7 which implies that it is associated with A041091.
For a similar statement about sequences with recurrence a(n) = k*a(n-1) + a(n-2) but with a(0) = 2, and a(-1) = 0, see A086902; a sequence that is associated with A041090.
For more information follow the Khovanova link and see A087130, A140455 and A178765.
(End)
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 7'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,...,7} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
From Michael A. Allen, Feb 21 2023: (Start)
Also called the 7-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 7 kinds of squares available. (End)

Crossrefs

Row n=7 of A073133, A172236 and A352361.
Cf. A099367 (squares).

Programs

Formula

a(3n) = A041091(5n), a(3n+1) = A041091(5n+3), a(3n+2) = 2*A041091(5n+4).
G.f.: 1/(1 - 7x - x^2).
a(n) = U(n, 7*i/2)*(-i)^n with i^2=-1 and Chebyshev's U(n, x/2) = S(n, x) polynomials. See A049310.
a(n) = F(n, 7), the n-th Fibonacci polynomial evaluated at x=7. - T. D. Noe, Jan 19 2006
From Sergio Falcon, Sep 24 2007: (Start)
a(n) = (sigma^n - (-sigma)^(-n))/(sqrt(53)) with sigma = (7+sqrt(53))/2;
a(n) = Sum_{i=0..floor((n-1)/2)} binomial(n-1-i,i)*7^(n-1-2i). (End)
a(n) = ((7 + sqrt(53))^n - (7 - sqrt(53))^n)/(2^n*sqrt(53)). Offset 1. a(3)=50. - Al Hakanson (hawkuu(AT)gmail.com), Jan 17 2009
From Johannes W. Meijer, Jun 12 2010: (Start)
a(2n+1) = 7*A097836(n), a(2n) = A097838(n).
Lim_{k->oo} a(n+k)/a(k) = (A086902(n) + A054413(n-1)*sqrt(53))/2.
Lim_{n->oo} A086902(n)/A054413(n-1) = sqrt(53).
(End)
Sum_{n>=0} (-1)^n/(a(n)*a(n+1)) = (sqrt(53)-7)/2. - Vladimir Shevelev, Feb 23 2013
From Kai Wang, Feb 24 2020: (Start)
Sum_{m>=0} 1/(a(m)*a(m+2)) = 1/49.
Sum_{m>=0} 1/(a(2*m)*a(2*m+2)) = (sqrt(53)-7)/14.
In general, for sequences with recurrence f(n)= k*f(n-1)+f(n-2) and f(0)=1,
Sum_{m>=0} 1/(f(m)*f(m+2)) = 1/(k^2).
Sum_{m>=0} 1/(f(2*m)*f(2*m+2)) = (sqrt(k^2+4) - k)/(2*k). (End)
E.g.f.: (1/53)*exp(7*x/2)*(53*cosh(sqrt(53)*x/2) + 7*sqrt(53)*sinh(sqrt(53)*x/2)). - Stefano Spezia, Feb 26 2020
G.f.: x/(1 - 7*x - x^2) = Sum_{n >= 0} x^(n+1) *( Product_{k = 1..n} (m*k + 7 - m + x)/(1 + m*k*x) ) for arbitrary m (a telescoping series). - Peter Bala, May 08 2024

Extensions

Formula corrected by Johannes W. Meijer, May 30 2010, Jun 02 2010
Extended by T. D. Noe, May 23 2011

A086902 a(n) = 7*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 7.

Original entry on oeis.org

2, 7, 51, 364, 2599, 18557, 132498, 946043, 6754799, 48229636, 344362251, 2458765393, 17555720002, 125348805407, 894997357851, 6390330310364, 45627309530399, 325781497023157, 2326097788692498, 16608466017870643, 118585359913786999, 846705985414379636
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Sep 18 2003

Keywords

Comments

a(n+1)/a(n) converges to (7+sqrt(53))/2 = 7.14005... = A176439.
Lim a(n)/a(n+1) as n approaches infinity = 0.1400549... = 2/(7+sqrt(53)) = (sqrt(53)-7)/2 = 1/A176439 = A176439 - 7.
From Johannes W. Meijer, Jun 12 2010: (Start)
In general sequences with recurrence a(n) = k*a(n-1)+a(n-2) with a(0)=2 and a(1)=k [and a(-1)=0] have generating function (2-k*x)/(1-k*x-x^2). If k is odd (k>=3) they satisfy a(3n+1) = b(5n), a(3n+2)=b(5*n+3), a(3n+3)=2*b(5n+4) where b(n) is the sequence of numerators of continued fraction convergents to sqrt(k^2+4). [If k is even then a(n)/2, for n>=1, is the sequence of numerators of continued fraction convergents to sqrt(k^2/4+1).]
For the sequence given above k=7 which implies that it is associated with A041090.
For a similar statement about sequences with recurrence a(n) = k*a(n-1)+a(n-2) but with a(0)=1 [and a(-1)=0] see A054413; a sequence that is associated with A041091.
For more information follow the Khovanova link and see A087130, A140455 and A178765.
(End)

Examples

			a(4) = 7*a(3) + a(2) = 7*364 + 51 = 2599.
		

Crossrefs

Cf. A000032 (k=1), A006497 (k=3), A087130 (k=5), A086902 (k=7), A087798 (k=9), A001946 (k=11), A088316 (k=13), A090301 (k=15), A090306 (k=17). - Johannes W. Meijer, Jun 12 2010

Programs

  • Magma
    I:=[2,7]; [n le 2 select I[n] else 7*Self(n-1)+Self(n-2): n in [1..30]]; // Vincenzo Librandi, Sep 19 2016
  • Mathematica
    RecurrenceTable[{a[0] == 2, a[1] == 7, a[n] == 7 a[n-1] + a[n-2]}, a, {n, 30}] (* Vincenzo Librandi, Sep 19 2016 *)
    LinearRecurrence[{7,1},{2,7},30] (* Harvey P. Dale, May 25 2023 *)
  • PARI
    a(n)=([0,1; 1,7]^n*[2;7])[1,1] \\ Charles R Greathouse IV, Apr 06 2016
    

Formula

a(n) = ((7+sqrt(53))/2)^n + ((7-sqrt(53))/2)^n.
E.g.f. : 2exp(7x/2)cosh(sqrt(53)x/2); a(n)=2^(1-n)sum{k=0..floor(n/2), C(n, 2k)53^k7^(n-2k)}. a(n)=2T(n, 7i/2)(-i)^n with T(n, x) Chebyshev's polynomials of the first kind (see A053120) and i^2=-1. - Paul Barry, Nov 15 2003
G.f.: (2-7x)/(1-7x-x^2). - Philippe Deléham, Nov 16 2008
From Johannes W. Meijer, Jun 12 2010: (Start)
a(2n+1) = 7*A097837(n), a(2n) = A099368(n).
a(3n+1) = A041090(5n), a(3n+2) = A041090(5*n+3), a(3n+3) = 2*A041090(5n+4).
Limit(a(n+k)/a(k), k=infinity) = (A086902(n) + A054413(n-1)*sqrt(53))/2.
Limit(A086902(n)/A054413(n-1), n=infinity) = sqrt(53). (End)

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

Views

Author

Keywords

Crossrefs

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

Programs

  • Magma
    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
  • Mathematica
    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 *)

Formula

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

Extensions

More terms from Vincenzo Librandi, Dec 10 2013

A041047 Denominators of continued fraction convergents to sqrt(29).

Original entry on oeis.org

1, 2, 3, 5, 13, 135, 283, 418, 701, 1820, 18901, 39622, 58523, 98145, 254813, 2646275, 5547363, 8193638, 13741001, 35675640, 370497401, 776670442, 1147167843, 1923838285, 4994844413, 51872282415
Offset: 0

Views

Author

Keywords

Comments

The terms of this sequence can be constructed with the terms of sequence A052918.
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. - Johannes W. Meijer, Jun 12 2010

Crossrefs

Programs

  • Magma
    I:=[1, 2, 3, 5, 13, 135, 283, 418, 701, 1820]; [n le 10 select I[n] else 140*Self(n-5)+Self(n-10): n in [1..50]]; // Vincenzo Librandi, Dec 10 2013
  • Mathematica
    Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[29],n]]],{n,1,50}] (* Vladimir Joseph Stephan Orlovsky, Mar 18 2011 *)
    Denominator[Convergents[Sqrt[29], 30]] (* Vincenzo Librandi, Dec 10 2013 *)

Formula

a(5*n) = A052918(3*n), a(5*n+1) = (A052918(3*n+1) - A052918(3*n))/2, a(5*n+2) = (A052918(3*n+1) + A052918(3*n))/2, a(5*n+3) = A052918(3*n+1) and a(5*n+4) = A052918(3*n+2)/2. - Johannes W. Meijer, Jun 12 2010
G.f.: (1 + 2*x + 3*x^2 + 5*x^3 + 13*x^4 - 5*x^5 + 3*x^6 - 2*x^7 + x^8)/(1 - 140*x^5 - x^10). - Peter J. C. Moses, Jul 29 2013
a(n) = 140*a(n-5) + a(n-10). - Vincenzo Librandi, Dec 10 2013

A041090 Numerators of continued fraction convergents to sqrt(53).

Original entry on oeis.org

7, 22, 29, 51, 182, 2599, 7979, 10578, 18557, 66249, 946043, 2904378, 3850421, 6754799, 24114818, 344362251, 1057201571, 1401563822, 2458765393, 8777860001, 125348805407, 384824276222, 510173081629, 894997357851, 3195165155182, 45627309530399
Offset: 0

Views

Author

Keywords

Comments

The terms of this sequence can be constructed with the terms of sequence A086902. For the terms of the periodical 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

Crossrefs

Programs

  • Mathematica
    Numerator[Convergents[Sqrt[53],30]] (* Harvey P. Dale, Sep 24 2013 *)
    CoefficientList[Series[-(x^9 - 7 x^8 + 22 x^7 - 29 x^6 + 51 x^5 + 182 x^4 + 51 x^3 + 29 x^2 + 22 x + 7)/(x^10 + 364 x^5 - 1), {x, 0, 40}], x] (* Vincenzo Librandi, Sep 27 2013 *)

Formula

a(5*n) = A086902(3*n+1), a(5*n+1) = (A086902(3*n+2)-A086902(3*n+1))/2, a(5*n+2) = (A086902(3*n+2)+A086902(3*n+1))/2, a(5*n+3) = A086902(3*n+2) and a(5*n+4) = A086902(3*n+3)/2. - Johannes W. Meijer, Jun 12 2010
G.f.: -(x^9-7*x^8+22*x^7-29*x^6+51*x^5+182*x^4+51*x^3+29*x^2+22*x+7) / (x^10+364*x^5-1). - Colin Barker, Sep 26 2013

Extensions

More terms from 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

Views

Author

Keywords

Comments

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)

Crossrefs

Programs

  • Magma
    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
  • Mathematica
    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 *)

Formula

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

Views

Author

Keywords

Comments

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

Crossrefs

Programs

  • Magma
    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
  • Mathematica
    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 *)

Formula

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

Views

Author

Keywords

Comments

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

Crossrefs

Programs

  • Magma
    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
  • Mathematica
    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 *)

Formula

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

Views

Author

Keywords

Comments

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

Crossrefs

Programs

  • Magma
    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
  • Mathematica
    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 *)

Formula

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

Views

Author

Keywords

Comments

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

Crossrefs

Programs

  • Magma
    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
  • Mathematica
    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 *)

Formula

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
Showing 1-10 of 10 results.