cp's OEIS Frontend

A175375 Number of integer triples (x,y,z) satisfying x^4+y^4+z^4=n, -n <= x,y,z <= n.

%I A175375 #11 May 01 2019 14:00:21
%S A175375 1,6,12,8,0,0,0,0,0,0,0,0,0,0,0,0,6,24,24,0,0,0,0,0,0,0,0,0,0,0,0,0,
%T A175375 12,24,0,0,0,0,0,0,0,0,0,0,0,0,0,0,8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
%U A175375 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,6,24,24,0,0,0,0,0,0,0,0,0,0,0,0,0,24,48,0,0
%N A175375 Number of integer triples (x,y,z) satisfying x^4+y^4+z^4=n, -n <= x,y,z <= n.
%C A175375 4th power variant of A004018 and A175361.
%H A175375 Robert Israel, <a href="/A175375/b175375.txt">Table of n, a(n) for n = 0..10000</a>
%F A175375 G.f.: ( 1 + 2*Sum_{j>0} x^(j^4) )^3.
%e A175375 a(0) = 1 counts (x,y,z) = (0,0,0). a(3) = 8 counts (x,y,z) = (-1,-1,-1), (-1,-1,1), (-1,1,-1), (-1,1,1), (1,-1,-1), (1,-1,1), (1,1,-1) and (1,1,1). a(17) = 24 counts triples where one of x, y and z is 0, one is +-1 and the third +-2.
%p A175375 N:= 200: # to get a(0)..a(N)
%p A175375 A:= Array(0..N):
%p A175375 for i from 0 while i^4 <= N do
%p A175375   if i=0 then ai:= 1 else ai:= 2 fi;
%p A175375   for j from 0 while i^4 + j^4 <= N do
%p A175375     if j=0 then aj:= 1 else aj:= 2 fi;
%p A175375     for k from 0 do
%p A175375       v:= i^4 + j^4 + k^4;
%p A175375       if v > N then break fi;
%p A175375       if k = 0 then ak:= 1 else ak:= 2 fi;
%p A175375       A[v]:= A[v] + ai*aj*ak;
%p A175375 od od od:
%p A175375 seq(A[i],i=0..N); # _Robert Israel_, May 01 2019
%Y A175375 Cf. A004018, A175361.
%K A175375 nonn,look
%O A175375 0,2
%A A175375 _R. J. Mathar_, Apr 24 2010