login
Differences of two positive 4th powers.
12

%I #14 Sep 28 2018 17:58:32

%S 0,15,65,80,175,240,255,369,544,609,624,671,1040,1105,1215,1280,1295,

%T 1695,1776,2145,2320,2385,2400,2465,2800,3439,3471,3840,4015,4080,

%U 4095,4160,4641,5265,5904,5936,6095,6305,6480,6545,6560,7599,7825,8080,8704

%N Differences of two positive 4th powers.

%C If n belongs to this sequence then so does n*m^4 for any positive integer m. Primitive elements (i.e., not of the form n*m^4 for m>1) are listed in A147858.

%C There is no square in this sequence except 0. - _Altug Alkan_, Apr 08 2016

%H Robert Israel, <a href="/A147857/b147857.txt">Table of n, a(n) for n = 1..10000</a>

%p N:= 10^4: # to get all terms <= N

%p Res:= {0}:

%p for a from 1 to floor(sqrt(N-2)) do

%p if a^4 > N then bmin:= ceil((a^4-N)^(1/4)) else bmin:= 1 fi;

%p Res:= Res union {seq(a^4-b^4, b=bmin..a-1)}

%p od:

%p sort(convert(Res,list)); # _Robert Israel_, Sep 28 2018

%t nn = 50; Take[Union @@ Map[Differences, Union@ Map[Sort@ # &, Tuples[Range[Ceiling[nn/3]], {2}]]^4], nn] (* _Michael De Vlieger_, Apr 09 2016 *)

%Y Cf. A147854, A147856, A147858

%K nonn

%O 1,2

%A _Max Alekseyev_, Nov 15 2008, Nov 19 2008

%E Offset changed by _Robert Israel_, Sep 28 2018