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 28 results. Next

A298212 Smallest n such that A060645(a(n)) = 0 (mod n), i.e., x=A023039(a(n)) and y=A060645(a(n)) is the fundamental solution of the Pell equation x^2 - 5*(n*y)^2 = 1.

Original entry on oeis.org

1, 1, 2, 1, 5, 2, 4, 2, 2, 5, 5, 2, 7, 4, 10, 4, 3, 2, 3, 5, 4, 5, 4, 2, 25, 7, 6, 4, 7, 10, 5, 8, 10, 3, 20, 2, 19, 3, 14, 10, 10, 4, 22, 5, 10, 4, 8, 4, 28, 25, 6, 7, 9, 6, 5, 4, 6, 7, 29, 10, 5, 5, 4, 16, 35, 10, 34, 3, 4, 20, 35, 2, 37, 19, 50
Offset: 1

Views

Author

A.H.M. Smeets, Jan 15 2018

Keywords

Comments

The fundamental solution of the Pell equation x^2 - 5*(n*y)^2 = 1, is the smallest solution of x^2 - 5*y^2 = 1 satisfying y = 0 (mod n).

References

  • Michael J. Jacobson, Jr. and Hugh C. Williams, Solving the Pell Equation, Springer, 2009, pages 1-17.

Crossrefs

Programs

  • Mathematica
    b[n_] := b[n] = Switch[n, 0, 0, 1, 4, _, 18 b[n - 1] - b[n - 2]];
    a[n_] := For[k = 1, True, k++, If[Mod[b[k], n] == 0, Return[k]]];
    a /@ Range[100] (* Jean-François Alcover, Nov 16 2019 *)
  • Python
    xf, yf = 9, 4
    x, n = 2*xf, 0
    while n < 20000:
        n = n+1
        y1, y0, i = 0, yf, 1
        while y0%n != 0:
            y1, y0, i = y0, x*y0-y1, i+1
        print(n, i)

Formula

a(n) <= n.
a(A000351(n)) = A000351(n).
A023039(a(n)) = A002350(5*n^2).
A060645(a(n)) = A002349(5*n^2).
if n | m then a(n) | a(m).
a(5^m) = 5^m for m>=0.
In general: if p is prime and p = 1 (mod 4) then: a(n) = n iff n = p^m, for m>=0.

A005248 Bisection of Lucas numbers: a(n) = L(2*n) = A000032(2*n).

Original entry on oeis.org

2, 3, 7, 18, 47, 123, 322, 843, 2207, 5778, 15127, 39603, 103682, 271443, 710647, 1860498, 4870847, 12752043, 33385282, 87403803, 228826127, 599074578, 1568397607, 4106118243, 10749957122, 28143753123, 73681302247, 192900153618, 505019158607, 1322157322203
Offset: 0

Views

Author

Keywords

Comments

Drop initial 2; then iterates of A050411 do not diverge for these starting values. - David W. Wilson
All nonnegative integer solutions of Pell equation a(n)^2 - 5*b(n)^2 = +4 together with b(n)=A001906(n), n>=0. - Wolfdieter Lang, Aug 31 2004
a(n+1) = B^(n)AB(1), n>=0, with compositions of Wythoff's complementary A(n):=A000201(n) and B(n)=A001950(n) sequences. See the W. Lang link under A135817 for the Wythoff representation of numbers (with A as 1 and B as 0 and the argument 1 omitted). E.g., 3=`10`, 7=`010`, 18=`0010`, 47=`00010`, ..., in Wythoff code. a(0) = 2 = B(1) in Wythoff code.
Output of Tesler's formula (as well as that of Lu and Wu) for the number of perfect matchings of an m X n Möbius band where m and n are both even specializes to this sequence for m=2. - Sarah-Marie Belcastro, Jul 04 2009
Numbers having two 1's in their base-phi representation. - Robert G. Wilson v, Sep 13 2010
Pisano period lengths: 1, 3, 4, 3, 2, 12, 8, 6, 12, 6, 5, 12, 14, 24, 4, 12, 18, 12, 9, 6, ... - R. J. Mathar, Aug 10 2012
From Wolfdieter Lang, Feb 18 2013: (Start)
a(n) is also one half of the total number of round trips, each of length 2*n, on the graph P_4 (o-o-o-o) (the simple path with 4 points (vertices) and 3 lines (or edges)). See the array and triangle A198632 for the general case of the graph P_N (there N is n and the length is l=2*k).
O.g.f. for w(4,l) (with zeros for odd l): y*(d/dy)S(4,y)/S(4,y) with y=1/x and Chebyshev S-polynomials (coefficients A049310). See also A198632 for a rewritten form. One half of this o.g.f. for x -> sqrt(x) produces the g.f. (2-3x)/(1-3x+x^2) given below. (End)
Solutions (x, y) = (a(n), a(n+1)) satisfying x^2 + y^2 = 3xy - 5. - Michel Lagneau, Feb 01 2014
Except for the first term, positive values of x (or y) satisfying x^2 - 7xy + y^2 + 45 = 0. - Colin Barker, Feb 16 2014
Except for the first term, positive values of x (or y) satisfying x^2 - 18xy + y^2 + 320 = 0. - Colin Barker, Feb 16 2014
a(n) are the numbers such that a(n)^2-2 are Lucas numbers. - Michel Lagneau, Jul 22 2014
All sequences of this form, b(n+1) = 3*b(n) - b(n-1), regardless of initial values, which includes this sequence, yield this sequence as follows: a(n) = (b(j+n) + b(j-n))/b(j), for any j, except where b(j) = 0. Also note formula below relating this a(n) to all sequences of the form G(n+1) = G(n) + G(n-1). - Richard R. Forberg, Nov 18 2014
A non-simple continued fraction expansion for F(2n*(k+1))/F(2nk) k>=1 is a(n) + (-1)/(a(n) + (-1)/(a(n) + ... + (-1)/a(n))) where a(n) appears exactly k times (F(n) denotes the n-th Fibonacci number). E.g., F(16)/F(12) equals 7 + (-1)/(7 + (-1)/7). Furthermore, these a(n) are exactly the positive integers k such that the non-simple infinite continued fraction k + (-1)/(k + (-1)/(k + (-1)/(k + ...))) belongs to Q(sqrt(5)). Compare to Benoit Cloitre and Thomas Baruchel's comments at A002878. - Greg Dresden, Aug 13 2019
For n >= 1, a(n) is the number of cyclic up-down words of length 2*n over an alphabet of size 3. - Sela Fried, Apr 08 2025

Examples

			G.f. = 2 + 3*x + 7*x^2 + 18*x^3 + 47*x^4 + 123*x^5 + 322*x^6 + 843*x^7 + ... - _Michael Somos_, Aug 11 2009
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Richard P. Stanley, Enumerative combinatorics, Vol. 2. Volume 62 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 1999.

Crossrefs

Cf. A000032, A002878 (odd-indexed Lucas numbers), A001906 (Chebyshev S(n-1, 3)), a(n) = sqrt(4+5*A001906(n)^2), A228842.
a(n) = A005592(n)+1 = A004146(n)+2 = A065034(n)-1.
First differences of A002878. Pairwise sums of A001519. First row of array A103997.
Cf. A153415, A201157. Also Lucas(k*n): A000032 (k = 1), A014448 (k = 3), A056854 (k = 4), A001946 (k = 5), A087215 (k = 6), A087281 (k = 7), A087265 (k = 8), A087287 (k = 9), A065705 (k = 10), A089772 (k = 11), A089775 (k = 12).

Programs

  • Haskell
    a005248 n = a005248_list !! n
    a005248_list = zipWith (+) (tail a001519_list) a001519_list
    -- Reinhard Zumkeller, Jan 11 2012
  • Magma
    [Lucas(2*n) : n in [0..100]]; // Vincenzo Librandi, Apr 14 2011
    
  • Maple
    a:= n-> (<<2|3>>. <<3|1>, <-1|0>>^n)[1$2]: seq(a(n), n=0..30); # Alois P. Heinz, Jul 31 2008
    with(combinat): seq(5*fibonacci(n)^2+2*(-1)^n, n= 0..26);
  • Mathematica
    a[0] = 2; a[1] = 3; a[n_] := 3a[n - 1] - a[n - 2]; Table[ a[n], {n, 0, 27}] (* Robert G. Wilson v, Jan 30 2004 *)
    Fibonacci[1 + 2n] + 1/2 (-Fibonacci[2n] + LucasL[2n]) (* Tesler. Sarah-Marie Belcastro, Jul 04 2009 *)
    LinearRecurrence[{3, -1}, {2, 3}, 50] (* Sture Sjöstedt, Nov 27 2011 *)
    LucasL[Range[0,60,2]] (* Harvey P. Dale, Sep 30 2014 *)
  • PARI
    {a(n) = fibonacci(2*n + 1) + fibonacci(2*n - 1)}; /* Michael Somos, Jun 23 2002 */
    
  • PARI
    {a(n) = 2 * subst( poltchebi(n), x, 3/2)}; /* Michael Somos, Jun 28 2003 */
    
  • Sage
    [lucas_number2(n,3,1) for n in range(37)] # Zerinvary Lajos, Jun 25 2008
    

Formula

a(n) = Fibonacci(2*n-1) + Fibonacci(2*n+1).
G.f.: (2-3*x)/(1-3*x+x^2). - Simon Plouffe in his 1992 dissertation.
a(n) = S(n, 3) - S(n-2, 3) = 2*T(n, 3/2) with S(n-1, 3) = A001906(n) and S(-2, x) = -1. U(n, x)=S(n, 2*x) and T(n, x) are Chebyshev's U- and T-polynomials.
a(n) = a(k)*a(n - k) - a(n - 2k) for all k, i.e., a(n) = 2*a(n) - a(n) = 3*a(n - 1) - a(n - 2) = 7*a(n - 2) - a(n - 4) = 18*a(n - 3) - a(n - 6) = 47*a(n - 4) - a(n - 8) etc., a(2n) = a(n)^2 - 2. - Henry Bottomley, May 08 2001
a(n) = A060924(n-1, 0) = 3*A001906(n) - 2*A001906(n-1), n >= 1. - Wolfdieter Lang, Apr 26 2001
a(n) ~ phi^(2*n) where phi=(1+sqrt(5))/2. - Joe Keane (jgk(AT)jgk.org), May 15 2002
a(0)=2, a(1)=3, a(n) = 3*a(n-1) - a(n-2) = a(-n). - Michael Somos, Jun 28 2003
a(n) = phi^(2*n) + phi^(-2*n) where phi=(sqrt(5)+1)/2, the golden ratio. E.g., a(4)=47 because phi^(8) + phi^(-8) = 47. - Dennis P. Walsh, Jul 24 2003
With interpolated zeros, trace(A^n)/4, where A is the adjacency matrix of path graph P_4. Binomial transform is then A049680. - Paul Barry, Apr 24 2004
a(n) = (floor((3+sqrt(5))^n) + 1)/2^n. - Lekraj Beedassy, Oct 22 2004
a(n) = ((3-sqrt(5))^n + (3+sqrt(5))^n)/2^n (Note: substituting the number 1 for 3 in the last equation gives A000204, substituting 5 for 3 gives A020876). - Creighton Dement, Apr 19 2005
a(n) = (1/(n+1/2))*Sum_{k=0..n} B(2k)*L(2n+1-2k)*binomial(2n+1, 2k) where B(2k) is the (2k)-th Bernoulli number. - Benoit Cloitre, Nov 02 2005
a(n) = term (1,1) in the 1 X 2 matrix [2,3] . [3,1; -1,0]^n. - Alois P. Heinz, Jul 31 2008
a(n) = 2*cosh(2*n*psi), where psi=log((1+sqrt(5))/2). - Al Hakanson, Mar 21 2009
From Sarah-Marie Belcastro, Jul 04 2009: (Start)
a(n) - (a(n) - F(2n))/2 - F(2n+1) = 0. (Tesler)
Product_{r=1..n} (1 + 4*(sin((4r-1)*Pi/(4n)))^2). (Lu/Wu) (End)
a(n) = Fibonacci(2n+6) mod Fibonacci(2n+2), n > 1. - Gary Detlefs, Nov 22 2010
a(n) = 5*Fibonacci(n)^2 + 2*(-1)^n. - Gary Detlefs, Nov 22 2010
a(n) = A033888(n)/A001906(n), n > 0. - Gary Detlefs, Dec 26 2010
a(n) = 2^(2*n) * Sum_{k=1..2} (cos(k*Pi/5))^(2*n). - L. Edson Jeffery, Jan 21 2012
From Peter Bala, Jan 04 2013: (Start)
Let F(x) = Product_{n>=0} (1 + x^(4*n+1))/(1 + x^(4*n+3)). Let alpha = 1/2*(3 - sqrt(5)). This sequence gives the simple continued fraction expansion of 1 + F(alpha) = 2.31829 56058 81914 31334 ... = 2 + 1/(3 + 1/(7 + 1/(18 + ...))).
Also F(-alpha) = 0.64985 97768 07374 32950 has the continued fraction representation 1 - 1/(3 - 1/(7 - 1/(18 - ...))) and the simple continued fraction expansion 1/(1 + 1/((3-2) + 1/(1 + 1/((7-2) + 1/(1 + 1/((18-2) + 1/(1 + ...))))))).
F(alpha)*F(-alpha) has the simple continued fraction expansion 1/(1 + 1/((3^2-4) + 1/(1 + 1/((7^2-4) + 1/(1 + 1/((18^2-4) + 1/(1 + ...))))))).
Added Oct 13 2019: 1/2 + 1/2*F(alpha)/F(-alpha) = 1.5142923542... has the simple continued fraction expansion 1 + 1/((3 - 2) + 1/(1 + 1/((18 - 2) + 1/(1 + 1/(123 - 2) + 1/(1 + ...))))). (End)
G.f.: (W(0)+6)/(5*x), where W(k) = 5*x*k + x - 6 + 6*x*(5*k-9)/W(k+1) (continued fraction). - Sergei N. Gladkovskii, Aug 19 2013
Sum_{n >= 1} 1/( a(n) - 5/a(n) ) = 1. Compare with A001906, A002878 and A023039. - Peter Bala, Nov 29 2013
0 = a(n) * a(n+2) - a(n+1)^2 - 5 for all n in Z. - Michael Somos, Aug 24 2014
a(n) = (G(j+2n) + G(j-2n))/G(j), for n >= 0 and any j, positive or negative, except where G(j) = 0, and for any sequence of the form G(n+1) = G(n) + G(n-1) with any initial values for G(0), G(1), including non-integer values. G(n) includes Lucas, Fibonacci. Compare with A081067 for odd number offsets from j. - Richard R. Forberg, Nov 16 2014
a(n) = [x^n] ( (1 + 3*x + sqrt(1 + 6*x + 5*x^2))/2 )^n for n >= 1. - Peter Bala, Jun 23 2015
From J. M. Bergot, Oct 28 2015: (Start)
For n>0, a(n) = F(n-1) * L(n) + F(2*n+1) - (-1)^n with F(k) = A000045(k).
For n>1, a(n) = F(n+1) * L(n) + F(2*n-1) - (-1)^n.
For n>2, a(n) = 5*F(2*n-3) + 2*L(n-3) * L(n) + 8*(-1)^n. (End)
For n>1, a(n) = L(n-2)*L(n+2) -7*(-1)^n. - J. M. Bergot, Feb 10 2016
a(n) = 6*F(n-1)*L(n-1) - F(2*n-6) with F(n)=A000045(n) and L(n)=A000032(n). - J. M. Bergot, Apr 21 2017
a(n) = F(2*n) + 2*F(n-1)*L(n) with F(n)=A000045(n) and L(n)=A000032(n). - J. M. Bergot, May 01 2017
E.g.f.: exp(4*x/(1+sqrt(5))^2) + exp((1/4)*(1+sqrt(5))^2*x). - Stefano Spezia, Aug 13 2019
From Peter Bala, Oct 14 2019: (Start)
a(n) = F(2*n+2) - F(2*n-2) = A001906(n+1) - A001906(n-1).
a(n) = trace(M^n), where M is the 2 X 2 matrix [0, 1; 1, 1]^2 = [1, 1; 1, 2].
Consequently the Gauss congruences hold: a(n*p^k) = a(n*p^(k-1)) ( mod p^k ) for all prime p and positive integers n and k. See Zarelua and also Stanley (Ch. 5, Ex. 5.2(a) and its solution).
Sum_{n >= 1} (-1)^(n+1)/( a(n) + 1/a(n) ) = 1/5.
Sum_{n >= 1} (-1)^(n+1)/( a(n) + 3/(a(n) + 2/(a(n))) ) = 1/6.
Sum_{n >= 1} (-1)^(n+1)/( a(n) + 9/(a(n) + 4/(a(n) + 1/(a(n)))) ) = 1/9.
x*exp(Sum_{n >= 1} a(n)*x^/n) = x + 3*x^2 + 8*x^3 + 21*x^4 + ... is the o.g.f. for A001906. (End)
a(n) = n + 2 + Sum_{k=1..n-1} k*a(n-k). - Yu Xiao, May 30 2020
Sum_{n>=1} 1/a(n) = A153415. - Amiram Eldar, Nov 11 2020
Sum_{n>=0} 1/(a(n) + 3) = (2*sqrt(5) + 1)/10 (André-Jeannin, 1991). - Amiram Eldar, Jan 23 2022
a(n) = 2*cosh(2*n*arccsch(2)) = 2*cosh(2*n*asinh(1/2)). - Peter Luschny, May 25 2022
a(n) = (5/2)*(Sum_{k=-n..n} binomial(2*n, n+5*k)) - (1/2)*4^n. - Greg Dresden, Jan 05 2023
a(n) = Sum_{k>=0} Lucas(2*n*k)/(Lucas(2*n)^(k+1)). - Diego Rattaggi, Jan 12 2025

Extensions

Additional comments from Michael Somos, Jun 23 2001

A059100 a(n) = n^2 + 2.

Original entry on oeis.org

2, 3, 6, 11, 18, 27, 38, 51, 66, 83, 102, 123, 146, 171, 198, 227, 258, 291, 326, 363, 402, 443, 486, 531, 578, 627, 678, 731, 786, 843, 902, 963, 1026, 1091, 1158, 1227, 1298, 1371, 1446, 1523, 1602, 1683, 1766, 1851, 1938, 2027, 2118, 2211, 2306, 2403, 2502, 2603
Offset: 0

Views

Author

Henry Bottomley, Feb 13 2001

Keywords

Comments

Let s(n) = Sum_{k>=1} 1/n^(2^k). Then I conjecture that the maximum element in the continued fraction for s(n) is n^2 + 2. - Benoit Cloitre, Aug 15 2002
Binomial transformation yields A081908, with A081908(0)=1 dropped. - R. J. Mathar, Oct 05 2008
1/a(n) = R(n)/r with R(n) the n-th radius of the Pappus chain of the symmetric arbelos with semicircle radii r, r1 = r/2 = r2. See the MathWorld link for Pappus chain (there are two of them, a left and a right one. In this case these two chains are congruent). - Wolfdieter Lang, Mar 01 2013
a(n) is the number of election results for an election with n+2 candidates, say C1, C2, ..., and C(n+2), and with only two voters (each casting a single vote) that have C1 and C2 receiving the same number of votes. See link below. - Dennis P. Walsh, May 08 2013
This sequence gives the set of values such that for sequences b(k+1) = a(n)*b(k) - b(k-1), with initial values b(0) = 2, b(1) = a(n), all such sequences are invariant under this transformation: b(k) = (b(j+k) + b(j-k))/b(j), except where b(j) = 0, for all integer values of j and k, including negative values. Examples are: at n=0, b(k) = 2 for all k; at n=1, b(k) = A005248; at n=2, b(k) = 2*A001541; at n=3, b(k)= A057076; at n=4, b(k) = 2*A023039. This b(k) family are also the transformation results for all related b'(k) (i.e., those with different initial values) including non-integer values. Further, these b(k) are also the bisections of the transformations of sequences of the form G(k+1) = n * G(k) + G(k-1), and those bisections are invariant for all initial values of g(0) and g(1), including non-integer values. For n = 1 this g(k) family includes Fibonacci and Lucas, where the invariant bisection is b(k) = A005248. The applicable bisection for this transformation of g(k) is for the odd values of k, and applies for all n. Also see A000032 for a related family of sequences. - Richard R. Forberg, Nov 22 2014
Also the number of maximum matchings in the n-gear graph. - Eric W. Weisstein, Dec 31 2017
Also the Wiener index of the n-dipyramidal graph. - Eric W. Weisstein, Jun 14 2018
Numbers of the form n^2+2 have no factors that are congruent to 7 (mod 8). - Gordon E. Michaels, Sep 12 2019
For n >= 1, the continued fraction expansion of sqrt(a(n)) is [n; {n, 2n}]. - Magus K. Chu, Sep 10 2022

Examples

			For n = 2, a(2) = 6 since there are 6 election results in a 4-candidate, 2-voter election that have candidates c1 and c2 tied. Letting <i,j> denote voter 1 voting for candidate i and voter 2 voting for candidate j, the six election results are <1,2>, <2,1>, <3,3>, <3,4>, <4,3>, and <4,4>. - _Dennis P. Walsh_, May 08 2013
		

Crossrefs

Apart from initial terms, same as A010000.
2nd row/column of A295707.

Programs

Formula

G.f.: (2 - 3*x + 3*x^2)/(1 - x)^3. - R. J. Mathar, Oct 05 2008
a(n) = ((n - 2)^2 + 2*(n + 1)^2)/3. - Reinhard Zumkeller, Feb 13 2009
a(n) = A000196(A156798(n) - A000290(n)). - Reinhard Zumkeller, Feb 16 2009
a(n) = 2*n + a(n-1) - 1 with a(0) = 2. - Vincenzo Librandi, Aug 07 2010
a(n+3) = (A166464(n+5) - A166464(n))/20. - Paul Curtz, Nov 07 2012
From Paul Curtz, Nov 07 2012: (Start)
a(3*n) mod 9 = 2.
a(3*n+1) = 3*A056109(n).
a(3*n+2) = 3*A056105(n+1). (End)
Sum_{n >= 1} 1/a(n) = Pi * coth(sqrt(2)*Pi) / 2^(3/2) - 1/4. - Vaclav Kotesovec, May 01 2018
From Amiram Eldar, Jan 29 2021: (Start)
Sum_{n>=0} (-1)^n/a(n) = (1 + sqrt(2)*Pi*(csch(sqrt(2)*Pi)))/4.
Product_{n>=0} (1 + 1/a(n)) = sqrt(3/2)*csch(sqrt(2)*Pi)*sinh(sqrt(3)*Pi).
Product_{n>=0} (1 - 1/a(n)) = csch(sqrt(2)*Pi)*sinh(Pi)/sqrt(2). (End)
E.g.f.: exp(x)*(2 + x + x^2). - Stefano Spezia, Aug 07 2024

A049660 a(n) = Fibonacci(6*n)/8.

Original entry on oeis.org

0, 1, 18, 323, 5796, 104005, 1866294, 33489287, 600940872, 10783446409, 193501094490, 3472236254411, 62306751484908, 1118049290473933, 20062580477045886, 360008399296352015, 6460088606857290384, 115921586524134874897, 2080128468827570457762
Offset: 0

Views

Author

Keywords

Comments

For n >= 2, a(n) equals the permanent of the (n-1) X (n-1) tridiagonal matrix with 18's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
For n >= 2, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,17}. - Milan Janjic, Jan 25 2015
10*a(n)^2 = Tri(4)*S(n-1, 18)^2 is the triangular number Tri((T(n, 9) - 1)/2), with Tri, S and T given in A000217, A049310 and A053120. This is instance k = 4 of the k-family of identities given in a comment on A001109. - Wolfdieter Lang, Feb 01 2016
Possible solutions for y in Pell equation x^2 - 80*y^2 = 1. The values for x are given in A023039. - Herbert Kociemba, Jun 05 2022

Examples

			a(3) = F(6 * 3) / 8 = F(18) / 8 = 2584 / 8 = 323. - _Indranil Ghosh_, Feb 06 2017
		

Crossrefs

Column m=6 of array A028412.
Partial sums of A007805.

Programs

Formula

G.f.: x/(1 - 18*x + x^2).
a(n) = A134492(n)/8.
a(n) ~ (1/40)*sqrt(5)*(sqrt(5) + 2)^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
For all terms k of the sequence, 80*k^2 + 1 is a square. Limit_{n->oo} a(n)/a(n-1) = 8*phi + 5 = 9 + 4*sqrt(5). - Gregory V. Richardson, Oct 14 2002
a(n) = S(n-1, 18) with S(n, x) := U(n, x/2), Chebyshev's polynomials of the second kind. S(-1, x) := 0. See A049310.
a(n) = (((9 + 4*sqrt(5))^n - (9 - 4*sqrt(5))^n))/(8*sqrt(5)).
a(n) = sqrt((A023039(n)^2 - 1)/80) (cf. Richardson comment).
a(n) = 18*a(n-1) - a(n-2). - Gregory V. Richardson, Oct 14 2002
a(n) = A001076(2n)/4.
a(n) = 17*(a(n-1) + a(n-2)) - a(n-3) = 19*(a(n-1) - a(n-2)) + a(n-3). - Mohamed Bouhamida, May 26 2007
a(n+1) = Sum_{k=0..n} A101950(n,k)*17^k. - Philippe Deléham, Feb 10 2012
Product_{n>=1} (1 + 1/a(n)) = (1/2)*(2 + sqrt(5)). - Peter Bala, Dec 23 2012
Product_{n>=2} (1 - 1/a(n)) = (2/9)*(2 + sqrt(5)). - Peter Bala, Dec 23 2012
a(n) = (1/32)*(F(6*n + 3) - F(6*n - 3)).
Sum_{n>=1} 1/(4*a(n) + 1/(4*a(n))) = 1/4. Compare with A001906 and A049670. - Peter Bala, Nov 29 2013
From Peter Bala, Apr 02 2015: (Start)
Sum_{n >= 1} a(n)*x^(2*n) = -G(x)*G(-x), where G(x) = Sum_{n >= 1} A001076(n)*x^n.
1 + 4*Sum_{n >= 1} a(n)*x^(2*n) = (1 + F(x))*(1 + F(-x)) = (1 + 2*x*G(x))*(1 - 2*x*G(-x)), where F(x) = Sum_{n >= 1} Fibonacci(3*n + 3)*x^n.
1 + 7*Sum_{n >= 1} a(n)*x^(2*n) = (1 + G(x))*(1 + G(-x)) = (1 + 7*G(x))*(1 + 7*G(-x)).
1 + 12*Sum_{n >= 1} a(n)*x^(2*n) = (1 + 2*G(x))*(1 + 2*G(-x)) = (1 + 6*G(x))*(1 + 6*G(-x)) = (1 + A(x))*(1 + A(-x)), where A(x) = Sum_{n >= 1} Fibonacci(3*n)*x^n is the o.g.f for A014445.
1 + 15*Sum_{n >= 1} a(n)*x^(2*n) = (1 + 5*G(x))*(1 + 5*G(-x)) = (1 + 3*G(x))*(1 + 3*G(-x)) = H(x)*H(-x), where H(x) = Sum_{n >= 0} A155179(n)*x^n.
1 + 16*Sum_{n >= 1} a(n)*x^(2*n) = (1 + 4*G(x))*(1 + 4*G(-x)) = (1 + 2* Sum_{n >= 1} Fibonacci(3*n - 1)*x^n)*(1 + 2* Sum_{n >= 1} Fibonacci(3*n - 1)*(-x)^n) = (1 + 2* Sum_{n >= 1} Fibonacci(3*n + 1)*x^n)*(1 + 2* Sum_{n >= 1} Fibonacci(3*n + 1)*(-x)^n).
1 + 20*Sum_{n >= 1} a(n)*x^(2*n) = (1 + Sum_{n >= 1} Lucas(3*n)*x^n)*(1 + Sum_{n >= 1} Lucas(3*n)*(-x)^n).
1 - 5*Sum_{n >= 1} a(n)*x^(2*n) = (1 + Sum_{n >= 1} A001077(n+1)*x^n)*(1 + Sum_{n >= 1} A001077(n+1)*(-x)^n).
1 - 9*Sum_{n >= 1} a(n)*x^(2*n) = (1 - G(x))*(1 - G(-x)) = (1 + 9*G(x))*(1 + 9*G(-x)).
1 - 16*Sum_{n >= 1} a(n)*x^(2*n) = (1 + 2*Sum_{n >= 1} A099843(n)*x^n)*(1 + 2*Sum_{n >= 1} A099843(n)*(-x)^n).
1 - 20*Sum_{n >= 1} a(n)*x^(2*n) = (1 - 2*G(x))*(1 - 2*G(-x)) = (1 + 10*G(x))*(1 + 10*G(-x)).
(End)

Extensions

Chebyshev and other comments from Wolfdieter Lang, Nov 08 2002

A001077 Numerators of continued fraction convergents to sqrt(5).

Original entry on oeis.org

1, 2, 9, 38, 161, 682, 2889, 12238, 51841, 219602, 930249, 3940598, 16692641, 70711162, 299537289, 1268860318, 5374978561, 22768774562, 96450076809, 408569081798, 1730726404001, 7331474697802, 31056625195209
Offset: 0

Views

Author

Keywords

Comments

a(2*n+1) with b(2*n+1) := A001076(2*n+1), n >= 0, give all (positive integer) solutions to Pell equation a^2 - 5*b^2 = -1.
a(2*n) with b(2*n) := A001076(2*n), n >= 1, give all (positive integer) solutions to Pell equation a^2 - 5*b^2 = +1 (see Emerson reference).
Bisection: a(2*n) = T(n,9) = A023039(n), n >= 0 and a(2*n+1) = 2*S(2*n, 2*sqrt(5)) = A075796(n+1), n >= 0, with T(n,x), resp. S(n,x), Chebyshev's polynomials of the first, resp. second kind. See A053120, resp. A049310.
From Greg Dresden, May 21 2023: (Start)
For n >= 2, 8*a(n) is the number of ways to tile this T-shaped figure of length n-1 with four colors of squares and one color of domino; shown here is the figure of length 5 (corresponding to n=6), and it has 8*a(6) = 23112 different tilings.
_
|| _
|||_|||
|_|
(End)

Examples

			1  2  9  38  161  (A001077)
-, -, -, --, ---, ...
0  1  4  17   72  (A001076)
1 + 2*x + 9*x^2 + 38*x^3 + 161*x^4 + 682*x^5 + 2889*x^6 + 12238*x^7 + ... - _Michael Somos_, Aug 11 2009
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • V. Thébault, Les Récréations Mathématiques, Gauthier-Villars, Paris, 1952, p. 282.

Crossrefs

Programs

  • Magma
    I:=[1, 2]; [n le 2 select I[n] else 4*Self(n-1) + Self(n-2): n in [1..30]]; // G. C. Greubel, Dec 19 2017
  • Maple
    A001077:=(-1+2*z)/(-1+4*z+z**2); # conjectured by Simon Plouffe in his 1992 dissertation
    with(combinat): a:=n->fibonacci(n+1, 4)-2*fibonacci(n, 4): seq(a(n), n=0..30); # Zerinvary Lajos, Apr 04 2008
  • Mathematica
    LinearRecurrence[{4, 1}, {1, 2}, 30]
    Join[{1},Numerator[Convergents[Sqrt[5],30]]] (* Harvey P. Dale, Mar 23 2016 *)
    CoefficientList[Series[(1-2*x)/(1-4*x-x^2), {x, 0, 30}], x] (* G. C. Greubel, Dec 19 2017 *)
    LucasL[3*Range[0,30]]/2 (* Rigoberto Florez, Apr 03 2019 *)
    a[ n_] := LucasL[n, 4]/2; (* Michael Somos, Nov 02 2021 *)
  • PARI
    {a(n) = fibonacci(3*n) / 2 + fibonacci(3*n - 1)}; /* Michael Somos, Aug 11 2009 */
    
  • PARI
    a(n)=if(n<2,n+1,my(t=4);for(i=1,n-2,t=4+1/t);numerator(2+1/t)) \\ Charles R Greathouse IV, Dec 05 2011
    
  • PARI
    x='x+O('x^30); Vec((1-2*x)/(1-4*x-x^2)) \\ G. C. Greubel, Dec 19 2017
    
  • Sage
    [lucas_number2(n,4,-1)/2 for n in range(0, 30)] # Zerinvary Lajos, May 14 2009
    

Formula

G.f.: (1-2*x)/(1-4*x-x^2).
a(n) = 4*a(n-1) + a(n-2), a(0)=1, a(1)=2.
a(n) = ((2 + sqrt(5))^n + (2 - sqrt(5))^n)/2.
a(n) = A014448(n)/2.
Limit_{n->infinity} a(n)/a(n-1) = phi^3 = 2 + sqrt(5). - Gregory V. Richardson, Oct 13 2002
a(n) = ((-i)^n)*T(n, 2*i), with T(n, x) Chebyshev's polynomials of the first kind A053120 and i^2 = -1.
Binomial transform of A084057. - Paul Barry, May 10 2003
E.g.f.: exp(2x)cosh(sqrt(5)x). - Paul Barry, May 10 2003
a(n) = Sum_{k=0..floor(n/2)} binomial(n, 2k)*5^k*2^(n-2k). - Paul Barry, Nov 15 2003
a(n) = 4*a(n-1) + a(n-2) when n > 2; a(1) = 1, a(2) = 2. - Alex Vinokur (alexvn(AT)barak-online.net), Oct 25 2004
a(n) = A001076(n+1) - 2*A001076(n) = A097924(n) - A015448(n+1); a(n+1) = A097924(n) + 2*A001076(n) = A097924(n) + 2(A048876(n) - A048875(n)). - Creighton Dement, Mar 19 2005
a(n) = F(3*n)/2 + F(3*n-1) where F() = Fibonacci numbers A000045. - Gerald McGarvey, Apr 28 2007
a(n) = A000032(3*n)/2.
For n >= 1: a(n) = (1/2)*Fibonacci(6*n)/Fibonacci(3*n) and a(n) = integer part of (2 + sqrt(5))^n. - Artur Jasinski, Nov 28 2011
a(n) = Sum_{k=0..n} A201730(n,k)*4^k. - Philippe Deléham, Dec 06 2011
a(n) = A001076(n) + A015448(n). - R. J. Mathar, Jul 06 2012
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(5*k-4)/(x*(5*k+1) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 27 2013
a(n) is the (1,1)-entry of the matrix W^n with W=[2, sqrt(5); sqrt(5), 2]. - Carmine Suriano, Mar 21 2014
From Rigoberto Florez, Apr 03 2019: (Start)
a(n) = A099919(n) + A049651(n) if n > 0.
a(n) = 1 + Sum_{k=0..n-1} L(3*k + 1) if n >= 0, L(n) = n-th Lucas number (A000032). (End)
From Christopher Hohl, Aug 22 2021: (Start)
For n >= 2, a(2n-1) = A079962(6n-9) + A079962(6n-3).
For n >= 1, a(2n) = sqrt(20*A079962(6n-3)^2 + 1). (End)
a(n) = Sum_{k=0..n-2} A168561(n-2,k)*4^k + 2 * Sum_{k=0..n-1} A168561(n-1,k)*4^k, n>0. - R. J. Mathar, Feb 14 2024
a(n) = 4^n*Sum_{k=0..n} A374439(n, k)*(-1/4)^k. - Peter Luschny, Jul 26 2024
From Peter Bala, Jul 08 2025: (Start)
The following series telescope:
Sum_{n >= 1} 1/(a(n) + 5*(-1)^(n+1)/a(n)) = 3/8, since 1/(a(n) + 5*(-1)^(n+1)/a(n)) = b(n) - b(n+1), where b(n) = (1/4) * (a(n) + a(n-1)) / (a(n)*a(n-1)).
Sum_{n >= 1} (-1)^(n+1)/(a(n) + 5*(-1)^(n+1)/a(n)) = 1/8, since 1/(a(n) + 5*(-1)^(n+1)/a(n)) = c(n) + c(n+1), where c(n) = (1/4) * (a(n) - a(n-1)) / (a(n)*a(n-1)). (End)

Extensions

Chebyshev comments from Wolfdieter Lang, Jan 10 2003

A001566 a(0) = 3; thereafter, a(n) = a(n-1)^2 - 2.

Original entry on oeis.org

3, 7, 47, 2207, 4870847, 23725150497407, 562882766124611619513723647, 316837008400094222150776738483768236006420971486980607
Offset: 0

Views

Author

Keywords

Comments

Expansion of 1/phi: 1/phi = (1-1/3)*(1-1/((3-1)*7))*(1-1/(((3-1)*7-1)*47))*(1-1/((((3-1)*7-1)*47-1)*2207))... (phi being the golden ration (1+sqrt(5))/2). - Thomas Baruchel, Nov 06 2003
An infinite coprime sequence defined by recursion. - Michael Somos, Mar 14 2004
Starting with 7, the terms end with 7,47,07,47,07,..., of the form 8a+7 where a = 0,1,55,121771,... Conjecture: Every a is squarefree, every other a is divisible by 55, the a's are a subset of A046194, the heptagonal triangular numbers (the first, 2nd, 3rd, 6th, 11th, ?, ... terms). - Gerald McGarvey, Aug 08 2004
Also the reduced numerator of the convergents to sqrt(5) using Newton's recursion x = (5/x+x)/2. - Cino Hilliard, Sep 28 2008
The subsequence of primes begins a(n) for n = 0, 1, 2, 3. - Jonathan Vos Post, Feb 26 2011
We have Sum_{n=0..N} a(n)^2 = 2*(N+1) + Sum_{n=1..N+1} a(n), Sum_{n=0..N} a(n)^4 = 5*(Sum_{n=1..N+1} a(n)) + a(N+1)^2 + 6*N -3, etc. which is very interesting with respect to the fact that a(n) = Lucas(2^(n+1)); see W. Webb's problem in Witula-Slota's paper. - Roman Witula, Nov 02 2012
From Peter Bala, Nov 11 2012: (Start)
The present sequence corresponds to the case x = 3 of the following general remarks.
The recurrence a(n+1) = a(n)^2 - 2 with initial condition a(0) = x > 2 has the solution a(n) = ((x + sqrt(x^2 - 4))/2)^(2^n) + ((x - sqrt(x^2 - 4))/2)^(2^n).
We have the product expansion sqrt(x + 2)/sqrt(x - 2) = Product_{n>=0} (1 + 2/a(n)) (essentially due to Euler - see Mendes-France and van der Poorten). Another expansion is sqrt(x^2 - 4)/(x + 1) = Product_{n>=0} (1 - 1/a(n)), which follows by iterating the identity sqrt(x^2 - 4)/(x + 1) = (1 - 1/x)*sqrt(y^2 - 4)/(y + 1), where y = x^2 - 2.
The sequence b(n) := a(n) - 1 satisfies b(n+1) = b(n)^2 + 2*b(n) - 2. Cases currently in the database are A145502 through A145510. The sequence c(n) := a(n)/2 satisfies c(n+1) = 2*c(n)^2 - 1. Cases currently in the database are A002812, A001601, A005828, A084764 and A084765.
(End)
E. Lucas in Section XIX of "The Theory of Simply Periodic Numerical Functions" (page 56 of English translation) equation "(127) (1-sqrt(5))/2 = -1/1 + 1/3 + 1/(3*7) + 1/(3*7*47) + 1/(3*7*47*2207) + ..." - Michael Somos, Oct 11 2022
Let b(n) = a(n) - 3. The sequence {b(n)} appears to be a strong divisibility sequence, that is, gcd(b(n),b(m)) = b(gcd(n,m)) for n, m >= 1. - Peter Bala, Dec 08 2022
The number of digits of a(n) is given by A094057(n+1). - Hans J. H. Tuenter, Jul 29 2025

Examples

			From _Cino Hilliard_, Sep 28 2008: (Start)
Init x=1;
x = (5/1 + 1)/2 = 3/1;
x = (5/3 + 3)/2 = 7/3;
x = ((5/7)/3 + 7/3)/2 = 47/21;
x = ((5/47)/21 + 47/21)/2 = 2207/987;
(2207/987)^2 = 5.000004106... (End)
		

References

  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 1, p. 397.
  • E.-B. Escott, Note #1741, L'Intermédiaire des Mathématiciens, 8 (1901), page 13. - N. J. A. Sloane, Mar 02 2022
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 223.
  • Édouard Lucas, Nouveaux théorèmes d'arithmétique supérieure, Comptes Rend., 83 (1876), 1286-1288.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 7.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Lucas numbers (A000032) with subscripts that are powers of 2 greater than 1 (Herbert S. Wilf). Cf. A000045.
Cf. A003010 (starting with 4), A003423 (starting with 6), A003487 (starting with 5).
Cf. A058635. - Artur Jasinski, Oct 05 2008

Programs

  • Maple
    a:= n-> simplify(2*ChebyshevT(2^n, 3/2), 'ChebyshevT'):
    seq(a(n), n=0..8);
  • Mathematica
    NestList[#^2-2&,3,10] (* Harvey P. Dale, Dec 17 2014 *)
    Table[LucasL[2^n], {n, 1, 8}] (* Amiram Eldar, Oct 22 2020 *)
  • Maxima
    a[0]:3$
    a[n]:=a[n-1]^2-2$
    A001566(n):=a[n]$
    makelist(A001566(n),n,0,7); /* Martin Ettl, Nov 12 2012 */
  • PARI
    {a(n) = if( n<1, 3*(n==0), a(n-1)^2 - 2)}; /* Michael Somos, Mar 14 2004 */
    
  • PARI
    g(n,p) = x=1;for(j=1,p,x=(n/x+x)/2;print1(numerator(x)","));
    g(5,8) \\ Cino Hilliard, Sep 28 2008
    
  • PARI
    {a(n) = my(w = quadgen(5)); if( n<0, 0, n++; imag( (2*w - 1) * w^2^n ))}; /* Michael Somos, Nov 30 2014 */
    
  • PARI
    {a(n) = my(y = x^2-x-1); if( n<0, 0, n++; for(i=1, n, y = polgraeffe(y)); -polcoeff(y, 1))}; /* Michael Somos, Nov 30 2014 */
    

Formula

a(n) = Fibonacci(2^(n+2))/Fibonacci(2^(n+1)) = A058635(n+2)/A058635(n+1). - Len Smiley, May 08 2000, and Artur Jasinski, Oct 05 2008
a(n) = ceiling(c^(2^n)) where c = (3+sqrt(5))/2 = tau^2 is the largest root of x^2-3*x+1=0. - Benoit Cloitre, Dec 03 2002
a(n) = round(G^(2^n)) where G is the golden ratio (A001622). - Artur Jasinski, Sep 22 2008
a(n) = (G^(2^(n+1))-(1-G)^(2^(n+1)))/((G^(2^n))-(1-G)^(2^n)) = G^(2^n)+(1-G)^(2^n) = G^(2^n)+(-G)^(-2^n) where G is the golden ratio. - Artur Jasinski, Oct 05 2008
a(n) = 2*cosh(2^(n+1)*arccosh(sqrt(5)/2)). - Artur Jasinski, Oct 09 2008
a(n) = Fibonacci(2^(n+1)-1) + Fibonacci(2^(n+1)+1). (3-sqrt(5))/2 = 1/3 + 1/(3*7) + 1/(3*7*47) + 1/(3*7*47*2207) + ... (E. Lucas). - Philippe Deléham, Apr 21 2009
a(n)*(a(n+1)-1)/2 = A023039(2^n). - M. F. Hasler, Sep 27 2009
For n >= 1, a(n) = 2 + Product_{i=0..n-1} (a(i) + 2). - Vladimir Shevelev, Nov 28 2010
a(n) = 2*T(2^n,3/2) where T(n,x) is the Chebyshev polynomial of the first kind. - Leonid Bedratyuk, Mar 17 2011
From Peter Bala, Oct 31 2012: (Start)
Engel expansion of 1/2*(3 - sqrt(5)). Thus 1/2*(3 - sqrt(5)) = 1/3 + 1/(3*7) + 1/(3*7*47) + ... as noted above by Deleham. See Liardet and Stambul.
sqrt(5)/4 = Product_{n>=0} (1 - 1/a(n)).
sqrt(5) = Product_{n>=0} (1 + 2/a(n)). (End)
a(n) - 1 = A145502(n+1). - Peter Bala, Nov 11 2012
a(n) == 2 (mod 9), for n > 1. - Ivan N. Ianakiev, Dec 25 2013
From Amiram Eldar, Oct 22 2020: (Start)
a(n) = A000032(2^(n+1)).
Sum_{k>=0} 1/a(k) = -1 + A338304. (End)
a(n) = (A000045(m+2^(n+2))+A000045(m))/A000045(m+2^(n+1)) for any m>=0. - Alexander Burstein, Apr 10 2021
a(n) = 2*cos(2^n*arccos(3/2)). - Peter Luschny, Oct 12 2022
a(n) == -1 ( mod 2^(n+2) ). - Peter Bala, Nov 07 2022
a(n) = 5*Fibonacci(2^n)^2+2 = 5*A058635(n)^2+2, for n>0. - Jianglin Luo, Sep 21 2023
Sum_{n>=0} a(n)/Fibonacci(2^(n+2)) = A094874 (Sanford, 2016). - Amiram Eldar, Mar 01 2024

A008310 Triangle of coefficients of Chebyshev polynomials T_n(x).

Original entry on oeis.org

1, 1, -1, 2, -3, 4, 1, -8, 8, 5, -20, 16, -1, 18, -48, 32, -7, 56, -112, 64, 1, -32, 160, -256, 128, 9, -120, 432, -576, 256, -1, 50, -400, 1120, -1280, 512, -11, 220, -1232, 2816, -2816, 1024, 1, -72, 840, -3584, 6912, -6144, 2048, 13, -364, 2912, -9984, 16640, -13312, 4096
Offset: 0

Views

Author

Keywords

Comments

The row length sequence of this irregular array is A008619(n), n >= 0. Even or odd powers appear in increasing order starting with 1 or x for even or odd row numbers n, respectively. This is the standard triangle A053120 with 0 deleted. - Wolfdieter Lang, Aug 02 2014
Let T* denote the triangle obtained by replacing each number in this triangle by its absolute value. Then T* gives the coefficients for cos(nx) as a polynomial in cos x. - Clark Kimberling, Aug 04 2024

Examples

			Rows are: (1), (1), (-1,2), (-3,4), (1,-8,8), (5,-20,16) etc., since if c = cos(x): cos(0x) = 1, cos(1x) = 1c; cos(2x) = -1+2c^2; cos(3x) = -3c+4c^3, cos(4x) = 1-8c^2+8c^4, cos(5x) = 5c-20c^3+16c^5, etc.
From _Wolfdieter Lang_, Aug 02 2014: (Start)
This irregular triangle a(n,k) begins:
  n\k   0    1     2      3      4      5      6      7 ...
  0:    1
  1:    1
  2:   -1    2
  3:   -3    4
  4:    1   -8     8
  5:    5  -20    16
  6:   -1   18   -48     32
  7:   -7   56  -112     64
  8:    1  -32   160   -256    128
  9:    9 -120   432   -576    256
 10:   -1   50  -400   1120  -1280    512
 11:  -11  220 -1232   2816  -2816   1024
 12:    1  -72   840  -3584   6912  -6144   2048
 13:   13 -364  2912  -9984  16640 -13312   4096
 14:   -1   98 -1568   9408 -26880  39424 -28672   8192
 15:  -15  560 -6048  28800 -70400  92160 -61440  16384
  ...
T(4,x) = 1 - 8*x^2 + 8*x^4, T(5,x) = 5*x - 20*x^3 +16*x^5.
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • E. A. Guilleman, Synthesis of Passive Networks, Wiley, 1957, p. 593.
  • Yaroslav Zolotaryuk, J. Chris Eilbeck, "Analytical approach to the Davydov-Scott theory with on-site potential", Physical Review B 63, p543402, Jan. 2001. The authors write, "Since the algebra of these is 'hyperbolic', contrary to the usual Chebyshev polynomials defined on the interval 0 <= x <= 1, we call the set of functions (21) the hyperbolic Chebyshev polynomials." (This refers to the triangle T* described in Comments.)

Crossrefs

A039991 is a row reversed version, but has zeros which enable the triangle to be seen. Columns/diagonals are A011782, A001792, A001793, A001794, A006974, A006975, A006976 etc.
Reflection of A028297. Cf. A008312, A053112.
Row sums are one. Polynomial evaluations include A001075 (x=2), A001541 (x=3), A001091, A001079, A023038, A011943, A001081, A023039, A001085, A077422, A077424, A097308, A097310, A068203.
Cf. A053120.

Programs

  • Maple
    A008310 := proc(n,m) local x ; coeftayl(simplify(ChebyshevT(n,x),'ChebyshevT'),x=0,m) ; end: i := 0 : for n from 0 to 100 do for m from n mod 2 to n by 2 do printf("%d %d ",i,A008310(n,m)) ; i := i+1 ; od ; od ; # R. J. Mathar, Apr 20 2007
    # second Maple program:
    b:= proc(n) b(n):= `if`(n<2, 1, expand(2*b(n-1)-x*b(n-2))) end:
    T:= n-> (p-> (d-> seq(coeff(p, x, d-i), i=0..d))(degree(p)))(b(n)):
    seq(T(n), n=0..15);  # Alois P. Heinz, Sep 04 2019
  • Mathematica
    Flatten[{1, Table[CoefficientList[ChebyshevT[n, x], x], {n, 1, 13}]}]//DeleteCases[#, 0, Infinity]& (* or *) Flatten[{1, Table[Table[((-1)^k*2^(n-2 k-1)*n*Binomial[n-k, k])/(n-k), {k, Floor[n/2], 0, -1}], {n, 1, 13}]}] (* Eugeniy Sokol, Sep 04 2019 *)

Formula

a(n,m) = 2^(m-1) * n * (-1)^((n-m)/2) * ((n+m)/2-1)! / (((n-m)/2)! * m!) if n>0. - R. J. Mathar, Apr 20 2007
From Paul Weisenhorn, Oct 02 2019: (Start)
T_n(x) = 2*x*T_(n-1)(x) - T_(n-2)(x), T_0(x) = 1, T_1(x) = x.
T_n(x) = ((x+sqrt(x^2-1))^n + (x-sqrt(x^2-1))^n)/2. (End)
From Peter Bala, Aug 15 2022: (Start)
T(n,x) = [z^n] ( z*x + sqrt(1 + z^2*(x^2 - 1)) )^n.
Sum_{k = 0..2*n} binomial(2*n,k)*T(k,x) = (2^n)*(1 + x)^n*T(n,x).
exp( Sum_{n >= 1} T(n,x)*t^n/n ) = Sum_{n >= 0} P(n,x)*t^n, where P(n,x) denotes the n-th Legendre polynomial. (End)

Extensions

Additional comments and more terms from Henry Bottomley, Dec 13 2000
Edited: Corrected Cf. A039991 statement. Cf. A053120 added. - Wolfdieter Lang, Aug 06 2014

A087215 Lucas(6*n): a(n) = 18*a(n-1) - a(n-2), starting with a(0) = 2 and a(1) = 18.

Original entry on oeis.org

2, 18, 322, 5778, 103682, 1860498, 33385282, 599074578, 10749957122, 192900153618, 3461452808002, 62113250390418, 1114577054219522, 20000273725560978, 358890350005878082, 6440026026380244498, 115561578124838522882, 2073668380220713167378
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Oct 19 2003

Keywords

Comments

a(n+1)/a(n) converges to 9 + sqrt(80) = 17.9442719... a(0)/a(1) = 2/18; a(1)/a(2) = 18/322; a(2)/a(3) = 322/5778; a(3)/a(4) = 5778/103682; etc.
Lim_{n -> oo} a(n)/a(n+1) = 0.05572809000084... = 1/(9 + sqrt(80)) = 9 - sqrt(80).
From Peter Bala, Oct 13 2019: (Start)
Let F(x) = Product_{n >= 0} (1 + x^(4*n+1))/(1 + x^(4*n+3)). Let Phi = (1/2)*(sqrt(5) - 1). This sequence gives the partial denominators in the simple continued fraction expansion of the number F(Phi^6) = 1.0555459720... = 1 + 1/(18 + 1/(322 + 1/(5778 + ...))).
Also F(-Phi^6) = 0.9444348576... has the continued fraction representation 1 - 1/(18 - 1/(322 - 1/(5788 - ...))) and the simple continued fraction expansion 1/(1 + 1/((18 - 2) + 1/(1 + 1/((322 - 2) + 1/(1 + 1/((5788 - 2) + 1/(1 + ...))))))).
F(Phi^6)*F(-Phi^6) = 0.9968944099... has the simple continued fraction expansion 1/(1 + 1/((18^2 - 4) + 1/(1 + 1/((322^2 - 4) + 1/(1 + 1/((5788^2 - 4) + 1/(1 + ...))))))).
1/2 + (1/2)*F(Phi^6)/F(-Phi^6) = 1.0588241282... has the simple continued fraction expansion 1 + 1/((18 - 2) + 1/(1 + 1/((5778 - 2) + 1/(1 + 1/(1860498 - 2) + 1/(1 + ...))))). (End)

Examples

			a(4) = 103682 = 18*a(3) - a(2) = 18*5778 - 322 = (9 + sqrt(80))^4 + (9 - sqrt(80))^4 = 103681.99999035512... + 0.00000964487... = 103682.
		

References

  • J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 91.
  • R. P. Stanley. Enumerative combinatorics. Vol. 2, volume 62 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 1999.

Crossrefs

Cf. A074919.
Row 2 * 2 of array A188645.
Cf. Lucas(k*n): A000032 (k = 1), A005248 (k = 2), A014448 (k = 3), A056854 (k = 4), A001946 (k = 5), A087281 (k = 7), A087265 (k = 8), A087287 (k = 9), A065705 (k = 10), A089772 (k = 11), A089775 (k = 12).

Programs

  • Magma
    [ Lucas(6*n) : n in [0..100]]; // Vincenzo Librandi, Apr 14 2011
    
  • Mathematica
    a[0] = 2; a[1] = 18; a[n_] := 18a[n - 1] - a[n - 2]; Table[ a[n], {n, 0, 15}] (* Robert G. Wilson v, Jan 30 2004 *)
    Table[LucasL[6n], {n, 0, 18}]  (* or *) CoefficientList[Series[2*(1 - 9*x)/(1 - 18*x + x^2), {x, 0, 17}], x] (* Indranil Ghosh, Mar 15 2017 *)
  • PARI
    Vec(2*(1-9*x)/(1-18*x+x^2) + O(x^20)) \\ Colin Barker, Jan 24 2016
    
  • PARI
    a(n) = if(n<2, 17^n + 1, 18*a(n - 1) - a(n - 2));
    for(n=0, 17, print1(a(n),", ")) \\ Indranil Ghosh, Mar 15 2017

Formula

a(n) = A000032(6*n).
a(n) = 18*a(n-1) - a(n-2), starting with a(0) = 2 and a(1) = 18.
a(n) = (9 + sqrt(80))^n + (9 - sqrt(80))^n.
G.f.: 2*(1-9*x)/(1-18*x+x^2). - Philippe Deléham, Nov 17 2008
a(n) = 2*A023039(n). - R. J. Mathar, Oct 22 2010
From Peter Bala, Oct 13 2019: (Start)
a(n) = F(6*n+6)/F(6) - F(6*n-6)/F(6) = A049660(n+1) - A049660(n-1).
a(n) = trace(M^n), where M is the 2 X 2 matrix [0, 1; 1, 1]^6 = [5, 8; 8, 13].
Consequently the Gauss congruences hold: a(n*p^k) = a(n*p^(k-1)) ( mod p^k ) for all prime p and positive integers n and k. See Zarelua and also Stanley (Ch. 5, Ex. 5.2(a) and its solution).
16*Sum_{n >= 1} 1/(a(n) - 20/a(n)) = 1: (20 = Lucas(6) + 2 and 16 = Lucas(6) - 2)
20*Sum_{n >= 1} (-1)^(n+1)/(a(n) + 16/a(n)) = 1.
Series acceleration formulas for sum of reciprocals:
Sum_{n >= 1} 1/a(n) = 1/16 - 20*Sum_{n >= 1} 1/(a(n)*(a(n)^2 - 20)).
Sum_{n >= 1} (-1)^(n+1)/a(n) = 1/20 + 16*Sum_{n >= 1} (-1)^(n+1)/(a(n)*(a(n)^2 + 16)).
Sum_{n >= 1} 1/a(n) = ( (theta_3(9-4*sqrt(5)))^2 - 1 )/4 and
Sum_{n >= 1} (-1)^(n+1)/a(n) = ( 1 - (theta_3(4*sqrt(5)-9))^2 )/4,
where theta_3(x) = 1 + 2*Sum_{n >= 1} x^(n^2) (see A000122). Cf. A153415 and A003499.
x*exp(Sum_{n >= 1} a(n)*x^/n) = x + 18*x^2 + 323*x^3 + ... is the o.g.f. for A049660. (End)
E.g.f.: 2*exp(9*x)*cosh(4*sqrt(5)*x). - Stefano Spezia, Oct 18 2019
a(n) = L(2n-1)^2 * F(2n+1) + L(2n+1)^2 * F(2n-1), where F(n) = A000045(n) and L(n) = A000032(n). - Diego Rattaggi, Nov 12 2020
From Peter Bala, Apr 16 2025: (Start)
a(n) = Lucas(2*n)^3 - 3*Lucas(2*n) = 2*T(3, (1/2)*Lucas(2*n)), where T(k, x) denotes the k-th Chebyshev polynomial of the first kind; more generally, for k >= 0, Lucas(2*k*n) = 2*T(k, Lucas(2*n)/2).
Sum_{n >= 1} 1/a(n) = (1/4) * (theta_3(9 - sqrt(80))^2 - 1) and
Sum_{n >= 1} (-1)^(n+1)/a(n) = (1/4) * (1 - theta_3(sqrt(80) - 9)^2), where theta_3(x) = 1 + 2*Sum_{n >= 1} x^(n^2) (see A000122). See Borwein and Borwein, Proposition 3.5 (i), p. 91. Cf. A153415 and A003499. (End)

A056771 a(n) = a(-n) = 34*a(n-1) - a(n-2), and a(0)=1, a(1)=17.

Original entry on oeis.org

1, 17, 577, 19601, 665857, 22619537, 768398401, 26102926097, 886731088897, 30122754096401, 1023286908188737, 34761632124320657, 1180872205318713601, 40114893348711941777, 1362725501650887306817, 46292552162781456490001
Offset: 0

Views

Author

Henry Bottomley, Aug 16 2000

Keywords

Comments

The sequence satisfies the Pell equation a(n)^2 - 18 * A202299(n+1)^2 = 1. - Vincenzo Librandi, Dec 19 2011
Also numbers n such that n - 1 and 2*n + 2 are squares. - Arkadiusz Wesolowski, Mar 15 2015
And they, n - 1 and 2*n + 2, are the squares of A005319 and A003499. - Michel Marcus, Mar 15 2015
This sequence {a(n)} gives all the nonnegative integer solutions of the Pell equation a(n)^2 - 32*(3*A091761(n))^2 = +1. - Wolfdieter Lang, Mar 09 2019

Examples

			G.f. = 1 + 17*x + 577*x^2 + 19601*x^3 + 665857*x^4 + 22619537*x^5 + ...
		

Crossrefs

Cf. A001075, A001541, A001091, A001079, A023038, A011943, A001081, A023039, A001085 and note relationship with square triangular number sequences A001110 and A001109. A091761.
Row 3 of array A188644.

Programs

  • Magma
    I:=[1, 17]; [n le 2 select I[n] else 34*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Dec 18 2011
    
  • Mathematica
    LinearRecurrence[{34,-1},{1,17},30] (* Vincenzo Librandi, Dec 18 2011 *)
    a[ n_] := ChebyshevT[ 2 n, 3]; (* Michael Somos, May 28 2014 *)
  • Maxima
    makelist(expand(((17+sqrt(288))^n+(17-sqrt(288))^n))/2, n, 0, 15); /* Vincenzo Librandi, Dec 18 2011 */
    
  • PARI
    {a(n) = polchebyshev( n, 1, 17)}; /* Michael Somos, Apr 05 2019 */
  • Sage
    [lucas_number2(n,34,1)/2 for n in range(0,15)] # Zerinvary Lajos, Jun 27 2008
    

Formula

a(n) = (r^n + 1/r^n)/2 with r = 17 + sqrt(17^2-1).
a(n) = 16*A001110(n) + 1 = A001541(2n) = (4*A001109(n))^2 + 1 = 3*A001109(2n-1) - A001109(2n-2) = A001109(2n) - 3*A001109(2n-1).
a(n) = T(n, 17) = T(2*n, 3) with T(n, x) Chebyshev's polynomials of the first kind. See A053120. T(n, 3)= A001541(n).
G.f.: (1-17*x)/(1-34*x+x^2).
G.f.: (1 - 17*x / (1 - 288*x / (17 - x))). - Michael Somos, Apr 05 2019
a(n) = cosh(2n*arcsinh(sqrt(8))). - Herbert Kociemba, Apr 24 2008
a(n) = (a^n + b^n)/2 where a = 17 + 12*sqrt(2) and b = 17 - 12*sqrt(2); sqrt(a(n)-1)/4 = A001109(n). - James R. Buddenhagen, Dec 09 2011
a(-n) = a(n). - Michael Somos, May 28 2014
a(n) = sqrt(1 + 32*9*A091761(n)^2), n >= 0. See one of the Pell comments above. - Wolfdieter Lang, Mar 09 2019

Extensions

More terms from James Sellers, Sep 07 2000
Chebyshev comments from Wolfdieter Lang, Nov 29 2002

A188645 Array of ((k^n)+(k^(-n)))/2 where k=(sqrt(x^2+1)+x)^2 for integers x>=1.

Original entry on oeis.org

1, 3, 1, 17, 9, 1, 99, 161, 19, 1, 577, 2889, 721, 33, 1, 3363, 51841, 27379, 2177, 51, 1, 19601, 930249, 1039681, 143649, 5201, 73, 1, 114243, 16692641, 39480499, 9478657, 530451, 10657, 99, 1, 665857, 299537289, 1499219281, 625447713, 54100801, 1555849, 19601, 129, 1
Offset: 0

Views

Author

Charles L. Hohn, Apr 06 2011

Keywords

Comments

Conjecture: Given function f(x, y)=(sqrt(x^2+y)+x)^2; and constant k=f(x, y); then for all integers x>=1 and y=[+-]1, k may be irrational, but ((k^n)+(k^(-n)))/2 always produces integer sequences; y=1 results shown here; y=-1 results are A188644.
Also square array A(n,k), n >= 1, k >= 0, read by antidiagonals, where A(n,k) is Chebyshev polynomial of the first kind T_{k}(x), evaluated at x=2*n^2+1. - Seiichi Manyama, Jan 01 2019

Examples

			Square array begins:
     | 0    1       2          3             4
-----+---------------------------------------------
   1 | 1,   3,     17,        99,          577, ...
   2 | 1,   9,    161,      2889,        51841, ...
   3 | 1,  19,    721,     27379,      1039681, ...
   4 | 1,  33,   2177,    143649,      9478657, ...
   5 | 1,  51,   5201,    530451,     54100801, ...
   6 | 1,  73,  10657,   1555849,    227143297, ...
   7 | 1,  99,  19601,   3880899,    768398401, ...
   8 | 1, 129,  33281,   8586369,   2215249921, ...
   9 | 1, 163,  53137,  17322499,   5647081537, ...
  10 | 1, 201,  80801,  32481801,  13057603201, ...
  11 | 1, 243, 118097,  57394899,  27893802817, ...
  12 | 1, 289, 167041,  96549409,  55805391361, ...
  13 | 1, 339, 229841, 155831859, 105653770561, ...
  14 | 1, 393, 308897, 242792649, 190834713217, ...
  15 | 1, 451, 406801, 366934051, 330974107201, ...
  ...
		

Crossrefs

Row 1 is A001541, row 2 is A023039, row 3 is A078986, row 4 is A099370, row 5 is A099397, row 6 is A174747, row 8 is A176368, (row 1)*2 is A003499, (row 2)*2 is A087215.
Column 1 is A058331, (column 1)*2 is A005899.
A188644 (f(x, y) as above with y=-1).
Diagonal gives A173128.
Cf. A188647.

Programs

  • Mathematica
    max = 9; y = 1; t = Table[k = ((x^2 + y)^(1/2) + x)^2; ((k^n) + (k^(-n)))/2 // FullSimplify, {n, 0, max - 1}, {x, 1, max}]; Table[ t[[n - k + 1, k]], {n, 1, max}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jul 17 2013 *)

Formula

A(n,k) = (A188647(n,k-1) + A188647(n,k))/2.
A(n,k) = Sum_{j=0..k} binomial(2*k,2*j)*(n^2+1)^(k-j)*n^(2*j). - Seiichi Manyama, Jan 01 2019

Extensions

Edited and extended by Seiichi Manyama, Jan 01 2019
Showing 1-10 of 28 results. Next