OFFSET
0,2
COMMENTS
LINKS
Vladimir Pletser, Table of n, a(n) for n = 0..650
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, Using Pell equation solutions to find all triangular numbers multiple of other triangular numbers, 2022.
Index entries for linear recurrences with constant coefficients, signature (1,1154,-1154,-1,1).
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)
EXAMPLE
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.
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 *)
PROG
(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
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladimir Pletser, Aug 07 2020
STATUS
approved