A001077 -id:A001077 - cp's OEIS frontend

A180140 Eight rooks and one berserker on a 3 X 3 chessboard. G.f.: (1+x+x^2)/(1-3x-5x^2).

1, 4, 18, 74, 312, 1306, 5478, 22964, 96282, 403666, 1692408, 7095554, 29748702, 124723876, 522915138, 2192364794, 9191670072, 38536834186, 161568852918, 677390729684, 2840016453642, 11907003009346, 49921091296248

Offset: 0

Johannes W. Meijer, Aug 13 2010, Jun 15 2013

a(n) gives the number of n-move routes of a fairy chess piece starting in a given side square (m = 2, 4, 6,  8) on a 3 X 3 chessboard. This fairy chess piece behaves like a rook on the four side and four corner (m = 1, 3, 7, 9) squares but on the center square (m = 5) it goes berserk and turns into a berserker. For this sequence, the berserker can move to three of the side squares and three of the corners from the center.
The berserker is one of the Lewis chessmen which were discovered in 1831 on the Isle of Lewis. They are carved from walrus ivory in Scandinavian style of the 12th century. The pawns look like decorated tombstones. The pieces have all human representations with facial expressions varying from gloom to anger. Some of the rooks show men biting their shield in the manner of berserkers. According to Hooper and Whyld none looks happy.
Let A be the adjacency matrix of the graph G, where V(G) = {v1, v2, v3, v4, v5, v6, v7, v8, v9}. Then the (m, k) entry of A^n is the number of different vm-vk walks of length n in G, see the Chartrand reference. In the adjacency matrix A, see the Maple program, the A[1], A[3], A[7] and A[9] vectors represent the rook moves on the corner squares, the A[2], A[4], A[6] and A[8] vectors represent the rook moves on the side squares and the A[5] vector represents the moves of the berserker. On a 3 X 3 chessboard there are 2^9 = 512 ways a berserker could move from the center square (off the center the berserker behaves like a rook) so there are 512 different berserkers.
For the side squares the 512 berserker vectors lead to 42 different sequences, see the overview of berserker sequences. There are 16 berserker vectors that lead to the sequence given above. Their decimal [binary] values are: 111 [001 101 111] , 207 [011 001 111], 231 [011 100 111], 237 [011 101 101], 303 [100 101 111], 363 [101 101 011], 366 [101 101 110], 399 [110 001 111], 423 [110 100 111], 429 [110 101 101], 459 [111 001 011], 462 [111 001 110], 483 [111 100 011], 486 [111 100 110], 489 [111 101 001] and 492 [111 101 100]. These berserker vectors lead for the corner squares to sequence 4*A179606 (with leading term 1 added) and for the central square to sequence 6*A179606 (with leading term 1 added).
This sequence belongs to a family of sequences with GF(x)=(1+x-k*x^2)/(1-3*x+(k-4)*x^2), see A180142.

Gary Chartrand, Introductory Graph Theory, pp. 217-221, 1984.
David Hooper and Kenneth Whyld, The Oxford Companion to Chess, pp. 131, 225, 1992.

Indranil Ghosh, Table of n, a(n) for n = 0..1603
Dougie MacLean, The Lewis Chessmen: "Marching Mystery"
Johannes W. Meijer, The berserker sequences.
Wikipedia, Lewis Chessmen
Index entries for linear recurrences with constant coefficients, signature (3,5)

Cf. A180141 (corner squares) and A180147 (central square).

Cf. Berserker sequences side squares: 4*A007482 (with leading 1 added), A180144, A003500 (n>=1 and a(0)=1), A180142, A000302, A180140 (this sequence), 2*A001077 (n>=1 and a(0)=1), A180146, 4*A154964 (n>=1 and a(0)=1), 4*A123347 (with leading 1 added).

nmax:=22; m:=2; A[1]:=[0, 1, 1, 1, 0, 0, 1, 0, 0]: A[2]:=[1, 0, 1, 0, 1, 0, 0, 1, 0]: A[3]:= [1, 1, 0, 0, 0, 1, 0, 0, 1]: A[4]:= [1, 0, 0, 0, 1, 1, 1, 0, 0]: A[5]:=[0, 0, 1, 1, 0, 1, 1, 1, 1]: A[6]:=[0, 0, 1, 1, 1, 0, 0, 0, 1]: A[7]:=[1, 0, 0, 1, 0, 0, 0, 1, 1]: A[8]:=[0, 1, 0, 0, 1, 0, 1, 0, 1]: A[9]:=[0, 0, 1, 0, 0, 1, 1, 1, 0]: A:= Matrix([A[1], A[2], A[3], A[4], A[5], A[6], A[7], A[8], A[9]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m, k], k=1..9): od: seq(a(n), n=0..nmax);

CoefficientList[Series[(1+x+x^2)/(1-3*x-5*x^2), {x, 0, 22}],x] (* or *) LinearRecurrence[{3,5,0},{1,4,18},23] (* Indranil Ghosh, Mar 05 2017 *)

print(Vec((1 + x + x^2)/(1- 3*x - 5*x^2) + O(x^23))); \\ Indranil Ghosh, Mar 05 2017

G.f.: (1+x+x^2)/(1-3*x-5*x^2).
a(n) = 3*a(n-1) + 5*a(n-2) for n>=3  with a(0)=1, a(1)=4 and a(2)=18.
a(n) = ((22+54*A)*A^(-n-1) + (22+54*B)*B^(-n-1))/145  with A=(-3+sqrt(29))/10 and B=(-3-sqrt(29))/10 for n>=1 with a(0)=1.
5*a(n) = 2*( A015523(n) + 3*A015523(n+1)), n>0 - R. J. Mathar, May 11 2013

A201730 Triangle T(n,k), read by rows, given by (2,1/2,3/2,0,0,0,0,0,0,0,...) DELTA (0,1/2,-1/2,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 2, 0, 5, 1, 0, 14, 6, 0, 0, 41, 26, 1, 0, 0, 122, 100, 10, 0, 0, 0, 365, 363, 63, 1, 0, 0, 0, 1094, 1274, 322, 14, 0, 0, 0, 0, 3281, 4372, 1462, 116, 1, 0, 0, 0, 0, 9842, 14760, 6156, 744, 18, 0, 0, 0, 0, 0

Offset: 0

Philippe Deléham, Dec 04 2011

Riordan array ((1-2x)/(1-4x+3x^2),x^2/(1-4x+3x^2)).

A007318*A201701 as lower triangular matrices.

			Triangle begins:
1
2, 0
5, 1, 0
14, 6, 0, 0
41, 26, 1, 0, 0
122, 100, 10, 0, 0, 0
365, 363, 63, 1, 0, 0, 0

Cf. A007051 (1st column), A261064 (2nd column).

A201730 := proc(n,k)
    (1-2*x)/(1-4*x+(3-y)*x^2) ;
    coeftayl(%,y=0,k) ;
    coeftayl(%,x=0,n) ;
end proc:
seq(seq(A201730(n,k),k=0..n),n=0..12) ; # R. J. Mathar, Dec 06 2011

m = 13;
(* DELTA is defined in A084938 *)
DELTA[Join[{2, 1/2, 3/2}, Table[0, {m}]], Join[{0, 1/2, -1/2}, Table[0, {m}]], m] // Flatten (* Jean-François Alcover, Feb 19 2020 *)

G.f.: (1-2x)/(1-4x+(3-y)*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = A139011(n), A000079(n), A007051(n), A006012(n), A001075(n), A081294(n), A001077(n), A084059(n), A108851(n), A084128(n), A081340(n), A084132(n) for x = -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively.
Sum_{k, k>+0} T(n+k,k) = A081704(n) .
T(n,k) = 3*T(n-1,k)+ Sum_{j>0} T(n-1-j,k-1).
T(n,k) = 4*T(n-1,k)+ T(n-2,k-1) - 3*T(n-2,k) with T(0,0)=1, T(1,0)= 2, T(1,1) = 0 and T(n,k) = 0 if k<0 or if n

A048875 Generalized Pellian with second term of 6.

Original entry on oeis.org

1, 6, 25, 106, 449, 1902, 8057, 34130, 144577, 612438, 2594329, 10989754, 46553345, 197203134, 835365881, 3538666658, 14990032513, 63498796710, 268985219353, 1139439674122, 4826743915841, 20446415337486, 86612405265785, 366896036400626, 1554196550868289

Offset: 0

Barry E. Williams

			G.f. = 1 + 6*x + 25*x^2 + 106*x^3 + 449*x^4 + 1902*x^5 + 8057*x^6 + 34130*x^7 + ...

T. D. Noe, Table of n, a(n) for n = 0..200
M. Bicknell, A Primer on the Pell Sequence and related sequences, Fibonacci Quarterly, Vol. 13, No. 4, 1975, pp. 345-349.
L. Carlitz, R. Scoville and V. E. Hoggatt, Jr., Pellian Representations, Fib. Quart. Vol. 10, No. 5, (1972), pp. 449-488.
Tanya Khovanova, Recursive Sequences
A. K. Whitford, Binet's Formula Generalized, Fibonacci Quarterly, Vol. 15, No. 1, 1979, pp. 21, 24, 29.
Index entries for linear recurrences with constant coefficients, signature (4,1).

Cf. A015448, A001076, A001077, A033887, A134418, A097924, A374439.

with(combinat): a:=n->2*fibonacci(n-1,4)+fibonacci(n,4): seq(a(n), n=1..17); # Zerinvary Lajos, Apr 04 2008

LinearRecurrence[{4,1},{1,6},40] (* Harvey P. Dale, Nov 30 2011 *)
a[ n_] := (4 I ChebyshevT[ n + 1, -2 I] - 3 ChebyshevT[ n, -2 I]) I^n / 5; (* Michael Somos, Feb 23 2014 *)
a[ n_] := If[ n < 0, SeriesCoefficient[ (1 + 6 x) / (1 + 4 x - x^2), {x, 0, -n}], SeriesCoefficient[ (1 + 2 x) / (1 - 4 x - x^2), {x, 0, n}]]; (* Michael Somos, Feb 23 2014 *)

a[0]:1$ a[1]:6$ a[n]:=4*a[n-1]+a[n-2]$ makelist(a[n], n, 0, 30); /* Martin Ettl, Nov 03 2012 */

{a(n) = ( 4*I*polchebyshev( n+1, 1, -2*I) - 3*polchebyshev( n, 1, -2*I)) * I^n / 5}; /* Michael Somos, Feb 23 2014 */

{a(n) = if( n<0, polcoeff( (1 + 6*x) / (1 + 4*x - x^2) + x * O(x^-n), -n), polcoeff( (1 + 2*x) / (1 - 4*x - x^2) + x * O(x^n), n))}; \\ Michael Somos, Feb 23 2014

a(n) = ((4+sqrt(5))*(2+sqrt(5))^n - (4-sqrt(5))*(2-sqrt(5))^n)*sqrt(5)/2.
a(n) = 4a(n-1) + a(n-2); a(0)=1, a(1)=6.
Binomial transform of A134418: (1, 5, 14, 48, 152, ...). - Gary W. Adamson, Nov 23 2007
G.f.: (1+2*x)/(1-4*x-x^2). - Philippe Deléham, Nov 03 2008
a(-1 - n) = (-1)^n * A097924(n) for all n in Z. - Michael Somos, Feb 23 2014
a(n) = A001076(n+1) + 2*A001076(n). - R. J. Mathar, Sep 11 2019
a(n) = 4^n*Sum_{k=0..n} A374439(n, k)*(1/4)^k. - Peter Luschny, Jul 26 2024

Corrected by T. D. Noe, Nov 07 2006

A075796 Numbers k such that 5*k^2 + 5 is a square.

Original entry on oeis.org

2, 38, 682, 12238, 219602, 3940598, 70711162, 1268860318, 22768774562, 408569081798, 7331474697802, 131557975478638, 2360712083917682, 42361259535039638, 760141959546795802, 13640194012307284798, 244763350261984330562, 4392100110703410665318, 78813038642399407645162

Offset: 1

Gregory V. Richardson, Oct 13 2002

Bisection of A001077; a(n) = A001077(2*n-1). - Greg Dresden, Jun 08 2021
From Peter Bala, Aug 25 2022: (Start)
The aerated sequence (b(n))n>=1 = [2, 0, 38, 0, 682, 0, 1238, 0, ...] is a fourth-order linear divisibility sequence; that is, if n | m then b(n) | b(m). The sequence (1/2)*(b(n))n>=1 is the case P1 = 0, P2 = -16, Q = -1  of the 3-parameter family of divisibility sequences found by Williams and Guy. See A100047. (End)

A. H. Beiler, "The Pellian." Ch. 22 in Recreations in the Theory of Numbers: The Queen of Mathematics Entertains. Dover, New York, New York, pp. 248-268, 1966.
L. E. Dickson, History of the Theory of Numbers, Vol. II, Diophantine Analysis. AMS Chelsea Publishing, Providence, Rhode Island, 1999, pp. 341-400.
Peter G. L. Dirichlet, Lectures on Number Theory (History of Mathematics Source Series, V. 16); American Mathematical Society, Providence, Rhode Island, 1999, pp. 139-147.

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Tanya Khovanova, Recursive Sequences
J. J. O'Connor and E. F. Robertson, Pell's Equation
Eric Weisstein's World of Mathematics, Pell Equation.
H. C. Williams and R. K. Guy, Some fourth-order linear divisibility sequences, Intl. J. Number Theory 7 (5) (2011) 1255-1277.
Index entries for linear recurrences with constant coefficients, signature (18,-1).

Cf. A000290, A306380, A001077.

I:=[2,38]; [n le 2 select I[n] else 18*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Nov 30 2011

[Lucas(6*n-3)/2: n in [1..20]]; // G. C. Greubel, Feb 13 2019

with(combinat); A075796:=n->fibonacci(6*n+3)+fibonacci(6*n)/2; seq(A075796(n), n=1..50); # Wesley Ivan Hurt, Nov 29 2013

LinearRecurrence[{18, -1}, {2, 38}, 50] (* Sture Sjöstedt, Nov 29 2011; typo fixed by Vincenzo Librandi, Nov 30 2011 *)
LucasL[6*Range[20]-3]/2 (* G. C. Greubel, Feb 13 2019 *)
CoefficientList[Series[2*(1+x)/( 1-18*x+x^2 ), {x,0,20}],x] (* Stefano Spezia, Mar 02 2019 *)

vector(20, n, (fibonacci(6*n-2) + fibonacci(6*n-4))/2) \\ G. C. Greubel, Feb 13 2019

[(fibonacci(6*n-2) + fibonacci(6*n-4))/2 for n in (1..20)] # G. C. Greubel, Feb 13 2019

a(n) = (((9 + 4*sqrt(5))^n - (9 - 4*sqrt(5))^n) + ((9 + 4*sqrt(5))^(n-1) - (9 - 4*sqrt(5))^(n-1)))/(4*sqrt(5)).
a(n) = 18*a(n-1) - a(n-2).
a(n) = 2*A049629(n-1).
Limit_{n->oo} a(n)/a(n-1) = 8*phi + 1 = 9 + 4*sqrt(5).
a(n+1) = 9*a(n) + 4*sqrt(5)*sqrt((a(n)^2+1)). - Richard Choulet, Aug 30 2007
G.f.: 2*x*(1 + x)/(1 - 18*x + x^2). - Richard Choulet, Oct 09 2007
From Johannes W. Meijer, Jul 01 2010: (Start)
a(n) = A000045(6*n+3) + A000045(6*n)/2.
a(n) = 2*A167808(6*n+4) - A167808(6*n+6).
Limit_{k->oo} a(n+k)/a(k) = A023039(n)*A060645(n)*sqrt(5).
(End)
5*A007805(n)^2 - 1 = a(n+1)^2. - Sture Sjöstedt, Nov 29 2011
From Peter Bala, Nov 29 2013: (Start)
a(n) = Lucas(6*n - 3)/2.
Sum_{n >= 1} 1/(a(n) + 5/a(n)) = 1/4. Compare with A002878, A005248, A023039. (End)
Limit_{n->oo} a(n)/A007805(n-1) = sqrt(5). - A.H.M. Smeets, May 29 2017
E.g.f.: (exp((9 - 4*sqrt(5))*x)*(- 5 + 2*sqrt(5) + (5 + 2*sqrt(5))*exp(8*sqrt(5)*x)))/(2*sqrt(5)). - Stefano Spezia, Feb 13 2019
Sum_{n > 0} 1/a(n) = (1/log(9 - 4*sqrt(5)))*(- 17 - 38/sqrt(5))*sqrt(5*(9 - 4*sqrt(5)))*(- 9 + 4*sqrt(5))*(psi_{9 - 4*sqrt(5)}(1/2) - psi_{9 - 4*sqrt(5)}(1/2 - (I*Pi)/log(9 - 4*sqrt(5)))) approximately equal to 0.527868600269500798938265500122302016..., where psi_q(x) is the q-digamma function. - Stefano Spezia, Feb 25 2019
a(n) = sinh((6*n - 3)*arccsch(2)). - Peter Luschny, May 25 2022

A041008 Numerators of continued fraction convergents to sqrt(7).

Original entry on oeis.org

2, 3, 5, 8, 37, 45, 82, 127, 590, 717, 1307, 2024, 9403, 11427, 20830, 32257, 149858, 182115, 331973, 514088, 2388325, 2902413, 5290738, 8193151, 38063342, 46256493, 84319835, 130576328, 606625147, 737201475, 1343826622, 2081028097, 9667939010, 11748967107, 21416906117

Offset: 0

N. J. A. Sloane

Vincenzo Librandi, Table of n, a(n) for n = 0..200
C. Elsner, Series of Error Terms for Rational Approximations of Irrational Numbers, J. Int. Seq. 14 (2011) # 11.1.4.
C. Elsner, M. Stein, On Error Sum Functions Formed by Convergents of Real Numbers, J. Int. Seq. 14 (2011) # 11.8.6.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,16,0,0,0,-1).

Cf. A010465, A041009 (denominators), A266698 (quadrisection), A001081 (quadrisection).

Analog for other sqrt(m): A001333 (m=2), A002531 (m=3), A001077 (m=5), A041006 (m=6), A041010 (m=8), A005667 (m=10), A041014 (m=11), A041016 (m=12), ..., A042934 (m=999), A042936 (m=1000).

Table[Numerator[FromContinuedFraction[ContinuedFraction[Sqrt[7],n]]],{n,1,50}] (* Vladimir Joseph Stephan Orlovsky, Mar 16 2011 *)
Numerator[Convergents[Sqrt[7], 30]] (* Vincenzo Librandi, Oct 28 2013 *)
LinearRecurrence[{0,0,0,16,0,0,0,-1},{2,3,5,8,37,45,82,127},40] (* Harvey P. Dale, Jul 23 2021 *)

A041008=contfracpnqn(c=contfrac(sqrt(7)),#c)[1,][^-1] \\ Discard possibly incorrect last element. NB: a(n)=A041008[n+1]! For more terms use:
A041008(n)={n<#A041008|| A041008=extend(A041008, [4, 16; 8, -1], n\.8); A041008[n+1]}
extend(A,c,N)={for(n=#A+1, #A=Vec(A, N), A[n]=[A[n-i]|i<-c[,1]]*c[,2]); A} \\ (End)

G.f.: (2 + 3*x + 5*x^2 + 8*x^3 + 5*x^4 - 3*x^5 + 2*x^6 - x^7)/(1 - 16*x^4 + x^8).

A081459 Consider the mapping f(r) = (1/2)*(r + N/r) from rationals to rationals where N = 5. Starting with r = 2 and applying the mapping to each new (reduced) rational number gives 2, 9/4, 161/72, 51841/23184, ..., tending to N^(1/2). Sequence gives values of the numerators.

Original entry on oeis.org

2, 9, 161, 51841, 5374978561, 57780789062419261441, 6677239169351578707225356193679818792961, 89171045849445921581733341920411050611581102638589828325078491812417901966295041

Offset: 1

Amarnath Murthy, Mar 22 2003

Related sequence pairs (numerator, denominator) can be obtained by choosing N = 2, 3, 6 etc.

The sequence satisfies the Pell equation a(n+1)^2 - 5*A081460(n+1)^2 = 1. - Vincenzo Librandi, Dec 20 2011

Vincenzo Librandi, Table of n, a(n) for n = 1..11

Cf. A000032, A000129, A001333, A081460.

m:=8; f:=[ n eq 1 select 2 else (Self(n-1)+5/Self(n-1))/2: n in [1..m] ]; [ Numerator(f[n]): n in [1..m] ]; // Bruno Berselli, Dec 20 2011

k = 4; Table[Simplify[Expand[(1/2) (((k + Sqrt[k^2 + 4])/2)^(2^(n - 1)) + ((k - Sqrt[k^2 + 4])/2)^(2^(n - 1)))]], {n, 1, 6}] (* Artur Jasinski, Oct 12 2008 *)
aa = {}; k = 9; Do[AppendTo[aa, k]; k = 2 k^2 - 1, {n, 1, 5}]; aa (* Artur Jasinski, Oct 12 2008 *)

{r=2; N=5; for(n=1,8,a=numerator(r); b=denominator(r); print1(a,","); r=(1/2)*(r + N/r) )}

a(n) = a(n-1)^2 + 5*A081460(n-1)^2. - Mario Catalani (mario.catalani(AT)unito.it), May 21 2003
a(n) = (1/2)*(((4+2*sqrt(5))/2)^(2^(n-1)) + ((4-2*sqrt(5))/2)^(2^(n-1))). a(n+1) = 2*a(n)^2 - 1 for n > 1. - Artur Jasinski, Oct 12 2008
a(n) = A000032(3*2^(n-1))/2. - Ehren Metcalfe, Oct 05 2017
a(n) = A001077(2^(n-1)). - Robert FERREOL, Apr 16 2023
From  Peter Bala, Jun 22 2025: (Start)
Product_{n >= 1} (1 + 1/a(n)) = (3/4)*sqrt(5).
Product_{n >= 1} (1 - 1/(2*a(n))) = (6/19)*sqrt(5). See A002812. (End)

Edited and extended by Klaus Brockhaus and Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Apr 06 2003

A098648 Expansion of (1-3x)/(1 - 6x + 4*x^2).

Original entry on oeis.org

1, 3, 14, 72, 376, 1968, 10304, 53952, 282496, 1479168, 7745024, 40553472, 212340736, 1111830528, 5821620224, 30482399232, 159607914496, 835717890048, 4375875682304, 22912382533632, 119970792472576, 628175224700928, 3289168178315264, 17222308171087872

Offset: 0

Paul Barry, Sep 18 2004

Binomial transform of A001077. Second binomial transform of A084057. Third binomial transform of 1/(1-5*x^2). Let A=[1,1,1,1;3,1,-1,-3;3,-1,-1,3;1,-1,1,-1], the 4 X 4 Krawtchouk matrix. Then a(n)=trace((16(A*A`)^(-1))^n)/4.

Vincenzo Librandi, Table of n, a(n) for n = 0..300
Index entries for linear recurrences with constant coefficients, signature (6,-4).

Cf. A098647.

a[n_]:=(MatrixPower[{{5,1},{1,1}},n].{{2},{1}})[[2,1]]; Table[a[n],{n,0,40}] (* Vladimir Joseph Stephan Orlovsky, Feb 20 2010 *)
CoefficientList[Series[(1-3x)/(1-6x+4x^2),{x,0,30}],x] (* or *) LinearRecurrence[{6,-4},{1,3},31] (* Harvey P. Dale, Jun 06 2011 *)
Table[2^(n - 1) LucasL[2 n], {n, 0, 20}] (* Eric W. Weisstein, Mar 31 2017 *)

Vec((1-3*x)/(1 - 6*x + 4*x^2) + O(x^25)) \\ Jinyuan Wang, Jul 24 2021

E.g.f.: exp(3*x)*cosh(sqrt(5)*x).
a(n) = ((3-sqrt(5))^n + (3+sqrt(5))^n)/2.
a(n) = 2*(3*a(n-1) - 2*a(n-2)). - Lekraj Beedassy, Oct 22 2004
a(n) = A084326(n+1) - 3*A084326(n). - R. J. Mathar, Nov 10 2013
a(n) = 2^(n-1)*Lucas(2*n) = 2^(n-1)*A005248(n), n>0. - Eric W. Weisstein, Mar 31 2017

A041014 Numerators of continued fraction convergents to sqrt(11).

Original entry on oeis.org

3, 10, 63, 199, 1257, 3970, 25077, 79201, 500283, 1580050, 9980583, 31521799, 199111377, 628855930, 3972246957, 12545596801, 79245827763, 250283080090, 1580944308303, 4993116004999, 31539640338297

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,20,0,-1).

Cf. A010468, A041015 (denominators).

Analog for other sqrt(m): A001333 (m=2), A002531 (m=3), A001077 (m=5), A041006 (m=6), A041008 (m=7), A041010 (m=8), A005667 (m=10), A041016 (m=12), ..., A042936 (m=1000).

Table[Numerator[FromContinuedFraction[ContinuedFraction[Sqrt[11],n]]],{n,1,50}] (* Vladimir Joseph Stephan Orlovsky, Mar 16 2011 *)
Numerator[Convergents[Sqrt[11], 30]] (* Vincenzo Librandi, Oct 28 2013 *)

A041014=contfracpnqn(c=contfrac(sqrt(11)), #c)[1,][^-1] \\ Discard last element which may be incorrect. Use e.g. \p999 to get more terms, or extend as follows:
{A041014_upto(N,A=Vec(A041014,N))=for(n=#A041014+1,N, A[n]=20*A[n-2]-A[n-4]); A041014=A} \\ M. F. Hasler, Nov 01 2019

G.f.: (3 + 10*x + 3*x^2 - x^3)/(1 - 20*x^2 + x^4).

Previous Showing 11-20 of 53 results. Next

cp's OEIS Frontend

A041006 Numerators of continued fraction convergents to sqrt(6).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

Extensions

A048876 a(n) = 4*a(n-1) + a(n-2); a(0)=1, a(1)=7.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A180140 Eight rooks and one berserker on a 3 X 3 chessboard. G.f.: (1+x+x^2)/(1-3*x-5*x^2).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A201730 Triangle T(n,k), read by rows, given by (2,1/2,3/2,0,0,0,0,0,0,0,...) DELTA (0,1/2,-1/2,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

Formula

A048875 Generalized Pellian with second term of 6.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

PARI

Formula

Extensions

A075796 Numbers k such that 5*k^2 + 5 is a square.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Magma

Maple

Mathematica

PARI

A180140 Eight rooks and one berserker on a 3 X 3 chessboard. G.f.: (1+x+x^2)/(1-3x-5x^2).

A098648 Expansion of (1-3x)/(1 - 6x + 4*x^2).