|
|
A038723
|
|
a(n) = 6*a(n-1) - a(n-2), n >= 2, a(0)=1, a(1)=4.
|
|
15
|
|
|
1, 4, 23, 134, 781, 4552, 26531, 154634, 901273, 5253004, 30616751, 178447502, 1040068261, 6061962064, 35331704123, 205928262674, 1200237871921, 6995498968852, 40772755941191, 237641036678294, 1385073464128573
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
This sequence gives one half of all positive solutions y = y1 = a(n) of the first class of the Pell equation x^2 - 2*y^2 = -7. For the corresponding x=x1 terms see A054490(n). Therefore it also gives one fourth of all positive solutions x = x1 of the first class of the Pell equation x^2 - 2*y^2 = 14, with the y=y1 terms given by A054490. - Wolfdieter Lang, Feb 26 2015
|
|
REFERENCES
|
A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 122-125, 194-196.
|
|
LINKS
|
Seyed Hassan Alavi, Ashraf Daneshkhah, Cheryl E Praeger, Symmetries of biplanes, arXiv:2004.04535 [math.GR], 2020. See Lemma 7.9 p. 21.
|
|
FORMULA
|
a(n) = ((4+sqrt(2))/8)*(3+2*sqrt(2))^(n-1) + ((4-sqrt(2))/8)*(3-2*sqrt(2))^(n-1). - Antonio Alberto Olivares, Mar 29 2008
Sequence satisfies -7 = f(a(n), a(n+1)) where f(u, v) = u^2 + v^2 - 6*u*v. - Michael Somos, Sep 28 2008
G.f.: (1 - 2*x) / (1 - 6*x + x^2). a(n) = (7 + a(n-1)^2) / a(n-2). - Michael Somos, Sep 28 2008
a(n) = S(n,6) - 2*S(n-1, 6), n >= 0, with the Chebyshev polynomials S(n, x) (A049310) with S(-1, x) = 0 evaluated at x = 6. S(n, 6) = A001109(n-1). See the g.f. and the Pell comment above. - Wolfdieter Lang, Feb 26 2015
E.g.f.: exp(3*x)*(4*cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x))/4. - Stefano Spezia, Apr 30 2020
|
|
EXAMPLE
|
n = 2: A054490(2)^2 - 2*(2*a(2))^2 =
65^2 - 2*(2*23)^2 = -7,
|
|
MAPLE
|
a[0]:=1: a[1]:=4: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..20); # Zerinvary Lajos, Jul 26 2006
|
|
MATHEMATICA
|
LinearRecurrence[{6, -1}, {1, 4}, 30] (* Harvey P. Dale, Aug 06 2020 *)
|
|
PROG
|
(PARI) {a(n) = real((3 + 2*quadgen(8))^n * (1 + quadgen(8) / 4))} /* Michael Somos, Sep 28 2008 */
(PARI) {a(n) = polchebyshev(n, 1, 3) + polchebyshev(n-1, 2, 3)} /* Michael Somos, Sep 28 2008 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|