A046090 -id:A046090 - cp's OEIS frontend

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

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

A038761 a(n) = 6*a(n-1) - a(n-2), n >= 2, a(0)=1, a(1)=9.

Original entry on oeis.org

1, 9, 53, 309, 1801, 10497, 61181, 356589, 2078353, 12113529, 70602821, 411503397, 2398417561, 13979001969, 81475594253, 474874563549, 2767771787041, 16131756158697, 94022765165141, 548004834832149, 3194006243827753, 18616032628134369, 108502189524978461

Offset: 0

Barry E. Williams, May 02 2000

Bisection of A048654. - Lambert Klasen (lambert.klasen(AT)gmx.de), Nov 24 2004

This gives part of the (increasingly sorted) positive solutions y to the Pell equation x^2 - 2*y^2 = +7. For the x solutions see A038762. For the other part of solutions see A101386 and A253811. - Wolfdieter Lang, Feb 05 2015

			A038762(3)^2 - 2*a(4)^2 = 2547^2 - 2*1801^2 = +7. - _Wolfdieter Lang_, Feb 05 2015

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

Vincenzo Librandi, Table of n, a(n) for n = 0..400
M. J. DeLeon, Pell's Equation and Pell Number Triples, Fib. Quart., 14(1976), pp. 456-460.
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (6,-1).
Index entries for sequences related to Chebyshev polynomials.

Cf. A001109, A001541, A001542, A001652, A001653, A038762, A046090, A048654, A055997, A101386, A124124, A253811.

I:=[1, 9]; [n le 2 select I[n] else 6*Self(n-1)-Self(n-2): n in [1..40]]; // Vincenzo Librandi, Nov 16 2011

a[0]:=1: a[1]:=9: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..19); # Zerinvary Lajos, Jul 26 2006

LinearRecurrence[{6,-1},{1,9},40] (* Vincenzo Librandi, Nov 16 2011 *)

a(n)=([0,1; -1,6]^n*[1;9])[1,1] \\ Charles R Greathouse IV, Apr 08 2016

a(n) = (9*((3+2*sqrt(2))^n -(3-2*sqrt(2))^n)-((3+2*sqrt(2))^(n-1) - (3-2*sqrt(2))^(n-1)))/(4*sqrt(2)).
a(n) = sqrt(2*(A038762(n))^2-14)/2.
For n>1, a(n)-4a(n-1)=A001541(n)-A001542(n-2); e.g. 309-4*53=97=99-2. - Charlie Marion, Nov 12 2003
For n>0, a(n)=A046090(n)+A001653(n)+A001652(n-1)=A055997(n+1)+A001652(n-1); e.g., 309=120+169+20. - Charlie Marion, Oct 11 2006
G.f.: (1+3*x)/(1-6*x+x^2). - Philippe Deléham, Nov 03 2008
a(n) = third binomial transform of 1,6,8,48,64,384. - Al Hakanson (hawkuu(AT)gmail.com), Aug 15 2009
a(n)^2 + 2^2 = A124124(2*n+1)^2 + (A124124(2*n+1)+1)^2. - Hermann Stamm-Wilbrandt, Aug 31 2014
a(n) = irrational part of z(n) = (3 + sqrt(2))*(3 + 2*sqrt(2))^n, n >= 0. z(n) gives only part of the general positive solutions to the Pell equation x^2 - 2*y^2 = 7. See the Nagell reference in A038762 on how to find z(n), and a comment above. - Wolfdieter Lang, Feb 05 2015
a(n) = S(n, 6) + 3*S(n-1, 6), n >= 0, with the Chebyshev S-polynomials evaluated at x=6. See S(n-1, 6) = A001109(n). - Wolfdieter Lang, Mar 30 2015
E.g.f.: exp(3*x)*(2*cosh(2*sqrt(2)*x) + 3*sqrt(2)*sinh(2*sqrt(2)*x))/2. - Stefano Spezia, Mar 16 2024

Edited: Replaced the unspecific Pell comment. Moved a formula from the comment section to the formula section. - Wolfdieter Lang, Feb 05 2015

A075528 Triangular numbers that are half other triangular numbers.

Original entry on oeis.org

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

Offset: 0

Christian G. Bower, Sep 19 2002

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

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

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

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

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

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

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

Original entry on oeis.org

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

Offset: 0

N. J. A. Sloane

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

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

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

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

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

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

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

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

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

A008844 Squares of sequence A001653: y^2 such that x^2 - 2*y^2 = -1 for some x.

Original entry on oeis.org

1, 25, 841, 28561, 970225, 32959081, 1119638521, 38034750625, 1292061882721, 43892069261881, 1491038293021225, 50651409893459761, 1720656898084610641, 58451683124983302025, 1985636569351347658201, 67453191674820837076801, 2291422880374557112953025

Offset: 0

N. J. A. Sloane

Numbers simultaneously square and centered square. E.g., a(1)=25 because 25 is the fourth centered square number and the fifth square number. - Steven Schlicker, Apr 24 2007
Solutions to A007913(x)=A007913(2x-1). - Benoit Cloitre, Apr 07 2002
From Ant King, Nov 09 2011: (Start)
Indices of positive hexagonal numbers that are also perfect squares.
As n increases, this sequence is approximately geometric with common ratio r = lim_{n -> infinity} a(n)/a(n-1) = (1 + sqrt(2))^4 = 17 + 12 * sqrt(2).
(End)
Also indices of hexagonal numbers (A000384) which are also centered octagonal numbers (A016754). - Colin Barker, Jan 25 2015
Also positive integers x in the solutions to 4*x^2 - 8*y^2 - 2*x + 8*y - 2 = 0, the corresponding values of y being A253826. - Colin Barker, Jan 25 2015
Squares that are sum of two consecutive squares: y^2 = (k + 1)^2 + k^2 is equivalent to x^2 - 2*y^2 = -1 with x = 2*k + 1. - Jean-Christophe Hervé, Nov 11 2015
Squares in the main diagonal of the natural number array, A000027. - Clark Kimberling, Mar 12 2023

			From _Ravi Kumar Davala_, May 26 2013: (Start)
A001333(0)=1, A001333(4)=17, A001333(8)=577, A000129(0)=0, A000129(2)=2, A000129(4)=12, A000129(8)=408 so clearly
a(n+m)=A001333(4*m)*a(n)-(A000129(2*m))^2+A000129(4*m)*sqrt(2*a(n)^2-a(n)), with m=1,2 is true.
A002203(0)=2, A002203(4)=34, A002203(8)=1154 so clearly
a(n+m)=(1/2)*A002203(4*m)*a(n)-(A000129(2*m))^2+A000129(4*m)*sqrt(2*a(n)^2-a(n)) is true for m=1,2
a(n+1)*a(n-1) = (a(n)+4)^2 , with n=1 is 841*1=(25+4)^2, for n=2 , 28561*25=(841+4)^2.
(End)
1 = 1 + 0, 25 = 16 + 9, 841 = 29^2 = 21^2 + 20^2 = 441 + 400.

Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 256.

Muniru A Asiru, Table of n, a(n) for n = 0..100
Daniel C. Fielder, Special integer sequences controlled by three parameters, Fibonacci Quarterly 6, 1968, 64-70.
M. A. Gruber, Artemas Martin, A. H. Bell, J. H. Drummond, A. H. Holmes, and H. C. Wilkes, Problem 47, Amer. Math. Monthly, 4 (1897), 25-28.
Thomas Koshy, Products Involving Reciprocals of Gibonacci Polynomials, The Fibonacci Quarterly, Vol. 60, No. 1 (2022), pp. 15-24.
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Circular Segment, Forum Geometricorum, Vol. 18 (2018), 47-55.
Steven C. Schlicker, Numbers Simultaneously Polygonal and Centered Polygonal, Mathematics Magazine, Vol. 84, No. 5, December 2011
Eric Weisstein's World of Mathematics, Hexagonal Square Number.
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).

Cf. A000290, A001844, A007913, A000384, A046177, A016754, A253826.

a := [1, 25, 841];; for i in [4..10^2] do a[i] := 35*a[i-1] - 35*a[i-2] + a[i-3]; od; a;  # Muniru A Asiru, Jan 17 2018

I:=[1,25,841]; [n le 3 select I[n] else 35*Self(n-1)-35*Self(n-2)+Self(n-3): n in [1..20]]; // Vincenzo Librandi, Jan 20 2018

CP := n -> 1+1/2*4*(n^2-n): N:=10: u:=3: v:=1: x:=4: y:=1: k_pcp:=[1]: for i from 1 to N do tempx:=x; tempy:=y; x:=tempx*u+8*tempy*v: y:=tempx*v+tempy*u: s:=(y+1)/2: k_pcp:=[op(k_pcp),CP(s)]: end do: k_pcp; # Steven Schlicker, Apr 24 2007

LinearRecurrence[{35, -35, 1}, {1, 25, 841}, 15] (* Ant King, Nov 09 2011 *)
CoefficientList[Series[(1 - 10 x + x^2) / ((1 - x) (1 - 34 x + x^2)), {x, 0, 33}], x] (* Vincenzo Librandi, Jan 20 2018 *)

a(n)=if(n<0,0,sqr(subst(poltchebi(n+1)+poltchebi(n),x,3)/4))

vector(40, n, n--; (([5, 2; 2, 1]^n)[1, 1])^2) \\ Altug Alkan, Nov 11 2015

From Benoit Cloitre, Jan 19 2003: (Start)
a(n) = A078522(n) + 1.
a(n) = ceiling(A*B^n) where A = (3 + 2*sqrt(2))/8 and B = 17 + 12*sqrt(2). (End)
G.f.: (1-10x+x^2)/((1-x)(1-34x+x^2)).
a(n) = ceiling(A046176(n)/sqrt(2)). - Helge Robitzsch (hrobi(AT)math.uni-goettingen.de), Jul 28 2000
a(n+1) = 17*a(n) - 4 + 12*sqrt(2*a(n)^2 - a(n)). - Richard Choulet, Sep 14 2007
Define x(n) + y(n)*sqrt(8) = (4+sqrt(8))*(3+sqrt(8))^n, s(n) = (y(n)+1)/2; then a(n) = (1/2)*(2+4*(s(n)^2 - s(n))). - Steven Schlicker, Apr 24 2007
From Ant King, Nov 09 2011: (Start)
a(n) = 35*a(n-1) - 35*a(n-2) + a(n-3).
a(n) = 34*a(n-1) - a(n-2) - 8.
a(n) = 1/8 * ((1 + sqrt(2))^(4*n-2) + (1 - sqrt(2))^(4*n-2) + 2).
a(n) = ceiling((1/8) * (1 + sqrt(2))^(4*n-2)). (End)
From Ravi Kumar Davala, May 26 2013: (Start)
a(n+2) = 577*a(n) - 144 + 408*sqrt(2*a(n)^2 - a(n)).
a(n+m) = A001333(4*m)*a(n) - (A000129(2*m))^2 + A000129(4*m)*sqrt(2*a(n)^2 - a(n)).
a(n+m) = (1/2)*A002203(4*m)*a(n) - (A000129(2*m))^2 + A000129(4*m)*sqrt(2*a(n)^2 - a(n)).
a(n+1)*a(n-1) = (a(n)+4)^2. (End)
a(n) = A001652(n)^2 + A046090(n)^2. - César Aguilera, Jan 15 2018
Limit_{n -> infinity} a(n)/a(n-1) = A156164. - César Aguilera, Jan 28 2018
sqrt(2*a(n))-1 = A002315(n). - Ezhilarasu Velayutham, Apr 05 2019
4*a(n) = 1 +3*A077420(n). - R. J. Mathar, Mar 05 2024
Product_{n>=0} (1 + 4/a(n)) = 2*sqrt(2) + 3 (Koshy, 2022, section 3, p. 19). - Amiram Eldar, Jan 23 2025

Entry edited by N. J. A. Sloane, Sep 14 2007

A011900 a(n) = 6*a(n-1) - a(n-2) - 2 with a(0) = 1, a(1) = 3.

Original entry on oeis.org

1, 3, 15, 85, 493, 2871, 16731, 97513, 568345, 3312555, 19306983, 112529341, 655869061, 3822685023, 22280241075, 129858761425, 756872327473, 4411375203411, 25711378892991, 149856898154533, 873430010034205, 5090723162050695, 29670908962269963, 172934730611569081

Offset: 0

Mario Velucchi (mathchess(AT)velucchi.it)

Members of Diophantine pairs.
Solution to b*(b-1) = 2*a*(a-1) in natural numbers; a = a(n), b = b(n) = A046090(n).
Also the indices of centered octagonal numbers which are also centered square numbers. - Colin Barker, Jan 01 2015
Also positive integers y in the solutions to 4*x^2 - 8*y^2 - 4*x + 8*y = 0. - Colin Barker, Jan 01 2015
Also the number of perfect matchings on a triangular lattice of width 3 and length n. - Sergey Perepechko, Jul 11 2019

			G.f. = 1 + 3*x + 15x^2 + 85*x^3 + 493*x^4 + 2871*x^5 + 16731*x^6 + ... - _Michael Somos_, Feb 23 2019

Mario Velucchi "The Pell's equation ... an amusing application" in Mathematics and Informatics Quarterly, to appear 1997.

Colin Barker, Table of n, a(n) for n = 0..1000
H. J. Hindin, Stars, hexes, triangular numbers and Pythagorean triples, J. Rec. Math., 16 (1983/1984), 191-193. (Annotated scanned copy)
Giovanni Lucca, Circle Chains Inscribed in Symmetrical Lenses and Integer Sequences, Forum Geometricorum, Volume 16 (2016) 419-427.
S. Northshield, An Analogue of Stern's Sequence for Z[sqrt(2)], Journal of Integer Sequences, 18 (2015), #15.11.6.
S. N. Perepechko, Number of perfect matchings on triangular lattices of fixed width, DIMA'2015 slides. [see: page 12]
Index entries for linear recurrences with constant coefficients, signature (7,-7,1).

Cf. A000194, A001109, A001541, A001653, A002024, A006062, A011906, A029549, A046090, A053141, A075528, A156035.

I:=[1,3]; [n le 2 select I[n] else 6*Self(n-1) - Self(n-2) - 2: n in [1..40]]; // Vincenzo Librandi, Dec 05 2015

f:= gfun:-rectoproc({a(n)=6*a(n-1)-a(n-2)-2,a(0)=1,a(1)=3},a(n),remember):
seq(f(n),n=0..40); # Robert Israel, Dec 16 2015

a[0] = 1; a[1] = 3; a[n_]:= a[n]= 6 a[n-1] -a[n-2] -2; Table[a@ n, {n,0,40}] (* Michael De Vlieger, Dec 05 2015 *)
Table[(Fibonacci[2n + 1, 2] + 1)/2, {n, 0, 40}] (* Vladimir Reshetnikov, Sep 16 2016 *)
LinearRecurrence[{7,-7,1},{1,3,15},40] (* Harvey P. Dale, Feb 16 2017 *)
a[ n_] := (4 + ChebyshevT[n, 3] + ChebyshevT[n + 1, 3])/8; (* Michael Somos, Feb 23 2019 *)

Vec((1-4*x+x^2)/((1-x)*(1-6*x+x^2)) + O(x^50)) \\ Altug Alkan, Dec 06 2015

[(1+lucas_number1(2*n+1,2,-1))//2 for n in range(41)] # G. C. Greubel, Oct 17 2024

a(n) = (A001653(n+1) + 1)/2.
a(n) = (((1+sqrt(2))^(2*n-1) - (1-sqrt(2))^(2*n-1))/sqrt(8) + 1)/2.
a(n) = 7*(a(n-1) - a(n-2)) + a(n-3); a(1) = 1, a(2) = 3, a(3) = 15. Also 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. - Antonio Alberto Olivares, Dec 23 2003
Sqrt(2) = Sum_{n>=0} 1/a(n); a(n) = a(n-1) + floor(1/(sqrt(2) - Sum_{k=0..n-1} 1/a(k))) (n>0) with a(0)=1. - Paul D. Hanna, Jan 25 2004
For n>k, a(n+k) = A001541(n)*A001653(k) - A053141(n-k-1); e.g., 493 = 99*5 - 2. For n<=k, a(n+k)=A001541(n)*A001653(k) - A053141(k-n); e.g., 85 = 3*29 - 2. - Charlie Marion, Oct 18 2004
a(n+1) = 3*a(n) - 1 + sqrt(8*a(n)^2 - 8*a(n) + 1), a(1)=1. - Richard Choulet, Sep 18 2007
a(n+1) = a(n) * (a(n) + 2) / a(n-1) for n>=1 with a(0)=1 and a(1)=3. - Paul D. Hanna, Apr 08 2012
G.f.: (1 - 4*x + x^2)/((1-x)*(1 - 6*x + x^2)). - R. J. Mathar, Oct 26 2009
Sum_{k=a(n)..A001109(n+1)} k = a(n)*A001109(n+1) = A011906(n+1). Example n=2, 3+4+5+6=18, 3*6=18. - Paul Cleary, Dec 05 2015
a(n) = (sqrt(1+8*A001109(n+1)^2)+1)/2 - A001109(n+1). - Robert Israel, Dec 16 2015
a(n) = a(-1-n) for all n in Z. - Michael Somos, Feb 23 2019
E.g.f.: (2*exp(x) + exp(3*x)*(2*cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x)))/4. - Stefano Spezia, Mar 16 2024
a(n) = A053141(n) + 1 = A000194(A029549(n))+1 = A002024(A075528(n))+1. - Pontus von Brömssen, Sep 11 2024

More terms and comments from Wolfdieter Lang

A053142 a(n) = A053141(n)/2.

Original entry on oeis.org

0, 1, 7, 42, 246, 1435, 8365, 48756, 284172, 1656277, 9653491, 56264670, 327934530, 1911342511, 11140120537, 64929380712, 378436163736, 2205687601705, 12855689446495, 74928449077266, 436715005017102, 2545361581025347, 14835454481134981, 86467365305784540

Offset: 0

Wolfdieter Lang

Partial sums of A001109. - Barry Williams, May 03 2000.
Number m such that 16*m*(2*m+1)+1 is a square. - Bruno Berselli, Oct 19 2012
From Robert K. Moniot, Sep 21 2020: (Start)
Consecutive terms (a(n-1),a(n))=(u,v) give all points on the hyperbola u^2-u+v^2-v-6*u*v=0 in quadrant 1 with both coordinates an integer.
Let T(n) denote the n-th triangular number. If i, j are any two successive elements of the above sequence then (T(i-1) + T(j-1))/T(i+j-1) = 3/4.
(End)

G. C. Greubel, Table of n, a(n) for n = 0..1000
R. C. Alperin, A nonlinear recurrence and its relations to Chebyshev polynomials, Fib. Q., Vol. 58, No. 2 (2020), 140-142.
Ioana-Claudia Lazăr, Lucas sequences in t-uniform simplicial complexes, arXiv:1904.06555 [math.GR], 2019.
Index entries for sequences related to Chebyshev polynomials.
Index entries for linear recurrences with constant coefficients, signature (7,-7,1).

Cf. A000129, A001109, A001653, A053141, A001652, A046090, A049310.

Cf. A212336 for more sequences with g.f. of the type 1/(1-k*x+k*x^2-x^3).

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

Join[{a=0,b=1}, Table[c=6*b-a+1; a=b; b=c, {n,60}]] (* Vladimir Joseph Stephan Orlovsky, Jan 18 2011 *)
Table[(Fibonacci[2n + 1, 2] - 1)/4, {n, 0, 20}] (* Vladimir Reshetnikov, Sep 16 2016 *)
LinearRecurrence[{7, -7, 1}, {0, 1, 7}, 30] (* G. C. Greubel, Jul 15 2018 *)

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

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

a(n) = (A001653(n)-1)/4.
a(n) = 6*a(n-1)-a(n-2)+1, a(0)=0, a(1)=1.
G.f.: x/((1-x)*(1-6*x+x^2)).
From Paul Barry, Nov 14 2003: (Start)
a(n+1) = Sum_{k=0..n} S(k, 6) = Sum_{k=0..n} U(n, 3), Chebyshev polynomials of 2nd kind, A049310.
a(n+1) = (sqrt(2)-1)^(2*n)(5/8-7*sqrt(2)/16)+(sqrt(2)+1)^(2*n)*(7*sqrt(2)/16 + 5/8)-1/4. (End)
From Antonio Alberto Olivares, Jan 13 2004: (Start)
a(n) = 7*a(n-1)-7*a(n-2)+a(n-3).
a(n) = -(1/4) + (1-sqrt(2))/(-8*sqrt(2))*(3-2*sqrt(2))^n + (1+sqrt(2))/(8*sqrt(2))*(3+2*sqrt(2))^n. (End)
a(n) = Sum_{k=0..n} Sum_{j=0..2*k} (-1)^(j+1)*A000129(j)*A000129(2*k-j). Paul Barry, Oct 23 2009
a(2*k) = A001109(k)*(A001109(k) + A001109(k-1)) and a(2*k-1) = A001109(k)*(A001109(k) + A001109(k+1)). Kenneth J Ramsey, Sep 10 2010
a(n)*a(n-2) = a(n-1)*(a(n-1) - 1) for n>1. - Robert K. Moniot, Sep 21 2020
E.g.f.: (exp(3*x)*(2*cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x)) - 2*exp(x))/8. - Stefano Spezia, Mar 16 2024

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

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

Original entry on oeis.org

0, 4, 20, 120, 696, 4060, 23660, 137904, 803760, 4684660, 27304196, 159140520, 927538920, 5406093004, 31509019100, 183648021600, 1070379110496, 6238626641380, 36361380737780, 211929657785304, 1235216565974040, 7199369738058940, 41961001862379596, 244566641436218640

Offset: 0

N. J. A. Sloane

Related to Pythagorean triples: alternate terms of A001652 and A046090.
Even-valued legs of nearly isosceles right triangles: legs differ by 1. 0 is smaller leg of degenerate triangle with legs 0 and 1 and hypotenuse 1. - Charlie Marion, Nov 11 2003
The complete (nearly isosceles) primitive Pythagorean triple is given by {a(n), a(n)+(-1)^n, A001653(n)}. - Lekraj Beedassy, Feb 19 2004
Note also that A046092 is the even leg of this other class of nearly isosceles Pythagorean triangles {A005408(n), A046092(n), A001844(n)}, i.e., {2n+1, 2n(n+1), 2n(n+1)+1} where longer sides (viz. even leg and hypotenuse) are consecutive. - Lekraj Beedassy, Apr 22 2004
Union of even terms of A001652 and A046090. Sum of legs of primitive Pythagorean triangles is A002315(n) = 2*a(n) + (-1)^n. - Lekraj Beedassy, Apr 30 2004

			[1,0,1]*[1,2,2; 2,1,2; 2,2,3]^0 gives (degenerate) primitive Pythagorean triple [1, 0, 1], so a(0) = 0. [1,0,1]*[1,2,2; 2,1,2; 2,2,3]^7 gives primitive Pythagorean triple [137903, 137904, 195025] so a(7) = 137904.
G.f. = 4*x + 20*x^2 + 120*x^3 + 696*x^4 + 4060*x^5 + 23660*x^6 + ...

A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, pp. 122-125, 1964.
W. Sierpiński, Pythagorean triangles, Dover Publications, Inc., Mineola, NY, 2003, p. 17. MR2002669.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (5,5,-1).
Index entries for two-way infinite sequences

Cf. A000129, A001109, A001333, A001652, A002315, A029549, A046090.

Cf. A046727, A047235, A084158, A084159, A089499, A114336.

[4*Floor(((Sqrt(2)+1)^(2*n+1)-(Sqrt(2)-1)^(2*n+1)-2*(-1)^n) / 16): n in [0..30]]; // Vincenzo Librandi, Jul 29 2019

LinearRecurrence[{5,5,-1}, {0,4,20}, 25] (* Vincenzo Librandi, Jul 29 2019 *)

a(n)=n%2+(real((1+quadgen(8))^(2*n+1))-1)/2

a(n)=if(n<0,-a(-1-n),polcoeff(4*x/(1+x)/(1-6*x+x^2)+x*O(x^n),n))

[(lucas_number2(2*n+1,2,-1) -2*(-1)^n)/4 for n in range(41)] # G. C. Greubel, Feb 11 2023

a(n) = ((1+sqrt(2))^(2n+1) + (1-sqrt(2))^(2n+1) + 2*(-1)^(n+1))/4.
a(n) = A089499(n)*A089499(n+1).
a(n) = 4*A084158(n). - Lekraj Beedassy, Jul 16 2004
a(n) = ceiling((sqrt(2)+1)^(2*n+1) - (sqrt(2)-1)^(2*n+1) - 2*(-1)^n)/4. - Lambert Klasen (Lambert.Klasen(AT)gmx.net), Nov 12 2004
a(n) is the k-th entry among the complete near-isosceles primitive Pythagorean triple A114336(n), where k = (3*(2n-1) - (-1)^n)/2, i.e., a(n) = A114336(A047235(n)), for positive n. - Lekraj Beedassy, Jun 04 2006
a(n) = A046727(n) - (-1)^n = 2*A114620(n). - Lekraj Beedassy, Aug 14 2006
From George F. Johnson, Aug 29 2012: (Start)
2*a(n)*(a(n) + (-1)^n) + 1 = (A000129(2*n+1))^2;
n > 0, 2*a(n)*(a(n) + (-1)^n) + 1 = ((a(n+1) - a(n-1))/4)^2, a perfect square.
a(n+1) = (3*a(n) + 2*(-1)^n) + 2*sqrt(2*a(n)*(a(n) + (-1)^n)+ 1).
a(n-1) = (3*a(n) + 2*(-1)^n) - 2*sqrt(2*a(n)*(a(n) + (-1)^n)+ 1).
a(n+1) = 6*a(n) - a(n-1) + 4*(-1)^n.
a(n+1) = 5*a(n) + 5*a(n-1) - a(n-2).
a(n+1) *a(n-1) = a(n)*(a(n) + 4*(-1)^n).
a(n) = (sqrt(1 + 8*A029549(n)) - (-1)^n)/2.
a(n) = A002315(n) - A084159(n) = A084159(n) - (-1)^n.
a(n)  = A001652(n) + (1 - (-1)^n)/2 = A046090(n) - (1 + (-1)^n)/2.
Limit_{n->oo} a(n)/a(n-1) =  3 + 2*sqrt(2).
Limit_{n->oo} a(n)/a(n-2) = 17 + 12*sqrt(2).
Limit_{n->oo} a(n)/a(n-r) = (3 + 2*sqrt(2))^r.
Limit_{n->oo} a(n-r)/a(n) = (3 - 2*sqrt(2))^r. (End)
From G. C. Greubel, Feb 11 2023: (Start)
a(n) = (A001333(2*n+1) - 2*(-1)^n)/4.
a(n) = (1/2)*(A001109(n+1) + A001109(n) - (-1)^n). (End)
E.g.f.: exp(-x)*(exp(4*x)*(cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x)) - 1)/2. - Stefano Spezia, Aug 03 2024

cp's OEIS Frontend

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

A038761 a(n) = 6*a(n-1) - a(n-2), n >= 2, a(0)=1, a(1)=9.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A075528 Triangular numbers that are half other triangular numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A008844 Squares of sequence A001653: y^2 such that x^2 - 2*y^2 = -1 for some x.

Views

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.

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