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

A114525 Triangle of coefficients of the Lucas (w-)polynomials.

Original entry on oeis.org

2, 0, 1, 2, 0, 1, 0, 3, 0, 1, 2, 0, 4, 0, 1, 0, 5, 0, 5, 0, 1, 2, 0, 9, 0, 6, 0, 1, 0, 7, 0, 14, 0, 7, 0, 1, 2, 0, 16, 0, 20, 0, 8, 0, 1, 0, 9, 0, 30, 0, 27, 0, 9, 0, 1, 2, 0, 25, 0, 50, 0, 35, 0, 10, 0, 1, 0, 11, 0, 55, 0, 77, 0, 44, 0, 11, 0, 1, 2, 0, 36, 0, 105, 0, 112, 0, 54, 0, 12, 0, 1
Offset: 0

Views

Author

Eric W. Weisstein, Dec 06 2005

Keywords

Comments

Unsigned version of A108045.
The row reversed triangle is A162514. - Paolo Bonzini, Jun 23 2016

Examples

			2, x, 2 + x^2, 3*x + x^3, 2 + 4*x^2 + x^4, 5*x + 5*x^3 + x^5, ... give triangle
  n\k   0  1  2  3  4  5  6  7  8  9 10 ...
  0:    2
  1:    0  1
  2:    2  0  1
  3:    0  3  0  1
  4:    2  0  4  0  1
  5:    0  5  0  5  0  1
  6:    2  0  9  0  6  0  1
  7:    0  7  0 14  0  7  0  1
  8:    2  0 16  0 20  0  8  0  1
  9:    0  9  0 30  0 27  0  9  0  1
  10:   2  0 25  0 50  0 35  0 10  0  1
  n\k   0  1  2  3  4  5  6  7  8  9 10 ...
  .... reformatted by _Wolfdieter Lang_, Feb 10 2023
		

Crossrefs

Cf. A108045 (signed version).
Cf. Sequences L(n,x): A000032(x = 1), A002203 (x = 2), A006497 (x = 3), A014448 (x = 4), A087130 (x = 5), A085447 (x = 6), A086902 (x = 7), A086594 (x = 8), A087798 (x = 9), A086927 (x = 10), A001946 (x = 11), A086928 (x = 12), A088316 (x = 13), A090300 (x = 14), A090301 (x = 15), A090305 (x = 16), A090306 (x = 17), A090307 (x = 18), A090308 (x = 19), A090309 (x = 20), A090310 (x = 21), A090313 (x = 22), A090314 (x = 23), A090316 (x = 24), A087281 (x = 29), A087287 (x = 76), A089772 (x = 199).

Programs

  • Maple
    Lucas := proc(n,x)
        option remember;
        if  n=0 then
            2;
        elif n =1 then
            x ;
        else
            x*procname(n-1,x)+procname(n-2,x) ;
        end if;
    end proc:
    A114525 := proc(n,k)
        coeftayl(Lucas(n,x),x=0,k) ;
    end proc:
    seq(seq(A114525(n,k),k=0..n),n=0..12) ; # R. J. Mathar, Aug 16 2019
  • Mathematica
    row[n_] := CoefficientList[LucasL[n, x], x];
    Table[row[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Aug 11 2018 *)

Formula

From Peter Bala, Mar 18 2015: (Start)
The Lucas polynomials L(n,x) satisfy the recurrence L(n+1,x) = x*L(n,x) + L(n-1,x) with L(0,x) = 2 and L(1,x) = x.
O.g.f.: Sum_{n >= 0} L(n,x)*t^n = (2 - x*t)/(1 - t^2 - x*t) = 2 + x*t + (x^2 + 2)*t^2 + (3*x + x^3)*t^3 + ....
L(n,x) = trace( [ x, 1; 1, 0 ]^n ).
exp( Sum_{n >= 1} L(n,x)*t^n/n ) = Sum_{n >= 0} F(n+1,x)*t^n, where F(n,x) denotes the n-th Fibonacci polynomial. (see Appendix A3 in Johnson).
exp( Sum_{n >= 1} L(n,x)*L(2*n,x)*t^n/n ) = 1/( F(1,x)*F(2*x)*F(3,x) ) * Sum_{n >= 0} F(n+1,x)*F(n+2,x)*F(n+3,x)*t^n.
exp( Sum_{n >= 1} L(3*n,x)/L(n,x)*t^n/n ) = Sum_{n >= 0} L(2*n + 1,x)*t^n.
L(n,1) = Lucas(n) = A000032(n); L(n,4) = Lucas(3*n) = A014448(n); L(n,11) = Lucas(5*n) = A001946(n); L(n,29) = Lucas(7*n) = A087281(n); L(n,76) = Lucas(9*n) = A087287(n); L(n,199) = Lucas(11*n) = A089772(n). The general result is L(n,Lucas(2*k + 1)) = Lucas((2*k + 1)*n). (End)
From Jeremy Dover, Jun 10 2016: (Start)
Read as a triangle T(n,k), n >= 0, n >= k >= 0, T(n,k) = (Binomial((n+k)/2,k) + Binomial((n+k-2)/2,k))*(1+(-1)^(n-k))/2.
T(n,k) = A046854(n-1,k-1) + A046854(n-1,k) + A046854(n-2,k) for even n+k with n+k > 0, assuming A046854(n,k) = 0 for n < 0, k < 0, k > n.
T(n,k) is the number of binary strings of length n with exactly k pairs of consecutive 0's and no pair of consecutive 1's, where the first and last bits are considered consecutive. (End)
From Peter Bala, Sep 03 2022: (Start)
L(n,x) = 2*(i)^n*T(n,-i*x/2), where i = sqrt(-1) and T(n,x) is the n-th Chebyshev polynomial of the first kind.
d/dx(L(n,x)) = n*F(n,x), where F(n,x) denotes the n-th Fibonacci polynomial.
Let P_n(x,y) = (L(n,x) - L(n,y))/(x - y). Then {P_n(x,y): n >= 1} is a fourth-order linear divisibility sequence of polynomials in the ring Z[x,y]: if m divides n in Z then P_m(x,y) divides P_n(x,y) in Z[x,y].
P_n(1,1) = A045925(n); P_n(1,4) = A273622; P_n(2,2) = A093967(n).
L(2*n,x)^2 - L(2*n-1,x)*L(2*n+1,x) = x^2 + 4 for n >= 1.
Sum_{n >= 1} L(2*n,x)/( L(2*n-1,x) * L(2*n+1,x) ) = 1/x^2 and
Sum_{n >= 1} (-1)^(n+1)/( L(2*n,x) + x^2/L(2*n,x) ) = 1/(x^2 + 4), both valid for all nonzero real x. (End)
From Peter Bala, Nov 18 2022: (Start)
L(n,x) = Sum_{k = 0..floor(n/2)} (n/(n-k))*binomial(n-k,k)*x^(n-2*k) for n >= 1.
For odd m, L(n, L(m,x)) = L(n*m, x).
For integral x, the sequence {u(n)} := {L(n,x)} satisfies the Gauss congruences: u(m*p^r) == u(m*p^(r-1)) (mod p^r) for all positive integers m and r and all primes p.
Let p be an odd prime and let 0 <= k <= p - 1. Let alpha_k = the p-adic limit_{n -> oo} L(p^n,k). Then alpha_k is a well-defined p-adic-integer and the polynomial L(p,x) - x of degree p factorizes as L(p,x) - x = Product_{k = 0..p-1} (x - alpha_k). For example, L(5,x) - x = x^5 + 5*x^3 + 4*x = x*(x - A269591)*(x - A210850)*(x - A210851)*(x - A269592) in the ring of 5-adic integers. (End)
The formula for L(n,x) given in the first line of the preceding section, with L(0, x) = 2, is rewritten L(n, x) = Sum_{k = 0..floor(n/2)} A034807(n, k)*x^(n - 2*k). See the formula by Alexander Elkins in A034807. - Wolfdieter Lang, Feb 10 2023

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

Original entry on oeis.org

2, 15, 227, 3420, 51527, 776325, 11696402, 176222355, 2655031727, 40001698260, 602680505627, 9080209282665, 136805819745602, 2061167505466695, 31054318401746027, 467875943531657100, 7049193471376602527
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.066372... = 2/(15+sqrt(229)) = (sqrt(229)-15)/2.
Lim_{n-> infinity} a(n+1)/a(n) = 15.066372... = (15+sqrt(229))/2 = 2/(sqrt(229)-15).
For more information about this type of recurrence follow the Khovanova link and see A054413, A086902 and A178765. - Johannes W. Meijer, Jun 12 2010

Examples

			a(4) = 15*a(3) + a(2) = 15*3420 + 227 = ((15+sqrt(229))/2)^4 + ((15-sqrt(229))/2)^4 = 51526.9999805 + 0.0000194 = 51527.
		

Crossrefs

Lucas polynomials: A114525.
Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), this sequence (m=15), A090305 (m=16), A090306 (m=17), A090307 (m=18), A090308 (m=19), A090309 (m=20), A090310 (m=21), A090313 (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25), A087281 (m=29), A087287 (m=76), A089772 (m=199).

Programs

  • GAP
    m:=15;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 31 2019
  • Magma
    m:=15; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 31 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 15*I/2)), n = 0..20); # G. C. Greubel, Dec 31 2019
  • Mathematica
    LucasL[Range[20]-1, 15] (* G. C. Greubel, Dec 31 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 15*I/2) ) \\ G. C. Greubel, Dec 31 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 15*I/2) for n in (0..20)] # G. C. Greubel, Dec 31 2019
    

Formula

a(n) = 15*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 15.
a(n) = ((15+sqrt(229))/2)^n + ((15-sqrt(229))/2)^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5...
(a(n))^2 = a(2n) + 2 if n=2, 4, 6...
G.f.: (2-15*x)/(1-15*x-x^2). - Philippe Deléham, Nov 02 2008
Contribution from Johannes W. Meijer, Jun 12 2010: (Start)
Lim_{k-> infinity} a(n+k)/a(k) = (A090301(n) + A154597(n)*sqrt(229))/2.
Lim_{n-> infinity} A090301(n)/ A154597(n) = sqrt(229).
a(2n+1) = 15*A098246(n).
a(3n+1) = A041426(5n), a(3n+2) = A041426(5n+3), a(3n+3) = 2*A041426(5n+4).
(End)
a(n) = Lucas(n, 15) = 2*(-i)^n * ChebyshevT(n, 15*i/2). - G. C. Greubel, Dec 31 2019
E.g.f.: 2*exp(15*x/2)*cosh(sqrt(229)*x/2). - Stefano Spezia, Jan 01 2020

Extensions

More terms from Ray Chandler, Feb 14 2004

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

Original entry on oeis.org

2, 17, 291, 4964, 84679, 1444507, 24641298, 420346573, 7170533039, 122319408236, 2086600473051, 35594527450103, 607193567124802, 10357885168571737, 176691241432844331, 3014108989526925364, 51416544063390575519
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.058621... = 2/(17+sqrt(293)) = (sqrt(293)-17)/2.
Lim_{n-> infinity} a(n+1)/a(n) = 17.058621... = (17+sqrt(293))/2 = 2/(sqrt(293)-17).
For more information about this type of recurrence follow the Khovanova link and see A054413, A086902 and A178765. - Johannes W. Meijer, Jun 12 2010

Examples

			a(4) = 17*a(3) + a(2) = 17*4964 + 291=((17+sqrt(293))/2)^4 + ((17-sqrt(293))/2)^4 = 84678.999988190 + 0.000011809 = 84679.
		

Crossrefs

Cf. A005074.
Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), A090305 (m=16), this sequence (m=17), A090307 (m=18), A090308 (m=19), A090309 (m=20), A090310 (m=21), A090313 (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25).

Programs

  • GAP
    m:=17;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 30 2019
  • Magma
    m:=17; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 30 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 17*I/2)), n = 0..20); # G. C. Greubel, Dec 30 2019
  • Mathematica
    LinearRecurrence[{17,1},{2,17},30] (* Harvey P. Dale, Jan 24 2018 *)
    LucasL[Range[20]-1, 17] (* G. C. Greubel, Dec 30 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 17*I/2) ) \\ G. C. Greubel, Dec 30 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 17*I/2) for n in (0..20)] # G. C. Greubel, Dec 30 2019
    

Formula

a(n) = 17*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 17.
a(n) = ((17+sqrt(293))/2)^n + ((17-sqrt(293))/2)^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5, ...
(a(n))^2 = a(2n) + 2 if n=2, 4, 6, ...
G.f.: (2-17*x)/(1-17*x-x^2). - Philippe Deléham, Nov 02 2008
From Johannes W. Meijer, Jun 12 2010: (Start)
a(2n+1) = 17*A098249(n).
a(3n+1) = A041550(5n), a(3n+2) = A041550(5n+3), a(3n+3) = 2*A041550(5n+4).
Lim_{k-> infinity} a(n+k)/a(k) = (A090306(n) + A178765(n)*sqrt(293))/2.
Lim_{n-> infinity} A090306(n)/A178765(n) = sqrt(293). (End)
a(n) = Lucas(n, 17) = 2*(-i)^n * ChebyshevT(n, 17*i/2). - G. C. Greubel, Dec 30 2019
E.g.f.: 2*exp(17*x/2)*cosh(sqrt(293)*x/2). - Stefano Spezia, Dec 31 2019

Extensions

More terms from Ray Chandler, Feb 14 2004

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

Original entry on oeis.org

2, 20, 402, 8060, 161602, 3240100, 64963602, 1302512140, 26115206402, 523606640180, 10498248010002, 210488566840220, 4220269584814402, 84615880263128260, 1696537874847379602, 34015373377210720300
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.0498756... = 1/(10+sqrt(101)) = (sqrt(101)-10).
Lim_{n-> infinity} a(n+1)/a(n) = 20.0498756... = (10+sqrt(101)) = 1/(sqrt(101)-10).

Examples

			a(4) = 20*a(3) + a(2) = 20*8060 + 402 = (10+sqrt(101))^4 + (10-sqrt(101))^4 = 161601.999993811 + 0.000006188 = 161602.
		

Crossrefs

Cf. A002116.
Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), A090305 (m=16), A090306 (m=17), A090307 (m=18), A090308 (m=19), this sequence (m=20), A090310 (m=21), A090313 (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25).

Programs

  • GAP
    m:=20;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 30 2019
  • Magma
    m:=20; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 30 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 10*I)), n = 0..20); # G. C. Greubel, Dec 30 2019
  • Mathematica
    LinearRecurrence[{20,1},{2,20},20] (* Harvey P. Dale, Nov 19 2015 *)
    LucasL[Range[20]-1,20] (* G. C. Greubel, Dec 30 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 10*I) ) \\ G. C. Greubel, Dec 30 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 10*I) for n in (0..20)] # G. C. Greubel, Dec 30 2019
    

Formula

a(n) = 20*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 20.
a(n) = (10 + sqrt(101))^n + (10 - sqrt(101))^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5, ... .
(a(n))^2 = a(2n) + 2 if n=2, 4, 6, ... .
G.f.: (2-20*x)/(1-20*x-x^2). - Philippe Deléham, Nov 02 2008
a(n) = Lucas(n, 20) = 2*(-i)^n * ChebyshevT(n, 10*i). - G. C. Greubel, Dec 30 2019

Extensions

More terms from Ray Chandler, Feb 14 2004

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

Original entry on oeis.org

2, 21, 443, 9324, 196247, 4130511, 86936978, 1829807049, 38512885007, 810600392196, 17061121121123, 359094143935779, 7558038143772482, 159077895163157901, 3348193836570088403, 70471148463135014364
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.0475115... = 2/(21+sqrt(445)) = (sqrt(445)-21)/2.
Lim_{n-> infinity} a(n+1)/a(n) = 21.0475115... = (21+sqrt(445))/2 = 2/(sqrt(445)-21).
a(2) = 443 divides a(14) = 3348193836570088403. Does this relate to the sequence being the (21,1)-weighted Fibonacci sequence with seed (2,21) and both 14 and 21 being multiples of 7? Primes in this sequence include: a(0) = 2, a(2) = 443, a(4) = 196247 Semiprimes in this sequence include: a(8) = 38512885007 = 97967 * 393121, a(14) = 3348193836570088403 = 443 * 7557999631083721. - Jonathan Vos Post, Feb 10 2005

Examples

			a(4) = 21*a(3) + a(2) = 21*9324 + 443 = ((21+sqrt(445))/2)^4 + ((21-sqrt(445))/2)^4 = 196246.9999949043 + 0.0000050956 = 196247.
		

Crossrefs

Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), A090305 (m=16), A090306 (m=17), A090307 (m=18), A090308 (m=19), A090309 (m=20), this sequence (m=21), A090313 (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25).

Programs

  • GAP
    m:=21;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 30 2019
  • Magma
    m:=21; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 30 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 21*I/2)), n = 0..20); # G. C. Greubel, Dec 30 2019
  • Mathematica
    LinearRecurrence[{21,1},{2,21},40] (* or *) CoefficientList[ Series[ (2-21x)/(1-21x-x^2),{x,0,40}],x]  (* Harvey P. Dale, Apr 24 2011 *)
    LucasL[Range[20]-1,21] (* G. C. Greubel, Dec 30 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 21*I/2) ) \\ G. C. Greubel, Dec 30 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 21*I/2) for n in (0..20)] # G. C. Greubel, Dec 30 2019
    

Formula

a(n) = 21*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 21.
a(n) = ((21+sqrt(445))/2)^n + ((21-sqrt(445))/2)^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5... .
(a(n))^2 = a(2n) + 2 if n=2, 4, 6... .
G.f.: (2-21*x)/(1-21*x-x^2). - Philippe Deléham, Nov 02 2008
a(n) = Lucas(n, 21) = 2*(-i)^n * ChebyshevT(n, 21*i/2). - G. C. Greubel, Dec 30 2019

Extensions

More terms from Ray Chandler, Feb 14 2004

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

Original entry on oeis.org

2, 16, 258, 4144, 66562, 1069136, 17172738, 275832944, 4430499842, 71163830416, 1143051786498, 18359992414384, 294902930416642, 4736806879080656, 76083812995707138, 1222077814810394864, 19629328849962024962
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.0622577... = 1/(8+sqrt(65)) = (sqrt(65)-8).
Lim_{n-> infinity} a(n+1)/a(n) = 16.0622577... = (8+sqrt(65)) = 1/(sqrt(65)-8).

Examples

			a(4) = 16*a(3) + a(2) = 16*4144 + 258 = (8+sqrt(65))^4 + (8-sqrt(65))^4 = 66561.99998497... + 0.00001502... = 66562.
		

Crossrefs

Lucas polynomials: A114525.
Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), this sequence (m=16), A090306 (m=17), A090307 (m=18), A090308 (m=19), A090309 (m=20), A090310 (m=21), A090313 (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25), A087281 (m=29), A087287 (m=76), A089772 (m=199).

Programs

  • GAP
    m:=16;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 31 2019
  • Magma
    m:=16; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 31 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 8*I)), n = 0..20); # G. C. Greubel, Dec 31 2019
  • Mathematica
    LinearRecurrence[{16,1},{2,16},40] (* or *) With[{c=Sqrt[65]}, Simplify/@ Table[(c-8)((8+c)^n-(8-c)^n (129+16c)),{n,20}]] (* Harvey P. Dale, Oct 31 2011 *)
    LucasL[Range[20]-1, 16] (* G. C. Greubel, Dec 31 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 8*I) ) \\ G. C. Greubel, Dec 31 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 8*I) for n in (0..20)] # G. C. Greubel, Dec 31 2019
    

Formula

a(n) = 16*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 16.
a(n) = (8+sqrt(65))^n + (8-sqrt(65))^n.
a(n)^2 = a(2n) - 2 if n = 1, 3, 5, ...;
a(n)^2 = a(2n) + 2 if n = 2, 4, 6, ....
G.f.: (2-16*x)/(1-16*x-x^2). - Philippe Deléham, Nov 02 2008
a(n) = Lucas(n, 16) = 2*(-i)^n * ChebyshevT(n, 8*i). - G. C. Greubel, Dec 31 2019
E.g.f.: 2*exp(8*x)*cosh(sqrt(65)*x). - Stefano Spezia, Jan 01 2020

Extensions

More terms from Ray Chandler, Feb 14 2004

A090307 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, 326, 5886, 106274, 1918818, 34644998, 625528782, 11294163074, 203920464114, 3681862517126, 66477445772382, 1200275886420002, 21671443401332418, 391286257110403526, 7064824071388595886
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.0553851... = 1/(9+sqrt(82)) = (sqrt(82)-9).
Lim_{n-> infinity} a(n+1)/a(n) = 18.0553851... = (9+sqrt(82)) = 1/(sqrt(82)-9).

Examples

			a(4) = 18*a(3) + a(2) = 18*5886 + 326 = (9+sqrt(82))^4 + (9-sqrt(82))^4 = 106273.9999905903 + 0.000009406 = 106274.
		

Crossrefs

Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), A090305 (m=16), A090306 (m=17), this sequence (m=18), A090308 (m=19), A090309 (m=20), A090310 (m=21), A090313 (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25).

Programs

  • GAP
    m:=18;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 30 2019
  • Magma
    m:=18; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 30 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 9*I)), n = 0..20); # G. C. Greubel, Dec 30 2019
  • Mathematica
    LinearRecurrence[{18,1},{2,18},25] (* or *) CoefficientList[ Series[ (2-18x)/(1-18x-x^2),{x,0,25}],x] (* Harvey P. Dale, Apr 22 2011 *)
    LucasL[Range[20]-1, 18] (* G. C. Greubel, Dec 30 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 9*I) ) \\ G. C. Greubel, Dec 30 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 9*I) for n in (0..20)] # G. C. Greubel, Dec 30 2019
    

Formula

a(n) = 18*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 18.
a(n) = (9+sqrt(82))^n + (9-sqrt(82))^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5, ...
(a(n))^2 = a(2n) + 2 if n=2, 4, 6, ...
G.f.: (2-18*x)/(1-18*x-x^2). - Philippe Deléham, Nov 02 2008
a(n) = Lucas(n, 18) = 2*(-i)^n * ChebyshevT(n, 9*i). - G. C. Greubel, Dec 30 2019
E.g.f.: 2*exp(9*x)*cosh(sqrt(82)*x). - Stefano Spezia, Dec 31 2019

Extensions

More terms from Ray Chandler, Feb 14 2004

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

Original entry on oeis.org

2, 19, 363, 6916, 131767, 2510489, 47831058, 911300591, 17362542287, 330799604044, 6302555019123, 120079344967381, 2287810109399362, 43588471423555259, 830468767156949283, 15822495047405591636
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.052486... = 2/(19+sqrt(365)) = (sqrt(365)-19)/2.
Lim_{n-> infinity} a(n+1)/a(n) = 19.052486... = (19+sqrt(365))/2 = 2/(sqrt(365)-19).

Examples

			a(4) = 19*a(3) + a(2) = 19*6916 + 363 = ((19+sqrt(365))/2)^4 + ((19-sqrt(365))/2)^4 = 131766.9999924108 + 0.0000075891 = 131767.
		

Crossrefs

Cf. A049270.
Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), A090305 (m=16), A090306 (m=17), A090307 (m=18), this sequence (m=19), A090309 (m=20), A090310 (m=21), A090313 (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25).

Programs

  • GAP
    m:=19;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 30 2019
  • Magma
    m:=19; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 30 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 19*I/2)), n = 0..20); # G. C. Greubel, Dec 30 2019
  • Mathematica
    LucasL[Range[20]-1,20] (* G. C. Greubel, Dec 30 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 19*I/2) ) \\ G. C. Greubel, Dec 30 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 19*I/2) for n in (0..20)] # G. C. Greubel, Dec 30 2019
    

Formula

a(n) = 19*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 19.
a(n) = ((19+sqrt(365))/2)^n + ((19-sqrt(365))/2)^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5....
(a(n))^2 = a(2n) + 2 if n=2, 4, 6....
G.f.: (2-19*x)/(1-19*x-x^2). - Philippe Deléham, Nov 02 2008
a(n) = Lucas(n, 19) = 2*(-i)^n * ChebyshevT(n, 19*i/2). - G. C. Greubel, Dec 30 2019

Extensions

More terms from Ray Chandler, Feb 14 2004

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

Original entry on oeis.org

2, 22, 486, 10714, 236194, 5206982, 114789798, 2530582538, 55787605634, 1229857906486, 27112661548326, 597708411969658, 13176697724880802, 290485058359347302, 6403847981630521446, 141175140654230819114
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n-> infinity} a(n)/a(n+1) = 0.045361... = 1/(11+sqrt(122)) = (sqrt(122)-11).
Lim_{n-> infinity} a(n+1)/a(n) = 22.045361... = (11+sqrt(122)) = 1/(sqrt(122)-11).

Examples

			a(4) = 236194 = 22*a(3) + a(2) = 22*10714 + 486 = (11 + sqrt(122))^4 + (11 - sqrt(122))^4 = 236193.999995766 + 0.000004233 = 236194.
		

Crossrefs

Cf. A079219.
Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), A090305 (m=16), A090306 (m=17), A090307 (m=18), A090308 (m=19), A090309 (m=20), A090310 (m=21), this sequence (m=22), A090314 (m=23), A090316 (m=24), A330767 (m=25).

Programs

  • GAP
    m:=22;; a:=[2,m];; for n in [3..20] do a[n]:=m*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 30 2019
  • Magma
    m:=22; I:=[2,m]; [n le 2 select I[n] else m*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 30 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 11*I)), n = 0..20); # G. C. Greubel, Dec 30 2019
  • Mathematica
    LucasL[Range[20]-1,22] (* G. C. Greubel, Dec 29 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 11*I) ) \\ G. C. Greubel, Dec 30 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 11*I) for n in (0..20)] # G. C. Greubel, Dec 30 2019
    

Formula

a(n) = 22*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 22.
a(n) = (11+sqrt(122))^n + (11-sqrt(122))^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5...,
(a(n))^2 = a(2n) + 2 if n=2, 4, 6....
G.f.: (2-22*x)/(1-22*x-x^2). - Philippe Deléham, Nov 02 2008
a(n) = Lucas(n, 22) = 2*(-i)^n * ChebyshevT(n, 11*i). - G. C. Greubel, Dec 30 2019

Extensions

More terms from Ray Chandler, Feb 14 2004

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

Original entry on oeis.org

2, 23, 531, 12236, 281959, 6497293, 149719698, 3450050347, 79500877679, 1831970236964, 42214816327851, 972772745777537, 22415987969211202, 516540496037635183, 11902847396834820411, 274282030623238504636, 6320389551731320427039, 145643241720443608326533, 3356114949121934311937298
Offset: 0

Views

Author

Nikolay V. Kosinov (kosinov(AT)unitron.com.ua), Jan 25 2004

Keywords

Comments

Lim_{n -> infinity} a(n)/a(n+1) = 0.04339638... = 2/(23+sqrt(533)) = (sqrt(533)-23)/2.
Lim_{n -> infinity} a(n+1)/a(n) = 23.04339638... = (23+sqrt(533))/2 = 2/(sqrt(533) - 23).

Examples

			a(4) = 281959 = 23*a(3) + a(2) = 23*12236 + 531 = ((23 + sqrt(533))/2)^4 + ((23 - sqrt(533))/2)^4 = 281958.999996453 + 0.000003546 = 281959.
		

Crossrefs

Lucas polynomials Lucas(n,m): A000032 (m=1), A002203 (m=2), A006497 (m=3), A014448 (m=4), A087130 (m=5), A085447 (m=6), A086902 (m=7), A086594 (m=8), A087798 (m=9), A086927 (m=10), A001946 (m=11), A086928 (m=12), A088316 (m=13), A090300 (m=14), A090301 (m=15), A090305 (m=16), A090306 (m=17), A090307 (m=18), A090308 (m=19), A090309 (m=20), A090310 (m=21), A090313 (m=22), this sequence (m=23), A090316 (m=24), A330767 (m=25).

Programs

  • GAP
    a:=[2,23];; for n in [3..20] do a[n]:=23*a[n-1]+a[n-2]; od; a; # G. C. Greubel, Dec 29 2019
  • Magma
    I:=[2,23]; [n le 2 select I[n] else 23*Self(n-1) +Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 29 2019
    
  • Maple
    seq(simplify(2*(-I)^n*ChebyshevT(n, 23*I/2)), n = 0..20); # G. C. Greubel, Dec 29 2019
  • Mathematica
    LinearRecurrence[{23,1},{2,23},20] (* Harvey P. Dale, Jul 11 2014 *)
    LucasL[Range[20]-1,23] (* G. C. Greubel, Dec 29 2019 *)
  • PARI
    vector(21, n, 2*(-I)^(n-1)*polchebyshev(n-1, 1, 23*I/2) ) \\ G. C. Greubel, Dec 29 2019
    
  • Sage
    [2*(-I)^n*chebyshev_T(n, 23*I/2) for n in (0..20)] # G. C. Greubel, Dec 29 2019
    

Formula

a(n) = 23*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 23.
a(n) = ((23 + sqrt(533))/2)^n + ((23 - sqrt(533))/2)^n.
(a(n))^2 = a(2n) - 2 if n=1, 3, 5....
(a(n))^2 = a(2n) + 2 if n=2, 4, 6....
G.f.: (2-23*x)/(1-23*x-x^2). - Philippe Deléham, Nov 02 2008
a(n) = Lucas(n, 23) = 2*(-i)^n * ChebyshevT(n, 23*i/2). - G. C. Greubel, Dec 29 2019

Extensions

More terms from Ray Chandler, Feb 14 2004
Terms a(16) onward added by G. C. Greubel, Dec 29 2019
Showing 1-10 of 12 results. Next