A082405 -id:A082405 - cp's OEIS frontend

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

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

A041041 Denominators of continued fraction convergents to sqrt(26).

Original entry on oeis.org

1, 10, 101, 1020, 10301, 104030, 1050601, 10610040, 107151001, 1082120050, 10928351501, 110365635060, 1114584702101, 11256212656070, 113676711262801, 1148023325284080, 11593909964103601, 117087122966320090, 1182465139627304501, 11941738519239365100

Offset: 0

N. J. A. Sloane

Generalized Fibonacci sequence.
Sqrt(26) = 10/2 + 10/101 + 10/(101*10301) + 10/(10301*1050601) + ... - Gary W. Adamson, Jun 13 2008
For positive n, a(n) equals the permanent of the n X n tridiagonal matrix with 10's along the main diagonal and 1's along the superdiagonal and the subdiagonal. - John M. Campbell, Jul 08 2011
a(n) equals the number of words of length n on alphabet {0, 1, ..., 10} avoiding runs of zeros of odd lengths. - Milan Janjic, Jan 28 2015
From Bruno Berselli, May 03 2018: (Start)
Numbers k for which m*k^2 + (-1)^k is a perfect square:
m =  2: 0, 1, 2, 5, 12, 29, 70, 169, ...   (A000129);
m =  3: 0, 4, 56, 780, 10864, 151316, ...  (4*A007655);
m =  5: 0, 1, 4, 17, 72, 305, 1292, ...    (A001076);
m =  6: 0, 2, 20, 198, 1960, 19402, ...    (A001078);
m =  7: 0, 48, 12192, 3096720, ...         (2*A175672);
m =  8: 0, 6, 204, 6930, 235416, ...       (A082405);
m = 10: 0, 1, 6, 37, 228, 1405, 8658, ...  (A005668);
m = 11: 0, 60, 23880, 9504180, ...         [°];
m = 12: 0, 2, 28, 390, 5432, 75658, ...    (A011944);
m = 13: 0, 5, 180, 6485, 233640, ...       (5*A041613);
m = 14: 0, 4, 120, 3596, 107760, ...       (A068204);
m = 15: 0, 8, 496, 30744, 1905632, ...     [°];
m = 17: 0, 1, 8, 65, 528, 4289, 34840, ... (A041025);
m = 18: 0, 4, 136, 4620, 156944, ...       (A202299);
m = 19: 0, 13260, 1532829480, ...          [°];
m = 20: 0, 2, 36, 646, 11592, 208010, ...  (A207832);
m = 21: 0, 12, 1320, 145188, ...           (A174745);
m = 22: 0, 42, 16548, 6519870, ...         (A174766);
m = 23: 0, 240, 552480, 1271808720, ...    [°];
m = 24: 0, 10, 980, 96030, 9409960, ...    (A168520);
m = 26: 0, 1, 10, 101, 1020, 10301, ...    (this sequence);
m = 27: 0, 260, 702520, 1898208780, ...    [°];
m = 28: 0, 24, 6096, 1548360, ...          (A175672);
m = 29: 0, 13, 1820, 254813, 35675640, ... [°];
m = 30: 0, 2, 44, 966, 21208, 465610, ...  (2*A077421), etc.
[°] apparently without related sequences in the OEIS.
(End)
From Michael A. Allen, Mar 12 2023: (Start)
Also called the 10-metallonacci sequence; the g.f. 1/(1-k*x-x^2) gives the k-metallonacci sequence.
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using unit squares and dominoes (with dimensions 2 X 1) if there are 10 kinds of squares available. (End)

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Michael A. Allen and Kenneth Edwards, Fence tiling derived identities involving the metallonacci numbers squared or cubed, Fib. Q. 60:5 (2022) 5-17.
Sergio Falcón and Ángel Plaza, On the Fibonacci k-numbers, Chaos, Solitons & Fractals 2007; 32(5): 1615-24.
Sergio Falcón and Ángel Plaza, The k-Fibonacci sequence and the Pascal 2-triangle Chaos, Solitons & Fractals 2007; 33(1): 38-49.
Sergio Falcón and Ángel Plaza, On k-Fibonacci sequences and polynomials and their derivatives, Chaos, Solitons & Fractals (2007).
Milan Janjic, On Linear Recurrence Equations Arising from Compositions of Positive Integers, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.7.
Tanya Khovanova, Recursive Sequences
Kai Wang, On k-Fibonacci Sequences And Infinite Series List of Results and Examples, 2020.
Index entries for sequences related to Chebyshev polynomials.
Index entries for linear recurrences with constant coefficients, signature (10,1).

Cf. A099374 (squares).
Cf. A041040.
Cf. A000045, A000129, A006190, A001076, A052918, A005668, A054413, A041025, A099371, A243399.
Row n=10 of A073133, A172236 and A352361.

I:=[1,10]; [n le 2 select I[n] else 10*Self(n-1) + Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 24 2018

seq(combinat:-fibonacci(n+1, 10), n=0..19); # Peter Luschny, May 04 2018

Denominator[Convergents[Sqrt[26], 30]] (* Vincenzo Librandi, Dec 10 2013 *)
LinearRecurrence[{10,1}, {1,10}, 30] (* G. C. Greubel, Jan 24 2018 *)

x='x+O('x^30); Vec(1/(1-10*x-x^2)) \\ G. C. Greubel, Jan 24 2018

[lucas_number1(n,10,-1) for n in range(1, 19)] # Zerinvary Lajos, Apr 26 2009

G.f.: 1/(1 - 10*x - x^2).
a(n) = 10*a(n-1) + a(n-2), n>=1; a(-1):=0, a(0)=1.
a(n) = S(n, 10*i)*(-i)^n where i^2:=-1 and S(n, x)=U(n, x/2) Chebyshev's polynomials of the second kind. See A049310.
a(n) = (ap^(n+1) - am^(n+1))/(ap-am) with ap = 5+sqrt(26), am = -1/ap = 5-sqrt(26).
a(n) = F(n+1, 10), the (n+1)-th Fibonacci polynomial evaluated at x=10. - T. D. Noe, Jan 19 2006
a(n) = Sum_{i=0..floor(n/2)} binomial(n-i,i)*10^(n-2*i). - Sergio Falcon, Sep 24 2007

Extended by T. D. Noe, May 23 2011

A278310 Numbers m such that T(m) + 3*T(m+1) is a square, where T = A000217.

Original entry on oeis.org

3, 143, 4899, 166463, 5654883, 192099599, 6525731523, 221682772223, 7530688524099, 255821727047183, 8690408031080163, 295218051329678399, 10028723337177985443, 340681375412721826703, 11573138040695364122499, 393146012008229658338303, 13355391270239113019379843

Offset: 1

Bruno Berselli, Nov 17 2016

Equivalently, both m+1 and 2*m+3 are squares for nonnegative m.
Corresponding triangular numbers T(m): 6, 10296, 12002550, 13855048416, 15988853699286, 18451128064030200, 21292585958400815526, ...
Square roots of T(m) + 3*T(m+1) are listed by A082405 (after 0).
Negative values of m for which T(m) + 3*T(m+1) is a square: -1, -2, -26, -842, -28562, -970226, -32959082, ...

			3 is in the sequence because T(3) + 3*T(4) = 6 + 3*10 = 6^2.
For n=5 is a(5) = 5654883, therefore floor(sqrt(5654883)) = 2377 = A182189(5) - 2 = 2379 - 2.

Colin Barker, Table of n, a(n) for n = 1..650
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).

Subsequence of A000466.
Cf. A000217, A001109, A029546, A046176, A082405, A156164, A182189.
Cf. A278438: numbers m such that T(m) + 2*T(m+1) is a square.
Cf. A078522: numbers m such that 3*T(m) + T(m+1) is a square.
Cf. similar sequences with closed form ((1 + sqrt(2))^(4*r) + (1 - sqrt(2))^(4*r))/8 + k/4: A084703 (k=-1), A076218 (k=3), this sequence (k=-5).

Iv:=[3,143]; [n le 2 select Iv[n] else 34*Self(n-1)-Self(n-2)+40: n in [1..20]];

P:=proc(q) local n; for n from 3 to q do if type(sqrt(2*n^2+5*n+3),integer) then print(n); fi; od; end: P(10^9); # Paolo P. Lava, Nov 18 2016

Table[((1 + Sqrt[2])^(4 n) + (1 - Sqrt[2])^(4 n))/8 - 5/4, {n, 1, 20}]
RecurrenceTable[{a[1] == 3, a[2] == 143, a[n] == 34 a[n - 1] - a[n - 2] + 40}, a, {n, 1, 20}]
LinearRecurrence[{35, -35, 1}, {3, 143, 4899}, 50] (* G. C. Greubel, Nov 20 2016 *)

Vec(x*(3 + 38*x - x^2)/((1 - x)*(1 - 34*x + x^2)) + O(x^50)) \\ G. C. Greubel, Nov 20 2016

def A278310():
    a, b = 3, 143
    yield a
    while True:
        yield b
        a, b = b, 34*b - a + 40
a = A278310(); print([next(a) for  in range(18)]) # _Peter Luschny, Nov 18 2016

O.g.f.: x*(3 + 38*x - x^2)/((1 - x)*(1 - 34*x + x^2)).
E.g.f.: (exp((1-sqrt(2))^4*x) + exp((1+sqrt(2))^4*x) - 10*exp(x))/8 + 1.
a(n) = 35*a(n-1) - 35*a(n-2) + a(n-3) for n>3.
a(n) = 34*a(n-1) - a(n-2) + 40 for n>2.
a(n) = a(-n) = ((1 + sqrt(2))^(4*n) + (1 - sqrt(2))^(4*n))/8 - 5/4.
a(n) = 4*A001109(n)^2 - 1.
a(n) = -A029546(n) + 38*A029546(n-1) + 3*A029546(n-2) for n>1.
Lim_{n -> infinity} a(n)/a(n-1) = A156164.
Floor(sqrt(a(n))) = A182189(n) - 2.
a(n) - a(n-1) = 4*A046176(n) for n>1.

A273626 A fourth-order divisibility sequence: a(n) = (1/14)(Pell(4n) + Pell(2*n)).

Original entry on oeis.org

1, 30, 995, 33660, 1142629, 38810970, 1318402631, 44786716920, 1521429030985, 51683794848150, 1755727563817259, 59643053188493940, 2026108079758297261, 68828031652259981010, 2338126968060165944975, 79427488882178225107440, 2698196495024745460575889

Offset: 1

Peter Bala, May 31 2016

This is a divisibility sequence, that is, a(n) divides a(m) whenever n divides m. The sequence satisfies a linear recurrence of order 4.

A000129(n) = Pell(n) is the Lucas sequence U_n(2,-1). In general, if U(n) = U_n(P,Q) is the Lucas sequence with integer parameters P and Q then when Q = 1 or Q = -1 both U(4*n) + U(2*n) and U(4*n) - 2*U(2*n) are divisibility sequences of the fourth order. Cf. A127595, A215466 and A273627.

Colin Barker, Table of n, a(n) for n = 1..600
P. Bala, Lucas sequences and divisibility sequences
Wikipedia, Lucas Sequence
Index entries for linear recurrences with constant coefficients, signature (40,-206,40,-1).

Cf. A000129, A001109, A029547, A081555, A127595, A215466, A273627.

I:=[1,30,995,33660]; [n le 4 select I[n] else 40*Self(n-1)-206*Self(n-2)+40*Self(n-3)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Jun 04 2016

#A273626
A000129 := proc (n) option remember;
if n <= 1 then n else 2*A000129(n-1) + A000129(n-2) end if
end proc:
seq(1/14*(A000129(4*n) + A000129(2*n)), n = 1..20);

LinearRecurrence[{40,-206,40,-1},{1,30,995,33660},100] (* G. C. Greubel, Jun 02 2016 *)

a(n) = sqrt(2)/56*( (sqrt(2) + 1)^(4*n) - (sqrt(2) - 1)^(4*n) + (sqrt(2) + 1)^(2*n) - (sqrt(2) - 1)^(2*n) ).
a(n) = (A082405(n) + A001109(n))/7 .
a(n) = 1/14*Pell(2*n)*A081555(n).
a(n) = -a(-n).
a(n) = 40*a(n-1) - 206*a(n-2) + 40*a(n-3) - a(n-4) for n>4.
O.g.f.: x*(x^2 - 10*x + 1)/((x^2 - 6*x + 1)*(x^2 - 34*x + 1)).

cp's OEIS Frontend

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

A041041 Denominators of continued fraction convergents to sqrt(26).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A278310 Numbers m such that T(m) + 3*T(m+1) is a square, where T = A000217.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A273626 A fourth-order divisibility sequence: a(n) = (1/14)*(Pell(4*n) + Pell(2*n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Formula

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

A273626 A fourth-order divisibility sequence: a(n) = (1/14)(Pell(4n) + Pell(2*n)).