The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A056161 Solutions (value of x) of Diophantine equation 2*x^2 + 3*x + 2 = r^2. 3
2, 7, 94, 263, 3218, 8959, 109342, 304367, 3714434, 10339543, 126181438, 351240119, 4286454482, 11931824527, 145613270974, 405330793823, 4946564758658, 13769315165479, 168037588523422, 467751384832487, 5708331445037714, 15889777769139103, 193915231542758878 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,1
COMMENTS
The same equation also has negative solutions x=-c(n), where c would be the sequence {1,2,17,46,553,1538,18761,52222,...} with the corresponding values of r being {1,2,23,64,781,2174,26531,73852,...}. Moreover, replacing x with x+K, one obtains the Diophantine equation 2*x^2+(4*K+3)*x+(2*K^2+3*K+2)=r^2. Since K can be any integer (for example K=-1, giving 2*x^2-x+1=r^2), this amounts to an infinite family of Diophantine equations with closely related solutions. For example, if the present equation has a solution pair {a(n), A055979(n)}, the one with x replaced by x+K will have a solution {a(n)-K, A055979(n)}.
LINKS
Seon-Hong Kim and Kenneth B. Stolarsky, Translations and extensions of the Nicomachus identity, arXiv:2306.17402 [math.NT], 2023.
FORMULA
a(n) = floor(A055979(n)/sqrt(2)).
G.f.: (x^4 + x^3 - 19*x^2 - 5*x - 2) / (x^5 - x^4 - 34*x^3 + 34*x^2 + x - 1). - Alois P. Heinz, Jun 03 2009
a(n) = a(n-1) + 34*a(n-2) - 34*a(n-3) - a(n-4) + a(n-5). - Colin Barker, May 17 2015
MAPLE
a:= n-> (Matrix([94, 7, 2, -1, -2]). Matrix([[1, 1, 0, 0, 0], [34, 0, 1, 0, 0], [ -34, 0, 0, 1, 0], [ -1, 0, 0, 0, 1], [1, 0, 0, 0, 0]])^n)[1, 3]: seq(a(n), n=0..25); # Alois P. Heinz, Jun 03 2009
MATHEMATICA
CoefficientList[Series[(x^4 + x^3 - 19 x^2 - 5 x - 2)/(x^5 - x^4 - 34 x^3 + 34 x^2 + x - 1), {x, 0, 22}], x] (* Michael De Vlieger, Jan 09 2016 *)
LinearRecurrence[{1, 34, -34, -1, 1}, {2, 7, 94, 263, 3218}, 30] (* Vincenzo Librandi, Jan 10 2016 *)
PROG
(PARI) Vec((x^4+x^3-19*x^2-5*x-2)/((x-1)*(x^2-6*x+1)*(x^2+6*x+1)) + O(x^100)) \\ Colin Barker, May 17 2015
(Magma) I:=[2, 7, 94, 263, 3218]; [n le 5 select I[n] else Self(n-1)+34*Self(n-2)-34*Self(n-3)-Self(n-4)+Self(n-5): n in [1..30]]; // Vincenzo Librandi, Jan 10 2016
CROSSREFS
Cf. A055979.
Sequence in context: A350754 A224439 A304722 * A076740 A112290 A072059
KEYWORD
nonn,easy
AUTHOR
Helge Robitzsch (hrobi(AT)math.uni-goettingen.de), Aug 01 2000
EXTENSIONS
More terms from Alois P. Heinz, Jun 03 2009
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 28 19:24 EDT 2024. Contains 372919 sequences. (Running on oeis4.)