A001108 -id:A001108 - cp's OEIS frontend

A001542 a(n) = 6*a(n-1) - a(n-2) for n > 1, a(0)=0 and a(1)=2.

0, 2, 12, 70, 408, 2378, 13860, 80782, 470832, 2744210, 15994428, 93222358, 543339720, 3166815962, 18457556052, 107578520350, 627013566048, 3654502875938, 21300003689580, 124145519261542, 723573111879672

Offset: 0

N. J. A. Sloane

Consider the equation core(x) = core(2x+1) where core(x) is the smallest number such that x*core(x) is a square: solutions are given by a(n)^2, n > 0. - Benoit Cloitre, Apr 06 2002
Terms > 0 give numbers k which are solutions to the inequality |round(sqrt(2)*k)/k - sqrt(2)| < 1/(2*sqrt(2)*k^2). - Benoit Cloitre, Feb 06 2006
Also numbers m such that A125650(6*m^2) is an even perfect square, where A124650(m) is a numerator of m*(m+3)/(4*(m+1)*(m+2)) = Sum_{k=1..m} 1/(k*(k+1)*(k+2)). Sequence A033581 is a bisection of A125651. - Alexander Adamchuk, Nov 30 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 = {a(n)}. - Clark Kimberling, Aug 26 2008
Even Pell numbers. - Omar E. Pol, Dec 10 2008
Numbers k such that 2*k^2+1 is a square. - Vladimir Joseph Stephan Orlovsky, Feb 19 2010
These are the integer square roots of the Half-Squares, A007590(k), which occur at values of k given by A001541. Also the numbers produced by adding m + sqrt(floor(m^2/2) + 1) when m is in A002315. See array in A227972. - Richard R. Forberg, Aug 31 2013
A001541(n)/a(n) is the closest rational approximation of sqrt(2) with a denominator not larger than a(n), and 2*a(n)/A001541(n) is the closest rational approximation of sqrt(2) with a numerator not larger than 2*a(n). These rational approximations together with those obtained from the sequences A001653 and A002315 give a complete set of closest rational approximations of sqrt(2) with restricted numerator as well as denominator. - A.H.M. Smeets, May 28 2017
Conjecture: Numbers k such that c/m < k for all natural a^2 + b^2 = c^2 (Pythagorean triples), a < b < c and a+b+c = m. Numbers which correspondingly minimize c/m are A002939. - Lorraine Lee, Jan 31 2020
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=a(n), b=A005319(n), c=a(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0 < n. - Michael Somos, Jun 26 2022

			G.f. = 2*x + 12*x^2 + 70*x^3 + 408*x^4 + 2378*x^5 + 13860*x^6 + ...

Jay Kappraff, Beyond Measure, A Guided Tour Through Nature, Myth and Number, World Scientific, 2002; pp. 480-481.
Thomas Koshy, Fibonacci and Lucas Numbers with Applications, 2001, Wiley, pp. 77-79.
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..100
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.
Hacène Belbachir, Soumeya Merwa Tebtoub, 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.
A. J. C. Cunningham, Binomial Factorisations, Vols. 1-9, Hodgson, London, 1923-1929. See Vol. 1, page xxxv.
S. Falcon, Relationships between Some k-Fibonacci Sequences, Applied Mathematics, 2014, 5, 2226-2234.
R. J. Hetherington, Letter to N. J. A. Sloane, Oct 26 1974.
J. M. Katri and D. R. Byrkit, Problem E1976, Amer. Math. Monthly, 75 (1968), 683-684.
Tanya Khovanova, Recursive Sequences.
E. Kilic, Y. T. Ulutas, and N. Omur, A Formula for the Generating Functions of Powers of Horadam's Sequence with Two Additional Parameters, J. Int. Seq. 14 (2011) #11.5.6, Table 3, k=1.
D. H. Lehmer, On the multiple solutions of the Pell equation, Annals Math., 30 (1928), 66-72.
D. H. Lehmer, On the multiple solutions of the Pell equation (annotated scanned copy).
Mathematical Reflections, Solution to Problem O271, Issue 5, 2013, p 22.
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.
Mark A. Shattuck, Tiling proofs of some formulas for the Pell numbers of odd index, Integers, 9 (2009), 53-64.
R. A. Sulanke, Moments, Narayana numbers and the cut and paste for lattice paths.
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.
Index entries for linear recurrences with constant coefficients, signature (6,-1).

Bisection of Pell numbers A000129: {a(n)} and A001653(n+1), n >= 0.

Cf. A001108, A001353, A001541, A001835, A003499, A007805, A007913, A115598, A125650, A125651, A125652.

a:=[0,2];; for n in [3..20] do a[n]:=6*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019

a001542 n = a001542_list !! n
a001542_list =
   0 : 2 : zipWith (-) (map (6 *) $ tail a001542_list) a001542_list
-- Reinhard Zumkeller, Aug 14 2011

I:=[0,2]; [n le 2 select I[n] else 6*Self(n-1) -Self(n-2): n in [1..20]]; // G. C. Greubel, Dec 23 2019

A001542:=2*z/(1-6*z+z**2); # conjectured by Simon Plouffe in his 1992 dissertation
seq(combinat:-fibonacci(2*n, 2), n = 0..20); # Peter Luschny, Jun 28 2018

LinearRecurrence[{6, -1}, {0, 2}, 30] (* Harvey P. Dale, Jun 11 2011 *)
Fibonacci[2*Range[0,20], 2] (* G. C. Greubel, Dec 23 2019 *)
Table[2 ChebyshevU[-1 + n, 3], {n, 0, 20}] (* Herbert Kociemba, Jun 05 2022 *)

a[0]:0$
a[1]:2$
a[n]:=6*a[n-1]-a[n-2]$
A001542(n):=a[n]$
makelist(A001542(x),x,0,30); /* Martin Ettl, Nov 03 2012 */

{a(n) = imag( (3 + 2*quadgen(8))^n )}; /* Michael Somos, Jan 20 2017 */

vector(21, n, 2*polchebyshev(n-1, 2, 33) ) \\ G. C. Greubel, Dec 23 2019

l=[0, 2]
for n in range(2, 51): l+=[6*l[n - 1] - l[n - 2], ]
print(l) # Indranil Ghosh, Jun 06 2017

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

a(n) = 2*A001109(n).
a(n) = ((3+2*sqrt(2))^n - (3-2*sqrt(2))^n) / (2*sqrt(2)).
G.f.: 2*x/(1-6*x+x^2).
a(n) = sqrt(2*(A001541(n))^2 - 2)/2. - Barry E. Williams, May 07 2000
a(n) = (C^(2n) - C^(-2n))/sqrt(8) where C = sqrt(2) + 1. - Gary W. Adamson, May 11 2003
For all terms x of the sequence, 2*x^2 + 1 is a square. Limit_{n->oo} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 10 2002
For n > 0: a(n) = A001652(n) + A046090(n) - A001653(n); e.g., 70 = 119 + 120 - 169. Also a(n) = A001652(n - 1) + A046090(n - 1) + A001653(n - 1); e.g., 70 = 20 + 21 + 29. Also a(n)^2 + 1 = A001653(n - 1)*A001653(n); e.g., 12^2 + 1 = 145 = 5*29. Also a(n + 1)^2 = A084703(n + 1) = A001652(n)*A001652(n + 1) + A046090(n)*A046090(n + 1). - Charlie Marion, Jul 01 2003
a(n) = ((1+sqrt(2))^(2*n) - (1-sqrt(2))^(2*n))/(2*sqrt(2)). - Antonio Alberto Olivares, Dec 24 2003
2*A001541(k)*A001653(n)*A001653(n+k) = A001653(n)^2 + A001653(n+k)^2 + a(k)^2; e.g., 2*3*5*29 = 5^2 + 29^2 + 2^2; 2*99*29*5741 = 29^2 + 5741^2 + 70^2. - Charlie Marion, Oct 12 2007
a(n) = sinh(2*n*arcsinh(1))/sqrt(2). - Herbert Kociemba, Apr 24 2008
For n > 0, a(n) = A001653(n) + A002315(n-1). - Richard R. Forberg, Aug 31 2013
a(n) = 3*a(n-1) + 2*A001541(n-1); e.g., a(4) = 70 = 3*12 + 2*17. - Zak Seidov, Dec 19 2013
a(n)^2 + 1^2 = A115598(n)^2 + (A115598(n)+1)^2. - Hermann Stamm-Wilbrandt, Jul 27 2014
Sum {n >= 1} 1/( a(n) + 1/a(n) ) = 1/2. - _Peter Bala, Mar 25 2015
E.g.f.: exp(3*x)*sinh(2*sqrt(2)*x)/sqrt(2). - Ilya Gutkovskiy, Dec 07 2016
A007814(a(n)) = A001511(n). See Mathematical Reflections link. - Michel Marcus, Jan 06 2017
a(n) = -a(-n) for all n in Z. - Michael Somos, Jan 20 2017
From A.H.M. Smeets, May 28 2017: (Start)
A051009(n) = a(2^(n-2)).
a(2n) = 2*a(2)*A001541(n).
A001541(n)/a(n) > sqrt(2) > 2*a(n)/A001541(n). (End)
a(A298210(n)) = A002349(2*n^2). - A.H.M. Smeets, Jan 25 2018
a(n) = A000129(n)*A002203(n). - Adam Mohamed, Jul 20 2024

A053141 a(0)=0, a(1)=2 then a(n) = a(n-2) + 2sqrt(8a(n-1)^2 + 8*a(n-1) + 1).

Original entry on oeis.org

0, 2, 14, 84, 492, 2870, 16730, 97512, 568344, 3312554, 19306982, 112529340, 655869060, 3822685022, 22280241074, 129858761424, 756872327472, 4411375203410, 25711378892990, 149856898154532, 873430010034204, 5090723162050694, 29670908962269962, 172934730611569080

Offset: 0

Wolfdieter Lang

Solution to b(b+1) = 2a(a+1) in natural numbers including 0; a = a(n), b = b(n) = A001652(n).
The solution of a special case of a binomial problem of H. Finner and K. Strassburger (strass(AT)godot.dfi.uni-duesseldorf.de).
Also the indices of triangular numbers that are half other triangular numbers [a of T(a) such that 2T(a)=T(b)]. The T(a)'s are in A075528, the T(b)'s are in A029549 and the b's are in A001652. - Bruce Corrigan (scentman(AT)myfamily.com), Oct 30 2002
Sequences A053141 (this entry), A016278, A077259, A077288 and A077398 are part of an infinite series of sequences. Each depends upon the polynomial p(n) = 4k*n^2 + 4k*n + 1, when 4k is not a perfect square. Equivalently, they each depend on the equation k*t(x)=t(z) where t(n) is the triangular number formula n(n+1)/2. The dependencies are these: they are the sequences of positive integers n such that p(n) is a perfect square and there exists a positive integer m such that k*t(n)=t(m). A053141 is for k=2, A016278 is for k=3, A077259 is for k=5. - Robert Phillips (bobanne(AT)bellsouth.net), Oct 11 2007, Nov 27 2007
Jason Holt observes that a pair drawn from a drawer with A053141(n)+1 red socks and A001652(n) - A053141(n) blue socks will as likely as not be matching reds: (A053141+1)*A053141/((A001652+1)*A001652) = 1/2, n>0. - Bill Gosper, Feb 07 2010
The values x(n)=A001652(n), y(n)=A046090(n) and z(n)=A001653(n) form a nearly isosceles Pythagorean triple since y(n)=x(n)+1 and x(n)^2 + y(n)^2 = z(n)^2; e.g., for n=2, 20^2 + 21^2 = 29^2. In a similar fashion, if we define b(n)=A011900(n) and c(n)=A001652(n), a(n), b(n) and c(n) form a nearly isosceles anti-Pythagorean triple since b(n)=a(n)+1 and a(n)^2 + b(n)^2 = c(n)^2 + c(n) + 1; i.e., the value a(n)^2 + b(n)^2 lies almost exactly between two perfect squares; e.g., 2^2 + 3^2 = 13 = 4^2 - 3 = 3^2 + 4; 14^2 + 15^2 = 421 = 21^2 - 20 = 20^2 + 21. - Charlie Marion, Jun 12 2009
Behera & Panda call these the balancers and A001109 are the balancing numbers. - Michel Marcus, Nov 07 2017

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Jeremiah Bartz, Bruce Dearden, and Joel Iiams, Classes of Gap Balancing Numbers, arXiv:1810.07895 [math.NT], 2018.
Jeremiah Bartz, Bruce Dearden, and Joel Iiams, Counting families of generalized balancing numbers, The Australasian Journal of Combinatorics (2020) Vol. 77, Part 3, 318-325.
A. Behera and G. K. Panda, On the Square Roots of Triangular Numbers, Fib. Quart., 37 (1999), pp. 98-105.
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.
P. Catarino, H. Campos, and P. Vasco, On some identities for balancing and cobalancing numbers, Annales Mathematicae et Informaticae, 45 (2015) pp. 11-24.
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.
aBa Mbirika, Janee Schrader, and Jürgen Spilker, Pell and Associated Pell Braid Sequences as GCDs of Sums of k Consecutive Pell, Balancing, and Related Numbers, J. Int. Seq. (2023) Vol. 26, Art. 23.6.4.
J. S. Myers, R. Schroeppel, S. R. Shannon, N. J. A. Sloane, and P. Zimmermann, Three Cousins of Recaman's Sequence, arXiv:2004:14000 [math.NT], April 2020.
G. K. Panda, Sequence balancing and cobalancing numbers, Fib. Q., Vol. 45, No. 3 (2007), 265-271. See p. 266.
Michael Penn, (co) balancing numbers, YouTube video, 2022.
Robert Phillips, Polynomials of the form 1+4ke+4ke^2, 2008.
Robert Phillips, A triangular number result, 2009.
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.
Burkard Polster, Nice merging together, Mathologer video (2015).
B. Polster and M. Ross, Marching in squares, arXiv preprint arXiv:1503.04658 [math.HO], 2015.
A. Tekcan, M. Tayat, and M. 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.
Index entries for linear recurrences with constant coefficients, signature (7,-7,1).

Cf. A000129, A001108, A001109, A001652, A001653.
Cf. A011900, A029549, A053142, A075528, A103200.
Partial sums of A001542.
Cf. A000194, A001541, A002024, A016278, A046090, A055997, A077259, A077288, A077398, A084068.

a053141 n = a053141_list !! n
a053141_list = 0 : 2 : map (+ 2)
   (zipWith (-) (map (* 6) (tail a053141_list)) a053141_list)
-- Reinhard Zumkeller, Jan 10 2012

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

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

Join[{a=0,b=1}, Table[c=6*b-a+1; a=b; b=c, {n,60}]]*2 (* Vladimir Joseph Stephan Orlovsky, Jan 18 2011 *)
a[n_] := Floor[1/8*(2+Sqrt[2])*(3+2*Sqrt[2])^n]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Nov 28 2013 *)
Table[(Fibonacci[2n + 1, 2] - 1)/2, {n, 0, 20}] (* Vladimir Reshetnikov, Sep 16 2016 *)

concat(0,Vec(2/(1-x)/(1-6*x+x^2)+O(x^30))) \\ Charles R Greathouse IV, May 14 2012

{x=1+sqrt(2); y=1-sqrt(2); P(n) = (x^n - y^n)/(x-y)};
a(n) = round((P(2*n+1) - 1)/2);
for(n=0, 30, print1(a(n), ", ")) \\ G. C. Greubel, Jul 15 2018

[(lucas_number1(2*n+1, 2, -1)-1)/2 for n in range(30)] # G. C. Greubel, Apr 27 2020

a(n) = (A001653(n)-1)/2 = 2*A053142(n) = A011900(n)-1. [Corrected by Pontus von Brömssen, Sep 11 2024]
a(n) = 6*a(n-1) - a(n-2) + 2, a(0) = 0, a(1) = 2.
G.f.: 2*x/((1-x)*(1-6*x+x^2)).
Let c(n) = A001109(n). Then a(n+1) = a(n)+2*c(n+1), a(0)=0. This gives a generating function (same as existing g.f.) leading to a closed form: a(n) = (1/8)*(-4+(2+sqrt(2))*(3+2*sqrt(2))^n + (2-sqrt(2))*(3-2*sqrt(2))^n). - Bruce Corrigan (scentman(AT)myfamily.com), Oct 30 2002
a(n) = 2*Sum_{k = 0..n} A001109(k). - Mario Catalani (mario.catalani(AT)unito.it), Mar 22 2003
For n>=1, a(n) = 2*Sum_{k=0..n-1} (n-k)*A001653(k). - Charlie Marion, Jul 01 2003
For n and j >= 1, A001109(j+1)*A001652(n) - A001109(j)*A001652(n-1) + a(j) = A001652(n+j). - Charlie Marion, Jul 07 2003
From Antonio Alberto Olivares, Jan 13 2004: (Start)
a(n) = 7*a(n-1) - 7*a(n-2) + a(n-3).
a(n) = -(1/2) - (1-sqrt(2))/(4*sqrt(2))*(3-2*sqrt(2))^n + (1+sqrt(2))/(4*sqrt(2))*(3+2*sqrt(2))^n. (End)
a(n) = sqrt(2)*cosh((2*n+1)*log(1+sqrt(2)))/4 - 1/2 = (sqrt(1+4*A029549)-1)/2. - Bill Gosper, Feb 07 2010 [typo corrected by Vaclav Kotesovec, Feb 05 2016]
a(n+1) + A055997(n+1) = A001541(n+1) + A001109(n+1). - Creighton Dement, Sep 16 2004
From Charlie Marion, Oct 18 2004: (Start)
For n>k, a(n-k-1) = A001541(n)*A001653(k)-A011900(n+k); e.g., 2 = 99*5 - 493.
For n<=k, a(k-n) = A001541(n)*A001653(k) - A011900(n+k); e.g., 2 = 3*29 - 85 + 2. (End)
a(n) = A084068(n)*A084068(n+1). - Kenneth J Ramsey, Aug 16 2007
Let G(n,m) = (2*m+1)*a(n)+ m and H(n,m) = (2*m+1)*b(n)+m where b(n) is from the sequence A001652 and let T(a) = a*(a+1)/2. Then T(G(n,m)) + T(m) = 2*T(H(n,m)). - Kenneth J Ramsey, Aug 16 2007
Let S(n) equal the average of two adjacent terms of G(n,m) as defined immediately above and B(n) be one half the difference of the same adjacent terms. Then for T(i) = triangular number i*(i+1)/2, T(S(n)) - T(m) = B(n)^2 (setting m = 0 gives the square triangular numbers). - Kenneth J Ramsey, Aug 16 2007
a(n) = A001108(n+1) - A001109(n+1). - Dylan Hamilton, Nov 25 2010
a(n) = (a(n-1)*(a(n-1) - 2))/a(n-2) for n > 2. - Vladimir Pletser, Apr 08 2020
a(n) = (ChebyshevU(n, 3) - ChebyshevU(n-1, 3) - 1)/2 = (Pell(2*n+1) - 1)/2. - G. C. Greubel, Apr 27 2020
E.g.f.: (exp(3*x)*(2*cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x)) - 2*exp(x))/4. - Stefano Spezia, Mar 16 2024
a(n) = A000194(A029549(n)) = A002024(A075528(n)). - Pontus von Brömssen, Sep 11 2024

Name corrected by Zak Seidov, Apr 11 2011

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

A281871 Number T(n,k) of k-element subsets of [n] having a square element sum; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 1, 0, 1, 2, 2, 2, 0, 0, 1, 2, 3, 3, 2, 1, 0, 1, 2, 4, 5, 5, 2, 1, 0, 1, 2, 5, 8, 8, 6, 3, 0, 1, 1, 3, 6, 11, 14, 13, 7, 4, 1, 0, 1, 3, 7, 15, 23, 24, 19, 10, 3, 1, 0, 1, 3, 8, 20, 34, 43, 39, 25, 13, 3, 1, 0, 1, 3, 9, 26, 49, 71, 74, 60, 34, 14, 5, 0, 0

Offset: 0

Alois P. Heinz, Jan 31 2017

			T(7,0) = 1: {}.
T(7,1) = 2: {1}, {4}.
T(7,2) = 4: {1,3}, {2,7}, {3,6}, {4,5}.
T(7,3) = 5: {1,2,6}, {1,3,5}, {2,3,4}, {3,6,7}, {4,5,7}.
T(7,4) = 5: {1,2,6,7}, {1,3,5,7}, {1,4,5,6}, {2,3,4,7}, {2,3,5,6}.
T(7,5) = 2: {1,2,3,4,6}, {3,4,5,6,7}.
T(7,6) = 1: {1,2,4,5,6,7}.
T(7,7) = 0.
T(8,8) = 1: {1,2,3,4,5,6,7,8}.
Triangle T(n,k) begins:
  1;
  1, 1;
  1, 1, 0;
  1, 1, 1,  0;
  1, 2, 1,  1,  0;
  1, 2, 2,  2,  0,  0;
  1, 2, 3,  3,  2,  1,  0;
  1, 2, 4,  5,  5,  2,  1,  0;
  1, 2, 5,  8,  8,  6,  3,  0,  1;
  1, 3, 6, 11, 14, 13,  7,  4,  1,  0;
  1, 3, 7, 15, 23, 24, 19, 10,  3,  1, 0;
  1, 3, 8, 20, 34, 43, 39, 25, 13,  3, 1, 0;
  1, 3, 9, 26, 49, 71, 74, 60, 34, 14, 5, 0, 0;
  ...

Alois P. Heinz, Rows n = 0..200, flattened

Columns k=0-10 give: A000012, A000196, A176615, A281706, A281864, A281865, A281866, A281867, A281868, A281869, A281870.
Main diagonal is characteristic function of A001108.
Diagonals T(n+k,n) for k=2-10 give: A281965, A281966, A281967, A281968, A281969, A281970, A281971, A281972, A281973.
Row sums give A126024.
T(2n,n) gives A281872.
Cf. A000217, A000290, A007318, A214857, A214858, A278339, A281994, A284249, A377572.

b:= proc(n, s) option remember; expand(`if`(n=0,
      `if`(issqr(s), 1, 0), b(n-1, s)+x*b(n-1, s+n)))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n, 0)):
seq(T(n), n=0..16);

b[n_, s_] := b[n, s] = Expand[If[n == 0, If[IntegerQ @ Sqrt[s], 1, 0], b[n - 1, s] + x*b[n - 1, s + n]]];
T[n_] := Function[p, Table[Coefficient[p, x, i], {i, 0, n}]][b[n, 0]];
Table[T[n], {n, 0, 10}] // Flatten (* Jean-François Alcover, Jun 03 2018, from Maple *)

T(n,n) = 1 for n in { A001108 }, T(n,n) = 0 otherwise.
T(n,n-1) = 1 for n in { A214857 }, T(n,n-1) = 0 for n in { A214858 }.
Sum_{k=0..n} k * T(n,k) = A377572(n).

A063440 Number of divisors of n-th triangular number.

Original entry on oeis.org

1, 2, 4, 4, 4, 4, 6, 9, 6, 4, 8, 8, 4, 8, 16, 8, 6, 6, 8, 16, 8, 4, 12, 18, 6, 8, 16, 8, 8, 8, 10, 20, 8, 8, 24, 12, 4, 8, 24, 12, 8, 8, 8, 24, 12, 4, 16, 24, 9, 12, 16, 8, 8, 16, 24, 24, 8, 4, 16, 16, 4, 12, 36, 24, 16, 8, 8, 16, 16, 8, 18, 18, 4, 12, 24, 16, 16, 8, 16, 40, 10, 4, 16

Offset: 1

Henry Bottomley, Jul 24 2001

a(n) = 4 iff either n is in A005383 or n/2 is in A005384.
a(n) is odd iff n is in A001108.
a(n) = 6 if either n = 18 or n = q^2 where q is in A048161 or n = 2 q^2 - 1 where q is in A106483. - Robert Israel, Oct 26 2015
From Bernard Schott, Aug 29 2020: (Start)
a(n-1) is the number of solutions in positive integers (x, y, z) to the simultaneous equations (x + y - z = n, x^2 + y^2 - z^2 = n) for n > 1. See the British Mathematical Olympiad link. In this case, one always has z > x and z > y.
For n = 12 as in the Olympiad problem, the a(11) = 8 solutions are (13,78,79), (14,45,47), (15,34,37), (18,23,29), (23,18,29), (34,15,37), (45,14,47), (78,13,79). (End)

			a(6) = 4 since 1+2+3+4+5+6 = 21 has four divisors {1,3,7,21}.

Steve Dinh, The Hard Mathematical Olympiad Problems And Their Solutions, AuthorHouse, 2011, Problem 2 of the British Mathematical Olympiad 2007, page 28.

Ray Chandler, Table of n, a(n) for n = 1..10000
British Mathematical Olympiad 2007/2008, Round 1, Problem 2.
Index to sequences related to Olympiads.

Cf. A000005, A000217.

Cf. A001108, A005383, A005384, A048161, A060778, A081978 (greedy inverse), A106483, A101755 (indices of records), A101756 (records).

seq(numtheory:-tau(n*(n+1)/2), n=1..100); # Robert Israel, Oct 26 2015

DivisorSigma[0,#]&/@Accumulate[Range[90]] (* Harvey P. Dale, Apr 15 2019 *)

for (n=1, 10000, write("b063440.txt", n, " ", numdiv(n*(n + 1)/2)) ) \\ Harry J. Smith, Aug 21 2009

a(n)=factorback(apply(numdiv,if(n%2,[n,(n+1)/2],[n/2,n+1]))) \\ Charles R Greathouse IV, Dec 27 2014

vector(100, n, numdiv(n*(n+1)/2)) \\ Altug Alkan, Oct 26 2015

a(n) = A000005(A000217(n)).
From Robert Israel, Oct 26 2015: (Start)
a(2k) = A000005(k)*A000005(2k+1).
a(2k+1) = A000005(2k+1)*A000005(k+1).
gcd(a(2k), a(2k+1)) = A000005(2k+1) * A060778(k). (End)

A055997 Numbers k such that k*(k - 1)/2 is a square.

Original entry on oeis.org

1, 2, 9, 50, 289, 1682, 9801, 57122, 332929, 1940450, 11309769, 65918162, 384199201, 2239277042, 13051463049, 76069501250, 443365544449, 2584123765442, 15061377048201, 87784138523762, 511643454094369, 2982076586042450, 17380816062160329, 101302819786919522

Offset: 1

Barry E. Williams, Jun 14 2000

Numbers k such that (k-th triangular number - k) is a square.
Gives solutions to A007913(2x)=A007913(x-1). - Benoit Cloitre, Apr 07 2002
Number of closed walks of length 2k on the grid graph P_2 X P_3. - Mitch Harris, Mar 06 2004
If x = A001109(n - 1), y = a(n) and z = x^2 + y, then x^4 + y^3 = z^2. - Bruno Berselli, Aug 24 2010
The product of any term a(n) with an even successor a(n + 2k) is always a square number. The product of any term a(n) with an odd successor a(n + 2k + 1) is always twice a square number. - Bradley Klee & Bill Gosper, Jul 22 2015
It appears that dividing even terms by two and taking the square root gives sequence A079496. - Bradley Klee, Jul 25 2015
The bisections of this sequence are a(2n - 1) = A055792(n) and a(2n) = A088920(n). - Bernard Schott, Apr 19 2020

A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, p. 193.
P. Tauvel, Exercices d'Algèbre Générale et d'Arithmétique, Dunod, 2004, Exercice 35 pages 346-347.

Colin Barker, Table of n, a(n) for n = 1..100
Dario Alpern, a^4+b^3=c^2.
Yurii S. Bystryk, Vitalii L. Denysenko, and Volodymyr I. Ostryk, Lune and Lens Sequences, ResearchGate preprint, 2024. See pp. 46, 56.
Phil Lafer, Discovering the square-triangular numbers, Fib. Quart., 9 (1971), 93-105.
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Circular Segment, Forum Geometricorum, Vol. 18 (2018), 47-55.
Kenneth Ramsey, Generalized Proof re Square Triangular Numbers, message 62 in Triangular_and_Fibonacci_Numbers Yahoo group, Oct 10, 2011.
Index entries for linear recurrences with constant coefficients, signature (7,-7,1).

Cf. A007913, A001541.
A001109(n-1) = sqrt{[(a(n))^2 - (a(n))]/2}.
a(n) = A001108(n-1)+1.
A001110(n-1)=a(n)*(a(n)-1)/2.
Cf. A001652, A001653, A046090.
Identical to A115599, but with additional leading term.

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

A:= gfun:-rectoproc({a(n) = 6*a(n-1)-a(n-2)-2, a(1) = 1, a(2) = 2}, a(n), remember):
map(A,[$1..100]); # Robert Israel, Jul 22 2015

Table[ 1/4*(2 + (3 - 2*Sqrt[2])^k + (3 + 2*Sqrt[2])^k ) // Simplify, {k, 0, 20}] (* Jean-François Alcover, Mar 06 2013 *)
CoefficientList[Series[(1 - 5 x + 2 x^2) / ((1 - x) (1 - 6 x + x^2)), {x, 0, 40}], x] (* Vincenzo Librandi, Mar 20 2015 *)
(1 + ChebyshevT[#, 3])/2 & /@ Range[0, 20] (* Bill Gosper, Jul 20 2015 *)
a[1]=1;a[2]=2;a[n_]:=(a[n-1]+1)^2/a[n-2];a/@Range[25] (* Bradley Klee, Jul 25 2015 *)
LinearRecurrence[{7,-7,1},{1,2,9},30] (* Harvey P. Dale, Dec 06 2015 *)

Vec((1-5*x+2*x^2)/((1-x)*(1-6*x+x^2))+O(x^66)) /* Joerg Arndt, Mar 06 2013 */

t(n)=(1+sqrt(2))^(n-1);
for(k=1,24,print1(round((1/4)*(t(k)^2 + t(k)^(-2) + 2)),", ")) \\ Hugo Pfoertner, Nov 29 2019

a(n) = (1 + polchebyshev(n-1, 1, 3))/2; \\ Michel Marcus, Apr 21 2020

a(n) = 6*a(n - 1) - a(n - 2) - 2; n >= 3, a(1) = 1, a(2) = 2.
G.f.: x*(1 - 5*x + 2*x^2)/((1 - x)*(1 - 6*x + x^2)).
a(n) - 1 + sqrt(2*a(n)*(a(n) - 1)) = A001652(n - 1). - Charlie Marion, Jul 21 2003; corrected by Michel Marcus, Apr 20 2020
a(n) = IF(mod(n; 2)=0; (((1 - sqrt(2))^n + (1 + sqrt(2))^n)/2)^2; 2*((((1 - sqrt(2))^(n + 1) + (1 + sqrt(2))^(n + 1)) - (((1 - sqrt(2))^n + (1 + sqrt(2))^n)))/4)^2). The odd-indexed terms are a(2n + 1) = [A001333(2n)]^2; the even-indexed terms are a(2n) = [A001333(2n - 1)]^2 + 1 = 2*[A001653(n)]^2. - Antonio Alberto Olivares, Jan 31 2004; corrected by Bernard Schott, Apr 20 2020
A053141(n + 1) + a(n + 1) = A001541(n + 1) + A001109(n + 1). - Creighton Dement, Sep 16 2004
a(n) = (1/2) + (1/4)*(3+2*sqrt(2))^(n-1) + (1/4)*(3-2*sqrt(2))^(n-1). - Antonio Alberto Olivares, Feb 21 2006; corrected by Michel Marcus, Apr 20 2020
a(n) = A001653(n)-A001652(n-1). - Charlie Marion, Apr 10 2006; corrected by Michel Marcus, Apr 20 2020
a(2k) = A001541(k)^2. - Alexander Adamchuk, Nov 24 2006
a(n) = 2*A001653(m)*A011900(n-m-1) +A002315(m)*A001652(n-m-1) - A001108(m) with mA001653(m)*A011900(m-n) - A002315(m)*A046090(m-n) - A001108(m). See Link to Generalized Proof re Square Triangular Numbers. - Kenneth J Ramsey, Oct 13 2011
a(n) = +7*a(n-1) -7*a(n-2) +1*a(n-3). - Joerg Arndt, Mar 06 2013
a(n) * a(n+2) = (A001108(n)-A001652(n)+3*A046090(n))^2. - Robert Israel, Jul 23 2015
sqrt(a(n+1)*a(n-1)) = a(n)+1 - Bradley Klee & Bill Gosper, Jul 25 2015
a(n) = 1 + sum{k=0..n-2} A002315(k). - David Pasino, Jul 09 2016; corrected by Michel Marcus, Apr 20 2020
E.g.f.: (2*exp(x) + exp((3-2*sqrt(2))*x) + exp((3+2*sqrt(2))*x))/4. - Ilya Gutkovskiy, Jul 09 2016
sqrt(a(n)*(a(n)-1)/2) = A001542(n)/2. - David Pasino, Jul 09 2016
Limit_{n -> infinity} a(n)/a(n-1) = A156035. - César Aguilera, Apr 07 2018
a(n) = (1/4)*(t^2 + t^(-2) + 2), where t = (1+sqrt(2))^(n-1). - Ridouane Oudra, Nov 29 2019
sqrt(a(n)) + sqrt(a(n) - 1) = (1 + sqrt(2))^(n - 1). - Ridouane Oudra, Nov 29 2019
sqrt(a(n)) - sqrt(a(n) - 1) = (-1 + sqrt(2))^(n - 1). - Bernard Schott, Apr 18 2020

A072221 a(n) = 6*a(n-1) - a(n-2) + 2, with a(0)=1, a(1)=4.

Original entry on oeis.org

1, 4, 25, 148, 865, 5044, 29401, 171364, 998785, 5821348, 33929305, 197754484, 1152597601, 6717831124, 39154389145, 228208503748, 1330096633345, 7752371296324, 45184131144601, 263352415571284, 1534930362283105, 8946229758127348, 52142448186480985

Offset: 0

Lekraj Beedassy, Jul 04 2002

The product of three consecutive triangular numbers with middle term A000217(m) where m is in this sequence is a square.
k is in this sequence iff the triangle with sides 3,k,k+1 has integer area. Equivalently, numbers k such that 2*(k+2)*(k-1) is a square. - James R. Buddenhagen, Oct 19 2008
Triangular numbers that are equal to a square plus one have this sequence as indices. For example, the 25th triangular number is 25*26/2 = 325 = 18^2 + 1. - Tanya Khovanova and Alexey Radul, Aug 08 2009
The triangle with sides 3, a(n), a(n)+1 has area A075848(n) if n>=0. - Michael Somos, Dec 25 2018
Compare with A016064 for integers m with triangles with sides 4, m, m+2 and integer area. - Michael Somos, May 11 2019

			From _Michael Somos_, Dec 25 2018: (Start)
For n=1, the triangle (3, 4, 5) has area 6 = A075848(1).
For n=2, the triangle (3, 25, 26) has area 36 = A075848(2). (End)

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Sylvester Robins, Certain Series of Integral, Rational, Scalene Triangles, The American Mathematical Monthly, Vol. 1, No. 1 (Jan., 1894), pp. 13-14. See Example I.
Index entries for linear recurrences with constant coefficients, signature (7,-7,1).

Cf. A000217, A001108, A001541, A016064, A075848.

a072221 n = a072221_list !! n
a072221_list = 1 : 4 : (map (+ 2) $
   zipWith (-) (map (* 6) $ tail a072221_list) a072221_list)
-- Reinhard Zumkeller, Apr 27 2012

a[n_] := a[n] = 6a[n - 1] - a[n - 2] + 2; a[0] = 1; a[1] = 4; Table[ a[n], {n, 0, 20}]
LinearRecurrence[{7, -7, 1}, {1, 4, 25}, 25] (* T. D. Noe, Dec 09 2013 *)
a[ n_] := (3 ChebyshevT[ n, 3] - 1) / 2; (* Michael Somos, Dec 25 2018 *)

{a(n) = (3 * polchebyshev( n, 1, 3) - 1) / 2}; /* Michael Somos, Dec 25 2018 */

a(n) = (3*A001541(n) - 1)/2.
a(n) = 3*A001108(n) + 1. - David Scheers, Dec 25 2006
From Franz Vrabec, Aug 21 2006: (Start)
a(n) = -1/2 + (3/4)*((3+sqrt(8))^n + (3-sqrt(8))^n) for n >= 0.
a(n) = floor((3/4)*(3+sqrt(8))^n) for n > 0. (End)
G.f.: (1-3x+4x^2)/((1-x)(1-6x+x^2)). - R. J. Mathar, Sep 09 2008
a(n) = a(-n) for all n in Z. - Michael Somos, Dec 25 2018

Edited by Robert G. Wilson v, Jul 08 2002

A077444 Numbers k such that (k^2 + 4)/2 is a square.

Original entry on oeis.org

2, 14, 82, 478, 2786, 16238, 94642, 551614, 3215042, 18738638, 109216786, 636562078, 3710155682, 21624372014, 126036076402, 734592086398, 4281516441986, 24954506565518, 145445522951122, 847718631141214, 4940866263896162, 28797478952235758, 167844007449518386

Offset: 1

Gregory V. Richardson, Nov 09 2002

The equation "(k^2 + 4)/2 is a square" is a version of the generalized Pell Equation x^2 - D*y^2 = C where x^2 - 2*y^2 = -4.
Sequence of all positive integers k such that continued fraction [k,k,k,k,k,k,...] belongs to Q(sqrt(2)). - Thomas Baruchel, Sep 15 2003
Equivalently, 2*n^2 + 8 is a square.
Numbers n such that (ceiling(sqrt(n*n/2)))^2 = 2 + n^2/2. - Ctibor O. Zizka, Nov 09 2009
The continued fraction [a(n);a(n),a(n),...] = (1 + sqrt(2))^(2*n-1). - Thomas Ordowski, Jun 07 2013
a((p+1)/2) == 2 (mod p) where p is an odd prime. - Altug Alkan, Mar 17 2016

A. H. Beiler, "The Pellian." Ch. 22 in Recreations in the Theory of Numbers: The Queen of Mathematics Entertains. Dover, New York, New York, pp. 248-268, 1966.
L. E. Dickson, History of the Theory of Numbers, Vol. II, Diophantine Analysis. AMS Chelsea Publishing, Providence, Rhode Island, 1999, pp. 341-400.
Peter G. L. Dirichlet, Lectures on Number Theory (History of Mathematics Source Series, V. 16); American Mathematical Society, Providence, Rhode Island, 1999, pp. 139-147.

Amiram Eldar, Table of n, a(n) for n = 1..1306
Sergio Falcon, Relationships between Some k-Fibonacci Sequences, Applied Mathematics, 2014, 5, 2226-2234.
Tanya Khovanova, Recursive Sequences
Prabha Sivaraman Nair and Rejikumar Karunakaran, On k-Fibonacci Brousseau Sums, J. Int. Seq. (2024) Art. No. 24.6.4. See p. 3.
J. J. O'Connor and E. F. Robertson, Pell's Equation. [Broken link]
Eric Weisstein's World of Mathematics, Pell Equation.
Eric Weisstein's World of Mathematics, NSW Number.
Index entries for linear recurrences with constant coefficients, signature (6,-1).

(A077445(n))^2 - 2*a(n) = 8.
First differences of A001541.
Pairwise sums of A001542.
Bisection of A002203 and A080039.
Cf. A001653.

[n: n in [0..10^8] | IsSquare((n^2 + 4) div 2)]; // Vincenzo Librandi, Jun 20 2015

LinearRecurrence[{6,-1},{2,14},30] (* Harvey P. Dale, Jul 25 2018 *)

for(n=1,20,q=(1+sqrt(2))^(2*n-1);print1(contfrac(q)[1],", ")) \\ Derek Orr, Jun 18 2015

Vec(2*x*(1+x)/(1-6*x+x^2) + O(x^100)) \\ Altug Alkan, Mar 17 2016

a(n) = (((3 + 2*sqrt(2))^n - (3 - 2*sqrt(2))^n) + ((3 + 2*sqrt(2))^(n-1) - (3 - 2*sqrt(2))^(n-1))) / (2*sqrt(2)).
a(n) = 2*A002315(n-1).
Recurrence: a(n) = 6*a(n-1) - a(n-2), starting 2, 14.
Offset 0, with a=3+2*sqrt(2), b=3-2*sqrt(2): a(n) = a^((2n+1)/2) - b^((2n+1)/2). a(n) = 2*(A001109(n+1) + A001109(n)) = (A003499(n+1) - A003499(n))/2 = 2*sqrt(A001108(2n+1)) = sqrt(A003499(2n+1)-2). - Mario Catalani (mario.catalani(AT)unito.it), Mar 31 2003
Limit_{n->oo} a(n)/a(n-1) = 5.82842712474619009760... = 3 + 2*sqrt(2). See A156035.
From R. J. Mathar, Nov 16 2007: (Start)
G.f.: 2*x*(1+x)/(1-6*x+x^2).
a(n) = 2*(7*A001109(n) - A001109(n+1)). (End)
a(n) = (1+sqrt(2))^(2*n-1) - (1+sqrt(2))^(1-2*n). - Gerson Washiski Barbosa, Sep 19 2010
a(n) = floor((1 + sqrt(2))^(2*n-1)). - Thomas Ordowski, Jun 07 2013
a(n) = sqrt(2*A075870(n)^2-4). - Derek Orr, Jun 18 2015
a(n) = 2*sqrt((2*A001653(n)^2)-1). - César Aguilera, Jul 13 2023
E.g.f.: 2*(1 + exp(3*x)*(sqrt(2)*sinh(2*sqrt(2)*x) - cosh(2*sqrt(2)*x))). - Stefano Spezia, Aug 27 2024

A143608 A005319 and A002315 interleaved.

Original entry on oeis.org

0, 1, 4, 7, 24, 41, 140, 239, 816, 1393, 4756, 8119, 27720, 47321, 161564, 275807, 941664, 1607521, 5488420, 9369319, 31988856, 54608393, 186444716, 318281039, 1086679440, 1855077841, 6333631924, 10812186007, 36915112104, 63018038201, 215157040700

Offset: 0

Originally submitted by Clark Kimberling, Aug 27 2008. Merged with an essentially identical sequence submitted by Kenneth J Ramsey, Jun 01 2012, by N. J. A. Sloane, Aug 02 2012

Also, numerators of the lower principal and intermediate convergents to 2^(1/2). The lower principal and intermediate convergents to 2^(1/2), beginning with 1/1, 4/3, 7/5, 24/17, 41/29, form a strictly increasing sequence; essentially, numerators=A143608 and denominators=A079496.
Sequence a(n) such that a(2*n) = sqrt(2*A001108(2*n)) and a(2*n+1) = sqrt(A001108(2*n+1)).
For n > 0, a(n) divides A******(k+1,n+1)-A******(k,n+1) where A****** is any one of A182431, A182439, A182440, A182441 and k is any nonnegative integer.
If p is a prime of the form 8*r +/- 3 then a(p+1) == 0 (mod p); if p is a prime of the form 8*r +/- 1 then a(p-1) == 0 (mod p).
Numbers n such that sqrt(floor(n^2/2 + 1)) is an integer. The integer square roots are given by A079496. - Richard R. Forberg, Aug 01 2013
From Peter Bala, Mar 23 2018: (Start)
Define a binary operation o on the real numbers by x o y = x*sqrt(1 + y^2) + y*sqrt(1 + x^2). The operation o is commutative and associative with identity 0. Then we have
a(2*n + 1) = 1 o 1 o ... o 1 (2*n + 1 terms) and
a(2*n) = sqrt(2)*(1 o 1 o ... o 1) (2*n terms). Cf. A084068.
This is a fourth-order divisibility sequence. Indeed, a(2*n) = sqrt(2)*U(2*n) and a(2*n+1) = U(2*n+1), where U(n) is the Lehmer sequence [Lehmer, 1930] defined by the recurrence U(n) = 2*sqrt(2)*U(n-1) - U(n-2) with U(0) = 0 and U(1) = 1. The solution to the recurrence is U(n) = (1/2)*( (sqrt(2) + 1)^n - (sqrt(2) - 1)^n ). (End)

Serge Lang, Introduction to Diophantine Approximations, Addison-Wesley, New York, 1966.

Colin Barker, Table of n, a(n) for n = 0..1000
John M. Campbell, An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences, arXiv preprint arXiv:1105.3399 [math.GM], 2011.
Creighton Kenneth Dement, Comments on A143608 and A143609
Vincent Granville, Successive records in mathematical sequences: surprising result, Mathematics Stack Exchange, 2019.
Clark Kimberling, Best lower and upper approximations to irrational numbers, Elem. Math. vol. 52 iss. 3 (1997) 122-126.
D. H. Lehmer, An extended theory of Lucas' functions, Annals of Mathematics, Second Series, Vol. 31, No. 3 (Jul., 1930), pp. 419-448.
Eric Weisstein's World of Mathematics, Lehmer Number
Index entries for linear recurrences with constant coefficients, signature (0,6,0,-1).

Cf. A002315, A005319, A182431, A001108, A049629, A084068.

I:=[0,1,4,7]; [n le 4 select I[n] else 6*Self(n-2) - Self(n-4): n in [1..30]]; // G. C. Greubel, Mar 27 2018

A143608 := proc(n)
    option remember;
    if n <= 3 then
        op(n+1,[0,1,4,7]) ;
    else
        6*procname(n-2)-procname(n-4) ;
    end if;
end proc: # R. J. Mathar, Jul 22 2012

a = -4; b = -1; Reap[While[b<2000000000, t = 4*b-a; Sow[t]; a=b; b=t; t = 2*b-a; Sow[t]; a=b; b=t]][[2,1]]
CoefficientList[Series[x*(1 + 4*x + x^2)/(1 - 6*x^2 + x^4), {x, 0, 30}], x] (* Wesley Ivan Hurt, Aug 24 2014 *)
LinearRecurrence[{0, 6, 0, -1}, {0, 1, 4, 7}, 31] (* Jean-François Alcover, Sep 21 2017 *)

a(n)=([0,1,0,0;0,0,1,0;0,0,0,1;-1,0,6,0]^n*[0;1;4;7])[1,1] \\ Charles R Greathouse IV, Jun 11 2015

concat(0, Vec(x*(1+4*x+x^2)/((1+2*x-x^2)*(1-2*x-x^2)) + O(x^50))) \\ Colin Barker, Mar 27 2016

a(2*n) = (a(2*n - 1) + a(2*n + 1))/2.
a(2*n + 1) = (a(2*n) + a(2*n + 2))/4.
a(2*n) = 4*A001109(n).
a(2*n + 1) = 4*A001109(n) + A001541(n).
From Colin Barker, Jun 29 2012: (Start)
a(n) = 6*a(n-2) - a(n-4).
G.f.: x*(1 + 4*x + x^2)/((1 + 2*x - x^2)*(1 - 2*x - x^2)) = x*(1 + 4*x + x^2)/(1 - 6*x^2 + x^4). (End)
2*a(n) = A078057(n) - A123335(n-1). - R. J. Mathar, Jul 04 2012
a(2n) = A005319(n); a(2n+1) = A002315(n). - R. J. Mathar, Jul 17 2009
a(n)*a(n+1) + 1 = A001653(n+1). - Charlie Marion, Dec 11 2012
a(n) = (((-2 - sqrt(2) + (-1)^n * (-2+sqrt(2))) * ((-1+sqrt(2))^n - (1+sqrt(2))^n)))/(4*sqrt(2)). - Colin Barker, Mar 27 2016
a(n) = A084068(n) - A079496(n). - César Aguilera, Feb 14 2023

A012132 Numbers z such that x(x+1) + y(y+1) = z*(z+1) is solvable in positive integers x,y.

Original entry on oeis.org

3, 6, 8, 10, 11, 13, 15, 16, 18, 20, 21, 23, 26, 27, 28, 31, 33, 36, 37, 38, 40, 41, 43, 44, 45, 46, 48, 49, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 62, 63, 64, 66, 67, 68, 71, 73, 74, 75, 76, 77, 78, 80, 81, 83, 86, 88, 89, 91, 92, 93

Offset: 1

Sander van Rijnswou (sander(AT)win.tue.nl)

nonn

Theorem (Sierpinski, 1963): n is a term iff n^2+(n+1)^2 is a composite number. - N. J. A. Sloane, Feb 29 2020
For n > 1, A047219 is a subset of this sequence. This is because n^2 + (n+1)^2 is divisible by 5 if n is (1 or 3) mod 5 (also see A027861). - Dmitry Kamenetsky, Sep 02 2008
From Hermann Stamm-Wilbrandt, Sep 10 2014: (Start)
For n > 0, A212160 is a subset of this sequence (n^2 + (n+1)^2 is divisible by 13 if n == (2 or 10) (mod 13)).
For n >= 0, A212161 is a subset of this sequence (n^2 + (n+1)^2 is divisible by 17 if n == (6 or 10) (mod 17)).
The above are for divisibility by 5, 13, 17; notation (1,3,5), (2,10,13), (6,10,17). Divisibility by p for a and p-a-1; notation (a,p-a-1,p). These are the next tuples: (8,20,29), (15,21,37), (4,36,41), (11,41,53), ... . The corresponding sequences are a subset of this sequence (8,20,37,49,66,78,... for (8,20,29)). These sequences have no entries in the OEIS yet. For any prime of the form 4*k+1 there is exactly one of these tuples/sequences.
For n > 1, A000217 (triangular numbers) is a subset of this sequence (3,6,10,15,...); z=A000217(n), y=z-1, x=n.
For n > 0, A001652 is a subset of this sequence; z=A001652(n), x=y=A053141(n).
For n > 1, A001108(=A115598) is a subset of this sequence; z=A001108(n), x=A076708(n), y=x+1.
For n > 0, A124124(2*n+1)(=A098790(2*n)) is a subset of this sequence (6,37,218,...); z=A124124(2*n+1), x=a(n)-1, y=a(n)+1, a(m) = 6*a(m-1) - a(m-2) + 2, a(0)=0, a(1)=4.
(End)

Aviezri S. Fraenkel, Diophantine equations involving generalized triangular and tetrahedral numbers, pp. 99-114 of A. O. L. Atkin and B. J. Birch, editors, Computers in Number Theory. Academic Press, NY, 1971.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
H. Finner and K. Strassburger, Increasing sample sizes do not necessarily increase the power of UMPU-tests for 2 X 2-tables, Metrika, 54, 77-91, (2001).
Heiko Harborth, Fermat-like binomial equations, Applications of Fibonacci numbers, Proc. 2nd Int. Conf., San Jose/Ca., August 1986, 1-5 (1988).
W. Sierpinski, On triangular numbers which are sums of two smaller triangular numbers, (Polish), Wiadom. Mat. (2) 7 (1963), 27-28. See MR0182602.

Complement of A027861. - Michael Somos, Jun 08 2000
Cf. A000217, A047219, A027861, A166080.
Cf. also A212160, A212161, A001652, A001108, A115598, A124124, A098790.

Select[Range[100], !PrimeQ[#^2 + (#+1)^2]& ] (* Jean-François Alcover, Jan 17 2013, after Michael Somos *)

More terms and references from Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Feb 09 2000

cp's OEIS Frontend

A001542 a(n) = 6*a(n-1) - a(n-2) for n > 1, a(0)=0 and a(1)=2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

Maxima

PARI

PARI

Python

Sage

Formula

A053141 a(0)=0, a(1)=2 then a(n) = a(n-2) + 2*sqrt(8*a(n-1)^2 + 8*a(n-1) + 1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

Extensions

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

A281871 Number T(n,k) of k-element subsets of [n] having a square element sum; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A063440 Number of divisors of n-th triangular number.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

A053141 a(0)=0, a(1)=2 then a(n) = a(n-2) + 2sqrt(8a(n-1)^2 + 8*a(n-1) + 1).

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

A012132 Numbers z such that x(x+1) + y(y+1) = z*(z+1) is solvable in positive integers x,y.