A154245 -id:A154245 - cp's OEIS frontend

A190958 a(n) = 2a(n-1) - 10a(n-2), with a(0) = 0, a(1) = 1.

0, 1, 2, -6, -32, -4, 312, 664, -1792, -10224, -2528, 97184, 219648, -532544, -3261568, -1197696, 30220288, 72417536, -157367808, -1038910976, -504143872, 9380822016, 23803082752, -46202054656, -330434936832, -198849327104, 2906650714112, 7801794699264

Offset: 0

Vladimir Joseph Stephan Orlovsky, May 24 2011

For the difference equation a(n) = c*a(n-1) - d*a(n-2), with a(0) = 0, a(1) = 1, the solution is a(n) = d^((n-1)/2) * ChebyshevU(n-1, c/(2*sqrt(d))) and has the alternate form a(n) = ( ((c + sqrt(c^2 - 4*d))/2)^n - ((c - sqrt(c^2 - 4*d))/2)^n )/sqrt(c^2 - 4*d). In the case c^2 = 4*d then the solution is a(n) = n*d^((n-1)/2). The generating function is x/(1 - c*x + d^2) and the exponential generating function takes the form (2/sqrt(c^2 - 4*d))*exp(c*x/2)*sinh(sqrt(c^2 - 4*d)*x/2) for c^2 > 4*d, (2/sqrt(4*d - c^2))*exp(c*x/2)*sin(sqrt(4*d - c^2)*x/2) for 4*d > c^2, and x*exp(sqrt(d)*x) if c^2 = 4*d. - G. C. Greubel, Jun 10 2022

Vincenzo Librandi, Table of n, a(n) for n = 0..190
Index entries for linear recurrences with constant coefficients, signature (2,-10).

Sequences of the form a(n) = c*a(n-1) - d*a(n-2), with a(0)=0, a(1)=1:
c/d...1.......2.......3.......4.......5.......6.......7.......8.......9......10
.A128834,A107920,A106852,A106853,A106854,A145934,A145976,A145978,A146078,A146080
.A000027,A009545,A088137,A088138,A045873,A088139,A089181,A090591,A127357,A190958
.A001906,A000225,A057083,A049072,A190959,A190960,A190961,A190962,A190963,A190964
.A001353,A007070,A003462,A001787,A099456,A190965,A168175,A190966,A190967,A190968
.A004254,A107839,A116415,A002450,A030191,A001047,A099450,A190969,A190970,A190971
.A001109,A154244,A138395,A084326,A003463,A030192,A081179,A006516,A027471,A106392
.A004187,A186446,A190972,A190973,A190974,A003464,A030240,A152268,A099459,A016127
.A001090,A190975,A190976,A099156,A190977,A190978,A023000,A057084,A154245,A154235
.A018913,A190979,A190980,A190981,A190982,A190983,A190984,A023001,A057085,A178869
A004189,A190869,A190985,A190986,A190987,A190988,A190989,A190990,A002452,A057086

I:=[0,1]; [n le 2 select I[n] else 2*Self(n-1)-10*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Sep 17 2011

Mathematica
```
LinearRecurrence[{2,-10}, {0,1}, 50]
```

a(n)=([0,1; -10,2]^n*[0;1])[1,1] \\ Charles R Greathouse IV, Apr 08 2016

[lucas_number1(n,2,10) for n in (0..50)] # G. C. Greubel, Jun 10 2022

G.f.: x / ( 1 - 2*x + 10*x^2 ). - R. J. Mathar, Jun 01 2011
E.g.f.: (1/3)*exp(x)*sin(3*x). - Franck Maminirina Ramaharo, Nov 13 2018
a(n) = 10^((n-1)/2) * ChebyshevU(n-1, 1/sqrt(10)). - G. C. Greubel, Jun 10 2022
a(n) = (1/3)*10^(n/2)*sin(n*arctan(3)) = Sum_{k=0..floor(n/2)} (-1)^k*3^(2*k)*binomial(n,2*k+1). - Gerry Martens, Oct 15 2022

A083099 a(n) = 2a(n-1) + 6a(n-2), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 2, 10, 32, 124, 440, 1624, 5888, 21520, 78368, 285856, 1041920, 3798976, 13849472, 50492800, 184082432, 671121664, 2446737920, 8920205824, 32520839168, 118562913280, 432250861568, 1575879202816, 5745263575040

Offset: 0

Mario Catalani (mario.catalani(AT)unito.it), Apr 22 2003

a(n+1) = a(n) + A083098(n+1). A083098(n+1)/a(n) converges to sqrt(7).
The same sequence may be obtained by the following process. Starting a priori with the fraction 1/1, the denominators of fractions built according to the rule: add top and bottom to get the new bottom, add top and 7 times the bottom to get the new top. The limit of the sequence of fractions is sqrt(7). - Cino Hilliard, Sep 25 2005
Pisano period lengths: 1, 1, 2, 1, 12, 2, 7, 1, 6, 12, 60, 2,168, 7, 12, 1,288, 6, 18, 12, ... - R. J. Mathar, Aug 10 2012
a(n) is divisible by 2^ceiling(n/2), see formula below. - Ralf Stephan, Dec 24 2013
Connect the center of a regular hexagon with side length 1 with its six vertices. a(n) is the number of paths of length n from the center to any of its vertices. Number of paths of length n from the center to itself is 6*a(n-1). - Jianing Song, Apr 20 2019

John Derbyshire, Prime Obsession, Joseph Henry Press, April 2004, see p. 16.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Project Euler, Problem 752, sequence beta(n).
Index entries for linear recurrences with constant coefficients, signature (2,6).

The following sequences (and others) belong to the same family: A000129, A001333, A002532, A002533, A002605, A015518, A015519, A026150, A046717, A063727, A083098, A083099, A083100, A084057.

Cf. A154245, A291008.

[n le 2 select n-1 else 2*Self(n-1) + 6*Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 24 2018

A083099 := proc(n)
    option remember;
    if n <= 1 then
        n;
    else
        2*procname(n-1)+6*procname(n-2) ;
    end if;
end proc: # R. J. Mathar, Sep 23 2016

CoefficientList[Series[x/(1-2x-6x^2), {x, 0, 25}], x] (* Adapted for offset 0 by Vincenzo Librandi, Feb 07 2014 *)
Expand[Table[((1 + Sqrt[7])^n - (1 - Sqrt[7])^n)7/(14Sqrt[7]), {n, 0, 25}]] (* Zerinvary Lajos, Mar 22 2007 *)
LinearRecurrence[{2,6}, {0,1}, 25] (* Sture Sjöstedt, Dec 06 2011 *)

a(n)=([0,1; 6,2]^n*[0;1])[1,1] \\ Charles R Greathouse IV, May 10 2016

my(x='x+O('x^30)); concat([0], Vec(x/(1-2*x-6*x^2))) \\ G. C. Greubel, Jan 24 2018

[lucas_number1(n,2,-6) for n in range(0, 25)] # Zerinvary Lajos, Apr 22 2009

A083099=BinaryRecurrenceSequence(2,6,0,1)
[A083099(n) for n in range(41)] # G. C. Greubel, Jun 01 2023

G.f.: x/(1 - 2*x - 6*x^2).
From Paul Barry, Sep 29 2004: (Start)
E.g.f.: (d/dx)(exp(x)*sinh(sqrt(7)*x)/sqrt(7));
a(n-1) = Sum_{k=0..n} binomial(n, 2k+1)*7^k. (End)
Simplified formula: a(n) = ((1 + sqrt(7))^n - (1 - sqrt(7))^n)/sqrt(28). - Al Hakanson (hawkuu(AT)gmail.com), Jan 05 2009
G.f.: G(0)*x/(2*(1-x)), where G(k) = 1 + 1/(1 - x*(7*k-1)/(x*(7*k+6) - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 26 2013
a(2n) = 2^n * A154245(n), a(2n+1) = 2^n * (5*A154245(n) - 9*A154245(n-1)). - Ralf Stephan, Dec 24 2013
a(n) = Sum_{k=1,3,5,...<=n} binomial(n,k)*7^((k-1)/2). - Vladimir Shevelev, Feb 06 2014
a(n) = i^(n-1)*6^((n-1)/2)*ChebyshevU(n-1, -i/sqrt(6)). - G. C. Greubel, Jun 01 2023

A094432 a(n) = rightmost term in M^n * [1 0 0]. M = the 3 X 3 stiffness matrix [1 -1 0 / -1 4 -3 / 0 -3 3].

Original entry on oeis.org

0, 3, 24, 165, 1104, 7347, 48840, 324597, 2157216, 14336355, 95275896, 633179973, 4207956720, 27965034003, 185848661544, 1235103986325, 8208193936704, 54549615616707, 362523179503320, 2409238895476197, 16011202548279696

Offset: 1

Gary W. Adamson, May 02 2004

A094431(n) = left term in M^n * [1 0 0]. A stiffness matrix in Hooke's Law governs the force on nodes of stretched or compressed springs (refer to A094431). a(n)/a(n-1) tends to 4 + sqrt(7) = 6.6457513...; a(n)/A094431(n) tends to 2 + sqrt(7). A stiffness matrix is symmetric.

			a(4) = 165 since M^4 * [1 0 0] = [38 -203 165].

Carl D. Meyer, "Matrix Analysis and Applied Linear Algebra", SIAM, 2000, pp. 86.-87.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (8,-9).

Cf. A094431, A154245.

Table[(MatrixPower[{{1, -1, 0}, {-1, 4, -3}, {0, -3, 3}}, n].{1, 0, 0})[[3]], {n, 21}] (* Robert G. Wilson v, May 08 2004 *)

a(n) = (3/(2*sqrt(7)))*((4+sqrt(7))^(n-1)-(4-sqrt(7))^(n-1)). For n>1, a(n) = 3*A154245(n-1). - Francesco Daddi, Aug 02 2011

G.f.: 3*x^2/(1-8*x+9*x^2). - Bruno Berselli, Aug 03 2011

More terms from Robert G. Wilson v, May 08 2004

cp's OEIS Frontend

A190958 a(n) = 2*a(n-1) - 10*a(n-2), with a(0) = 0, a(1) = 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

SageMath

Formula

A083099 a(n) = 2*a(n-1) + 6*a(n-2), a(0) = 0, a(1) = 1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Sage

SageMath

Formula

A094432 a(n) = rightmost term in M^n * [1 0 0]. M = the 3 X 3 stiffness matrix [1 -1 0 / -1 4 -3 / 0 -3 3].

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A190958 a(n) = 2a(n-1) - 10a(n-2), with a(0) = 0, a(1) = 1.

A083099 a(n) = 2a(n-1) + 6a(n-2), a(0) = 0, a(1) = 1.