cp's OEIS Frontend

Maple

N:= 100: # to get a(1)..a(N)
A:= Array(0..N):
for i from 0 while i^4 <= N do
  if i=0 then ai:= 1 else ai:= 2 fi;
  for j from 0 while i^4 + j^4 <= N do
    if j=0 then aj:= 1 else aj:= 2 fi;
    for k from 0 do
      v:= i^4 + j^4 + k^4;
      if v > N then break fi;
      if k = 0 then ak:= 1 else ak:= 2 fi;
      A[v]:= A[v] + ai*aj*ak;
od od od:
ListTools:-PartialSums(convert(A,list)); # Robert Israel, May 01 2019