login
Larger member of a coprime pair (x,y) which solves x^2 + y^2 = z^3 with positive x, y and z.
1

%I #17 Dec 04 2017 03:01:48

%S 11,46,52,117,142,198,236,286,415,488,524,549,621,666,835,873,908,970,

%T 1001,1199,1388,1432,1692,1757,1962,1964,1971,2035,2041,2366,2392,

%U 2630,2655,2681,2702,2815,2826,3195,3421,3544,3664,3715,4048,4070,4097,4356

%N Larger member of a coprime pair (x,y) which solves x^2 + y^2 = z^3 with positive x, y and z.

%C If x and y are coprime, so obviously are also (x,z) and (y,z).

%C The ordered values of the bases of the cubes, z, are a subsequence of (and conjecturally the same as) A008846.

%C For production purposes we advice to use the parametrized representations (see references).

%H R. J. Mathar, <a href="/A282095/b282095.txt">Table of n, a(n) for n = 1..1172</a>

%H Imin Chen, <a href="http://dx.doi.org/10.1090/S0025-5718-07-02083-2">On the equation s^2+y^(2p)=alpha^3</a>, Math. Comp. 77 (262) (2008) 1223-1227.

%H Sander R. Dahmen, <a href="http://arxiv.org/abs/1002.0020">A refined modular approach to the diophantine equation x^2+y^(2n)=z^3</a>, arXiv:1002.0020 [math.NT] (2010).

%F {y: x^2 + y^2 = z^3; gcd(x,y) = 1; 1 <= x <= y; x, y, z in N}

%e 2^2 + 11^2 = 5^3, so 11 is in the sequence.

%e 9^2 + 46^2 = 13^3, so 46 is in the sequence.

%e 47^2 + 52^2 = 17^3, so 52 is in the sequence.

%e 44^2 + 117^2 = 25^2, so 117 is in the sequence.

%p # slow version for demonstration only.

%p isA282095 := proc(y)

%p local x,z3 ;

%p for x from 1 to y do

%p if igcd(x,y) = 1 then

%p z3 := x^2+y^2 ;

%p if isA000578(z3) then

%p return true ;

%p end if;

%p end if;

%p end do:

%p return false ;

%p end proc:

%p for y from 1 do

%p if isA282095(y) then

%p printf("%d,\n",y) ;

%p end if;

%p end do:

%t okQ[y_] := Module[{x, z3}, For[x=1, x<y, x++, If[CoprimeQ[x, y], z3 = x^2 + y^2; If[IntegerQ[z3^(1/3)], Print[y]; Return[True]]]]; Return[False]];

%t Select[Range[5000], okQ] (* _Jean-François Alcover_, Dec 04 2017, after _R. J. Mathar_ *)

%Y Subsequence of A282093. Cf. A099533.

%K nonn

%O 1,1

%A _R. J. Mathar_, Feb 06 2017