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-9 of 9 results.

A077241 Combined Diophantine Chebyshev sequences A054488 and A077413.

Original entry on oeis.org

1, 2, 8, 13, 47, 76, 274, 443, 1597, 2582, 9308, 15049, 54251, 87712, 316198, 511223, 1842937, 2979626, 10741424, 17366533, 62605607, 101219572, 364892218, 589950899, 2126747701, 3438485822, 12395593988, 20040964033, 72246816227, 116807298376
Offset: 0

Views

Author

Wolfdieter Lang, Nov 08 2002

Keywords

Comments

-8*a(n)^2 + b(n)^2 = 17, with the companion sequence b(n)= A077242(n).
The number a > 0 belongs to the sequence A077241, if a^2 belongs to the sequence A034856. - Alzhekeyev Ascar M, Apr 27 2012
Numbers k such that k^2 + 2 is a triangular number (see A214838). - Alex Ratushnyak, Mar 07 2013

Examples

			8*a(2)^2 + 17 = 8*8^2+17 = 529 = 23^2 = A077242(2)^2.
		

Programs

  • Magma
    I:=[1,2,8,13]; [n le 4 select I[n] else 6*Self(n-2)-Self(n-4): n in [1..40]]; // Vincenzo Librandi, Feb 18 2014
  • Mathematica
    LinearRecurrence[{0, 6, 0, -1}, {1, 2, 8, 13}, 30] (* Bruno Berselli, Mar 10 2013 *)
    CoefficientList[Series[(1 + x) (1 + x + x^2)/(1 - 6 x^2 + x^4), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 18 2014 *)
  • Maxima
    makelist(expand((-1)^n*((4-5*sqrt(2))*(1-(-1)^n*sqrt(2))^(2*floor((n+1)/2))+(4+5*sqrt(2))*(1+(-1)^n*sqrt(2))^(2*floor((n+1)/2)))/8), n, 0, 30); /* Bruno Berselli, Mar 10 2013 */
    

Formula

a(2k) = A054488(k) and a(2k+1)= A077413(k) for k>=0.
G.f.: (1+x)*(1+x+x^2)/(1-6*x^2+x^4).
a(n) = (-1)^n*((4-5*sqrt(2))*(1-(-1)^n*sqrt(2))^(2*floor((n+1)/2))+(4+5*sqrt(2))*(1+(-1)^n*sqrt(2))^(2*floor((n+1)/2)))/8. [Bruno Berselli, Mar 10 2013]

A108931 a(2n) = -A106328(n), a(2n+1) = A054488(n).

Original entry on oeis.org

0, 1, -3, 8, -18, 47, -105, 274, -612, 1597, -3567, 9308, -20790, 54251, -121173, 316198, -706248, 1842937, -4116315, 10741424, -23991642, 62605607, -139833537, 364892218, -815009580, 2126747701, -4750223943, 12395593988, -27686334078, 72246816227
Offset: 0

Views

Author

Creighton Dement, Jul 20 2005

Keywords

Crossrefs

Programs

  • Maple
    seriestolist(series(x*(2*x-1)*(x-1)/((x^2+2*x-1)*(x^2-2*x-1)), x=0,30)); -or- Floretion Algebra Multiplication Program, FAMP Code: 1jescycrokseq[ + 'jj' + 'kk' + 'ij'] Roktype is set to: Y[sqa.Findk()] = Y[sqa.Findk()] + 1/2;

Formula

G.f. x*(2*x-1)*(x-1)/((x^2+2*x-1)*(x^2-2*x-1))

A001109 a(n)^2 is a triangular number: a(n) = 6*a(n-1) - a(n-2) with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 6, 35, 204, 1189, 6930, 40391, 235416, 1372105, 7997214, 46611179, 271669860, 1583407981, 9228778026, 53789260175, 313506783024, 1827251437969, 10650001844790, 62072759630771, 361786555939836, 2108646576008245, 12290092900109634, 71631910824649559, 417501372047787720
Offset: 0

Views

Author

Keywords

Comments

8*a(n)^2 + 1 = 8*A001110(n) + 1 = A055792(n+1) is a perfect square. - Gregory V. Richardson, Oct 05 2002
For n >= 2, A001108(n) gives exactly the positive integers m such that 1,2,...,m has a perfect median. The sequence of associated perfect medians is the present sequence. Let a_1,...,a_m be an (ordered) sequence of real numbers, then a term a_k is a perfect median if Sum_{j=1..k-1} a_j = Sum_{j=k+1..m} a_j. See Puzzle 1 in MSRI Emissary, Fall 2005. - Asher Auel, Jan 12 2006
(a(n), b(n)) where b(n) = A082291(n) are the integer solutions of the equation 2*binomial(b,a) = binomial(b+2,a). - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de); comment revised by Michael Somos, Apr 07 2003
This sequence gives the values of y in solutions of the Diophantine equation x^2 - 8y^2 = 1. It also gives the values of the product xy where (x,y) satisfies x^2 - 2y^2 = +-1, i.e., a(n) = A001333(n)*A000129(n). a(n) also gives the inradius r of primitive Pythagorean triangles having legs whose lengths are consecutive integers, with corresponding semiperimeter s = a(n+1) = {A001652(n) + A046090(n) + A001653(n)}/2 and area rs = A029549(n) = 6*A029546(n). - Lekraj Beedassy, Apr 23 2003 [edited by Jon E. Schoenfield, May 04 2014]
n such that 8*n^2 = floor(sqrt(8)*n*ceiling(sqrt(8)*n)). - Benoit Cloitre, May 10 2003
For n > 0, ratios a(n+1)/a(n) may be obtained as convergents to continued fraction expansion of 3+sqrt(8): either successive convergents of [6;-6] or odd convergents of [5;1, 4]. - Lekraj Beedassy, Sep 09 2003
a(n+1) + A053141(n) = A001108(n+1). Generating floretion: - 2'i + 2'j - 'k + i' + j' - k' + 2'ii' - 'jj' - 2'kk' + 'ij' + 'ik' + 'ji' + 'jk' - 2'kj' + 2e ("jes" series). - Creighton Dement, Dec 16 2004
Kekulé numbers for certain benzenoids (see the Cyvin-Gutman reference). - Emeric Deutsch, Jun 19 2005
Number of D steps on the line y=x in all Delannoy paths of length n (a Delannoy path of length n is a path from (0,0) to (n,n), consisting of steps E=(1,0), N=(0,1) and D=(1,1)). Example: a(2)=6 because in the 13 (=A001850(2)) Delannoy paths of length 2, namely (DD), (D)NE, (D)EN, NE(D), NENE, NEEN, NDE, NNEE, EN(D), ENNE, ENEN, EDN and EENN, we have altogether six D steps on the line y=x (shown between parentheses). - Emeric Deutsch, Jul 07 2005
Define a T-circle to be a first-quadrant circle with integral radius that is tangent to the x- and y-axes. Such a circle has coordinates equal to its radius. Let C(0) be the T-circle with radius 1. Then for n > 0, define C(n) to be the smallest T-circle that does not intersect C(n-1). C(n) has radius a(n+1). Cf. A001653. - Charlie Marion, Sep 14 2005
Numbers such that there is an m with t(n+m)=2t(m), where t(n) are the triangular numbers A000217. For instance, t(20)=2*t(14)=210, so 6 is in the sequence. - Floor van Lamoen, Oct 13 2005
One half the bisection of the Pell numbers (A000129). - Franklin T. Adams-Watters, Jan 08 2006
Pell trapezoids: for n > 0, a(n) = (A000129(n-1)+A000129(n+1))*A000129(n)/2; see also A084158. - Charlie Marion, Apr 01 2006
Tested for 2 < p < 27: If and only if 2^p - 1 (the Mersenne number M(p)) is prime then M(p) divides a(2^(p-1)). - Kenneth J Ramsey, May 16 2006
If 2^p - 1 is prime then M(p) divides a(2^(p-1)-1). - Kenneth J Ramsey, Jun 08 2006; comment corrected by Robert Israel, Mar 18 2007
If 8*n+5 and 8*n+7 are twin primes then their product divides a(4*n+3). - Kenneth J Ramsey, Jun 08 2006
If p is an odd prime, then if p == 1 or 7 (mod 8), then a((p-1)/2) == 0 (mod p) and a((p+1)/2) == 1 (mod p); if p == 3 or 5 (mod 8), then a((p-1)/2) == 1 (mod p) and a((p+1)/2) == 0 (mod p). Kenneth J Ramsey's comment about twin primes follows from this. - Robert Israel, Mar 18 2007
a(n)*(a(n+b) - a(b-2)) = (a(n+1)+1)*(a(n+b-1) - a(b-1)). This identity also applies to any series a(0) = 0 a(1) = 1 a(n) = b*a(n-1) - a(n-2). - Kenneth J Ramsey, Oct 17 2007
For n < 0, let a(n) = -a(-n). Then (a(n+j) + a(k+j)) * (a(n+b+k+j) - a(b-j-2)) = (a(n+j+1) + a(k+j+1)) * (a(n+b+k+j-1) - a(b-j-1)). - Charlie Marion, Mar 04 2011
Sequence gives y values of the Diophantine equation: 0+1+2+...+x = y^2. If (a,b) and (c,d) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with aMohamed Bouhamida, Aug 29 2009
If (p,q) and (r,s) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with p < r then r = 3*p+4*q+1 and s = 2*p+3*q+1. - Mohamed Bouhamida, Sep 02 2009
a(n)/A002315(n) converges to cos^2(Pi/8) (see A201488). - Gary Detlefs, Nov 25 2009
Binomial transform of A086347. - Johannes W. Meijer, Aug 01 2010
If x=a(n), y=A055997(n+1) and z = x^2+y, then x^4 + y^3 = z^2. - Bruno Berselli, Aug 24 2010
In general, if b(0)=1, b(1)=k and for n > 1, b(n) = 6*b(n-1) - b(n-2), then
for n > 0, b(n) = a(n)*k-a(n-1); e.g.,
for k=2, when b(n) = A038725(n), 2 = 1*2 - 0, 11 = 6*2 - 1, 64 = 35*2 - 6, 373 = 204*2 - 35;
for k=3, when b(n) = A001541(n), 3 = 1*3 - 0, 17 = 6*3 - 1; 99 = 35*3 - 6; 577 = 204*3 - 35;
for k=4, when b(n) = A038723(n), 4 = 1*4 - 0, 23 = 6*4 - 1; 134 = 35*4 - 6; 781 = 204*4 - 35;
for k=5, when b(n) = A001653(n), 5 = 1*5 - 0, 29 = 6*5 - 1; 169 = 35*5 - 6; 985 = 204*5 - 35.
- Charlie Marion, Dec 08 2010
See a Wolfdieter Lang comment on A001653 on a sequence of (u,v) values for Pythagorean triples (x,y,z) with x=|u^2-v^2|, y=2*u*v and z=u^2+v^2, with u odd and v even, generated from (u(0)=1,v(0)=2), the triple (3,4,5), by a substitution rule given there. The present a(n) appears there as b(n). The corresponding generated triangles have catheti differing by one length unit. - Wolfdieter Lang, Mar 06 2012
a(n)*a(n+2k) + a(k)^2 and a(n)*a(n+2k+1) + a(k)*a(k+1) are triangular numbers. Generalizes description of sequence. - Charlie Marion, Dec 03 2012
a(n)*a(n+2k) + a(k)^2 is the triangular square A001110(n+k). a(n)*a(n+2k+1) + a(k)*a(k+1) is the triangular oblong A029549(n+k). - Charlie Marion, Dec 05 2012
From Richard R. Forberg, Aug 30 2013: (Start)
The squares of a(n) are the result of applying triangular arithmetic to the squares, using A001333 as the "guide" on what integers to square, as follows:
a(2n)^2 = A001333(2n)^2 * (A001333(2n)^2 - 1)/2;
a(2n+1)^2 = A001333(2n+1)^2 * (A001333(2n+1)^2 + 1)/2. (End)
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,5}. - Milan Janjic, Jan 25 2015
Panda and Rout call these "balancing numbers" and note that the period of the sequence modulo a prime p is the same as that modulo p^2 when p = 13, 31, 1546463. But these are precisely the p in A238736 such that p^2 divides A000129(p - (2/p)), where (2/p) is a Jacobi symbol. In light of the above observation by Franklin T. Adams-Watters that the present sequence is one half the bisection of the Pell numbers, i.e., a(n) = A000129(2*n)/2, it follows immediately that modulo a fixed prime p, or any power thereof, the period of a(n) is half that of A000129(n). - John Blythe Dobson, Mar 06 2015
The triangular number = square number identity Tri((T(n, 3) - 1)/2) = S(n-1, 6)^2 with Tri, T, and S given in A000217, A053120 and A049310, is the special case k = 1 of the k-family of identities Tri((T(n, 2*k+1) - 1)/2) = Tri(k)*S(n-1, 2*(2*k+1))^2, k >= 0, n >= 0, with S(-1, x) = 0. For k=2 see A108741(n) for S(n-1, 10)^2. This identity boils down to the identities S(n-1, 2*x)^2 = (T(2*n, x) - 1)/(2*(x^2-1)) and 2*T(n, x)^2 - 1 = T(2*n, x) with x = 2*k+1. - Wolfdieter Lang, Feb 01 2016
a(2)=6 is perfect. For n=2*k, k > 0, k not equal to 1, a(n) is a multiple of a(2) and since every multiple (beyond 1) of a perfect number is abundant, then a(n) is abundant. sigma(a(4)) = 504 > 408 = 2*a(4). For n=2*k+1, k > 0, a(n) mod 10 = A000012(n), so a(n) is odd. If a(n) is a prime number, it is deficient; otherwise a(n) has one or two distinct prime factors and is therefore deficient again. So for n=2k+1, k > 0, a(n) is deficient. sigma(a(5)) = 1260 < 2378 = 2*a(5). - Muniru A Asiru, Apr 14 2016
Behera & Panda call these the balancing numbers, and A001541 are the balancers. - Michel Marcus, Nov 07 2017
In general, a second-order linear recurrence with constant coefficients having a signature of (c,d) will be duplicated by a third-order recurrence having a signature of (x,c^2-c*x+d,-d*x+c*d). The formulas of Olivares and Bouhamida in the formula section which have signatures of (7,-7,1) and (5,5,-1), respectively, are specific instances of this general rule for x = 7 and x = 5. - Gary Detlefs, Jan 29 2021
Note that 6 is the largest triangular number in the sequence, because it is proved that 8 and 9 are the largest perfect powers which are consecutive (Catalan's conjecture). 0 and 1 are also in the sequence because they are also perfect powers and 0*1/2 = 0^2 and 8*9/2 = (2*3)^2. - Metin Sariyar, Jul 15 2021

Examples

			G.f. = x + 6*x^2 + 35*x^3 + 204*x^4 + 1189*x^5 + 6930*x^6 + 40391*x^7 + ...
6 is in the sequence since 6^2 = 36 is a triangular number: 36 = 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8. - _Michael B. Porter_, Jul 02 2016
		

References

  • Julio R. Bastida, Quadratic properties of a linearly recurrent sequence. Proceedings of the Tenth Southeastern Conference on Combinatorics, Graph Theory and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1979), pp. 163--166, Congress. Numer., XXIII-XXIV, Utilitas Math., Winnipeg, Man., 1979. MR0561042 (81e:10009) - From N. J. A. Sloane, May 30 2012
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, pp. 193, 197.
  • D. M. Burton, The History of Mathematics, McGraw Hill, (1991), p. 213.
  • 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. 2, p. 10.
  • P. Franklin, E. F. Beckenbach, H. S. M Coxeter, N. H. McCoy, K. Menger, and J. L. Synge, Rings And Ideals, No 8, The Carus Mathematical Monographs, The Mathematical Association of America, (1967), pp. 144-146.
  • A. Patra, G. K. Panda, and T. Khemaratchatakumthorn. "Exact divisibility by powers of the balancing and Lucas-balancing numbers." Fibonacci Quart., 59:1 (2021), 57-64; see B(n).
  • 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).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 257-258.
  • P.-F. Teilhet, Query 2376, L'Intermédiaire des Mathématiciens, 11 (1904), 138-139. - N. J. A. Sloane, Mar 08 2022

Crossrefs

Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), this sequence (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.

Programs

  • GAP
    a:=[0,1];; for n in [3..25] do a[n]:=6*a[n-1]-a[n-2]; od; a; # Muniru A Asiru, Dec 18 2018
  • Haskell
    a001109 n = a001109_list !! n :: Integer
    a001109_list = 0 : 1 : zipWith (-)
       (map (* 6) $ tail a001109_list) a001109_list
    -- Reinhard Zumkeller, Dec 17 2011
    
  • Magma
    [n le 2 select n-1 else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Jul 25 2015
    
  • Maple
    a[0]:=1: a[1]:=6: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n],n=0..26); # Emeric Deutsch
    with (combinat):seq(fibonacci(2*n,2)/2, n=0..20); # Zerinvary Lajos, Apr 20 2008
  • Mathematica
    Transpose[NestList[Flatten[{Rest[#],ListCorrelate[{-1,6},#]}]&, {0,1}, 30]][[1]]  (* Harvey P. Dale, Mar 23 2011 *)
    CoefficientList[Series[x/(1-6x+x^2),{x,0,30}],x]  (* Harvey P. Dale, Mar 23 2011 *)
    LinearRecurrence[{6, -1}, {0, 1}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 12 2012 *)
    a[ n_]:= ChebyshevU[n-1, 3]; (* Michael Somos, Sep 02 2012 *)
    Table[Fibonacci[2n, 2]/2, {n, 0, 20}] (* Vladimir Reshetnikov, Sep 16 2016 *)
    TrigExpand@Table[Sinh[2 n ArcCsch[1]]/(2 Sqrt[2]), {n, 0, 10}] (* Federico Provvedi, Feb 01 2021 *)
  • PARI
    {a(n) = imag((3 + quadgen(32))^n)}; /* Michael Somos, Apr 07 2003 */
    
  • PARI
    {a(n) = subst( poltchebi( abs(n+1)) - 3 * poltchebi( abs(n)), x, 3) / 8}; /* Michael Somos, Apr 07 2003 */
    
  • PARI
    {a(n) = polchebyshev( n-1, 2, 3)}; /* Michael Somos, Sep 02 2012 */
    
  • PARI
    is(n)=ispolygonal(n^2,3) \\ Charles R Greathouse IV, Nov 03 2016
    
  • Sage
    [lucas_number1(n,6,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
    
  • Sage
    [chebyshev_U(n-1,3) for n in (0..20)] # G. C. Greubel, Dec 23 2019
    

Formula

G.f.: x / (1 - 6*x + x^2). - Simon Plouffe in his 1992 dissertation.
a(n) = S(n-1, 6) = U(n-1, 3) with U(n, x) Chebyshev's polynomials of the second kind. S(-1, x) := 0. Cf. triangle A049310 for S(n, x).
a(n) = sqrt(A001110(n)).
a(n) = A001542(n)/2.
a(n) = sqrt((A001541(n)^2-1)/8) (cf. Richardson comment).
a(n) = 3*a(n-1) + sqrt(8*a(n-1)^2+1). - R. J. Mathar, Oct 09 2000
a(n) = A000129(n)*A001333(n) = A000129(n)*(A000129(n)+A000129(n-1)) = ceiling(A001108(n)/sqrt(2)). - Henry Bottomley, Apr 19 2000
a(n) ~ (1/8)*sqrt(2)*(sqrt(2) + 1)^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
Limit_{n->oo} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 05 2002
a(n) = ((3 + 2*sqrt(2))^n - (3 - 2*sqrt(2))^n) / (4*sqrt(2)). - Gregory V. Richardson, Oct 13 2002. Corrected for offset 0, and rewritten. - Wolfdieter Lang, Feb 10 2015
a(2*n) = a(n)*A003499(n). 4*a(n) = A005319(n). - Mario Catalani (mario.catalani(AT)unito.it), Mar 21 2003
a(n) = floor((3+2*sqrt(2))^n/(4*sqrt(2))). - Lekraj Beedassy, Apr 23 2003
a(-n) = -a(n). - Michael Somos, Apr 07 2003
For n >= 1, a(n) = Sum_{k=0..n-1} A001653(k). - Charlie Marion, Jul 01 2003
For n > 0, 4*a(2*n) = A001653(n)^2 - A001653(n-1)^2. - Charlie Marion, Jul 16 2003
For n > 0, a(n) = Sum_{k = 0..n-1}((2*k+1)*A001652(n-1-k)) + A000217(n). - Charlie Marion, Jul 18 2003
a(2*n+1) = a(n+1)^2 - a(n)^2. - Charlie Marion, Jan 12 2004
a(k)*a(2*n+k) = a(n+k)^2 - a(n)^2; e.g., 204*7997214 = 40391^2 - 35^2. - Charlie Marion, Jan 15 2004
For j < n+1, a(k+j)*a(2*n+k-j) - Sum_{i = 0..j-1} a(2*n-(2*i+1)) = a(n+k)^2 - a(n)^2. - Charlie Marion, Jan 18 2004
From Paul Barry, Feb 06 2004: (Start)
a(n) = A000129(2*n)/2;
a(n) = ((1+sqrt(2))^(2*n) - (1-sqrt(2))^(2*n))*sqrt(2)/8;
a(n) = Sum_{i=0..n} Sum_{j=0..n} A000129(i+j)*n!/(i!*j!*(n-i-j)!)/2. (End)
E.g.f.: exp(3*x)*sinh(2*sqrt(2)*x)/(2*sqrt(2)). - Paul Barry, Apr 21 2004
A053141(n+1) + A055997(n+1) = A001541(n+1) + a(n+1). - Creighton Dement, Sep 16 2004
a(n) = Sum_{k=0..n} binomial(2*n, 2*k+1)*2^(k-1). - Paul Barry, Oct 01 2004
a(n) = A001653(n+1) - A038723(n); (a(n)) = chuseq[J]( 'ii' + 'jj' + .5'kk' + 'ij' - 'ji' + 2.5e ), apart from initial term. - Creighton Dement, Nov 19 2004, modified by Davide Colazingari, Jun 24 2016
a(n+1) = Sum_{k=0..n} A001850(k)*A001850(n-k), self convolution of central Delannoy numbers. - Benoit Cloitre, Sep 28 2005
a(n) = 7*(a(n-1) - a(n-2)) + a(n-3), a(1) = 0, a(2) = 1, a(3) = 6, n > 3. Also a(n) = ( (1 + sqrt(2) )^(2*n) - (1 - sqrt(2) )^(2*n) ) / (4*sqrt(2)). - Antonio Alberto Olivares, Oct 23 2003
a(n) = 5*(a(n-1) + a(n-2)) - a(n-3). - Mohamed Bouhamida, Sep 20 2006
Define f(x,s) = s*x + sqrt((s^2-1)*x^2+1); f(0,s)=0. a(n) = f(a(n-1),3), see second formula. - Marcos Carreira, Dec 27 2006
The perfect median m(n) can be expressed in terms of the Pell numbers P() = A000129() by m(n) = P(n + 2) * (P(n + 2) + P(n + 1)) for n >= 0. - Winston A. Richards (ugu(AT)psu.edu), Jun 11 2007
For k = 0..n, a(2*n-k) - a(k) = 2*a(n-k)*A001541(n). Also, a(2*n+1-k) - a(k) = A002315(n-k)*A001653(n). - Charlie Marion, Jul 18 2007
[A001653(n), a(n)] = [1,4; 1,5]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(n) = Sum_{k=0..n-1} 4^k*binomial(n+k,2*k+1). - Paul Barry, Apr 20 2009
a(n+1)^2 - 6*a(n+1)*a(n) + a(n)^2 = 1. - Charlie Marion, Dec 14 2010
a(n) = A002315(m)*A011900(n-m-1) + A001653(m)*A001652(n-m-1) - a(m) = A002315(m)*A053141(n-m-1) + A001653(m)*A046090(n-m-1) + a(m) with m < n; otherwise a(n) = A002315(m)*A053141(m-n) - A001653(m)*A011900(m-n) + a(m) = A002315(m)*A053141(m-n) - A001653(m)*A046090(m-n) - a(m) = (A002315(n) - A001653(n))/2. - Kenneth J Ramsey, Oct 12 2011
16*a(n)^2 + 1 = A056771(n). - James R. Buddenhagen, Dec 09 2011
A010054(A000290(a(n))) = 1. - Reinhard Zumkeller, Dec 17 2011
In general, a(n+k)^2 - A003499(k)*a(n+k)*a(n) + a(n)^2 = a(k)^2. - Charlie Marion, Jan 11 2012
a(n+1) = Sum_{k=0..n} A101950(n,k)*5^k. - Philippe Deléham, Feb 10 2012
PSUM transform of a(n+1) is A053142. PSUMSIGN transform of a(n+1) is A084158. BINOMIAL transform of a(n+1) is A164591. BINOMIAL transform of A086347 is a(n+1). BINOMIAL transform of A057087(n-1). - Michael Somos, May 11 2012
a(n+k) = A001541(k)*a(n) + sqrt(A132592(k)*a(n)^2 + a(k)^2). Generalizes formula dated Oct 09 2000. - Charlie Marion, Nov 27 2012
a(n) + a(n+2*k) = A003499(k)*a(n+k); a(n) + a(n+2*k+1) = A001653(k+1)*A002315(n+k). - Charlie Marion, Nov 29 2012
From Peter Bala, Dec 23 2012: (Start)
Product_{n >= 1} (1 + 1/a(n)) = 1 + sqrt(2).
Product_{n >= 2} (1 - 1/a(n)) = (1/3)*(1 + sqrt(2)). (End)
G.f.: G(0)*x/(2-6*x), where G(k) = 1 + 1/(1 - x*(8*k-9)/( x*(8*k-1) - 3/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 12 2013
G.f.: H(0)*x/2, where H(k) = 1 + 1/( 1 - x*(6-x)/(x*(6-x) + 1/H(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Feb 18 2014
a(n) = (a(n-1)^2 - a(n-3)^2)/a(n-2) + a(n-4) for n > 3. - Patrick J. McNab, Jul 24 2015
a(n-k)*a(n+k) + a(k)^2 = a(n)^2, a(n+k) + a(n-k) = A003499(k)*a(n), for n >= k >= 0. - Alexander Samokrutov, Sep 30 2015
Dirichlet g.f.: (PolyLog(s,3+2*sqrt(2)) - PolyLog(s,3-2*sqrt(2)))/(4*sqrt(2)). - Ilya Gutkovskiy, Jun 27 2016
4*a(n)^2 - 1 = A278310(n) for n > 0. - Bruno Berselli, Nov 24 2016
From Klaus Purath, Jan 18 2020: (Start)
a(n) = (a(n-3) + a(n+3))/198.
a(n) = Sum_{i=1..n} A001653(i), n>=1.
a(n) = sinh( 2 * n * arccsch(1) ) / ( 2 * sqrt(2) ). - Federico Provvedi, Feb 01 2021
(End)
a(n) = A002965(2*n)*A002965(2*n+1). - Jon E. Schoenfield, Jan 08 2022
a(n) = A002965(4*n)/2. - Gerry Martens, Jul 14 2023
a(n) = Sum_{k = 0..n-1} (-1)^(n+k+1)*binomial(n+k, 2*k+1)*8^k. - Peter Bala, Jul 17 2023

Extensions

Additional comments from Wolfdieter Lang, Feb 10 2000
Duplication of a formula removed by Wolfdieter Lang, Feb 10 2015

A156035 Decimal expansion of 3 + 2*sqrt(2).

Original entry on oeis.org

5, 8, 2, 8, 4, 2, 7, 1, 2, 4, 7, 4, 6, 1, 9, 0, 0, 9, 7, 6, 0, 3, 3, 7, 7, 4, 4, 8, 4, 1, 9, 3, 9, 6, 1, 5, 7, 1, 3, 9, 3, 4, 3, 7, 5, 0, 7, 5, 3, 8, 9, 6, 1, 4, 6, 3, 5, 3, 3, 5, 9, 4, 7, 5, 9, 8, 1, 4, 6, 4, 9, 5, 6, 9, 2, 4, 2, 1, 4, 0, 7, 7, 7, 0, 0, 7, 7, 5, 0, 6, 8, 6, 5, 5, 2, 8, 3, 1, 4, 5, 4, 7, 0, 0, 2
Offset: 1

Views

Author

Klaus Brockhaus, Feb 02 2009

Keywords

Comments

Limit_{n -> oo} b(n+1)/b(n) = 3+2*sqrt(2) for b = A155464, A155465, A155466.
Limit_{n -> oo} b(n)/b(n-1) = 3+2*sqrt(2) for b = A001652, A001653, A002315, A156156, A156157, A156158. - Klaus Brockhaus, Sep 23 2009
From Richard R. Forberg, Aug 14 2013: (Start)
Ratios b(n+1)/b(n) for all sequences of the form b(n) = 6*b(n-1) - b(n-2), for any initial values of b(0) and b(1), converge to this ratio.
Ratios b(n+1)/b(n) for all sequences of the form b(n) = 5*b(n-1) + 5*b(n-2) + b(n-3), for all b(0), b(1) and b(2) also converge to 3 + 2*sqrt(2). For example see A084158 (Pell Triangles).
Ratios of alternating values, b(n+2)/b(n), for all sequences of the form b(n) = 2*b(n-1) + b(n-2), also converge to 3 + 2*sqrt(2). These include A000129 (Pell Numbers). Also see A014176. (End)
Let ABCD be a square inscribed in a circle. When P is the midpoint of the arc AB, then the ratio (PC*PD)/(PA*PB) is equal to 3+2*sqrt(2). See the Mathematical Reflections link. - Michel Marcus, Jan 10 2017
Limit of ratios of successive terms of A001652 when n-> infinity. - Harvey P. Dale, Jun 16 2017; improved by Bernard Schott, Feb 28 2022
A quadratic integer with minimal polynomial x^2 - 6x + 1. - Charles R Greathouse IV, Jul 11 2020
Ratio between radii of the large circumscribed circle R and the small internal circle r drawn on the Sangaku tablet at Isaniwa Jinjya shrine in Ehime Prefecture (pictures in links). - Bernard Schott, Feb 25 2022

Examples

			3 + 2*sqrt(2) = 5.828427124746190097603377448...
		

References

  • Diogo Queiros-Condé and Michel Feidt, Fractal and Trans-scale Nature of Entropy, Iste Press and Elsevier, 2018, page 45.

Crossrefs

Cf. A002193 (sqrt(2)), A090488, A010466, A014176.
Cf. A104178 (decimal expansion of log_10(3+2*sqrt(2))).
Cf. A242412 (sangaku).

Programs

Formula

Equals 1 + A090488 = 3 + A010466. - R. J. Mathar, Feb 19 2009
Equals exp(arccosh(3)), since arccosh(x) = log(x+sqrt(x^2-1)). - Stanislav Sykora, Nov 01 2013
Equals (1+sqrt(2))^2, that is, A014176^2. - Michel Marcus, May 08 2016
The periodic continued fraction is [5; [1, 4]]. - Stefano Spezia, Mar 17 2024

A054490 Expansion of (1+5*x)/(1-6*x+x^2).

Original entry on oeis.org

1, 11, 65, 379, 2209, 12875, 75041, 437371, 2549185, 14857739, 86597249, 504725755, 2941757281, 17145817931, 99933150305, 582453083899, 3394785353089, 19786259034635, 115322768854721, 672150354093691, 3917579355707425, 22833325780150859
Offset: 0

Views

Author

Barry E. Williams, May 04 2000

Keywords

Comments

A Pellian-related second-order recursive sequence.
Third binomial transform of 1,8,8,64,64,512. - Al Hakanson (hawkuu(AT)gmail.com), Aug 17 2009
Binomial transform of A164607. - R. J. Mathar, Oct 26 2011
Pisano period lengths: 1, 1, 4, 2, 6, 4, 3, 2, 12, 6, 12, 4, 14, 3, 12, 2, 8, 12, 20, 6, ... - R. J. Mathar, Aug 10 2012
From Wolfdieter Lang, Feb 26 2015: (Start)
This sequence gives all positive solutions x = x1 = a(n) of the first class of the (generalized) Pell equation x^2 - 2*y^2 = -7. For the corresponding y1 terms see 2*A038723(n). All positive solutions of the second class are given by (x2(n), y2(n)) = (A255236(n), A038725(n+1)), n >= 0. See (A254938(1), 2*A255232(1)) for the fundamental solution (1, 2) of the first class. See the Nagell reference, Theorem 111, p. 210, Theorem 110, p. 208, Theorem 108a, pp. 206-207.
This sequence also gives all positive solutions y = y1 of the first class of the Pell equation x^2 - 2*y^2 = 14. The corresponding solutions x1 are given in 4*A038723. This follows from the preceding comment. (End)
From Wolfdieter Lang, Mar 19 2015: (Start)
a(0) = -(2*A038761(0) - A038762(0)), a(n) = 2*A253811(n-1) + A101386(n-1), for n >= 1.
This follows from the general trivial fact that if X^2 - D*Y^2 = N (X, Y positive integers, D > 1, not a square, and N a non-vanishing integer) then x:= D*Y +/- X and y:= Y +/- X (correlated signs) satisfy x^2 - D*y^2 = -(D-1)*N. with integers x and y. Here D = 2 and N = 7. (End)

Examples

			n = 2: sqrt(8*23^2-7) = 65.
2*19 + 27  = 65. - _Wolfdieter Lang_, Mar 19 2015
		

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, N. Y., 1964, pp. 122-125, 194-196.
  • T. Nagell, Introduction to Number Theory, Chelsea Publishing Company, New York, 1964.

Crossrefs

Programs

  • GAP
    a:=[1,11];; for n in [3..30] do a[n]:=6*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Jan 20 2020
  • Magma
    I:=[1,11]; [n le 2 select I[n] else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Mar 20 2015
    
  • Maple
    a[0]:=1: a[1]:=11: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..30); # Zerinvary Lajos, Jul 26 2006
  • Mathematica
    CoefficientList[Series[(1+5x)/(1-6x+x^2), {x, 0, 30}], x] (* Vincenzo Librandi, Mar 20 2015 *)
    LinearRecurrence[{6, -1}, {1, 11}, 30] (* G. C. Greubel, Jul 26 2018 *)
  • PARI
    my(x='x+O('x^30)); Vec((1+5*x)/(1-6*x+x^2)) \\ G. C. Greubel, Jul 26 2018
    
  • Sage
    [lucas_number1(2*n+1,2,-1) + 3*lucas_number1(2*n,2,-1) for n in (0..30)] # G. C. Greubel, Jan 20 2020
    

Formula

a(n) = 6*a(n-1) - a(n-2) for n>1, a(0)=1, a(1)=11.
a(n) = sqrt(8*A038723(n)^2 - 7).
a(n) = (11*((3+2*sqrt(2))^n - (3-2*sqrt(2))^n) - ((3+2*sqrt(2))^(n-1) - (3-2*sqrt(2))^(n-1)))/(4*sqrt(2)).
a(n) = 11*S(n, 6) + 5*S(n-1, 6), n >= 0, with Chebyshev's polynomials S(n, x) (A049310) evaluated at x=6: S(n, 6) = A001109(n-1). See the g.f. and the Pell equation comments above. - Wolfdieter Lang, Feb 26 2015
a(n) = 2*A253811(n-1) + A101386(n-1), for n >= 1. See the Mar 19 2015 comment above. - Wolfdieter Lang, Mar 19 2015
From G. C. Greubel, Jan 20 2020: (Start)
a(n) = Pell(2*n+1) + 3*Pell(2*n).
a(n) = ChebyshevU(n,3) + 5*ChebyshevU(n-1,3).
E.g.f.: exp(3*x)*( cosh(2*sqrt(2)*x) + 2*sqrt(2)*sinh(2*sqrt(2)*x) ). (End)

Extensions

More terms from James Sellers, May 05 2000
More terms from Vincenzo Librandi, Mar 20 2015

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

Original entry on oeis.org

1, 10, 59, 344, 2005, 11686, 68111, 396980, 2313769, 13485634, 78600035, 458114576, 2670087421, 15562409950, 90704372279, 528663823724, 3081278570065, 17959007596666, 104672767009931, 610077594462920
Offset: 0

Views

Author

Barry E. Williams, May 04 2000

Keywords

Comments

Numbers n such that 8*n^2 + 41 is a square.
(x, y) = (a(n), a(n+1)) are solutions to x^2 + y^2 - 6*x*y = 41. - John O. Oladokun, Mar 17 2021

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, N. Y., 1964, pp. 122-125, 194-196.

Crossrefs

Programs

  • GAP
    a:=[1,10];; for n in [3..30] do a[n]:=6*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Jan 19 2020
  • Magma
    I:=[1,10]; [n le 2 select I[n] else 6*Self(n-1) - Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 19 2020
    
  • Maple
    a[0]:=1: a[1]:=10: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..19); # Zerinvary Lajos, Jul 26 2006
  • Mathematica
    Table[(LucasL[2*n+1, 2] + 3*Fibonacci[2*n, 2])/2, {n,0,30}] (* G. C. Greubel, Jan 19 2020 *)
    LinearRecurrence[{6,-1},{1,10},20] (* Harvey P. Dale, Jun 11 2024 *)
  • PARI
    vector(31, n, polchebyshev(n-1,2,3) +4*polchebyshev(n-2,2,3) ) \\ G. C. Greubel, Jan 19 2020
    
  • Sage
    [chebyshev_U(n,3) +4*chebyshev_U(n-1,3) for n in (0..30)] # G. C. Greubel, Jan 19 2020
    

Formula

a(n) = 6*a(n-1) - a(n-2), a(0)=1, a(1)=10.
a(n) = (10*((3+2*sqrt(2))^n - (3-2*sqrt(2))^n) - ((3+2*sqrt(2))^(n-1) - (3-2*sqrt(2))^(n-1)))/(4*sqrt(2)).
From G. C. Greubel, Jan 19 2020: (Start)
a(n) = ChebyshevU(n,3) + 4*ChebyshevU(n-1,3).
a(n) = (Pell(2*n+2) + 4*Pell(2*n))/2 = (Pell-Lucas(2*n+1) + 3*Pell(2*n))/2.
E.g.f.: exp(3*x)*( cosh(2*sqrt(2)*x) + 7*sinh(2*sqrt(2)*x)/(2*sqrt(2)) ). (End)

Extensions

More terms from James Sellers, May 05 2000

A077413 Bisection (odd part) of Chebyshev sequence with Diophantine property.

Original entry on oeis.org

2, 13, 76, 443, 2582, 15049, 87712, 511223, 2979626, 17366533, 101219572, 589950899, 3438485822, 20040964033, 116807298376, 680802826223, 3968009658962, 23127255127549, 134795521106332, 785645871510443, 4579079707956326, 26688832376227513, 155553914549408752
Offset: 0

Views

Author

Wolfdieter Lang, Nov 08 2002

Keywords

Comments

-8*a(n)^2 + b(n)^2 = 17, with the companion sequence b(n) = A077239(n).
The even part is A054488(n) with Diophantine companion A077240(n).

Examples

			8*a(1)^2 + 17 = 8*13^2+17 = 1369 = 37^2 = A077239(1)^2.
G.f. = 2 + 13*x + 76*x^2 + 443*x^3 + 2582*x^4 + ... - _Michael Somos_, Jul 30 2024
		

Crossrefs

Cf. A077241 (even and odd parts), A001109, A054488.

Programs

  • Magma
    I:=[2,13]; [n le 2 select I[n] else 6*Self(n-1) - Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 18 2018
  • Mathematica
    LinearRecurrence[{6,-1}, {2,13}, 30] (* or *) CoefficientList[Series[ (2+x)/(1-6*x+x^2), {x, 0, 50}], x] (* G. C. Greubel, Jan 18 2018 *)
    a[ n_] := 2*ChebyshevU[n, 3] + ChebyshevU[n-1, 3]; (* Michael Somos, Jul 30 2024 *)
  • PARI
    Vec((2+x)/(1-6*x+x^2) + O(x^30)) \\ Colin Barker, Jun 16 2015
    
  • PARI
    {a(n) = 2*polchebyshev(n, 2, 3) + polchebyshev(n-1, 2, 3)}; /* Michael Somos, Jul 30 2024 */
    
  • PARI
    {a(n) = my(w = 3 + quadgen(32)); imag(w^n + 2*w^(n+1))}; /* Michael Somos, Jul 30 2024 */
    

Formula

a(n) = 6*a(n-1) - a(n-2), a(-1)=-1, a(0)=2.
a(n) = 2*S(n, 6)+S(n-1, 6), with S(n, x) = U(n, x/2), Chebyshev polynomials of 2nd kind, A049310. S(n, 6) = A001109(n+1).
G.f.: (2+x)/(1-6*x+x^2).
a(n) = (((3-2*sqrt(2))^n*(-7+4*sqrt(2))+(3+2*sqrt(2))^n*(7+4*sqrt(2))))/(4*sqrt(2)). - Colin Barker, Oct 12 2015
a(n) = -A054488(-1-n) for all n in Z. - Michael Somos, Jul 30 2024

A077239 Bisection (odd part) of Chebyshev sequence with Diophantine property.

Original entry on oeis.org

7, 37, 215, 1253, 7303, 42565, 248087, 1445957, 8427655, 49119973, 286292183, 1668633125, 9725506567, 56684406277, 330380931095, 1925601180293, 11223226150663, 65413755723685, 381259308191447, 2222142093424997, 12951593252358535, 75487417420726213
Offset: 0

Views

Author

Wolfdieter Lang, Nov 08 2002

Keywords

Comments

a(n)^2 - 8*b(n)^2 = 17, with the companion sequence b(n)= A077413(n).
The even part is A077240(n) with Diophantine companion A054488(n).

Examples

			37 = a(1) = sqrt(8*A077413(1)^2 +17) = sqrt(8*13^2 + 17)= sqrt(1369) = 37.
		

Crossrefs

Cf. A077242 (even and odd parts).

Programs

  • Mathematica
    Table[2*ChebyshevT[n+1, 3] + ChebyshevT[n, 3], {n, 0, 19}]  (* Jean-François Alcover, Dec 19 2013 *)
  • PARI
    Vec((7-5*x)/(1-6*x+x^2) + O(x^40)) \\ Colin Barker, Oct 12 2015

Formula

a(n) = 6*a(n-1) - a(n-2), a(-1) := 5, a(0)=7.
a(n) = 2*T(n+1, 3)+T(n, 3), with T(n, x) Chebyshev's polynomials of the first kind, A053120. T(n, 3)= A001541(n).
G.f.: (7-5*x)/(1-6*x+x^2).
a(n) = (((3-2*sqrt(2))^n*(-8+7*sqrt(2))+(3+2*sqrt(2))^n*(8+7*sqrt(2))))/(2*sqrt(2)). - Colin Barker, Oct 12 2015

A077240 Bisection (even part) of Chebyshev sequence with Diophantine property.

Original entry on oeis.org

5, 23, 133, 775, 4517, 26327, 153445, 894343, 5212613, 30381335, 177075397, 1032071047, 6015350885, 35060034263, 204344854693, 1191009093895, 6941709708677, 40459249158167, 235813785240325, 1374423462283783, 8010726988462373, 46689938468490455
Offset: 0

Views

Author

Wolfdieter Lang, Nov 08 2002

Keywords

Comments

a(n)^2 - 8*b(n)^2 = 17, with the companion sequence b(n)= A054488(n).
The odd part is A077239(n) with Diophantine companion A077413(n).

Examples

			23 = a(1) = sqrt(8*A054488(1)^2 + 17) = sqrt(8*8^2 + 17)= sqrt(529) = 23.
		

Crossrefs

Cf. A077242 (even and odd parts).

Programs

  • Mathematica
    Table[ChebyshevT[n+1, 3] + 2*ChebyshevT[n, 3], {n, 0, 19}]  (* Jean-François Alcover, Dec 19 2013 *)
    LinearRecurrence[{6,-1},{5,23},30] (* Harvey P. Dale, Mar 29 2017 *)
  • PARI
    Vec((5-7*x)/(1-6*x+x^2) + O(x^40)) \\ Colin Barker, Oct 12 2015

Formula

a(n) = 6*a(n-1) - a(n-2), a(-1) = 7, a(0) = 5.
a(n) = T(n+1, 3)+2*T(n, 3), with T(n, x) Chebyshev's polynomials of the first kind, A053120. T(n, 3)= A001541(n).
G.f.: (5-7*x)/(1-6*x+x^2).
a(n) = (((3-2*sqrt(2))^n*(-4+5*sqrt(2))+(3+2*sqrt(2))^n*(4+5*sqrt(2))))/(2*sqrt(2)). - Colin Barker, Oct 12 2015
Showing 1-9 of 9 results.