login
Ordered z such that, for 0<x<y<z, the primitive quadruples (x,y,z,w) satisfy x^3 + y^3 + z^3 = w^3.
7

%I #18 Jul 17 2020 03:49:49

%S 5,8,17,18,21,22,27,33,37,37,40,41,44,49,53,54,57,61,64,65,66,69,69,

%T 70,72,74,75,78,79,79,79,84,85,86,86,87,89,90,92,96,97,97,97,99,101,

%U 102,102,104,105,108,114,116,118,121,122,123,124,124,128,131,136,136,137

%N Ordered z such that, for 0<x<y<z, the primitive quadruples (x,y,z,w) satisfy x^3 + y^3 + z^3 = w^3.

%C For corresponding values w see A096546.

%D Y. Perelman, Solutions to x^3 + y^3 + z^3 = u^3, Mathematics can be Fun, pp. 316-9 Mir Moscow 1985.

%H David A. Corneth, <a href="/A096545/b096545.txt">Table of n, a(n) for n = 1..13798</a> (terms corresponding to z <= 8000)

%H Fred Richman, <a href="http://math.fau.edu/Richman/cubes.htm">Sums of Three Cubes</a>

%e 21 and 22, for instance, are terms because we have: 18^3 + 19^3 + 21^3 = 28^3 and 4^3 + 17^3 + 22^3 = 25^3.

%t s[w_] := Solve[0 < x < y < z && x^3 + y^3 + z^3 == w^3 && GCD[x, y, z, w] == 1, {x, y, z}, Integers];

%t xyzw = Reap[For[w = 1, w <= 200, w++, sw = s[w]; If[sw != {}, Print[{x, y, z, w} /. sw; Sow[{x, y, z, w} /. sw ]]]]][[2, 1]] // Flatten[#, 1]&;

%t Sort[xyzw[[All, 3]]] (* _Jean-François Alcover_, Mar 06 2020 *)

%Y Primitive quadruples (x, y, z, w) = (A095868, A095867, A096545, A096546).

%K nonn

%O 1,1

%A _Lekraj Beedassy_, Jun 25 2004

%E Edited, corrected and extended by _Ray Chandler_, Jun 28 2004