%I #26 Sep 14 2020 17:06:20
%S 8,19,24,28,32,36,75,81,88,96,136,176,224,225,250,328,369,395,432,432,
%T 468,500,512,537,648,701,710,864,980,1000,1078,1089,1125,1216,1225,
%U 1296,1440,1536,1620,1734,1764,1792,1800,1944,1944,2000,2028,2048,2178,2304
%N Numbers whose fifth power is the sum of a positive square and cube: n^5 = x^2 + y^3 with x, y > 0.
%C If a number appears more than once, it can be expressed in as many ways as its occurrence.
%C Zeros are not allowed otherwise 8 would appear twice, since 8^5 = 0^2 + 32^3 = 104^2 + 28^3. - _Michel Marcus_, Feb 25 2013
%C Cubes are not allowed to be negative, otherwise there would be two 8 and two 32, since 8^5 = 192^2 - 16^3 and 32^5=36352^2 - 1088^3. - _Giovanni Resta_, Feb 25 2013
%C Sequence is infinite: if k is a term then k*m^6, for any m, is also a term. For example, {8, 19, 24, 28, 32, 36}*2^6 = {512, 1216, 1536, 1792, 2048, 2304}, or terms 1, 2, 3, 4, 5, and 6 multiplied by 2^6 = terms 23, 34, 39, 42, 48, and 50. - _Zak Seidov_, Jun 27 2013
%e For n=4 a(4)=28 since 28^5 = 3912^2 + 124^3.
%e Multiple representation can happen like 3528^5 = 714208320^2 + 331632^3 = 464679936^2 + 691488^3.
%e 432 is the smallest having multiple representation, since 432^5 = 3732480^2 + 10368^3 = 3359232^2 + 15552^3.
%p f:= proc(z) local z5, y,x,m;
%p z5:= z^5; m:= 0;
%p for y from 1 do
%p x:= z5 - y^3;
%p if x <= 0 then return z$m fi;
%p if issqr(x) then m:= m+1 fi;
%p od
%p end proc:
%p map(f, [$1..2500]); # _Robert Israel_, Nov 13 2016
%t Reap[Do[m5 = m^5; Do[x2 = m5 - y^3; If[x2 <= 0, Break[]]; If[IntegerQ@ Sqrt@ x2, Print[m]; Sow[m]], {y, m5}], {m, 2500}]][[2, 1]] (* _Jean-François Alcover_, Sep 14 2020, after PARI *)
%o (PARI) lista() = for (m=1, 2500, n = m^5; for (i=1, n, v = n - i^3; if (v <= 0, break); if (issquare(v), print1(m, ", ")))) \\ _Michel Marcus_, Feb 25 2013
%Y Cf. A227029.
%K nonn
%O 1,1
%A _Carmine Suriano_, May 20 2010
%E It seems that 1392 is wrong and should be omitted. - _Zak Seidov_, Jun 25 2010
%E Script also found that 1392 is wrong, so 1392 has been removed and 2304 has been appended instead. - _Michel Marcus_, Feb 27 2013
%E 1444 removed and 1734 inserted by _Robert Israel_, Nov 13 2016
|