A046729 -id:A046729 - cp's OEIS frontend

A002315 NSW numbers: a(n) = 6a(n-1) - a(n-2); also a(n)^2 - 2b(n)^2 = -1 with b(n) = A001653(n+1).

1, 7, 41, 239, 1393, 8119, 47321, 275807, 1607521, 9369319, 54608393, 318281039, 1855077841, 10812186007, 63018038201, 367296043199, 2140758220993, 12477253282759, 72722761475561, 423859315570607, 2470433131948081, 14398739476117879, 83922003724759193

Offset: 0

N. J. A. Sloane

Named after the Newman-Shanks-Williams reference.
Also numbers k such that A125650(3*k^2) is an odd perfect square. Such numbers 3*k^2 form a bisection of A125651. - Alexander Adamchuk, Nov 30 2006
For positive n, a(n) corresponds to the sum of legs of near-isosceles primitive Pythagorean triangles (with consecutive legs). - Lekraj Beedassy, Feb 06 2007
Also numbers m such that m^2 is a centered 16-gonal number; or a number of the form 8k(k+1)+1, where k = A053141(m) = {0, 2, 14, 84, 492, 2870, ...}. - Alexander Adamchuk, Apr 21 2007
The lower principal convergents to 2^(1/2), beginning with 1/1, 7/5, 41/29, 239/169, comprise a strictly increasing sequence; numerators=A002315 and denominators=A001653. - Clark Kimberling, Aug 27 2008
The upper intermediate convergents to 2^(1/2) beginning with 10/7, 58/41, 338/239, 1970/1393 form a strictly decreasing sequence; essentially, numerators=A075870, denominators=A002315. - Clark Kimberling, Aug 27 2008
General recurrence is a(n) = (a(1)-1)*a(n-1) - a(n-2), a(1) >= 4, lim_{n->oo} a(n) = x*(k*x+1)^n, k = (a(1)-3), x = (1+sqrt((a(1)+1)/(a(1)-3)))/2. Examples in OEIS: a(1)=4 gives A002878. a(1)=5 gives A001834. a(1)=6 gives A030221. a(1)=7 gives A002315. a(1)=8 gives A033890. a(1)=9 gives A057080. a(1)=10 gives A057081. - Ctibor O. Zizka, Sep 02 2008
Numbers k such that (ceiling(sqrt(k*k/2)))^2 = (1+k*k)/2. - Ctibor O. Zizka, Nov 09 2009
A001109(n)/a(n) converges to cos^2(Pi/8) = 1/2 + 2^(1/2)/4. - Gary Detlefs, Nov 25 2009
The values 2(a(n)^2+1) are all perfect squares, whose square root is given by A075870. - Neelesh Bodas (neelesh.bodas(AT)gmail.com), Aug 13 2010
a(n) represents all positive integers K for which 2(K^2+1) is a perfect square. - Neelesh Bodas (neelesh.bodas(AT)gmail.com), Aug 13 2010
For positive n, a(n) equals the permanent of the (2n) X (2n) tridiagonal matrix with sqrt(8)'s along the main diagonal, and i's along the superdiagonal and subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
Integers k such that A000217(k-2) + A000217(k-1) + A000217(k) + A000217(k+1) is a square (cf. A202391). - Max Alekseyev, Dec 19 2011
Integer square roots of floor(k^2/2 - 1) or A047838. - Richard R. Forberg, Aug 01 2013
Remark: x^2 - 2*y^2 = +2*k^2, with positive k, and X^2 - 2*Y^2 = +2 reduce to the present Pell equation a^2 - 2*b^2 = -1 with x = k*X = 2*k*b and y = k*Y = k*a. (After a proposed solution for k = 3 by Alexander Samokrutov.) - Wolfdieter Lang, Aug 21 2015
If p is an odd prime, a((p-1)/2) == 1 (mod p). - Altug Alkan, Mar 17 2016
a(n)^2 + 1 = 2*b(n)^2, with b(n) = A001653(n), is the necessary and sufficient condition for a(n) to be a number k for which the diagonal of a 1 X k rectangle is an integer multiple of the diagonal of a 1 X 1 square. If squares are laid out thus along one diagonal of a horizontal 1 X a(n) rectangle, from the lower left corner to the upper right, the number of squares is b(n), and there will always be a square whose top corner lies exactly within the top edge of the rectangle. Numbering the squares 1 to b(n) from left to right, the number of the one square that has a corner in the top edge of the rectangle is c(n) = (2*b(n) - a(n) + 1)/2, which is A055997(n). The horizontal component of the corner of the square in the edge of the rectangle is also an integer, namely d(n) = a(n) - b(n), which is A001542(n). - David Pasino, Jun 30 2016
(a(n)^2)-th triangular number is a square; a(n)^2 = A008843(n) is a subsequence of A001108. - Jaroslav Krizek, Aug 05 2016
a(n-1)/A001653(n) is the closest rational approximation of sqrt(2) with a numerator not larger than a(n-1). These rational approximations together with those obtained from the sequences A001541 and A001542 give a complete set of closest rational approximations of sqrt(2) with restricted numerator or denominator. a(n-1)/A001653(n) < sqrt(2). - A.H.M. Smeets, May 28 2017
Consider the quadrant of a circle with center (0,0) bounded by the positive x and y axes. Now consider, as the start of a series, the circle contained within this quadrant which kisses both axes and the outer bounding circle. Consider further a succession of circles, each kissing the x-axis, the outer bounding circle, and the previous circle in the series. See Holmes link. The center of the n-th circle in this series is ((A001653(n)*sqrt(2)-1)/a(n-1), (A001653(n)*sqrt(2)-1)/a(n-1)^2), the y-coordinate also being its radius. It follows that a(n-1) is the cotangent of the angle subtended at point (0,0) by the center of the n-th circle in the series with respect to the x-axis. - Graham Holmes, Aug 31 2019
There is a link between the two sequences present at the numerator and at the denominator of the fractions that give the coordinates of the center of the kissing circles. A001653 is the sequence of numbers k such that 2*k^2 - 1 is a square, and here, we have 2*A001653(n)^2 - 1 = a(n-1)^2. - Bernard Schott, Sep 02 2019
Let G be a sequence satisfying G(i) = 2*G(i-1) + G(i-2) for arbitrary integers i and without regard to the initial values of G. Then a(n) = (G(i+4*n+2) - G(i))/(2*G(i+2*n+1)) as long as G(i+2*n+1) != 0. - Klaus Purath, Mar 25 2021
All of the positive integer solutions of a*b+1=x^2, a*c+1=y^2, b*c+1=z^2, x+z=2*y, 0 < a < b < c are given by a=A001542(n), b=A005319(n), c=A001542(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0 < n. - Michael Somos, Jun 26 2022
3*a(n-1) is the n-th almost Lucas-cobalancing number of second type (see Tekcan and Erdem). - Stefano Spezia, Nov 26 2022
In Moret-Blanc (1881) on page 259 some solution of m^2 - 2n^2 = -1 are listed. The values of m give this sequence, and the values of n give A001653. - Michael Somos, Oct 25 2023
From Klaus Purath, May 11 2024: (Start)
For any two consecutive terms (a(n), a(n+1)) = (x,y): x^2 - 6xy + y^2 = 8 = A028884(1). In general, the following applies to all sequences (t) satisfying t(i) = 6t(i-1) - t(i-2) with t(0) = 1 and two consecutive terms (x,y): x^2 - 6xy + y^2 = A028884(t(1)-6). This includes and interprets the Feb 04 2014 comment on A001541 by Colin Barker as well as the Mar 17 2021 comment on A054489 by John O. Oladokun and the Sep 28 2008 formula on A038723 by Michael Somos. By analogy to this, for three consecutive terms (x,y,z) y^2 - xz = A028884(t(1)-6) always applies.
If (t) is a sequence satisfying t(k) = 7t(k-1) - 7t(k-2) + t(k-3) or t(k) = 6t(k-1) - t(k-2) without regard to initial values and including this sequence itself, then a(n) = (t(k+2n+1) - t(k))/(t(k+n+1) - t(k+n)) always applies, as long as t(k+n+1) - t(k+n) != 0 for integer k and n >= 0. (End)

			G.f. = 1 + 7*x + 41*x^2 + 239*x^3 + 1393*x^4 + 8119*x^5 + 17321*x^6 + ... - _Michael Somos_, Jun 26 2022

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)
A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 256.
Paulo Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 288.
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 247.
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).
P.-F. Teilhet, Reply to Query 2094, L'Intermédiaire des Mathématiciens, 10 (1903), 235-238.
P.-F. Teilhet, Query 2376, L'Intermédiaire des Mathématiciens, 11 (1904), 138-139. - N. J. A. Sloane, Mar 08 2022

Indranil Ghosh, Table of n, a(n) for n = 0..1303 (terms 0..200 from T. D. Noe)
Marco Abrate, Stefano Barbero, Umberto Cerruti, and Nadir Murru, Polynomial sequences on quadratic curves, Integers, Vol. 15, 2015, #A38.
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3, pp. 245-278, 2011. See Section 9.
E. Barcucci et al., A combinatorial interpretation of the recurrence f_{n+1} = 6 f_n - f_{n-1}, Discrete Math., 190 (1998), 235-240.
Elena Barcucci, Antonio Bernini, and Renzo Pinzani, A Gray code for a regular language, Semantic Sensor Networks Workshop 2018, CEUR Workshop Proceedings (2018) Vol. 2113.
Hacène Belbachir and Yassine Otmani, Quadrinomial-Like Versions for Wolstenholme, Morley and Glaisher Congruences, Integers (2023) Vol. 23.
J. Bonin, L. Shapiro and R. Simion, Some q-analogues of the Schroeder numbers arising from combinatorial statistics on lattice paths, H. Statistical Planning and Inference, 16, 1993, 35-55 (p. 50).
P. Catarino, H. Campos, and P. Vasco, On some identities for balancing and cobalancing numbers, Annales Mathematicae et Informaticae, 45 (2015) pp. 11-24.
Enrica Duchi, Andrea Frosini, Renzo Pinzani and Simone Rinaldi, A Note on Rational Succession Rules, J. Integer Seqs., Vol. 6, 2003.
Melissa Emory, The Diophantine equation X^4 + Y^4 = D^2 Z^4 in quadratic fields, INTEGERS 12 (2012), #A65. - _N. J. A. Sloane_, Feb 06 2013
S. Falcon, Relationships between Some k-Fibonacci Sequences, Applied Mathematics, 2014, 5, 2226-2234.
Alex Fink, Richard K. Guy, and Mark Krusemeyer, Partitions with parts occurring at most thrice, Contributions to Discrete Mathematics, Vol 3, No 2 (2008), pp. 76-114. See Section 13.
A. S. Fraenkel, On the recurrence f(m+1)= b(m)*f(m)-f(m-1) and applications, Discrete Mathematics 224 (2000), pp. 273-279.
A. S. Fraenkel, Recent results and questions in combinatorial game complexities, Theoretical Computer Science, vol. 249, no. 2 (2000), 265-288.
A. S. Fraenkel, Arrays, numeration systems and Frankenstein games, Theoret. Comput. Sci. 282 (2002), 271-284.
Bernard Frénicle de Bessy, Solutio duorum problematum circa numeros cubos et quadratos, (1657), page 9. Bibliothèque Nationale de Paris.
M. A. Gruber, Artemas Martin, A. H. Bell, J. H. Drummond, A. H. Holmes and H. C. Wilkes, Problem 47, Amer. Math. Monthly, 4 (1897), 25-28.
R. J. Hetherington, Letter to N. J. A. Sloane, Oct 26 1974.
Graham Holmes, Kissing circles and cotangents.
Tanya Khovanova, Recursive Sequences.
Ioana-Claudia Lazăr, Lucas sequences in t-uniform simplicial complexes, arXiv:1904.06555 [math.GR], 2019.
D. H. Lehmer, Lacunary recurrence formulas for the numbers of Bernoulli and Euler, Annals Math., 36 (1935), 637-649.
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Hyperbola, Forum Geometricorum (2019) Vol. 19, 11-16.
aBa Mbirika, Janeè Schrader, and Jürgen Spilker, Pell and associated Pell braid sequences as GCDs of sums of k consecutive Pell, balancing, and related numbers, arXiv:2301.05758 [math.NT], 2023. See also J. Int. Seq. (2023) Vol. 26, Art. 23.6.4.
Donatella Merlini and Renzo Sprugnoli, Arithmetic into geometric progressions through Riordan arrays, Discrete Mathematics 340.2 (2017): 160-174.
Claude Moret-Blanc, Questions Nouvelles d'Arithmétique Supérieure Proposées Par M. Edouard Lucas, Nouvelles annales de mathématiques 2^e serie, tome 20 (1881), 253-265.
Morris Newman, Daniel Shanks, and H. C. Williams, Simple groups of square order and an interesting sequence of primes, Acta Arith., 38 (1980/1981) 129-140.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
The Prime Glossary, NSW number.
S. F. Santana and J. L. Diaz-Barrero, Some properties of sums involving Pell numbers, Missouri Journal of Mathematical Sciences 18(1), 2006.
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1.
R. A. Sulanke, Bijective recurrences concerning Schroeder paths, Electron. J. Combin. 5 (1998), Research Paper 47, 11 pp.
R. A. Sulanke, Moments of generalized Motzkin paths, J. Integer Sequences, Vol. 3 (2000), #00.1.
Ahmet Tekcan and Alper Erdem, General Terms of All Almost Balancing Numbers of First and Second Type, arXiv:2211.08907 [math.NT], 2022.
Eric Weisstein's World of Mathematics, NSW Number.
Eric Weisstein's World of Mathematics, Centered Polygonal Number.
H. C. Williams and R. K. Guy, Some fourth-order linear divisibility sequences, Intl. J. Number Theory 7 (5) (2011) 1255-1277.
H. C. Williams and R. K. Guy, Some Monoapparitic Fourth Order Linear Divisibility Sequences, Integers, Volume 12A (2012) The John Selfridge Memorial Volume.
Index entries for sequences related to Chebyshev polynomials.
Index entries for two-way infinite sequences.
Index entries for linear recurrences with constant coefficients, signature (6,-1).

Bisection of A001333. Cf. A001109, A001653. A065513(n)=a(n)-1.
First differences of A001108 and A055997. Bisection of A084068 and A088014. Cf. A077444.
Cf. A125650, A125651, A125652.
Row sums of unsigned triangle A127675.
Cf. A053141, A075870. Cf. A000045, A002878, A004146, A026003, A100047, A119915, A192425, A088165 (prime subsequence), A057084 (binomial transform), A108051 (inverse binomial transform).
See comments in A301383.
Cf. similar sequences of the type (1/k)*sinh((2*n+1)*arcsinh(k)) listed in A097775.
Cf. A000129, A001653, A003499, A358682, A002203.

a002315 n = a002315_list !! n
a002315_list = 1 : 7 : zipWith (-) (map (* 6) (tail a002315_list)) a002315_list
-- Reinhard Zumkeller, Jan 10 2012

I:=[1,7]; [n le 2 select I[n] else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Mar 22 2015

A002315 := proc(n)
    option remember;
    if n = 0 then
        1 ;
    elif n = 1 then
        7;
    else
        6*procname(n-1)-procname(n-2) ;
    end if;
end proc: # Zerinvary Lajos, Jul 26 2006, modified R. J. Mathar, Apr 30 2017
a:=n->abs(Im(simplify(ChebyshevT(2*n+1,I)))):seq(a(n),n=0..20); # Leonid Bedratyuk, Dec 17 2017
# third Maple program:
a:= n-> (<<0|1>, <-1|6>>^n. <<1, 7>>)[1, 1]:
seq(a(n), n=0..22);  # Alois P. Heinz, Aug 25 2024

a[0] = 1; a[1] = 7; a[n_] := a[n] = 6a[n - 1] - a[n - 2]; Table[ a[n], {n, 0, 20}] (* Robert G. Wilson v, Jun 09 2004 *)
Transpose[NestList[Flatten[{Rest[#],ListCorrelate[{-1,6},#]}]&, {1,7},20]][[1]]  (* Harvey P. Dale, Mar 23 2011 *)
Table[ If[n>0, a=b; b=c; c=6b-a, b=-1; c=1], {n, 0, 20}] (* Jean-François Alcover, Oct 19 2012 *)
LinearRecurrence[{6, -1}, {1, 7}, 20] (* Bruno Berselli, Apr 03 2018 *)
a[ n_] := -I*(-1)^n*ChebyshevT[2*n + 1, I]; (* Michael Somos, Jun 26 2022 *)

{a(n) = subst(poltchebi(abs(n+1)) - poltchebi(abs(n)), x, 3)/2};

{a(n) = if(n<0, -a(-1-n), polsym(x^2-2*x-1, 2*n+1)[2*n+2]/2)};

{a(n) = my(w=3+quadgen(32)); imag((1+w)*w^n)};

for (i=1,10000,if(Mod(sigma(i^2+1,2),2)==1,print1(i,",")))

{a(n) = -I*(-1)^n*polchebyshev(2*n+1, 1, I)}; /* Michael Somos, Jun 26 2022 */

a(n) = (1/2)*((1+sqrt(2))^(2*n+1) + (1-sqrt(2))^(2*n+1)).
a(n) = A001109(n)+A001109(n+1).
a(n) = (1+sqrt(2))/2*(3+sqrt(8))^n+(1-sqrt(2))/2*(3-sqrt(8))^n. - Ralf Stephan, Feb 23 2003
a(n) = sqrt(2*(A001653(n+1))^2-1), n >= 0. [Pell equation a(n)^2 - 2*Pell(2*n+1)^2 = -1. - Wolfdieter Lang, Jul 11 2018]
G.f.: (1 + x)/(1 - 6*x + x^2). - Simon Plouffe in his 1992 dissertation
a(n) = S(n, 6)+S(n-1, 6) = S(2*n, sqrt(8)), S(n, x) = U(n, x/2) are Chebyshev's polynomials of the 2nd kind. Cf. A049310. S(n, 6)= A001109(n+1).
a(n) ~ (1/2)*(sqrt(2) + 1)^(2*n+1). - 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 06 2002
Let q(n, x) = Sum_{i=0..n} x^(n-i)*binomial(2*n-i, i); then (-1)^n*q(n, -8) = a(n). - Benoit Cloitre, Nov 10 2002
With a=3+2*sqrt(2), b=3-2*sqrt(2): a(n) = (a^((2n+1)/2)-b^((2n+1)/2))/2. a(n) = A077444(n)/2. - Mario Catalani (mario.catalani(AT)unito.it), Mar 31 2003
a(n) = Sum_{k=0..n} 2^k*binomial(2*n+1, 2*k). - Zoltan Zachar (zachar(AT)fellner.sulinet.hu), Oct 08 2003
Same as: i such that sigma(i^2+1, 2) mod 2 = 1. - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 26 2004
a(n) = L(n, -6)*(-1)^n, where L is defined as in A108299; see also A001653 for L(n, +6). - Reinhard Zumkeller, Jun 01 2005
a(n) = A001652(n)+A046090(n); e.g., 239=119+120. - Charlie Marion, Nov 20 2003
A001541(n)*a(n+k) = A001652(2n+k) + A001652(k)+1; e.g., 3*1393 = 4069 + 119 + 1; for k > 0, A001541(n+k)*a(n) = A001652(2n+k) - A001652(k-1); e.g., 99*7 = 696 - 3. - Charlie Marion, Mar 17 2003
a(n) = Jacobi_P(n,1/2,-1/2,3)/Jacobi_P(n,-1/2,1/2,1). - Paul Barry, Feb 03 2006
P_{2n}+P_{2n+1} where P_i are the Pell numbers (A000129). Also the square root of the partial sums of Pell numbers: P_{2n}+P_{2n+1} = sqrt(Sum_{i=0..4n+1} P_i) (Santana and Diaz-Barrero, 2006). - David Eppstein, Jan 28 2007
a(n) = 2*A001652(n) + 1 = 2*A046729(n) + (-1)^n. - Lekraj Beedassy, Feb 06 2007
a(n) = sqrt(A001108(2*n+1)). - Anton Vrba (antonvrba(AT)yahoo.com), Feb 14 2007
a(n) = sqrt(8*A053141(n)*(A053141(n) + 1) + 1). - Alexander Adamchuk, Apr 21 2007
a(n+1) = 3*a(n) + sqrt(8*a(n)^2 + 8), a(1)=1. - Richard Choulet, Sep 18 2007
a(n) = A001333(2*n+1). - Ctibor O. Zizka, Aug 13 2008
a(n) = third binomial transform of 1, 4, 8, 32, 64, 256, 512, ... . - Al Hakanson (hawkuu(AT)gmail.com), Aug 15 2009
a(n) = (-1)^(n-1)*(1/sqrt(-1))*cos((2*n - 1)*arcsin(sqrt(2)). - Artur Jasinski, Feb 17 2010 *WRONG*
a(n+k) = A001541(k)*a(n) + 4*A001109(k)*A001653(n); e.g., 8119 = 17*239 + 4*6*169. - Charlie Marion, Feb 04 2011
In general, a(n+k) = A001541(k)*a(n)) + sqrt(A001108(2k)*(a(n)^2+1)). See Sep 18 2007 entry above. - Charlie Marion, Dec 07 2011
a(n) = floor((1+sqrt(2))^(2n+1))/2. - Thomas Ordowski, Jun 12 2012
(a(2n-1) + a(2n) + 8)/(8*a(n)) = A001653(n). - Ignacio Larrosa Cañestro, Jan 02 2015
(a(2n) + a(2n-1))/a(n) = 2*sqrt(2)*( (1 + sqrt(2))^(4*n) - (1 - sqrt(2))^(4*n))/((1 + sqrt(2))^(2*n+1) + (1 - sqrt(2))^(2*n+1)). [This was my solution to problem 5325, School Science and Mathematics 114 (No. 8, Dec 2014).] - Henry Ricardo, Feb 05 2015
From Peter Bala, Mar 22 2015: (Start)
The aerated sequence (b(n))n>=1 = [1, 0, 7, 0, 41, 0, 239, 0, ...] is a fourth-order linear divisibility sequence; that is, if n | m then b(n) | b(m). It is the case P1 = 0, P2 = -4, Q = -1 of the 3-parameter family of divisibility sequences found by Williams and Guy. See A100047.
b(n) = 1/2*((-1)^n - 1)*Pell(n) + 1/2*(1 + (-1)^(n+1))*Pell(n+1). The o.g.f. is x*(1 + x^2)/(1 - 6*x^2 + x^4).
Exp( Sum_{n >= 1} 2*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 2*A026003(n-1)*x^n.
Exp( Sum_{n >= 1} (-2)*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 2*A026003(n-1)*(-x)^n.
Exp( Sum_{n >= 1} 4*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 4*Pell(n)*x^n.
Exp( Sum_{n >= 1} (-4)*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 4*Pell(n)*(-x)^n.
Exp( Sum_{n >= 1} 8*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 8*A119915(n)*x^n.
Exp( Sum_{n >= 1} (-8)*b(n)*x^n/n ) = 1 + Sum_{n >= 1} 8*A119915(n)*(-x)^n. Cf. A002878, A004146, A113224, and A192425. (End)
E.g.f.: (sqrt(2)*sinh(2*sqrt(2)*x) + cosh(2*sqrt(2)*x))*exp(3*x). - Ilya Gutkovskiy, Jun 30 2016
a(n) = Sum_{k=0..n} binomial(n,k) * 3^(n-k) * 2^k * 2^ceiling(k/2). - David Pasino, Jul 09 2016
a(n) = A001541(n) + 2*A001542(n). - A.H.M. Smeets, May 28 2017
a(n+1) = 3*a(n) + 4*b(n), b(n+1) = 2*a(n) + 3*b(n), with b(n)=A001653(n). - Zak Seidov, Jul 13 2017
a(n) = |Im(T(2n-1,i))|, i=sqrt(-1), T(n,x) is the Chebyshev polynomial of the first kind, Im is the imaginary part of a complex number, || is the absolute value. - Leonid Bedratyuk, Dec 17 2017
a(n) = sinh((2*n + 1)*arcsinh(1)). - Bruno Berselli, Apr 03 2018
a(n) = 5*a(n-1) + A003499(n-1), a(0) = 1. - Ivan N. Ianakiev, Aug 09 2019
From Klaus Purath, Mar 25 2021: (Start)
a(n) = A046090(2*n)/A001541(n).
a(n+1)*a(n+2) = a(n)*a(n+3) + 48.
a(n)^2 + a(n+1)^2 = 6*a(n)*a(n+1) + 8.
a(n+1)^2 = a(n)*a(n+2) + 8.
a(n+1) = a(n) + 2*A001541(n+1).
a(n) = 2*A046090(n) - 1. (End)
3*a(n-1) = sqrt(8*b(n)^2 + 8*b(n) - 7), where b(n) = A358682(n). - Stefano Spezia, Nov 26 2022
a(n) = -(-1)^n - 2 + Sum_{i=0..n} A002203(i)^2. - Adam Mohamed, Aug 22 2024
From Peter Bala, May 09 2025: (Start)
a(n) = Dir(n, 3), where Dir(n, x) denotes the n-th row polynomial of the triangle A244419.
For arbitrary x, a(n+x)^2 - 6*a(n+x)*a(n+x+1) + a(n+x+1)^2 = 8 with a(n) := (1/2)*((1+sqrt(2))^(2*n+1) + (1-sqrt(2))^(2*n+1)) as above. The particular case x = 0 is noted above,
a(n+1/2) = sqrt(2) * A001542(n+1).
Sum_{n >= 1} (-1)^(n+1)/(a(n) - 1/a(n)) = 1/8 (telescoping series: for n >= 1, 1/(a(n) - 1/a(n)) = 1/A081554(n) + 1/A081554(n+1)).
Product_{n >= 1} (a(n) + 1)/(a(n) - 1) = sqrt(2) (telescoping product: Product_{n = 1..k} ((a(n) + 1)/(a(n) - 1))^2 = 2*(1 - 1/A055997(k+2))). (End)

A001541 a(0) = 1, a(1) = 3; for n > 1, a(n) = 6*a(n-1) - a(n-2).

Original entry on oeis.org

1, 3, 17, 99, 577, 3363, 19601, 114243, 665857, 3880899, 22619537, 131836323, 768398401, 4478554083, 26102926097, 152139002499, 886731088897, 5168247530883, 30122754096401, 175568277047523, 1023286908188737, 5964153172084899, 34761632124320657

Offset: 0

N. J. A. Sloane

Chebyshev polynomials of the first kind evaluated at 3.
This sequence gives the values of x in solutions of the Diophantine equation x^2 - 8*y^2 = 1, the corresponding values of y are in A001109. For n > 0, the ratios a(n)/A001090(n) may be obtained as convergents to sqrt(8): either successive convergents of [3; -6] or odd convergents of [2; 1, 4]. - Lekraj Beedassy, Sep 09 2003 [edited by Jon E. Schoenfield, May 04 2014]
Also gives solutions to the equation x^2 - 1 = floor(x*r*floor(x/r)) where r = sqrt(8). - Benoit Cloitre, Feb 14 2004
Appears to give all solutions greater than 1 to the equation: x^2 = ceiling(x*r*floor(x/r)) where r = sqrt(2). - Benoit Cloitre, Feb 24 2004
This sequence give numbers n such that (n-1)*(n+1)/2 is a perfect square. Remark: (i-1)*(i+1)/2 = (i^2-1)/2 = -1 = i^2 with i = sqrt(-1) so i is also in the sequence. - Pierre CAMI, Apr 20 2005
a(n) is prime for n = {1, 2, 4, 8}. Prime a(n) are {3, 17, 577, 665857}, which belong to A001601(n). a(2k-1) is divisible by a(1) = 3. a(4k-2) is divisible by a(2) = 17. a(8k-4) is divisible by a(4) = 577. a(16k-8) is divisible by a(8) = 665857. - Alexander Adamchuk, Nov 24 2006
The upper principal convergents to 2^(1/2), beginning with 3/2, 17/12, 99/70, 577/408, comprise a strictly decreasing sequence; essentially, numerators=A001541 and denominators=A001542. - Clark Kimberling, Aug 26 2008
Also index of sequence A082532 for which A082532(n) = 1. - Carmine Suriano, Sep 07 2010
Numbers n such that sigma(n-1) and sigma(n+1) are both odd numbers. - Juri-Stepan Gerasimov, Mar 28 2011
Also, numbers such that floor(a(n)^2/2) is a square: base 2 analog of A031149, A204502, A204514, A204516, A204518, A204520, A004275, A001075. - M. F. Hasler, Jan 15 2012
Numbers such that 2n^2 - 2 is a square. Also integer square roots of the expression 2*n^2 + 1, at values of n given by A001542. Also see A228405 regarding 2n^2 -+ 2^k generally for k >= 0. - Richard R. Forberg, Aug 20 2013
Values of x (or y) in the solutions to x^2 - 6xy + y^2 + 8 = 0. - Colin Barker, Feb 04 2014
Panda and Ray call the numbers in this sequence the Lucas-balancing numbers C_n (see references and links).
Partial sums of X or X+1 of Pythagorean triples (X,X+1,Z). - Peter M. Chema, Feb 03 2017
a(n)/A001542(n) is the closest rational approximation to sqrt(2) with a numerator not larger than a(n), and 2*A001542(n)/a(n) is the closest rational approximation to sqrt(2) with a denominator not larger than a(n). These rational approximations together with those obtained from the sequences A001653 and A002315 give a complete set of closest rational approximations to sqrt(2) with restricted numerator or denominator. a(n)/A001542(n) > sqrt(2) > 2*A001542(n)/a(n). - A.H.M. Smeets, May 28 2017
x = a(n), y = A001542(n) are solutions of the Diophantine equation x^2 - 2y^2 = 1 (Pell equation). x = 2*A001542(n), y = a(n) are solutions of the Diophantine equation x^2 - 2y^2 = -2. Both together give the set of fractional approximations for sqrt(2) obtained from limited fractions obtained from continued fraction representation to sqrt(2). - A.H.M. Smeets, Jun 22 2017
a(n) is the radius of the n-th circle among the sequence of circles generated as follows: Starting with a unit circle centered at the origin, every subsequent circle touches the previous circle as well as the two limbs of hyperbola x^2 - y^2 = 1, and lies in the region y > 0. - Kaushal Agrawal, Nov 10 2018
All of the positive integer solutions of a*b+1=x^2, a*c+1=y^2, b*c+1=z^2, x+z=2*y, 0A001542(n), b=A005319(n), c=A001542(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0Michael Somos, Jun 26 2022

			99^2 + 99^2 = 140^2 + 2. - _Carmine Suriano_, Jan 05 2015
G.f. = 1 + 3*x + 17*x^2 + 99*x^3 + 577*x^4 + 3363*x^5 + 19601*x^6 + 114243*x^7 + ...

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)
J. W. L. Glaisher, On Eulerian numbers (formulas, residues, end-figures), with the values of the first twenty-seven, Quarterly Journal of Mathematics, vol. 45, 1914, pp. 1-51.
G. K. Panda, Some fascinating properties of balancing numbers, In Proc. of Eleventh Internat. Conference on Fibonacci Numbers and Their Applications, Cong. Numerantium 194 (2009), 185-189.
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

T. D. Noe, Table of n, a(n) for n = 0..200
I. Adler, Three Diophantine equations - Part II, Fib. Quart., 7 (1969), 181-193.
Christian Aebi and Grant Cairns, Lattice Equable Parallelograms, arXiv:2006.07566 [math.NT], 2020.
Jean-Paul Allouche, Zeta-regularization of arithmetic sequences, EPJ Web of Conferences (2020) Vol. 244, 01008.
Hacène Belbachir, Soumeya Merwa Tebtoub, and László Németh, Ellipse Chains and Associated Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.8.5.
H. Brocard, Notes élémentaires sur le problème de Peel, Nouvelle Correspondance Mathématique, 4 (1878), 161-169.
John M. Campbell, An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences, arXiv preprint arXiv:1105.3399 [math.GM], 2011.
P. Catarino, H. Campos, and P. Vasco, On some identities for balancing and cobalancing numbers, Annales Mathematicae et Informaticae, 45 (2015) pp. 11-24.
Kwang-Wu Chen and Yu-Ren Pan, Greatest Common Divisors of Shifted Horadam Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.5.8.
S. Falcon, Relationships between Some k-Fibonacci Sequences, Applied Mathematics, 2014, 5, 2226-2234.
Morgan Fiebig, aBa Mbirika, and Jürgen Spilker, Period patterns, entry points, and orders in the Lucas sequences: theory and applications, arXiv:2408.14632 [math.NT], 2024. See p. 5.
Robert Frontczak, A Note on Hybrid Convolutions Involving Balancing and Lucas-Balancing Numbers, Applied Mathematical Sciences, Vol. 12, 2018, No. 25, 1201-1208.
Robert Frontczak, Sums of Balancing and Lucas-Balancing Numbers with Binomial Coefficients, International Journal of Mathematical Analysis (2018) Vol. 12, No. 12, 585-594.
Robert Frontczak, Powers of Balancing Polynomials and Some Consequences for Fibonacci Sums, International Journal of Mathematical Analysis (2019) Vol. 13, No. 3, 109-115.
Robert Frontczak, Identities for generalized balancing numbers, Notes on Number Theory and Discrete Mathematics (2019) Vol. 25, No. 2, 169-180.
Robert Frontczak and Taras Goy, Additional close links between balancing and Lucas-balancing polynomials, arXiv:2007.14048 [math.NT], 2020.
Robert Frontczak and Taras Goy, More Fibonacci-Bernoulli relations with and without balancing polynomials, arXiv:2007.14618 [math.NT], 2020.
Robert Frontczak and Taras Goy, Lucas-Euler relations using balancing and Lucas-balancing polynomials, arXiv:2009.09409 [math.NT], 2020.
O. Khadir, K. Liptai, and L. Szalay, On the Shifted Product of Binary Recurrences, J. Int. Seq. 13 (2010), 10.6.1.
J. M. Katri and D. R. Byrkit, Problem E1976, Amer. Math. Monthly, 75 (1968), 683-684.
Tanya Khovanova, Recursive Sequences
D. H. Lehmer, A cotangent analogue of continued fractions, Duke Math. J., 4 (1935), 323-340.
D. H. Lehmer, A cotangent analogue of continued fractions, Duke Math. J., 4 (1935), 323-340. [Annotated scanned copy]
D. H. Lehmer, Lacunary recurrence formulas for the numbers of Bernoulli and Euler, Annals Math., 36 (1935), 637-649.
Dino Lorenzini and Z. Xiang, Integral points on variable separated curves, Preprint 2016.
aBa Mbirika, Janeè Schrader, and Jürgen Spilker, Pell and associated Pell braid sequences as GCDs of sums of k consecutive Pell, balancing, and related numbers, arXiv:2301.05758 [math.NT], 2023. See also J. Int. Seq. (2023) Vol. 26, Art. 23.6.4.
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 C(n).
Robert Phillips, Polynomials of the form 1+4ke+4ke^2, 2008.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Prasanta K. Ray, Curious congruences for balancing numbers, Int. J. Contemp. Math. Sci. 7 (2012), 881-889.
Soumeya M. Tebtoub, Hacène Belbachir, and László Németh, Integer sequences and ellipse chains inside a hyperbola, Proceedings of the 1st International Conference on Algebras, Graphs and Ordered Sets (ALGOS 2020), hal-02918958 [math.cs], 17-18.
Bibhu Prasad Tripathy and Bijan Kumar Patel, On balancing and Lucas-balancing numbers expressible as product of two k-Fibonacci numbers, arXiv:2508.12238 [math.NT], 2025.
N. J. Wildberger, Pell's equation without irrational numbers, J. Int. Seq. 13 (2010), 10.4.3, Section 4.
Index entries for two-way infinite sequences
Index entries for linear recurrences with constant coefficients, signature (6,-1).
Index entries for sequences related to Chebyshev polynomials.

Bisection of A001333. A003499(n) = 2a(n).
Cf. A055997 = numbers n such that n(n-1)/2 is a square.
Row 1 of array A188645.
Cf. A046090, A001109, A053142, A084130, A001601, A056771, A077420, A005319, A082532, A001542.
Cf. A055792 (terms squared), A132592.

a001541 n = a001541_list !! (n-1)
a001541_list =
1 : 3 : zipWith (-) (map (* 6) $ tail a001541_list) a001541_list
-- Reinhard Zumkeller, Oct 06 2011
(Scheme, with memoization-macro definec)
(definec (A001541 n) (cond ((zero? n) 1) ((= 1 n) 3) (else (- (* 6 (A001541 (- n 1))) (A001541 (- n 2))))))
;; Antti Karttunen, Oct 04 2016

[n: n in [1..10000000] |IsSquare(8*(n^2-1))]; // Vincenzo Librandi, Nov 18 2010

a[0]:=1: a[1]:=3: 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
A001541:=-(-1+3*z)/(1-6*z+z**2); # Simon Plouffe in his 1992 dissertation

Table[Simplify[(1/2) (3 + 2 Sqrt[2])^n + (1/2) (3 - 2 Sqrt[2])^n], {n, 0, 20}] (* Artur Jasinski, Feb 10 2010 *)
a[ n_] := If[n == 0, 1, With[{m = Abs @ n}, m Sum[4^i Binomial[m + i, 2 i]/(m + i), {i, 0, m}]]]; (* Michael Somos, Jul 11 2011 *)
a[ n_] := ChebyshevT[ n, 3]; (* Michael Somos, Jul 11 2011 *)
LinearRecurrence[{6, -1}, {1, 3}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 12 2012 *)

{a(n) = real((3 + quadgen(32))^n)}; /* Michael Somos, Apr 07 2003 */

{a(n) = subst( poltchebi( abs(n)), x, 3)}; /* Michael Somos, Apr 07 2003 */

{a(n) = if( n<0, a(-n), polsym(1 - 6*x + x^2, n) [n+1] / 2)}; /* Michael Somos, Apr 07 2003 */

{a(n) = polchebyshev( n, 1, 3)}; /* Michael Somos, Jul 11 2011 */

a(n)=([1,2,2;2,1,2;2,2,3]^n)[3,3] \\ Vim Wenders, Mar 28 2007

G.f.: (1-3*x)/(1-6*x+x^2). - Barry E. Williams and Wolfdieter Lang, May 05 2000
E.g.f.: exp(3*x)*cosh(2*sqrt(2)*x). Binomial transform of A084128. - Paul Barry, May 16 2003
From N. J. A. Sloane, May 16 2003: (Start)
a(n) = sqrt(8*((A001109(n))^2) + 1).
a(n) = T(n, 3), with Chebyshev's T-polynomials A053120. (End)
a(n) = ((3+2*sqrt(2))^n + (3-2*sqrt(2))^n)/2.
a(n) = cosh(2*n*arcsinh(1)). - Herbert Kociemba, Apr 24 2008
a(n) ~ (1/2)*(sqrt(2) + 1)^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
For all elements x of the sequence, 2*x^2 - 2 is a square. Limit_{n -> infinity} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 10 2002 [corrected by Peter Pein, Mar 09 2009]
a(n) = 3*A001109(n) - A001109(n-1), n >= 1. - Barry E. Williams and Wolfdieter Lang, May 05 2000
For n >= 1, a(n) = A001652(n) - A001652(n-1). - Charlie Marion, Jul 01 2003
From Paul Barry, Sep 18 2003: (Start)
a(n) = ((-1+sqrt(2))^n + (1+sqrt(2))^n + (1-sqrt(2))^n + (-1-sqrt(2))^n)/4 (with interpolated zeros).
E.g.f.: cosh(x)*cosh(sqrt(2)x) (with interpolated zeros). (End)
For n > 0, a(n)^2 + 1 = 2*A001653(n-1)*A001653(n). - Charlie Marion, Dec 21 2003
a(n)^2 + a(n+1)^2 = 2*(A001653(2*n+1) - A001652(2*n)). - Charlie Marion, Mar 17 2003
a(n) = Sum_{k >= 0} binomial(2*n, 2*k)*2^k = Sum_{k >= 0} A086645(n, k)*2^k. - Philippe Deléham, Feb 29 2004
a(n)*A002315(n+k) = A001652(2*n+k) + A001652(k) + 1; for k > 0, a(n+k)*A002315(n) = A001652(2*n+k) - A001652(k-1). - Charlie Marion, Mar 17 2003
For n > k, a(n)*A001653(k) = A011900(n+k) + A053141(n-k-1). For n <= k, a(n)*A001653(k) = A011900(n+k) + A053141(k-n). - Charlie Marion, Oct 18 2004
A053141(n+1) + A055997(n+1) = a(n+1) + A001109(n+1). - Creighton Dement, Sep 16 2004
a(n+1) - A001542(n+1) = A090390(n+1) - A046729(n) = A001653(n); a(n+1) - 4*A079291(n+1) = (-1)^(n+1). Formula generated by the floretion - .5'i + .5'j - .5i' + .5j' - 'ii' + 'jj' - 2'kk' + 'ij' + .5'ik' + 'ji' + .5'jk' + .5'ki' + .5'kj' + e. - Creighton Dement, Nov 16 2004
a(n) = sqrt( A055997(2*n) ). - Alexander Adamchuk, Nov 24 2006
a(2n) = A056771(n). a(2*n+1) = 3*A077420(n). - Alexander Adamchuk, Feb 01 2007
a(n) = (A000129(n)^2)*4 + (-1)^n. - Vim Wenders, Mar 28 2007
2*a(k)*A001653(n)*A001653(n+k) = A001653(n)^2 + A001653(n+k)^2 + A001542(k)^2. - Charlie Marion, Oct 12 2007
a(n) = A001333(2*n). - Ctibor O. Zizka, Aug 13 2008
A028982(a(n)-1) + 2 = A028982(a(n)+1). - Juri-Stepan Gerasimov, Mar 28 2011
a(n) = 2*A001108(n) + 1. - Paul Weisenhorn, Dec 17 2011
a(n) = sqrt(2*x^2 + 1) with x being A001542(n). - Zak Seidov, Jan 30 2013
a(2n) = 2*a(n)^2 - 1 = a(n)^2 + 2*A001542(n)^2. a(2*n+1) = 1 + 2*A002315(n)^2. - Steven J. Haker, Dec 04 2013
a(n) = 3*a(n-1) + 4*A001542(n-1); e.g., a(4) = 99 = 3*17 + 4*12. - Zak Seidov, Dec 19 2013
a(n) = cos(n * arccos(3)) = cosh(n * log(3 + 2*sqrt(2))). - Daniel Suteu, Jul 28 2016
From Ilya Gutkovskiy, Jul 28 2016: (Start)
Inverse binomial transform of A084130.
Exponential convolution of A000079 and A084058.
Sum_{n>=0} (-1)^n*a(n)/n! = cosh(2*sqrt(2))/exp(3) = 0.4226407909842764637... (End)
a(2*n+1) = 2*a(n)*a(n+1) - 3. - Timothy L. Tiffin, Oct 12 2016
a(n) = a(-n) for all n in Z. - Michael Somos, Jan 20 2017
a(2^n) = A001601(n+1). - A.H.M. Smeets, May 28 2017
a(A298210(n)) = A002350(2*n^2). - A.H.M. Smeets, Jan 25 2018
a(n) = S(n, 6) - 3*S(n-1, 6), for n >= 0, with S(n, 6) = A001109(n+1), (Chebyshev S of A049310). See the first comment and the formula a(n) = T(n, 3). - Wolfdieter Lang, Nov 22 2020
From Peter Bala, Dec 31 2021: (Start)
a(n) = [x^n] (3*x + sqrt(1 + 8*x^2))^n.
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) hold for all prime p and positive integers n and k.
O.g.f. A(x) = 1 + x*d/dx(log(B(x))), where B(x) = 1/sqrt(1 - 6*x + x^2) is the o.g.f. of A001850. (End)
From Peter Bala, Aug 17 2022: (Start)
Sum_{n >= 1} 1/(a(n) - 2/a(n)) = 1/2.
Sum_{n >= 1} (-1)^(n+1)/(a(n) + 1/a(n)) = 1/4.
Sum_{n >= 1} 1/(a(n)^2 - 2) = 1/2 - 1/sqrt(8). (End)
From  Peter Bala, Jun 23 2025: (Start)
Product_{n >= 0} (1 + 1/a(2^n)) = sqrt(2).
Product_{n >= 0} (1 - 1/(2*a(2^n))) = (4/7)*sqrt(2). See A002812. (End)

A029549 a(n + 3) = 35a(n + 2) - 35a(n + 1) + a(n), with a(0) = 0, a(1) = 6, a(2) = 210.

Original entry on oeis.org

0, 6, 210, 7140, 242556, 8239770, 279909630, 9508687656, 323015470680, 10973017315470, 372759573255306, 12662852473364940, 430164224521152660, 14612920781245825506, 496409142337836914550, 16863297918705209269200, 572855720093639278238256

Offset: 0

Don N. Page

Triangular numbers that are twice other triangular numbers. - Don N. Page
Triangular numbers that are also pronic numbers. These will be shown to have a Pythagorean connection in a paper in preparation. - Stuart M. Ellerstein (ellerstein(AT)aol.com), Mar 09 2002
In other words, triangular numbers which are products of two consecutive numbers. E.g., a(2) = 210: 210 is a triangular number which is the product of two consecutive numbers: 14 * 15. - Shyam Sunder Gupta, Oct 26 2002
Coefficients of the series giving the best rational approximations to sqrt(8). The partial sums of the series 3 - 1/a(1) - 1/a(2) - 1/a(3) - ... give the best rational approximations to sqrt(8) = 2 sqrt(2), which constitute every second convergent of the continued fraction. The corresponding continued fractions are [2; 1, 4, 1], [2; 1, 4, 1, 4, 1], [2; 1, 4, 1, 4, 1, 4, 1], [2; 1, 4, 1, 4, 1, 4, 1, 4, 1] and so forth. - Gene Ward Smith, Sep 30 2006
This sequence satisfy the same recurrence as A165518. - Ant King, Dec 13 2010
Intersection of A000217 and A002378.
This is the sequence of areas, x(n)*y(n)/2, of the ordered Pythagorean triples (x(n), y(n) = x(n) + 1,z(n)) with x(0) = 0, y(0) = 1, z(0) = 1, a(0) = 0 and x(1) = 3, y(1) = 4, z(1) = 5, a(1) = 6. - George F. Johnson, Aug 20 2012

Reinhard Zumkeller, Table of n, a(n) for n = 0..100
H. J. Hindin, Stars, hexes, triangular numbers and Pythagorean triples, J. Rec. Math., 16 (1983/1984), 191-193. (Annotated scanned copy)
Shyam Sunder Gupta Fascinating Triangular Numbers
Roger B. Nelson, Multi-Polygonal Numbers, Mathematics Magazine, Vol. 89, No. 3 (June 2016), pp. 159-164.
Vladimir Pletser, Recurrent Relations for Multiple of Triangular Numbers being Triangular Numbers, arXiv:2101.00998 [math.NT], 2021.
Vladimir Pletser, Closed Form Equations for Triangular Numbers Multiple of Other Triangular Numbers, arXiv:2102.12392 [math.GM], 2021.
Vladimir Pletser, Triangular Numbers Multiple of Triangular Numbers and Solutions of Pell Equations, arXiv:2102.13494 [math.NT], 2021.
Vladimir Pletser, Congruence Properties of Indices of Triangular Numbers Multiple of Other Triangular Numbers, arXiv:2103.03019 [math.GM], 2021.
Vladimir Pletser, Searching for multiple of triangular numbers being triangular numbers, 2021.
Vladimir Pletser, Using Pell equation solutions to find all triangular numbers multiple of other triangular numbers, 2021.
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).

Cf. A123478, A123479, A123480, A123482, A075528, A082405 (first differences).
Cf. A000129, A001108, A002203, A009111, A245031.
Cf. A000217, A001109, A001652, A002315, A002378, A011900, A029546, A046090, A046729, A053141, A084159, A096979, A157259, A165518.

List([0..20], n-> (Lucas(2,-1, 4*n+2)[2] -6)/32 ); # G. C. Greubel, Jan 13 2020

a029549 n = a029549_list !! n
a029549_list = [0,6,210] ++
   zipWith (+) a029549_list
               (map (* 35) $ tail delta)
   where delta = zipWith (-) (tail a029549_list) a029549_list
-- Reinhard Zumkeller, Sep 19 2011

(makelist(binom(n,2),n,1,999999),intersection(%%,2*%%)) /* Bill Gosper, Feb 07 2010 */

R:=PowerSeriesRing(Integers(), 25); [0] cat Coefficients(R!(6/(1-35*x+35*x^2-x^3))); // G. C. Greubel, Jul 15 2018

A029549 := proc(n)
    option remember;
    if n <= 1 then
        op(n+1,[0,6]) ;
    else
        34*procname(n-1)-procname(n-2)+6 ;
    end if;
end proc: # R. J. Mathar, Feb 05 2016

Table[Floor[(Sqrt[2] + 1)^(4n + 2)/32], {n, 0, 20} ] (* Original program from author, corrected by Ray Chandler, Jul 09 2015 *)
CoefficientList[Series[6/(1 - 35x + 35x^2 - x^3), {x, 0, 14}], x]
Intersection[#, 2#] &@ Table[Binomial[n, 2], {n, 999999}] (* Bill Gosper, Feb 07 2010 *)
LinearRecurrence[{35, -35, 1}, {0, 6, 210}, 20] (* Harvey P. Dale, Jun 06 2011 *)
(LucasL[4Range[20] - 2, 2] -6)/32 (* G. C. Greubel, Jan 13 2020 *)

concat(0,Vec(6/(1-35*x+35*x^2-x^3)+O(x^25))) \\ Charles R Greathouse IV, Jun 13 2013

[(lucas_number2(4*n+2, 2, -1) -6)/32 for n in (0..20)] # G. C. Greubel, Jan 13 2020

val triNums = (0 to 39999).map(n => (n * n + n)/2)
triNums.filter( % 2 == 0).filter(n => (triNums.contains(n/2))) // _Alonso del Arte, Jan 12 2020

G.f.: 6*x/(1 - 35*x + 35*x^2 - x^3) = 6*x /( (1-x)*(1 - 34*x + x^2) ).
a(n) = 6*A029546(n-1) = 2*A075528(n).
a(n) = -3/16 + ((3+2*sqrt(2))/32) *(17 + 12*sqrt(2))^n + ((3-2*sqrt(2))/32) *(17 - 12*sqrt(2))^n. - Gene Ward Smith, Sep 30 2006
From Bill Gosper, Feb 07 2010: (Start)
a(n) = (cosh((4*n + 2)*log(1 + sqrt(2))) - 3)/16.
a(n) = binomial(A001652(n) + 1, 2) = 2*binomial(A053141(n) + 1, 2). (End)
a(n) = binomial(A046090(n), 2) = A000217(A001652(n)). - Mitch Harris, Apr 19 2007, R. J. Mathar, Jun 26 2009
a(n) = ceiling((3 + 2*sqrt(2))^(2n + 1) - 6)/32 = floor((1/32) (1+sqrt(2))^(4n+2)). - Ant King, Dec 13 2010
Sum_{n >= 1} 1/a(n) = 3 - 2*sqrt(2) = A157259 - 4. - Ant King, Dec 13 2010
a(n) = a(n - 1) + A001109(2n). - Charlie Marion, Feb 10 2011
a(n+2) = 34*a(n + 1) - a(n) + 6. - Charlie Marion, Feb 11 2011
From George F. Johnson, Aug 20 2012: (Start)
a(n) = ((3 + 2*sqrt(2))^(2*n + 1) + (3 - 2*sqrt(2))^(2*n + 1) - 6)/32.
8*a(n) + 1 = (A002315(n))^2, 4*a(n) + 1 = (A000129(2*n + 1))^2, 32*a(n)^2 + 12*a(n) + 1 are perfect squares.
a(n + 1) = 17*a(n) + 3 + 3*sqrt((8*a(n) + 1)*(4*a(n) + 1)).
a(n - 1) = 17*a(n) + 3 - 3*sqrt((8*a(n) + 1)*(4*a(n) + 1)).
a(n - 1)*a(n + 1) = a(n)*(a(n) - 6), a(n) = A096979(2*n).
a(n) = (1/2)*A084159(n)*A046729(n) = (1/2)*A001652(n)*A046090(n).
Limit_{n->infinity} a(n)/a(n - 1) = 17 + 12*sqrt(2).
Limit_{n->infinity} a(n)/a(n - 2) = (17 + 12*sqrt(2))^2 = 577 + 408*sqrt(2).
Limit_{n->infinity} a(n)/a(n - r) = (17 + 12*sqrt(2))^r.
Limit_{n->infinity} a(n - r)/a(n) = (17 + 12*sqrt(2))^(-r) = (17 - 12*sqrt(2))^r. (End)
a(n) = 3 * T( b(n) ) + (2*b(n) + 1)*sqrt( T( b(n) ) ) where b(n) = A001108(n) (indices of the square triangular numbers), T(n) = A000217(n) (the n-th triangular number). - Dimitri Papadopoulos, Jul 07 2017
a(n) = (Pell(2*n + 1)^2 - 1)/4 = (Q(4*n + 2) - 6)/32, where Q(n) are the Pell-Lucas numbers (A002203). - G. C. Greubel, Jan 13 2020
a(n) = A002378(A011900(n)-1) = A002378(A053141(n)). - Pontus von Brömssen, Sep 11 2024

Additional comments from Christian G. Bower, Sep 19 2002; T. D. Noe, Nov 07 2006; and others

Edited by N. J. A. Sloane, Apr 18 2007, following suggestions from Andrew S. Plewe and Tanya Khovanova

A046090 Consider all Pythagorean triples (X,X+1,Z) ordered by increasing Z; sequence gives X+1 values.

Original entry on oeis.org

1, 4, 21, 120, 697, 4060, 23661, 137904, 803761, 4684660, 27304197, 159140520, 927538921, 5406093004, 31509019101, 183648021600, 1070379110497, 6238626641380, 36361380737781, 211929657785304, 1235216565974041, 7199369738058940, 41961001862379597, 244566641436218640

Offset: 0

Eric W. Weisstein

Solution to a*(a-1) = 2b*(b-1) in natural numbers: a = a(n), b = b(n) = A011900(n).
n such that n^2 = (1/2)*(n+floor(sqrt(2)*n*floor(sqrt(2)*n))). - Benoit Cloitre, Apr 15 2003
Place a(n) balls in an urn, of which b(n) = A011900(n) are red; draw 2 balls without replacement; 2*Probability(2 red balls) = Probability(2 balls); this is equivalent to the Pell equation A(n)^2-2*B(n)^2 = -1 with a(n) = (A(n)+1)/2; b(n) = (B(n)+1)/2; and the fundamental solution (7;5) and the solution (3;2) for the unit form. - Paul Weisenhorn, Aug 03 2010
Find base x in which repdigit yy has a square that is repdigit zzzz, corresponding to Diophantine equation zzzz_x = (yy_x)^2; then, solution z = a(n) with x = A002315(n) and y = A001653(n+1) for n >= 1 (see Maurice Protat reference). - Bernard Schott, Dec 21 2022

			For n=4: a(4)=697; b(4)=493; 2*binomial(493,2)=485112=binomial(697,2). - _Paul Weisenhorn_, Aug 03 2010

A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, pp. 122-125, 1964.
Maurice Protat, Des Olympiades à l'Agrégation, De zzzz_x = (yy_x)^2 à Pell-Fermat, Problème 23, pp. 52-54, Ellipses, Paris, 1997.

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
T. W. Forget and T. A. Larkin, Pythagorean triads of the form X, X+1, Z described by recurrence sequences, Fib. Quart., 6 (No. 3, 1968), 94-104.
L. J. Gerstein, Pythagorean triples and inner products, Math. Mag., 78 (2005), 205-213.
H. J. Hindin, Stars, hexes, triangular numbers and Pythagorean triples, J. Rec. Math., 16 (1983/1984), 191-193. (Annotated scanned copy)
Ron Knott, Pythagorean Triples and Online Calculators
S. Northshield, An Analogue of Stern's Sequence for Z[sqrt(2)], Journal of Integer Sequences, 18 (2015), #15.11.6.
Eric Weisstein's World of Mathematics, Pythagorean Triple
Index entries for two-way infinite sequences
Index entries for linear recurrences with constant coefficients, signature (7,-7,1).

Other 2 sides are A001652 and A001653.
Cf. A001009, A001541, A001542, A011900, A046729, A053141, A084159, A084703, A156035.
See comments in A301383.
Cf. A001109, A001653, A002024, A002315, A029549.

a046090 n = a046090_list !! n
a046090_list = 1 : 4 : map (subtract 2)
   (zipWith (-) (map (* 6) (tail a046090_list)) a046090_list)
-- Reinhard Zumkeller, Jan 10 2012

m:=30; R:=PowerSeriesRing(Integers(), m); Coefficients(R!((1-3*x)/((1-6*x+x^2)*(1-x)))); // G. C. Greubel, Jul 15 2018

Digits:=100: seq(round((1+(7+5*sqrt(2))*(3+2*sqrt(2))^(n-1))/2)/2, n=0..20); # Paul Weisenhorn, Aug 03 2010

Join[{1},#+1&/@With[{c=3+2Sqrt[2]},NestList[Floor[c #]+3&,3,20]]] (* Harvey P. Dale, Aug 19 2011 *)
LinearRecurrence[{7,-7,1},{1,4,21},25] (* Harvey P. Dale, Apr 13 2012 *)
a[n_] := (2-ChebyshevT[n, 3]+ChebyshevT[n+1, 3])/4; Array[a, 21, 0] (* Jean-François Alcover, Jul 10 2016, adapted from PARI *)

a(n)=(2-subst(poltchebi(abs(n))-poltchebi(abs(n+1)),x,3))/4

x='x+O('x^30); Vec((1-3*x)/((1-6*x+x^2)*(1-x))) \\ G. C. Greubel, Jul 15 2018

a(n) = (-1+sqrt(1+8*b(n)*(b(n)+1)))/2 with b(n) = A011900(n). [corrected by Michel Marcus, Dec 23 2022]
a(n) = 6*a(n-1) - a(n-2) - 2, n >= 2, a(0) = 1, a(1) = 4.
a(n) = (A(n+1) - 3*A(n) + 2)/4 with A(n) = A001653(n).
A001652(n) = -a(-1-n).
From Barry E. Williams, May 03 2000: (Start)
G.f.: (1-3*x)/((1-6*x+x^2)*(1-x)).
a(n) = partial sums of A001541(n). (End)
From Charlie Marion, Jul 01 2003: (Start)
A001652(n)*A001652(n+1) + a(n)*a(n+1) = A001542(n+1)^2 = A084703(n+1).
Let a(n) = A001652(n), b(n) = this sequence and c(n) = A001653(n). Then for k > j, c(i)*(c(k) - c(j)) = a(k+i) + ... + a(i+j+1) + a(k-i-1) + ... + a(j-i) + k - j. For n < 0, a(n) = -b(-n-1). Also a(n)*a(n+2k+1) + b(n)*b(n+2k+1) + c(n)*c(n+2k+1) = (a(n+k+1) - a(n+k))^2; a(n)*a(n+2k) + b(n)*b(n+2k) + c(n)*c(n+2k) = 2*c(n+k)^2. (End)
a(n) = 1/2 + ((1-2^(1/2))/4)*(3 - 2^(3/2))^n + ((1+2^(1/2))/4)*(3 + 2^(3/2))^n. - Antonio Alberto Olivares, Oct 13 2003
2*a(n) = 2*A084159(n) + 1 + (-1)^(n+1) = 2*A046729(n) + 1 - (-1)^(n+1). - Lekraj Beedassy, Jul 16 2004
a(n) = A001109(n+1) - A053141(n). - Manuel Valdivia, Apr 03 2010
From Paul Weisenhorn, Aug 03 2010: (Start)
a(n+1) = round((1+(7+5*sqrt(2))*(3+2*sqrt(2))^n)/2);
b(n+1) = round((2+(10+7*sqrt(2))*(3+2*sqrt(2))^n)/4) = A011900(n+1).
(End)
a(n)*(a(n)-1)/2 = b(n)*b(n+1) and 2*a(n) - 1 = b(n) + b(n+1), where b(n) = A001109. - Kenneth J Ramsey, Apr 24 2011
T(a(n)) = A011900(n)^2 + A001109(n), where T(n) is the n-th triangular number. See also A001653. - Charlie Marion, Apr 25 2011
a(0)=1, a(1)=4, a(2)=21, a(n) = 7*a(n-1) - 7*a(n-2) + a(n-3). - Harvey P. Dale, Apr 13 2012
Limit_{n->oo} a(n+1)/a(n) = 3 + 2*sqrt(2) = A156035. - Ilya Gutkovskiy, Jul 10 2016
a(n) = A001652(n)+1. - Dimitri Papadopoulos, Jul 06 2017
a(n) = (A002315(n) + 1)/2. - Bernard Schott, Dec 21 2022
E.g.f.: (exp(x) + exp(3*x)*(cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x)))/2. - Stefano Spezia, Mar 16 2024
a(n) = A002024(A029549(n))+1. - Pontus von Brömssen, Sep 11 2024

Additional comments from Wolfdieter Lang

Comment moved to A001653 by Claude Morin, Sep 22 2023

A075528 Triangular numbers that are half other triangular numbers.

Original entry on oeis.org

0, 3, 105, 3570, 121278, 4119885, 139954815, 4754343828, 161507735340, 5486508657735, 186379786627653, 6331426236682470, 215082112260576330, 7306460390622912753, 248204571168918457275, 8431648959352604634600, 286427860046819639119128

Offset: 0

Christian G. Bower, Sep 19 2002

This is the sequence of 1/2 the areas, x(n)*y(n)/2, of the ordered Pythagorean triples (x(n), y(n)=x(n)+1, z(n)) with x(0)=0, y(0)=1, z(0)=1, a(0)=0 and x(1)=3, y(1)=4, z(1)=5, a(1)=3. - George F. Johnson, Aug 24 2012

Colin Barker, Table of n, a(n) for n = 0..653
Martin V. Bonsangue, Gerald E. Gannon and Laura J. Pheifer, Misinterpretations can sometimes be a good thing, Math. Teacher, vol. 95, No. 6 (2002) pp. 446-449.
Harvey J. Hindin, Stars, hexes, triangular numbers and Pythagorean triples, J. Rec. Math., 16 (1983/1984), 191-193. (Annotated scanned copy)
Vladimir Pletser, Recurrent Relations for Multiple of Triangular Numbers being Triangular Numbers, arXiv:2101.00998 [math.NT], 2021.
Vladimir Pletser, Closed Form Equations for Triangular Numbers Multiple of Other Triangular Numbers, arXiv:2102.12392 [math.GM], 2021.
Vladimir Pletser, Triangular Numbers Multiple of Triangular Numbers and Solutions of Pell Equations, arXiv:2102.13494 [math.NT], 2021.
Vladimir Pletser, Congruence Properties of Indices of Triangular Numbers Multiple of Other Triangular Numbers, arXiv:2103.03019 [math.GM], 2021.
Vladimir Pletser, Searching for multiple of triangular numbers being triangular numbers, 2021.
Vladimir Pletser, Using Pell equation solutions to find all triangular numbers multiple of other triangular numbers, 2022.
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).

Cf. A000129, A000217, A001652, A002315, A029546, A029549, A046090, A046729, A053141, A084159, A096979.

CoefficientList[ Series[ 3x/(1 - 35 x + 35 x^2 - x^3), {x, 0, 15}], x] (* Robert G. Wilson v, Jun 24 2011 *)

concat(0, Vec(3*x/((1-x)*(1-34*x+x^2)) + O(x^20))) \\ Colin Barker, Jun 18 2015

a(n) = 3*A029546(n-1) = A029549(n)/2.
G.f.: 3*x/((1-x)*(1-34*x+x^2)).
From George F. Johnson, Aug 24 2012: (Start)
a(n) = ((3+2*sqrt(2))^(2*n+1) + (3-2*sqrt(2))^(2*n+1) - 6)/64.
8*a(n)+1 = A000129(2*n+1)^2.
16*a(n)+1 = A002315(n)^2.
128*a(n)^2 + 24*a(n) + 1 is a perfect square.
a(n+1) = 17*a(n) + 3/2 + 3*sqrt((8*a(n)+1)*(16*a(n)+1))/2.
a(n-1) = 17*a(n) + 3/2 - 3*sqrt((8*a(n)+1)*(16*a(n)+1))/2.
a(n-1)*a(n+1) = a(n)*(a(n)-3); a(n+1) = 34*a(n) - a(n-1) + 3.
a(n+1) = 35*a(n) - 35*a(n-1) + a(n-2); a(n) = A096979(2*n)/2.
a(n) = A084159(n)*A046729(n)/4 = A001652(n)*A046090(n)/4.
Lim_{n->infinity} a(n)/a(n-1) =  17 + 12*sqrt(2).
Lim_{n->infinity} a(n)/a(n-2) = (17 + 12*sqrt(2))^2 = 577 + 408*sqrt(2).
Lim_{n->infinity} a(n)/a(n-r) = (17 + 12*sqrt(2))^r.
Lim_{n->infinity} a(n-r)/a(n) = (17 + 12*sqrt(2))^(-r) = (17 - 12*sqrt(2))^r.
(End)
a(n) = 34*a(n-1) - a(n-2) + 3, n >= 2. - R. J. Mathar, Nov 07 2015
a(n) = A000217(A053141(n)). - R. J. Mathar, Aug 16 2019
a(n) = (a(n-1)*(a(n-1)-3))/a(n-2) for n > 2. - Vladimir Pletser, Apr 08 2020
Sum_{n>=1} 1/a(n) = 2*(3 - 2*sqrt(2)). - Amiram Eldar, Dec 04 2024

A005319 a(n) = 6*a(n-1) - a(n-2).

Original entry on oeis.org

0, 4, 24, 140, 816, 4756, 27720, 161564, 941664, 5488420, 31988856, 186444716, 1086679440, 6333631924, 36915112104, 215157040700, 1254027132096, 7309005751876, 42600007379160, 248291038523084, 1447146223759344, 8434586304032980

Offset: 0

N. J. A. Sloane

Solutions y of the equation 2x^2-y^2=2; the corresponding x values are given by A001541. - N-E. Fahssi, Feb 25 2008
The lower intermediate convergents to 2^(1/2) beginning with 4/3, 24/17, 140/99, 816/577, form a strictly increasing sequence; essentially, numerators=A005319 and denominators=A001541. - Clark Kimberling, Aug 26 2008
Numbers n such that (ceiling(sqrt(n*n/2)))^2 = 1 + n*n/2. - Ctibor O. Zizka, Nov 09 2009
All nonnegative solutions of the indefinite binary quadratic form X^2 + 4*X*Y -4*Y^2 of discriminant 32, representing -4 are (X(n), Y(n)) = (a(n), A001653(n+1)), for n >= 0. - Wolfdieter Lang, Jun 13 2018
Also the number of edge covers in the n-triangular snake graph. - Eric W. Weisstein, Jun 08 2019
All of the positive integer solutions of a*b+1=x^2, a*c+1=y^2, b*c+1=z^2, x+z=2*y, 0A001542(n), b=A005319(n), c=A001542(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0Michael Somos, Jun 26 2022
a(n) is the sum of 4*n consecutive powers of the silver ratio 1+sqrt(2), starting at (1+sqrt(2))^(-2*n) and ending at (1+sqrt(2))^(2*n-1). - Greg Dresden and Ruxin Sheng, Jul 25 2024

			G.f. = 4*x + 24*x^2 + 140*x^3 + 816*x^4 + 4756*x^5 + ... - _Michael Somos_, Jun 26 2022

P. de la Harpe, Topics in Geometric Group Theory, Univ. Chicago Press, 2000, p. 160, middle display.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Marius A. Burtea, Table of n, a(n) for n = 0..100
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
Hacène Belbachir, Soumeya Merwa Tebtoub, and László Németh, Ellipse Chains and Associated Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.8.5.
John M. Campbell, An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences, arXiv preprint arXiv:1105.3399 [math.GM], 2011.
Tanya Khovanova, Recursive Sequences
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Soumeya M. Tebtoub, Hacène Belbachir, and László Németh, Integer sequences and ellipse chains inside a hyperbola, Proceedings of the 1st International Conference on Algebras, Graphs and Ordered Sets (ALGOS 2020), hal-02918958 [math.cs], 17-18.
Eric Weisstein's World of Mathematics, Edge Cover
Eric Weisstein's World of Mathematics, Triangular Snake Graph
Index entries for linear recurrences with constant coefficients, signature (6,-1).

Cf. A000129, A001109, A001541, A001542, A001652, A001653, A046090, A046729, A090390.

a:=[0,4]; [n le 2 select a[n] else 6*Self(n-1) - Self(n-2):n in [1..22]]; // Marius A. Burtea, Sep 19 2019

LinearRecurrence[{6, -1}, {0, 4}, 22] (* Jean-François Alcover, Sep 26 2017 *)
Table[((3 + 2 Sqrt[2])^n - (3 - 2 Sqrt[2])^n)/Sqrt[2], {n, 20}] // Expand (* Eric W. Weisstein, Jun 08 2019 *)
CoefficientList[Series[(4 x)/(1 - 6 x + x^2), {x, 0, 20}], x] (* Eric W. Weisstein, Jun 08 2019 *)
a[ n_] := 4*ChebyshevU[n-1, 3]; (* Michael Somos, Jun 26 2022 *)

{a(n) = 4*polchebyshev(n-1, 2, 3)}; /* Michael Somos, Jun 26 2022 */

G.f.: 4*x / ( 1-6*x+x^2 ). - Simon Plouffe in his 1992 dissertation.
G.f. for signed version beginning with 1: (1+2*x+x^2)/(1+6*x+x^2).
For any term n of the sequence, 2*n^2 + 4 is a perfect square. Limit_{n->infinity} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 06 2002
a(n) = ((3+2*sqrt(2))^n - (3-2*sqrt(2))^n) / sqrt(2). - Gregory V. Richardson, Oct 06 2002
(-1)^(n+1) = A090390(n+1) + A001542(n+1) + A046729(n) - a(n) (conjectured). - Creighton Dement, Nov 17 2004
For n > 0, a(n) = A000129(n+1)^2 - A000129(n-1)^2; a(n) = A046090(n-1) + A001652(n); e.g., 816 = 120 + 696; a(n) = A001653(n) - A001653(n-1); e.g., 816 = 985 - 169. - Charlie Marion Jul 22 2005
a(n) = 4*A001109(n). - M. F. Hasler, Mar 2009
For n > 1, a(n) is the denominator of continued fraction [1,4,1,4,...,1,4] with (n-1) repetitions of 1,4. For the numerators, see A001653. - Greg Dresden, Sep 10 2019
1/a(n) - 1/a(n+1) = 1/(Pell(2*n+1) - 1/Pell(2*n+1)) for n >= 1, where Pell(n) = A000129(n). - Peter Bala, Aug 21 2022
E.g.f.: sqrt(2)*exp(3*x)*sinh(2*sqrt(2)*x). - Stefano Spezia, Nov 25 2022
a(n) = 2*A000129(2*n). - Tanya Khovanova and MIT PRIMES STEP senior group, Apr 17 2024

cp's OEIS Frontend

A089498 Duplicate of A046729.

Views

Author

Keywords

A002315 NSW numbers: a(n) = 6*a(n-1) - a(n-2); also a(n)^2 - 2*b(n)^2 = -1 with b(n) = A001653(n+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Formula

A001541 a(0) = 1, a(1) = 3; for n > 1, a(n) = 6*a(n-1) - a(n-2).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Formula

A029549 a(n + 3) = 35*a(n + 2) - 35*a(n + 1) + a(n), with a(0) = 0, a(1) = 6, a(2) = 210.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Haskell

Macsyma

Magma

Maple

Mathematica

PARI

Sage

Scala

Formula

Extensions

A046090 Consider all Pythagorean triples (X,X+1,Z) ordered by increasing Z; sequence gives X+1 values.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

A002315 NSW numbers: a(n) = 6a(n-1) - a(n-2); also a(n)^2 - 2b(n)^2 = -1 with b(n) = A001653(n+1).

A029549 a(n + 3) = 35a(n + 2) - 35a(n + 1) + a(n), with a(0) = 0, a(1) = 6, a(2) = 210.