login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A077398 First member of the Diophantine pair (m,k) that satisfies 7*(m^2+m) = k^2+k; a(n)=m. 6
0, 2, 5, 39, 87, 629, 1394, 10032, 22224, 159890, 354197, 2548215, 5644935, 40611557, 89964770, 647236704, 1433791392, 10315175714, 22850697509, 164395574727, 364177368759, 2620014019925, 5803987202642, 41755828744080, 92499617873520, 665473245885362 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

Equivalently, m such that 28*m*(m+1)+1 is a square. - Bruno Berselli, May 19 2014

LINKS

Colin Barker, Table of n, a(n) for n = 0..1000

Mohammad K. Azarian, Diophantine Pair, Problem B-881, Fibonacci Quarterly, Vol. 37, No. 3, August 1999, pp. 277-278; Solution to Problem B-881, Fibonacci Quarterly, Vol. 38, No. 2, May 2000, pp. 183-184.

Index entries for linear recurrences with constant coefficients, signature (1,16,-16,-1,1).

FORMULA

G.f.: x*(2+3*x+2*x^2)/((1-x)*(1-16*x^2+x^4)).

a(n) = 16*a(n-2) - a(n-4) + 7, n >= 3. [corrected by Vladimir Pletser, Feb 29 2020]

Let b(n) be A077397 then a(n+2) = 2*a(n+1) - a(n) + b(n) with a(0)=0 a(1)=2.

a(0)=0, a(1)=2; a(n+2) = (7 + 16*a(n) + 3*sqrt((1+28*a(n)+28*a(n)^2))/2. - Herbert Kociemba, May 12 2008

a(n) = a(n-1) + 16*a(n-2) - 16*a(n-3) - a(n-4) + a(n-5). - Wesley Ivan Hurt, Jul 24 2020

MAPLE

f := gfun:-rectoproc({a(-2) = 2, a(-1) = 0, a(0) = 0, a(1) = 2, a(n) = 16*a(n - 2) - a(n - 4) + 7}, a(n), remember): map(f, [$ (0 .. 40)])[]; # Vladimir Pletser, Jul 24 2020

MATHEMATICA

LinearRecurrence[{1, 16, -16, -1, 1}, {0, 2, 5, 39, 87}, 30] (* G. C. Greubel, Jan 18 2018 *)

PROG

(PARI) {a(n)=if(n<0, 0, polcoeff(x*(2+3*x+2*x^2)/((1-x)*(1-16*x^2+x^4)) + x*O(x^n), n))};

(MAGMA) I:=[0, 2, 5, 39, 87]; [n le 5 select I[n] else Self(n-1)+16*Self(n-2) -16*Self(n-3)-Self(n-4)+Self(n-5): n in [1..30]]; // G. C. Greubel, Jan 18 2018

CROSSREFS

Cf. A077397, A077399, A077400. The k values are in A077401.

Cf. A053141.

Sequence in context: A228837 A206155 A135378 * A255550 A067083 A183255

Adjacent sequences:  A077395 A077396 A077397 * A077399 A077400 A077401

KEYWORD

nonn,easy

AUTHOR

Bruce Corrigan (scentman(AT)myfamily.com), Nov 05 2002

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 23 21:27 EDT 2020. Contains 337315 sequences. (Running on oeis4.)