|
%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
|
