cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-20 of 38 results. Next

A084703 Squares k such that 2*k+1 is also a square.

Original entry on oeis.org

0, 4, 144, 4900, 166464, 5654884, 192099600, 6525731524, 221682772224, 7530688524100, 255821727047184, 8690408031080164, 295218051329678400, 10028723337177985444, 340681375412721826704, 11573138040695364122500, 393146012008229658338304, 13355391270239113019379844
Offset: 0

Views

Author

Amarnath Murthy, Jun 08 2003

Keywords

Comments

With the exception of 0, a subsequence of A075114. - R. J. Mathar, Dec 15 2008
Consequently, A014105(k) is a square if and only if k = a(n). - Bruno Berselli, Oct 14 2011
From M. F. Hasler, Jan 17 2012: (Start)
Bisection of A079291. The squares 2*k+1 are given in A055792.
A204576 is this sequence written in binary. (End)
a(n+1), n >= 0, is the perimeter squared (x(n) + y(n) + z(n))^2 of the ordered primitive Pythagorean triple (x(n), y(n) = x(n) + 1, z(n)). The first two terms are (x(0)=0, y(0)=1, z(0)=1), a(1) = 2^2, and (x(1)=3, y(1)=4, z(1)=5), a(2) = 12^2. - George F. Johnson, Nov 02 2012

Crossrefs

Cf. similar sequences with closed form ((1 + sqrt(2))^(4*r) + (1 - sqrt(2))^(4*r))/8 + k/4: this sequence (k=-1), A076218 (k=3), A278310 (k=-5).

Programs

  • Magma
    [4*Evaluate(ChebyshevU(n), 3)^2: n in [0..30]]; // G. C. Greubel, Aug 18 2022
    
  • Mathematica
    b[n_]:= b[n]= If[n<2, n, 34*b[n-1] -b[n-2] +2]; (* b=A001110 *)
    a[n_]:= 4*b[n]; Table[a[n], {n, 0, 30}]
    4*ChebyshevU[Range[-1,30], 3]^2 (* G. C. Greubel, Aug 18 2022 *)
  • SageMath
    [4*chebyshev_U(n-1, 3)^2 for n in (0..30)] # G. C. Greubel, Aug 18 2022

Formula

a(n) = 4*A001110(n) = A001542(n)^2.
a(n+1) = A001652(n)*A001652(n+1) + A046090(n)*A046090(n+1) = A001542(n+1)^2. - Charlie Marion, Jul 01 2003
a(n) = A001653(k+n)*A001653(k-n) - A001653(k)^2, for k >= n >= 0; e.g. 144 = 5741*5 - 169^2. - Charlie Marion, Jul 16 2003
G.f.: 4*x*(1+x)/((1-x)*(1-34*x+x^2)). - R. J. Mathar, Dec 15 2008
a(n) = A079291(2n). - M. F. Hasler, Jan 16 2012
From George F. Johnson, Nov 02 2012: (Start)
a(n) = ((17+12*sqrt(2))^n + (17-12*sqrt(2))^n - 2)/8.
a(n+1) = 17*a(n) + 4 + 12*sqrt(a(n)*(2*a(n) + 1)).
a(n-1) = 17*a(n) + 4 - 12*sqrt(a(n)*(2*a(n) + 1)).
a(n-1)*a(n+1) = (a(n) - 4)^2.
2*a(n) + 1 = (A001541(n))^2.
a(n+1) = 34*a(n) - a(n-1) + 8 for n>1, a(0)=0, a(1)=4.
a(n+1) = 35*a(n) - 35*a(n-1) + a(n-2) for n>0, a(0)=0, a(1)=4, a(2)=144.
a(n)*a(n+1) = (4*A029549(n))^2.
a(n+1) - a(n) = 4*A046176(n).
a(n) + a(n+1) = 4*(6*A029549(n) + 1).
a(n) = (2*A001333(n)*A000129(n))^2.
Limit_{n -> infinity} a(n)/a(n-r) = (17+12*sqrt(2))^r. (End)
Empirical: a(n) = A089928(4*n-2), for n > 0. - Alex Ratushnyak, Apr 12 2013
a(n) = 4*A001109(n)^2. - G. C. Greubel, Aug 18 2022
Product_{n>=2} (1 - 4/a(n)) = sqrt(2)/3 + 1/2 (Koshy, 2022, section 3, p. 19). - Amiram Eldar, Jan 23 2025

Extensions

Edited and extended by Robert G. Wilson v, Jun 15 2003

A029546 Expansion of 1/( (1-x)*(1-34*x+x^2) ).

Original entry on oeis.org

1, 35, 1190, 40426, 1373295, 46651605, 1584781276, 53835911780, 1828836219245, 62126595542551, 2110475412227490, 71694037420192110, 2435486796874304251, 82734857056306152425, 2810549653117534878200, 95475953348939879706376, 3243371864210838375138585
Offset: 0

Views

Author

Keywords

Comments

Numbers m such that r = 24*m+1 and 2*r-1 are both squares. - Bruno Berselli, Jul 17 2014

Crossrefs

Cf. A029547 (first differences), A245031 (see Comments line).

Programs

  • GAP
    List([0..20], n-> (Lucas(2,-1, 4*n+6)[2] -6)/192 ); # G. C. Greubel, Jan 13 2020
  • Magma
    I:=[1,35,1190]; [n le 3 select I[n] else 35*Self(n-1)-35*Self(n-2) +Self(n-3): n in [1..20]]; // Vincenzo Librandi, Nov 22 2011
    
  • Maple
    seq(coeff(series(1/( (1-x)*(1-34*x+x^2) ), x, n+1), x, n), n = 0..20); # G. C. Greubel, Jan 13 2020
  • Mathematica
    LinearRecurrence[{35,-35,1},{1,35,1190},20] (* Vincenzo Librandi, Nov 22 2011 *)
    Table[(Fibonacci[2*n+3, 2]^2 -1)/24, {n,0,20}] (* G. C. Greubel, Jan 13 2020 *)
  • PARI
    Vec(1/(1-35*x+35*x^2-x^3)+O(x^20)) \\ Charles R Greathouse IV, Sep 23 2012
    
  • Sage
    [(lucas_number2(4*n+6, 2,-1) -6)/192 for n in (0..20)] # G. C. Greubel, Jan 13 2020
    

Formula

a(n) = A029549(n+1)/6 = A075528(n+1)/3.
From Colin Barker, Mar 02 2016: (Start)
a(n) = 35*a(n-1) - 35*a(n-2) + a(n-3) for n > 2.
a(n) = ( (99-70*sqrt(2))*(17-12*sqrt(2))^n - 6 + (99+70*sqrt(2))*(17+12*sqrt(2))^n )/192. (End)
a(n) = ( Pell(2*n+3)^2 - 1)/24 = (Q(4*n+6) - 6)/192, where Q(n) = Pell-Lucas numbers. - G. C. Greubel, Jan 13 2020
Sum_{n>=0} 1/a(n) = 6*(3 - 2*sqrt(2)). - Amiram Eldar, Dec 04 2024

A123480 Coefficients of the series giving the best rational approximations to sqrt(3).

Original entry on oeis.org

4, 60, 840, 11704, 163020, 2270580, 31625104, 440480880, 6135107220, 85451020204, 1190179175640, 16577057438760, 230888624967004, 3215863692099300, 44791203064423200, 623860979209825504, 8689262505873133860, 121025814103014048540, 1685672134936323545704
Offset: 1

Views

Author

Gene Ward Smith, Sep 28 2006

Keywords

Comments

The partial sums of the series 2 - 1/a(1) - 1/a(2) - 1/a(3) - ... give the best rational approximations to sqrt(3), which constitute every second convergent of the continued fraction. The corresponding continued fractions are [1;1,2,1], [1;1,2,1,2,1], [1;1,2,1,2,1,2,1], [1;1,2,1,2,1,2,1,2,1] and so forth.

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[-4*x/((x - 1)*(x^2 - 14*x + 1)), {x, 0, 50}], x] (* G. C. Greubel, Oct 13 2017 *)
  • PARI
    my(x='x+O('x^50)); Vec(-4*x/((x-1)*(x^2-14*x+1))) \\ G. C. Greubel, Oct 13 2017

Formula

a(n+3) = 15*a(n+2) - 15*a(n+1) + a(n).
a(n) = -1/3 + (1/6 + 1/12*3^(1/2))*(7 + 4*3^(1/2))^n + (1/6 - 1/12*3^(1/2))*(7 - 4*3^(1/2))^n.
a(n) = 4*A076139(n) = 2*A217855(n) = 1/2*A045899(n) = 4/3*A076140(n). - Peter Bala, Dec 31 2012
G.f.: -4*x/((x-1)*(x^2-14*x+1)). - Colin Barker, Jan 20 2013
a(n) = A001353(n)*A001353(n+1). - Antonio Alberto Olivares, Apr 06 2020

A200993 Triangular numbers, T(m), that are two-thirds of another triangular number; T(m) such that 3*T(m) = 2*T(k) for some k.

Original entry on oeis.org

0, 10, 990, 97020, 9506980, 931587030, 91286021970, 8945098566040, 876528373449960, 85890835499530050, 8416425350580494950, 824723793521388975060, 80814515339745539060940, 7918997779501541438997070, 775980967875811315482651930, 76038215854050007375860892080
Offset: 0

Views

Author

Charlie Marion, Dec 15 2011

Keywords

Comments

For n>1, a(n) = 98*a(n-1) - a (n-2) + 10. In general, for m>0, let b(n) be those triangular numbers such that for some triangular number c(n), (m+1)*b(n) = m*c(n). Then b(0) = 0, b(1)= A014105(m) and for n>1, b(n) = 2*A069129(m+1)*b(n-1) - b(n-2) + A014105(m).
Further, c(0) = 0, c(1) = A000384(m+1) and for n>1, c(n) = 2*A069129(m+1)*c(n-1) - c(n-2) + A000384(m+1).

Examples

			3*0 = 2*0.
3*10 = 2*15.
3*990 = 2*1485.
3*97020 = 2*145530.
		

Crossrefs

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(10*x/((1-x)*(x^2-98*x+1)))); // G. C. Greubel, Jul 15 2018
  • Mathematica
    LinearRecurrence[{99,-99,1},{0,10,990},20] (* Harvey P. Dale, Feb 25 2018 *)
  • PARI
    concat(0, Vec(10*x/((1-x)*(1-98*x+x^2)) + O(x^40))) \\ Colin Barker, Mar 02 2016
    

Formula

G.f. 10*x / ((1-x)*(x^2-98*x+1)). - R. J. Mathar, Dec 20 2011
a(n) = 99*a(n-1)-99*a(n-2)+a(n-3) for n>2. - Colin Barker, Mar 02 2016
a(n) = (-10+(5-2*sqrt(6))*(49+20*sqrt(6))^(-n)+(5+2*sqrt(6))*(49+20*sqrt(6))^n)/96. - Colin Barker, Mar 07 2016

A336624 Triangular numbers that are one-eighth of other triangular numbers; T(t) such that 8*T(t)=T(u) for some u where T(k) is the k-th triangular number.

Original entry on oeis.org

0, 15, 66, 17391, 76245, 20069280, 87986745, 23159931810, 101536627566, 26726541239541, 117173180224500, 30842405430498585, 135217748442445515, 35592109140254127630, 156041164529401899891, 41073263105447832786516, 180071368649181350028780, 47398510031577658781511915
Offset: 0

Views

Author

Vladimir Pletser, Aug 07 2020

Keywords

Comments

The triangular numbers T(t) that are one-eighth of other triangular numbers T(u) : T(t)=T(u)/8. The t's are in A336623, the T(u)'s are in A336626 and the u's are in A336625.
Can be defined for negative n by setting a(n) = a(-1-n) for all n in Z.

Examples

			a(1)= 15 is a term because it is triangular and 8*15 = 120 is also triangular.
a(2) = 1154*a(0) - a(-2) + 81 = 0 - 15 + 81 = 66;
a(3) = 1154*a(1) - a(-1) + 81 = 1154*15 - 0 + 81 = 17391, etc.
		

Crossrefs

Programs

  • Maple
    f := gfun:-rectoproc({a(n) = 1154*a(n - 2) - a(n - 4) + 81, a(1) = 15, a(0) = 0, a(-1) = 0, a(-2) = 15}, a(n), remember): map(f, [$ (0 .. 40)])[]; #
  • Mathematica
    LinearRecurrence[{1, 1154, -1154, -1, 1}, {0, 15, 66, 17391, 76245}, 18] (* Amiram Eldar, Aug 08 2020 *)
    FullSimplify[Table[((Sqrt[2] + 1)^(4*n + 2)*(11 - 6*(-1)^n*Sqrt[2]) + (Sqrt[2] - 1)^(4*n + 2)*(11 + 6*(-1)^n*Sqrt[2]) - 18)/256, {n, 0, 17}]] (* Vaclav Kotesovec, Sep 08 2020 *)
    Select[Accumulate[Range[0, 10^6]]/8, OddQ[Sqrt[8 # + 1]] &] (* The program generates the first 8 terms of the sequence. *) (* Harvey P. Dale, Jan 15 2024 *)
  • PARI
    concat(0, Vec(3*x*(5 + 17*x + 5*x^2) / ((1 - x)*(1 - 34*x + x^2)*(1 + 34*x + x^2)) + O(x^40))) \\ Colin Barker, Aug 08 2020

Formula

a(n) = 1154*a(n-2) - a(n-4) + 81, for n>=2 with a(1)=15, a(0)=0, a(-1)=0, a(-2)=15.
a(n) = a(n-1) + 1154*a(n-2) - 1154*a(n-3) - a(n-4) + a(n-5), for n>=3 with a(2)=66, a(1)=15, a(0)=0, a(-1)=0, a(-2)=15.
a(n) = b(n)*(b(n)+1)/2 where b(n) is A336623(n).
G.f.: 3*x*(5 + 17*x + 5*x^2) / ((1 - x)*(1 - 34*x + x^2)*(1 + 34*x + x^2)). - Colin Barker, Aug 08 2020
a(n) = ((sqrt(2) + 1)^(4*n + 2) * (11 - 6*(-1)^n*sqrt(2)) + (sqrt(2) - 1)^(4*n + 2) * (11 + 6*(-1)^n*sqrt(2)) - 18)/256. - Vaclav Kotesovec, Sep 08 2020
From Vladimir Pletser, Feb 21 2021: (Start)
a(n) = ((11 - 6*sqrt(2))*(1 + sqrt(2))^(4n + 2) + (11 + 6*sqrt(2))*(1 - sqrt(2) )^(4n + 2) - 18) / 256 for even n.
a(n) = ((11 + 6*sqrt(2))*(1 + sqrt(2) )^(4n + 2) + (11 - 6*sqrt(2))*(1 - sqrt(2) )^(4n + 2) - 18) / 256 for odd n. (End)
128*a(n) = -9+33*A077420(n)-24*(-1)^n*A046176(n+1). - R. J. Mathar, May 05 2023

A336625 Indices of triangular numbers that are eight times other triangular numbers.

Original entry on oeis.org

0, 15, 32, 527, 1104, 17919, 37520, 608735, 1274592, 20679087, 43298624, 702480239, 1470878640, 23863649055, 49966575152, 810661587647, 1697392676544, 27538630330959, 57661384427360, 935502769664975, 1958789677853712, 31779555538278207, 66541187662598864, 1079569385531794079, 2260441590850507680
Offset: 1

Views

Author

Vladimir Pletser, Aug 13 2020

Keywords

Comments

Second member of the Diophantine pair (b(n), a(n)) that satisfies a(n)^2 + a(n) = 8*(b(n)^2 + b(n)) or T(a(n)) = 8*T(b(n)) where T(x) is the triangular number of x. The T(a)'s are in A336626, the T(b)'s are in A336624 and the b's are in A336623.
Can be defined for negative n by setting a(-n) = -a(n+1) - 1 for all n in Z.

Examples

			a(3) = 34*a(1) - a(-1) + 16 = 0 - (-16) + 16 = 32,
a(4) = 34*a(2) - a(0) + 16 = 34*15 - (-1) + 16 = 527, etc.
		

Crossrefs

Programs

  • Maple
    f := gfun:-rectoproc({a(n) = 34*a(n - 2) - a(n - 4) + 16, a(2) = 15, a(1) = 0, a(0) = -1, a(-1) = -16}, a(n), remember); map(f, [$ (0 .. 1000)]); #
  • Mathematica
    LinearRecurrence[{1, 34, -34, -1, 1}, {0, 15, 32, 527, 1104, 17919}, 29] (* Amiram Eldar, Aug 18 2020 *)
    FullSimplify[Table[((Sqrt[2] + 1)^(2*n + 1) * (3 - Sqrt[2]*(-1)^n) - (Sqrt[2] - 1)^(2*n + 1) * (3 + Sqrt[2]*(-1)^n) - 2)/4, {n, 0, 20}]] (* Vaclav Kotesovec, Sep 08 2020 *)
  • PARI
    concat(0, Vec(x*(15 + 17*x - 15*x^2 - x^3) / ((1 - x)*(1 - 6*x + x^2)*(1 + 6*x + x^2)) + O(x^22))) \\ Colin Barker, Aug 14 2020

Formula

a(n) = 34*a(n-2) - a(n-4) + 16, for n>=2 with a(2)=15, a(1)=0, a(0)=-1, a(-1)=-16.
a(n) = a(n-1) + 34*a(n-2) - 34*a(n-3) - a(n-4) + a(n-5), for n>=3 with a(3)=32, a(2)=15, a(1)=0, a(0)=-1, a(-1)=-16.
a(n) = (-1 + sqrt(8*b(n) + 1))/2, where b(n) is A336626(n).
G.f.: x^2*(15 + 17*x - 15*x^2 - x^3) / ((1 - x)*(1 - 6*x + x^2)*(1 + 6*x + x^2)). - Colin Barker, Aug 14 2020
a(n) = ((sqrt(2) + 1)^(2*n+1) * (3 - sqrt(2)*(-1)^n) - (sqrt(2) - 1)^(2*n+1) * (3 + sqrt(2)*(-1)^n) - 2)/4. - Vaclav Kotesovec, Sep 08 2020
From Vladimir Pletser, Feb 21 2021: (Start)
a(n) = ((3 - sqrt(2))*(1 + sqrt(2))^(2*n+1) + (3 + sqrt(2))*(1 - sqrt(2))^(2*n+1))/4 - 1/2 for even n.
a(n) = ((3 + sqrt(2))*(1 + sqrt(2))^(2*n+1) + (3 - sqrt(2))*(1 - sqrt(2))^(2*n+1))/4 - 1/2 for odd n. (End)

A200994 Triangular numbers, T(m), that are three-halves of another triangular number; T(m) such that 2*T(m) = 3*T(k) for some k.

Original entry on oeis.org

0, 15, 1485, 145530, 14260470, 1397380545, 136929032955, 13417647849060, 1314792560174940, 128836253249295075, 12624638025870742425, 1237085690282083462590, 121221773009618308591410, 11878496669252312158495605, 1163971451813716973223977895
Offset: 0

Views

Author

Charlie Marion, Feb 15 2012

Keywords

Comments

For n > 1, a(n) = 98*a(n-1) - a(n-2) + 15. In general, for m>0, let b(n) be those triangular numbers such that for some triangular number c(n), (m+1)*b(n) = m*c(n). Then b(0) = 0, b(1) = A014105(m) and for n > 1, b(n) = 2*A069129(m+1)*b(n-1) - b(n-2) + A014105(m). Further, c(0) = 0, c(1) = A000384(m+1) and for n>1, c(n) = 2*A069129(m+1)*c(n-1) - c(n-2) + A000384(m+1).

Examples

			2*0 = 3*0.
2*15 = 3*10.
2*1485 = 3*990.
2*145530 = 3*97020.
		

Crossrefs

Programs

  • Magma
    m:=25; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(15*x/((1-x)*(1-98*x+x^2)))); // G. C. Greubel, Jul 15 2018
  • Mathematica
    LinearRecurrence[{99, -99, 1}, {0, 15, 1485}, 20] (* T. D. Noe, Feb 15 2012 *)
  • PARI
    concat(0, Vec(15*x/((1-x)*(1-98*x+x^2)) + O(x^20))) \\ Colin Barker, Mar 02 2016
    

Formula

From Colin Barker, Mar 02 2016: (Start)
a(n) = 99*a(n-1) - 99*a(n-2) + a(n-3) for n>2.
G.f.: 15*x / ((1-x)*(1-98*x+x^2)). (End)
a(n) = (-10+(5-2*sqrt(6))*(49+20*sqrt(6))^(-n)+(5+2*sqrt(6))*(49+20*sqrt(6))^n)/64. - Colin Barker, Mar 03 2016

A336623 First member of the Diophantine pair (m, k) that satisfies 8*(m^2 + m) = k^2 + k; a(n) = m.

Original entry on oeis.org

0, 5, 11, 186, 390, 6335, 13265, 215220, 450636, 7311161, 15308375, 248364270, 520034130, 8437074035, 17665852061, 286612152936, 600118935960, 9736376125805, 20386377970595, 330750176124450, 692536732064286, 11235769612105511, 23525862512215145, 381685416635462940
Offset: 0

Views

Author

Vladimir Pletser, Aug 07 2020

Keywords

Comments

The indices of triangular numbers that are one-eighth of other triangular numbers [m of T(m) such that T(m)=T(k)/8]. The T(m)'s are in A336624, the T(k)'s are in A336626 and the k's are in A336625.
Also, nonnegative m such that 32*m^2 + 32*m + 1 is a square.
Can be defined for negative n by setting a(n) = a(-1-n) for all n in Z.

Examples

			a(2) = 34 a(0) - a(-2)+16=0 -5 +16 = 11 ; a(3) = 34 a(1) - a(-1)+16 = 34*5 -0 +16 = 186, etc.
		

Crossrefs

Programs

  • Maple
    f := gfun:-rectoproc({a(n) = 34*a(n - 2) - a(n - 4) + 16, a(1) = 5, a(0) = 0, a(-1) = 0,  a(-2) = 5}, a(n), remember); map(f, [$ (0 .. 50)]); #
  • Mathematica
    LinearRecurrence[{1, 34, -34, -1, 1}, {0, 5, 11, 186, 390}, 24] (* Amiram Eldar, Aug 08 2020 *)
    FullSimplify[Table[((3*Sqrt[2] - 2*(-1)^n)*(1 + Sqrt[2])^(2*n + 1) + (3*Sqrt[2] + 2*(-1)^n)*(Sqrt[2] - 1)^(2*n + 1) - 8)/16, {n, 0, 20}]] (* Vaclav Kotesovec, Sep 08 2020 *)
  • PARI
    concat(0, Vec(x*(5 + 6*x + 5*x^2) / ((1 - x)*(1 - 6*x + x^2)*(1 + 6*x + x^2)) + O(x^22))) \\ Colin Barker, Aug 08 2020

Formula

a(n) = 34 a(n-2) - a(n-4) + 16 for n>=2, with a(1)=5, a(0)=0, a(-1)=0, a(-2)=5.
a(n) = a(n-1) + 34 a(n-2) - 34 a(n-3) - a(n-4)+ a(n-5) for n>=3 with a(2)=11, a(1)=5, a(0)=0, a(-1)=0, a(-2)=5.
a(n) = (C+((-1)^n)*D)*A^n + (E+((-1)^n)*F)*B^n -1/2 with A = (sqrt(2) + 1)^2 ; B = (sqrt(2) - 1)^2 ; C = 3*(2 + sqrt(2))/16 ; D = -(1 + sqrt(2))/8 ; E = 3*(2 - sqrt(2))/16 ; F = (sqrt(2) - 1)/8 and n>=0.
a(n) = (-1 + sqrt(8*b(n) + 1))/2 where b(n) = A336624(n).
G.f.: x*(5 + 6*x + 5*x^2) / ((1 - x)*(1 - 6*x + x^2)*(1 + 6*x + x^2)). - Colin Barker, Aug 08 2020
a(n) = ((3*sqrt(2) - 2*(-1)^n) * (1 + sqrt(2))^(2*n + 1) + (3*sqrt(2) + 2*(-1)^n) * (sqrt(2) - 1)^(2*n + 1) - 8)/16. - Vaclav Kotesovec, Sep 08 2020
Comment from _Vladimir Pletser, Feb 21 2021: (Start)
a(n) = ((4 + sqrt(2))(1 + sqrt(2))^(2n) + (4 - sqrt(2))(1 - sqrt(2))^(2n))/16 - 1/2 for even n.
a(n) = ((8 + 5 sqrt(2))(1 + sqrt(2))^(2n) + (8 - 5 sqrt(2))(1 - sqrt(2))^(2n))/16 - 1/2 for odd n. (End)

A336626 Triangular numbers that are eight times another triangular number.

Original entry on oeis.org

0, 120, 528, 139128, 609960, 160554240, 703893960, 185279454480, 812293020528, 213812329916328, 937385441796000, 246739243443988680, 1081741987539564120, 284736873122033021040, 1248329316235215199128, 328586104843582662292128, 1440570949193450800230240, 379188080252621270252095320
Offset: 1

Views

Author

Vladimir Pletser, Oct 04 2020

Keywords

Comments

The triangular numbers T(t) that are eight times another triangular number T(u) : T(t) = 8*T(u). The t's are in A336625, the T(u)'s are in A336624 and the u's are in A336623.
Can be defined for negative n by setting a(n) = a(-1-n) for all n in Z.

Examples

			a(2) = 120 is a term because it is triangular and 120/8 = 15 is also triangular.
a(3) = 1154*a(1) - a(-1) + 648 = 0 - 120 + 648 = 528;
a(4) = 1154*a(2) - a(0) + 648 = 1154*120 - 0 + 648 = 139128, etc.
.
From _Peter Luschny_, Oct 19 2020: (Start)
Related sequences in context, as computed by the Julia function:
n   [A336623, A336624,        A336625,  A336626        ]
[0] [0,       0,              0,        0              ]
[1] [5,       15,             15,       120            ]
[2] [11,      66,             32,       528            ]
[3] [186,     17391,          527,      139128         ]
[4] [390,     76245,          1104,     609960         ]
[5] [6335,    20069280,       17919,    160554240      ]
[6] [13265,   87986745,       37520,    703893960      ]
[7] [215220,  23159931810,    608735,   185279454480   ]
[8] [450636,  101536627566,   1274592,  812293020528   ]
[9] [7311161, 26726541239541, 20679087, 213812329916328] (End)
		

Crossrefs

Programs

  • Julia
    function omnibus()
        println("[A336623, A336624, A336625, A336626]")
        println([0, 0, 0, 0])
        t, h = 1, 1
        for n in 1:999999999
            d, r = divrem(t, 8)
            if r == 0
                d2 = 2*d
                s = isqrt(d2)
                d2 == s * (s + 1) && println([s, d, n, t])
            end
            t, h = t + h + 1, h + 1
        end
    end
    omnibus() # Peter Luschny, Oct 19 2020
  • Maple
    f := gfun:-rectoproc({a(n) = 1154*a(n - 2) - a(n - 4) + 648, a(2) = 120, a(1) = 0, a(0) = 0, a(-1) = 120}, a(n), remember); map(f, [$ (1 .. 1000)])[]; #
  • Mathematica
    LinearRecurrence[{1, 1154, -1154, -1, 1}, {0, 120, 528, 139128, 609960}, 18]

Formula

a(n) = 8*A336624(n).
a(n) = 1154*a(n-2) - a(n-4) + 648, for n>=2 with a(2)=120, a(1)=0, a(0)=0, a(-1)=120.
a(n) = a(n-1) + 1154*a(n-2) - 1154*a(n-3) - a(n-4) + a(n-5), for n>=3 with a(3)=528, a(2)=120, a(1)=0, a(0)=0, a(-1)=120.
a(n) = ((10*sqrt(2))/17 + 15/17)*(17 + 12*sqrt(2))^n + (-(10*sqrt(2))/17 + 15/17)*(17 - 12*sqrt(2))^n + (-15/17 - (45*sqrt(2))/68)*(-17 - 12*sqrt(2))^n + (-15/17 + (45*sqrt(2))/68)*(-17 + 12*sqrt(2))^n - 27*(-4 + 3*sqrt(2))*sqrt(2)*(-1/(-17 + 12*sqrt(2)))^n/(1088*(-17 + 12*sqrt(2))) - 27*(4 + 3*sqrt(2))*sqrt(2)*(-1/(-17 - 12*sqrt(2)))^n/(1088*(-17 - 12*sqrt(2))) - 9/16 - 9*(-3 + 2*sqrt(2))*sqrt(2)*(-1/(17 - 12*sqrt(2)))^n/(272*(17 - 12*sqrt(2))) - 9*(3 + 2*sqrt(2))*sqrt(2)*(-1/(17 + 12*sqrt(2)))^n/(272*(17 + 12*sqrt(2))).
Let b(n) be A336625(n). Then a(n) = b(n)*(b(n)+1)/2.
G.f.: 24*x^2*(5 + 17*x + 5*x^2)/(1 - x - 1154*x^2 + 1154*x^3 + x^4 - x^5). - Stefano Spezia, Oct 05 2020
From Vladimir Pletser, Feb 21 2021: (Start)
a(n) = ((11*(1 + sqrt(2))^2 - (-1)^n*6*(4 + 3*sqrt(2)))*(1 + sqrt(2))^(4n) + (11*(1 - sqrt(2))^2 - (-1)^n*6*(4 - 3*sqrt(2)))*(1 - sqrt(2))^(4n))/32 - 9/16.
a(n) = ((1 + 2*sqrt(2))^2*(1 + sqrt(2))^(4n) + (1 - 2*sqrt(2))^2*(1 - sqrt(2))^(4n))/32 - 9/16 for even n.
a(n) = ((5 + 4*sqrt(2))^2*(1 + sqrt(2))^(4n) + (5 - 4*sqrt(2))^2*(1 - sqrt(2))^(4n))/32 - 9/16 for odd n. (End)

A098602 a(n) = A001652(n) * A046090(n).

Original entry on oeis.org

0, 12, 420, 14280, 485112, 16479540, 559819260, 19017375312, 646030941360, 21946034630940, 745519146510612, 25325704946729880, 860328449042305320, 29225841562491651012, 992818284675673829100, 33726595837410418538400, 1145711440187278556476512
Offset: 0

Views

Author

Charlie Marion, Oct 26 2004

Keywords

Comments

From Ron Knott, Nov 25 2013: (Start)
a(n) = 2*r*(r+1) which is also of form s(s+1) where the s is in A053141.
a(n) is an oblong number (A002378) which is twice another oblong number. (End)
2*a(n)+1 and 4*a(n)+1 are both square. - Paul Cleary, Jun 23 2014

Examples

			a(1) = 12 = 2(2*3) = 3*4, a(2) = 420 = 2(14*15) = 20*21.
		

Crossrefs

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(12*x/((1-x)*(x^2-34*x+1)))); // G. C. Greubel, Jul 15 2018
  • Mathematica
    2*Table[ Floor[(Sqrt[2] + 1)^(4n + 2)/32], {n, 0, 20} ] (* Ray Chandler, Nov 10 2004, copied incorrect program from A029549, revised Jul 09 2015 *)
    RecurrenceTable[{a[n+3] == 35 a[n+2] - 35 a[n+1] + a[n], a[1] == 0, a[2] == 12, a[3] == 420}, a, {n, 1, 10}] (* Ron Knott, Nov 25 2013 *)
    LinearRecurrence[{35, -35, 1}, {0, 12, 420}, 25] (* T. D. Noe, Nov 25 2013 *)
    Table[(LucasL[4*n+2, 2] - 6)/16, {n,0,30}] (* G. C. Greubel, Jul 15 2018 *)
  • PARI
    concat(0, Vec(12*x/((1-x)*(1-34*x+x^2)) + O(x^20))) \\ Colin Barker, Mar 02 2016
    
  • PARI
    {a=1+sqrt(2); b=1-sqrt(2); Q(n) = a^n + b^n};
    for(n=0, 30, print1(round((Q(4*n+2) - 6)/16), ", ")) \\ G. C. Greubel, Jul 15 2018
    

Formula

a(n) = 2*A029549(n) = 2*A001109(n)*A001109(n+1).
a(n) = (A001653(n)^2 - 1)/2.
a(n) = A053141(n)^2 + A011900(n)^2 - 1.
For n>0, a(n) = A053141(2n) - 2*A001109(n-1)^2.
For n>0, a(n) = 3*(A001542(n)^2 - A001542(n-1)^2).
For n>0, a(n) = A053141(2n-1) + 2*(A001653(2n-1) - A001109(n-1)^2).
a(n+1) + a(n) = 3*A001542(n+1)^2.
a(n+1) - a(n) = A001542(2*n).
a(n+1)*a(n) = 4*(A001109(n)^4 - A001109(n)^2) = 4*A001110(n)*(A001110(n) - 1).
From Ron Knott, Nov 25 2013: (Start)
a(n) = 35*a(n-1) - 35*a(n-2) + a(n-3).
G.f.: 12*x / ((1-x)*(x^2-34*x+1)). (End)
a(n) = (-6 + (3-2*sqrt(2))*(17+12*sqrt(2))^(-n)+(3+2*sqrt(2))*(17+12*sqrt(2))^n)/16. - Colin Barker, Mar 02 2016

Extensions

More terms from Ray Chandler, Nov 10 2004
Corrected by Bill Lam (bill_lam(AT)myrealbox.com), Feb 27 2006
Previous Showing 11-20 of 38 results. Next