A055792 -id:A055792 - cp's OEIS frontend

A204575 Squares such that [a(n)/2] is again a square (A055792), written in binary.

0, 1, 1001, 100100001, 10011001001001, 1010001010010000001, 101011001001001011001001, 10110111001100110101000100001, 1100001001111011011000010110001001, 110011100111010101001010101001000000001

Offset: 1

M. F. Hasler, Jan 16 2012

See also A031149=sqrt(A023110) (base 10), A204502=sqrt(A204503) (base 9), A204514=sqrt(A055872) (base 8), A204516=sqrt(A055859) (base 7), A204518=sqrt(A055851) (base 6), A204520=sqrt(A055812) (base 5), A004275=sqrt(A055808) (base 4), A001075=sqrt(A055793) (base 3), A001541=sqrt(A055792) (base 2).

A204576 Floor[A055792(n-1)/2]=A084703(n-2) (truncated squares), written in binary.

Original entry on oeis.org

0, 0, 100, 10010000, 1001100100100, 101000101001000000, 10101100100100101100100, 1011011100110011010100010000, 110000100111101101100001011000100, 11001110011101010100101010100100000000, 1101101100101100000000000111010111101000100

Offset: 1

M. F. Hasler, Jan 16 2012

A204575 with the last (binary) digit (necessarily = 1, except for a(1)=0) deleted.

Also: Squares, written in binary, such that appending a (binary) digit (necessarily = 1) yields another square (except for a(1)=0 which corresponds to A204575(1)=00, the only square which remains square when multiplied by 2).

David W. Wilson, Table of n, a(n) for n = 1..100

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

N. J. A. Sloane

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.
  See also A002315, A054488, A038761, A054489, A054490.
  - 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

			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

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

Indranil Ghosh, Table of n, a(n) for n = 0..1304 (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.
Irving Adler, Three Diophantine equations - Part II, Fib. Quart., 7 (1969), pp. 181-193.
Seyed Hassan Alavi, Ashraf Daneshkhah, and Cheryl E. Praeger, Symmetries of biplanes, arXiv:2004.04535 [math.GR], 2020. See v_n in Lemma 7.9 p. 21.
Jean-Paul Allouche, Zeta-regularization of arithmetic sequences, EPJ Web of Conferences (2020) Vol. 244, 01008.
Dario Alpern for Diophantine equation a^4+b^3=c^2.
Kasper Andersen, Lisa 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.
Francesca Arici and Jens Kaad, Gysin sequences and SU(2)-symmetries of C*-algebras, arXiv:2012.11186 [math.OA], 2020.
Muniru A. Asiru, All square chiliagonal numbers, International Journal of Mathematical Education in Science and Technology, Volume 47, 2016 - Issue 7.
Jeremiah Bartz, Bruce Dearden, and Joel Iiams, Classes of Gap Balancing Numbers, arXiv:1810.07895 [math.NT], 2018.
Jeremiah Bartz, Bruce Dearden, Joel Iiams, and Julia Peterson, Powers of Two as Sums of Two Balancing Numbers, Combinatorics, Graph Theory and Computing (SEICCGTC 2021) Springer Proc. Math. Stat., Vol 448, pp. 383-392. See p. 384.
Raymond A. Beauregard and Vladimir A. Dobrushkin, Powers of a Class of Generating Functions, Mathematics Magazine, Volume 89, Number 5, December 2016, pp. 359-363.
A. Behera and G. K. Panda, On the Square Roots of Triangular Numbers, Fib. Quart., 37 (1999), pp. 98-105.
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.
Elwyn Berlekamp and Joe P. Buhler, Puzzle Column, Emissary, MSRI Newsletter, Fall 2005. Problem 1, (6 MB).
Kisan Bhoi and Prasanta Kumar Ray, On the Diophantine equation Bn1+Bn2=2^a1+2^a2+2^a3, arXiv:2212.06372 [math.NT], 2022.
Daniel Birmajer, Juan B. Gil, and Michael D. Weiner, On the Enumeration of Restricted Words over a Finite Alphabet, J. Int. Seq. 19 (2016) # 16.1.3, example 12.
Alexander Bogomolny, There exist triangular numbers that are also squares
John C. Butcher, On Ramanujan, continued Fractions and an interesting number
Paula Catarino, Helena Campos, and Paulo Vasco, On some identities for balancing and cobalancing numbers, Annales Mathematicae et Informaticae, 45 (2015) pp. 11-24.
E. K. Çetinalp, N. Yilmaz, and Ö. Deveci, The balancing-like sequences in groups, Acta Univ. Apulensis Math. (2023) No. 73, 139-153. See p. 144.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (pp. 301, 302, P_{13}).
Mahadi Ddamulira, Repdigits as sums of three balancing numbers, Mathematica Slovaca, (2019) hal-02405969.
Tomislav Doslic, Planar polycyclic graphs and their Tutte polynomials, Journal of Mathematical Chemistry, Volume 51, Issue 6, 2013, pp. 1599-1607.
D. B. Eperson, Triangular numbers, Math. Gaz., 47 (1963), 236-237.
Leonhard Euler, De solutione problematum diophanteorum per numeros integros, Par. 19.
Sergio Falcon, Relationships between Some k-Fibonacci Sequences, Applied Mathematics, 2014, 5, 2226-2234.
Bernadette Faye, Florian Luca, and Pieter Moree, On the discriminator of Lucas sequences, arXiv:1708.03563 [math.NT], 2017.
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.
Rigoberto Flórez, Robinson A. Higuita, and Antara Mukherjee, Alternating Sums in the Hosoya Polynomial Triangle, Article 14.9.5 Journal of Integer Sequences, Vol. 17 (2014).
Aviezri S. Fraenkel, On the recurrence f(m+1)= b(m)*f(m)-f(m-1) and applications, Discrete Mathematics 224 (2000), pp. 273-279.
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 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.
Robert Frontczak and Kalika Prasad, Balancing polynomials, Fibonacci numbers and some new series for $\pi$, Mediterranean Journal of Mathematics (2023) Vol. 20, Article number: 207.
Bill Gosper, The Triangular Squares, 2014.
H. Harborth, Fermat-like binomial equations, Applications of Fibonacci numbers, Proc. 2nd Int. Conf., San Jose/Ca., August 1986, 1-5 (1988).
Brian Hayes, Calculemus!, American Scientist, 96 (Sep-Oct 2008), 362-366.
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Michael A. Jones, Proof Without Words: The Square of a Balancing Number Is a Triangular Number, The College Mathematics Journal, Vol. 43, No. 3 (May 2012), p. 212.
Vedat Kabasakal and Fatih Yılmaz, On balancing numbers and application to coding theory, Ch. 16, Recent Developments in Mathematics, Pegem Akademi (Turkiye, 2024), 200-209. See p. 201.
Refik Keskin and Olcay Karaatli, Some New Properties of Balancing Numbers and Square Triangular Numbers, Journal of Integer Sequences, Vol. 15 (2012), Article #12.1.4.
Omar Khadir, Kalman Liptai, and Laszlo Szalay, On the Shifted Product of Binary Recurrences, J. Int. Seq. 13 (2010), 10.6.1.
Tanya Khovanova, Recursive Sequences
Phil Lafer, Discovering the square-triangular numbers, Fib. Quart., 9 (1971), 93-105.
Ioana-Claudia Lazăr, Lucas sequences in t-uniform simplicial complexes, arXiv:1904.06555 [math.GR], 2019.
Kalman Liptai, Fibonacci Balancing Numbers, Fib. Quart. 42 (4) (2004) 330-340.
Madras College, St Andrews, Square Triangular Numbers
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.
Roger B. Nelson, Multi-Polygonal Numbers, Mathematics Magazine, Vol. 89, No. 3 (June 2016), pp. 159-164.
G. K. Panda, Sequence balancing and cobalancing numbers, Fib. Q., Vol. 45, No. 3 (2007), 265-271. See p. 266.
G. K. Panda and S. S. Rout, Periodicity of Balancing Numbers, Acta Mathematica Hungarica 143 (2014), 274-286.
G. K. Panda and Ravi Kumar Davala, Perfect Balancing Numbers, Fibonacci Quart. 53 (2015), no. 3, 261-264.
Ashish Kumar Pandey and B. K. Sharma, On Inequalities Related to a Generalized Euler Totient Function and Lucas Sequences, J. Int. Seq. (2023) Vol. 26, Art. 23.8.6.
Poo-Sung Park, Ramanujan's Continued Fraction for a Puzzle, College Mathematics Journal, 2005, 363-365.
Michael Penn, Balancing Numbers, Youtube video, 2020.
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
B. Polster and M. Ross, Marching in squares, arXiv preprint arXiv:1503.04658 [math.HO], 2015.
Kalika Prasad, Munesh Kumari, Ritanjali Mohanty, and Hrishikesh Mahato, Spinor Algebra of k-Balancing and k-Lucas-Balancing Numbers, Journal of Algebra and Its Applications (2024).
Kalika Prasad, Munesh Kumari, and Jagmohan Tanti, Octonions and hyperbolic octonions with the k-balancing and k-Lucas balancing numbers, The Journal of Analysis, (2024), Vol. 32, No. 3, 1281-1296.
Kalika Prasad, Munesh Kumari, and Jagmohan Tanti, Generalized k-balancing and k-Lucas balancing numbers and associated polynomials, Kyungpook Mathematical Journal (2023), Vol. 63, No. 4, 539-550.
Kalika Prasad, Hrishikesh Mahato, and Munesh Kumari, Some properties of r-circulant matrices with k-balancing and k-Lucas balancing numbers, Boletín de la Sociedad Matemática Mexicana (2023), Vol. 29, No. 2, Article no. 44.
Helmut Prodinger, How to sum powers of balancing numbers efficiently, arXiv:2008.03916 [math.NT], 2020.
Rajesh Ram, Triangle Numbers that are Perfect Squares
K. J. Ramsey, Relation of Mersenne Primes To Square Triangular Numbers [edited by K. J. Ramsey, May 14 2011]
Kenneth Ramsay and Andras Erszegi, Relation of Square Triangular Numbers To Mersenne Primes, digest of 4 messages in Triangular_and_Fibonacci_Numbers Yahoo Group, May 15 - Jun 28, 2006.
Kenneth Ramsey, Generalized Proof re Square Triangular Numbers
Kenneth Ramsey, Generalized Proof re Square Triangular Numbers, digest of 2 messages in Triangular_and_Fibonacci_Numbers Yahoo group, May 27, 2005 - Oct 10, 2011.
Salah E. Rihane, Bernadette Faye, Florian Luca, and Alain Togbe, An exponential Diophantine equation related to the difference between powers of two consecutive Balancing numbers, arXiv:1811.03015 [math.NT], 2018.
A. Sandhya, Puzzle 4: A problem Srinivasa Ramanujan, the famous 20th century Indian Mathematician Solved
Sci.math Newsgroup, Square numbers which are triangular
Sci.math Newsgroup, Square numbers which are triangular [Cached copy]
R. A. Sulanke, Moments, Narayana numbers and the cut and paste for lattice paths.
R. A. Sulanke, Bijective recurrences concerning Schroeder paths, Electron. J. Combin. 5 (1998), Research Paper 47, 11 pp.
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.
Ahmet Tekcan, Merve Tayat, and Meltem E. Ozbek, The diophantine equation 8x^2-y^2+8x(1+t)+(2t+1)^2=0 and t-balancing numbers, ISRN Combinatorics, Volume 2014, Article ID 897834, 5 pages.
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.
Eric Weisstein's World of Mathematics, Binomial coefficient.
Eric Weisstein's World of Mathematics, Square Triangular Number.
Eric Weisstein's World of Mathematics, Triangular Number.
Wikipedia, Triangular square number
Rick Young, Relevant quotation from biography of Ramanujan
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).

Cf. A000217, A000290, A001108, A001542, A001653, A001850, A002315, A002965, A278310.
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.

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

a001109 n = a001109_list !! n :: Integer
a001109_list = 0 : 1 : zipWith (-)
   (map (* 6) $ tail a001109_list) a001109_list
-- Reinhard Zumkeller, Dec 17 2011

[n le 2 select n-1 else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Jul 25 2015

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

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 *)

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

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

{a(n) = polchebyshev( n-1, 2, 3)}; /* Michael Somos, Sep 02 2012 */

is(n)=ispolygonal(n^2,3) \\ Charles R Greathouse IV, Nov 03 2016

[lucas_number1(n,6,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008

[chebyshev_U(n-1,3) for n in (0..20)] # G. C. Greubel, Dec 23 2019

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

Additional comments from Wolfdieter Lang, Feb 10 2000

Duplication of a formula removed by Wolfdieter Lang, Feb 10 2015

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)

A023110 Squares which remain squares when the last digit is removed.

Original entry on oeis.org

0, 1, 4, 9, 16, 49, 169, 256, 361, 1444, 3249, 18496, 64009, 237169, 364816, 519841, 2079364, 4678569, 26666896, 92294449, 341991049, 526060096, 749609641, 2998438564, 6746486769, 38453641216, 133088524969, 493150849009, 758578289296, 1080936581761

Offset: 1

David W. Wilson

This A023110 = A031149^2 is the base 10 version of A001541^2 = A055792 (base 2), A001075^2 = A055793 (base 3), A004275^2 = A055808 (base 4), A204520^2 = A055812 (base 5), A204518^2 = A055851 (base 6), A204516^2 = A055859 (base 7), A204514^2 = A055872 (base 8) and A204502^2 = A204503 (base 9). - M. F. Hasler, Sep 28 2014

For the first 4 terms the square has only one digit. It is understood that deleting this digit yields 0. - Colin Barker, Dec 31 2017

R. K. Guy, Neg and Reg, preprint, Jan 2012.

Jon E. Schoenfield, Table of n, a(n) for n = 1..70 (terms 1..38 from David W. Wilson, terms 39..40 from Robert G. Wilson v, terms 41..67 from Dmitry Petukhov)
M. F. Hasler, Truncated squares, OEIS wiki, Jan 16 2012
Joshua Stucky, Pell's Equation and Truncated Squares, Number Theory Seminar, Kansas State University, Feb 19 2018.
Index to sequences related to truncating digits of squares.

Cf. A023111.
Cf. A031150, A053784, A031149, A055792, A055793, A055808, A055812, A055851, A055859, A055872.
Cf. A001541, A001075, A004275, A204520, A204518, A204516, A204514, A204502, A204503.

count:= 1: A[1]:= 0:
for n from 0 while count < 35 do
  for t in [1,4,6,9] do
    if issqr(10*n^2+t) then
       count:= count+1;
       A[count]:= 10*n^2+t;
    fi
  od
od:
seq(A[i],i=1..count); # Robert Israel, Sep 28 2014

fQ[n_] := IntegerQ@ Sqrt@ Quotient[n^2, 10]; Select[ Range@ 1000000, fQ]^2 (* Robert G. Wilson v, Jan 15 2011 *)

for(n=0,1e7, issquare(n^2\10) & print1(n^2",")) \\  M. F. Hasler, Jan 16 2012

Appears to satisfy a(n)=1444*a(n-7)+a(n-14)-76*sqrt(a(n-7)*a(n-14)) for n >= 16. For n = 15, 14, 13, ... this would require a(1) = 16, a(0) = 49, a(-1) = 169, ... - Henry Bottomley, May 08 2001; edited by Robert Israel, Sep 28 2014
a(n) = A031149(n)^2. - M. F. Hasler, Sep 28 2014
Conjectures from Colin Barker, Dec 31 2017: (Start)
G.f.: x^2*(1 + 4*x + 9*x^2 + 16*x^3 + 49*x^4 + 169*x^5 + 256*x^6 - 1082*x^7 - 4328*x^8 - 9738*x^9 - 4592*x^10 - 6698*x^11 - 6698*x^12 - 4592*x^13 + 361*x^14 + 1444*x^15 + 3249*x^16 + 256*x^17 + 169*x^18 + 49*x^19 + 16*x^20) / ((1 - x)*(1 + x + x^2 + x^3 + x^4 + x^5 + x^6)*(1 - 1442*x^7 + x^14)).
a(n) = 1443*a(n-7) - 1443*a(n-14) + a(n-21) for n>22.
(End)

More terms from M. F. Hasler, Jan 16 2012

A055793 Numbers k such that k and floor[k/3] are both squares; i.e., squares which remain squares when written in base 3 and last digit is removed.

Original entry on oeis.org

0, 1, 4, 49, 676, 9409, 131044, 1825201, 25421764, 354079489, 4931691076, 68689595569, 956722646884, 13325427460801, 185599261804324, 2585064237799729, 36005300067391876, 501489136705686529, 6984842613812219524, 97286307456665386801, 1355023461779503195684, 18873042157456379352769

Offset: 1

Henry Bottomley, Jul 14 2000

Or, squares of the form 3k^2+1.

See A023110, A204503, A204512, A204517, A204519, A055812, A055808 and A055792 for the analog in other bases.

			a(3) = 49 because 49 = 7^2 = 1211 base 3 and 121 base 3 = 16 = 4^2.

Vincenzo Librandi, Table of n, a(n) for n = 1..800
Tom C. Brown and Peter J Shiue, Squares of second-order linear recurrence sequences, Fib. Quart., 33 (1994), 352-356.
Yurii S. Bystryk, Vitalii L. Denysenko, and Volodymyr I. Ostryk, Lune and Lens Sequences, ResearchGate preprint, 2024. See pp. 48, 56.
M. F. Hasler, Truncated squares, OEIS wiki, Jan 16 2012
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Circular Segment, Forum Geometricorum, Vol. 18 (2018), 47-55.
Index to sequences related to truncating digits of squares.
Index entries for linear recurrences with constant coefficients, signature (15,-15,1).

Cf. A001075, A023110, A098301.

Cf. also A023110, A204503, A204512, A204517, A204519, A055812, A055808 and A055792 for the analog in other bases.

I:=[0, 1, 4]; [n le 3 select I[n] else 14*Self(n-1) - Self(n-2) - 6: n in [1..30]]; // Vincenzo Librandi, Jan 27 2013

A055793 := proc(n) coeftayl(x*(1-11*x+4*x^2)/((1-x)*(1-14*x+x^2)), x=0, n); end proc: seq(A055793(n), n=0..20); # Wesley Ivan Hurt, Sep 28 2014

CoefficientList[Series[x*(1 - 11*x + 4*x^2)/((1 - x)*(1 - 14*x + x^2)), {x, 0, 20}], x] (* Wesley Ivan Hurt, Sep 28 2014 *)
LinearRecurrence[{15,-15,1},{0,1,4,49},40] (* Harvey P. Dale, Jun 19 2021 *)

sq3nsqplus1(n) = { for(x=1,n, y = 3*x*x+1; \ print1(y" ") if(issquare(y),print1(y" ")) ) }

a(n) = 3*A098301(n-2)+1. - R. J. Mathar, Jun 11 2009
a(n) = 14*a(n-1)-a(n-2)-6, with a(0)=1, a(1)=4. (See Brown and Shiue)
a(n) = (A001075(n-2))^2. - Johannes Boot Dec 16 2011, corrected by M. F. Hasler, Jan 15 2012
G.f.: x*(1 - 11*x + 4*x^2)/((1 - x)*(1 - 14*x + x^2)). - M. F. Hasler, Jan 15 2012

More terms from Cino Hilliard, Mar 01 2003

A132592 X-values of solutions to the equation X(X + 1) - 8Y^2 = 0.

Original entry on oeis.org

0, 8, 288, 9800, 332928, 11309768, 384199200, 13051463048, 443365544448, 15061377048200, 511643454094368, 17380816062160328, 590436102659356800, 20057446674355970888, 681362750825443653408, 23146276081390728245000, 786292024016459316676608, 26710782540478226038759688

Offset: 0

Mohamed Bouhamida, Nov 14 2007

Equivalently, numbers k such that both k/2 and k+1 are squares. - Karl-Heinz Hofmann, Sep 20 2022

Equivalently, numbers k such that the k-dimensional volume and total (k-1)-dimensional volume are equal, with side length being a positive integer, for all regular polyhedra constructible in k dimensions. - Matt Moir, Jul 09 2024

Seiichi Manyama, Table of n, a(n) for n = 0..500
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).

Intersection between A132411 and A001105.
Cf. A007654.
Cf. A001541, A058331, A001079, A037270, A055792, A071253, A108741, A132592, A146311, A146312, A146313, A173115, A173116, A173121.

I:=[0,8,288]; [n le 3 select I[n] else 35*Self(n-1)-35*Self(n-2)+ Self(n-3): n in [1..30]]; // Vincenzo Librandi, Dec 24 2018

Table[Round[N[Sinh[2 n ArcCosh[Sqrt[2]]]^2, 100]], {n, 0, 20}] (* Artur Jasinski, Feb 10 2010 *)
LinearRecurrence[{35, -35, 1}, {0, 8, 288}, 30] (* Vincenzo Librandi, Dec 24 2018 *)

A132592 = [0, 8]
for n in range(2, 18): A132592.append(34 * A132592[-1] - A132592[-2] + 16)
print(A132592) # Karl-Heinz Hofmann, Sep 20 2022

a(0)=0, a(1)=8 and a(n) = 34*a(n-1) - a(n-2) + 16.
a(n) = (A056771(n) - 1)/2. - Max Alekseyev, Nov 13 2009
a(n) = sinh(2*n*arccosh(sqrt(2))^2) (n=0,1,2,3,...). - Artur Jasinski, Feb 10 2010
G.f.: -8*x*(x+1)/((x-1)*(x^2-34*x+1)). - Colin Barker, Oct 24 2012
a(n) = A055792(n+1)-1 = A001541(n)^2 - 1. - Antti Karttunen, Oct 03 2016

More terms from Max Alekseyev, Nov 13 2009

A031149 Numbers whose square with its last digit deleted is also a square.

Original entry on oeis.org

0, 1, 2, 3, 4, 7, 13, 16, 19, 38, 57, 136, 253, 487, 604, 721, 1442, 2163, 5164, 9607, 18493, 22936, 27379, 54758, 82137, 196096, 364813, 702247, 870964, 1039681, 2079362, 3119043, 7446484, 13853287, 26666893, 33073696, 39480499, 78960998, 118441497, 282770296

Offset: 1

Patrick De Geest

Square root of A023110(n).
For the first 4 terms, the square has only one digit, but in analogy to the sequences in other bases (A204502, A204512,  A204514, A204516, A204518, A204520, A004275, A055793, A055792), it is understood that deleting this digit yields 0.
From Robert Israel, Feb 16 2016: (Start)
Solutions x to x^2 = 10*y^2 + j, j in {0,1,4,6,9}, in increasing order of x.
j=0 occurs only for x=0.
Let M be the 2 X 2 matrix [19, 60; 6, 19].
Solutions of x^2 = 10*y^2 + 1 are (x,y)^T = M^k (1,0)^T for k >= 0.
Solutions of x^2 = 10*y^2 + 4 are (x,y)^T = M^k (2,0)^T for k >= 0.
Solutions of x^2 = 10*y^2 + 6 are (x,y)^T = M^k (4,1)^T and M^k (16,5)^T for k >= 0.
Solutions of x^2 = 10*y^2 + 9 are (x,y)^T = M^k (3,0)^T, M^k (7,2)^T, M^k (13,4)^T for k >= 0.
Since (1,0)^T <= (2,0)^T <= (3,0)^T <= (4,1)^T <= (7,2)^T <= (13,4)^T <= (16,5)^T <= (19,6)^T = M (1,0)^T (element-wise) and M has positive entries, we see that the terms always occur in this order, for successive k.
The eigenvalues of M are 19 + 6*sqrt(10) and 19 - 6*sqrt(10).
From this follow my formulas below and the G.f. (End)

			364813^2 = 133088524969, 115364^2 = 13308852496.

R. K. Guy, Neg and Reg, preprint, Jan 2012. [From N. J. A. Sloane, Jan 12 2012]

Dmitry Petukhov and Robert Israel, Table of n, a(n) for n = 1..4400 (n = 1..67 from Dmitry Petukhov)
M. F. Hasler, Truncated squares, OEIS wiki, Jan 16 2012
Index to sequences related to truncating digits of squares.

Cf. A000290, A023110, A030686, A030687, A031150.

for i from 1 to 150000 do if (floor(sqrt(10 * i^2 + 9)) > floor(sqrt(10 * i^2))) then print(floor(sqrt(10 * i^2 + 9))) end if end do;

fQ[n_] := IntegerQ@ Sqrt@ Quotient[n^2, 10]; Select[ Range[ 0, 40000000], fQ] (* Harvey P. Dale, Jun 15 2011 *) (* modified by Robert G. Wilson v, Jan 16 2012 *)

s=[]; for(n=0, 1e7, if(issquare(n^2\10), s=concat(s,n))); s \\ Colin Barker, Jan 17 2014; typo fixed by Zak Seidov, Jan 31 2014

Appears to satisfy: a(n)=38a(n-7)-a(n-14) which would require a(-k) to look like 3, 2, 1, 4, 7, 13, 16, 57, 38, 19, 136, ... for k>0. - Henry Bottomley, May 08 2001
Empirical g.f.: x^2*(1 + 2*x + 3*x^2 + 4*x^3 + 7*x^4 + 13*x^5 + 16*x^6 - 19*x^7 - 38*x^8 - 57*x^9 - 16*x^10 - 13*x^11 - 7*x^12 - 4*x^13) / ((1 - 38*x^7 + x^14)). - Colin Barker, Jan 17 2014
a(n) = 38*a(n-7) - a(n-14) for n>15 (conjectured). - Colin Barker, Dec 31 2017
With e1 = 19 + 6*sqrt(10) and e2 = 19 - 6*sqrt(10),
   a(2+7k) = (e1^k + e2^k)/2,
   a(3+7k) = e1^k + e2^k,
   a(4+7k) = (3/2) (e1^k + e2^k),
   a(5+7k) = (2+sqrt(10)/2) e1^k + (2-sqrt(10)/2) e2^k,
   a(6+7k) = (7/2+sqrt(10)) e1^k + (7/2-sqrt(10)) e2^k,
   a(7+7k) = (13/2+2 sqrt(10)) e1^k + (13/2-2 sqrt(10)) e2^k,
a(8+7k) = (8+5 sqrt(10)/2) e1^k + (8-5 sqrt(10)/2) e2^k. - Robert Israel, Feb 16 2016

4 initial terms added by M. F. Hasler, Jan 15 2012

a(40) from Robert G. Wilson v, Jan 15 2012

A055812 a(n) and floor(a(n)/5) are both squares; i.e., squares which remain squares when written in base 5 and last digit is removed.

Original entry on oeis.org

0, 1, 4, 9, 49, 81, 324, 2209, 15129, 25921, 103684, 710649, 4870849, 8346321, 33385284, 228826129, 1568397609, 2687489281, 10749957124, 73681302249, 505019158609, 865363202001, 3461452808004

Offset: 1

Henry Bottomley, Jul 14 2000

For the first 3 terms, the above "base 5" interpretation is questionable, since they have only 1 digit in base 5. It is understood that dropping this digit yields 0. - M. F. Hasler, Jan 15 2012

			a(4) = 49 because 49 = 7^2 = 144 base 5 and 14 base 5 = 9 = 3^2.

For analogs in other bases see A055792 (base 2), A055793 (base 3), A055808 (base 4), A055851 (base 6), A204517 (base 7), A204512 (base 8), A204503 (base 9) and A023110 (base 10).

Squares of A204520. The square roots of floor[a(n)/5] are given in A204521.

b=5;for(n=1,2e9,issquare(n^2\b) && print1(n^2,","))  \\ M. F. Hasler, Jan 15 2012

Empirical g.f.: -x^2*(9*x^11 +49*x^10 +324*x^9 +81*x^8 -698*x^7 -698*x^6 -968*x^5 -242*x^4 +49*x^3 +9*x^2 +4*x +1) / ((x -1)*(x +1)*(x^2 -4*x -1)*(x^2 +1)*(x^2 +4*x -1)*(x^4 +18*x^2 +1)). - Colin Barker, Sep 15 2014

More terms added and offset changed to 1 by M. F. Hasler, Jan 15 2012

A204502 Numbers such that floor[a(n)^2 / 9] is a square.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63, 66, 69, 72, 75, 78, 81, 84, 87, 90, 93, 96, 99, 102, 105, 108, 111, 114, 117, 120, 123, 126, 129, 132, 135, 138, 141, 144, 147, 150, 153, 156, 159, 162, 165, 168, 171, 174, 177

Offset: 1

M. F. Hasler, Jan 15 2012

nonn

Or, numbers n such that n^2, with its last base-9 digit dropped, is again a square. (Except maybe for the 3 initial terms whose square has only 1 digit in base 9.)

M. F. Hasler, Truncated squares, OEIS wiki, Jan 16 2012
Index to sequences related to truncating digits of squares.

The squares are in A204503, the squares with last base-9 digit dropped in A204504, and the square roots of the latter in A028310.

Cf. A031149=sqrt(A023110) (base 10), A204514=sqrt(A055872) (base 8), A204516=sqrt(A055859) (base 7), A204518=sqrt(A055851) (base 6), A204520=sqrt(A055812) (base 5), A004275=sqrt(A055808) (base 4), A001075=sqrt(A055793) (base 3), A001541=sqrt(A055792) (base 2).

Select[Range[0,200],IntegerQ[Sqrt[Floor[#^2/9]]]&] (* Harvey P. Dale, May 05 2018 *)

b=9;for(n=0,200,issquare(n^2\b) & print1(n","))

Conjecture: a(n) = 3*n-12 for n>5. G.f.: x^2*(x^2+x+1)*(x^3-x+1)/(x-1)^2. [Colin Barker, Nov 23 2012]

cp's OEIS Frontend

A204575 Squares such that [a(n)/2] is again a square (A055792), written in binary.

Views

Author

Keywords

Crossrefs

A204576 Floor[A055792(n-1)/2]=A084703(n-2) (truncated squares), written in binary.

Views

Author

Keywords

Comments

Links

Crossrefs

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Sage

Sage

Formula

Extensions

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

A023110 Squares which remain squares when the last digit is removed.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A055793 Numbers k such that k and floor[k/3] are both squares; i.e., squares which remain squares when written in base 3 and last digit is removed.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

A132592 X-values of solutions to the equation X(X + 1) - 8Y^2 = 0.