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.

Showing 1-1 of 1 results.

A336449 Values z of primitive solutions (x, y, z) to the Diophantine equation x^3 + y^3 + 2*z^3 = 2*4^6.

Original entry on oeis.org

1, 9, 16, 25, 49, 81, 121, 169, 225, 289, -356, 361, 441, 529, 625, 729, 841, -948, 961, 1045, 1089, 1225, 1369, 1521, 1681, -1715, 1849, 1876, 2025, 2209, 2401, 2601, 2809, 3025, 3249, 3481, -3587, 3721, 3969, 4225, 4489, 4761, 5041, 5329, 5625, 5769, 5929
Offset: 1

Views

Author

XU Pingya, Aug 08 2020

Keywords

Comments

Terms are arranged in order of increasing absolute value (if equal, the negative number comes first).
Let x = a^(2*k) - (a^k)*t - t^2, y = a^(2*k) + (a^k)*t - t^2, z = t^2; then x^3 + y^3 + 2*z^3 = 2*a^(6*k). When a = 4, k = 1, t = 2*n + 1; (x, y, z) are primitive solutions of equation. Thus, terms of A016754 are terms of the sequence.

Examples

			(-15)^3 + (-27)^3 + 2*25^3 = 11^3 + (-29)^3 + 2*25^3 = 8192, 25 is a term.
(-65)^3 + (449)^3 + 2*(-356)^3 = 8192, -356 is a term.
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, D5.

Crossrefs

Programs

  • Mathematica
    Clear[t]
    t = {};
    Do[y = (8192 - x^3 - 2z^3)^(1/3) /. (-1)^(1/3) -> -1;
    If[Abs@x <= Abs@y && IntegerQ[y] && GCD[x, y, z] == 1, AppendTo[t, z]], {z, -5929, 5929}, {x, -Round[(Abs[8192 - 2z^3]/3)^(1/2)], Round[(Abs[8192 - 2z^3]/3)^(1/2)]}]
    u = Union@t;
    v = Table[(-1)^n*Floor[(n + 1)/2], {n, 0, 12000}];
    Select[v, MemberQ[u, #] &]
Showing 1-1 of 1 results.