OFFSET
0,2
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = a(n-1) + A295848(n) for n > 0.
MAPLE
N:= 100:
V:= Vector(N):
for x from 0 to floor(sqrt(N/3)) do
for y from x to floor(sqrt((N-x^2)/2)) do
for z from y to floor(sqrt(N-x^2-y^2)) do
if igcd(x, y, z) = 1 then
r:= x^2 + y^2 + z^2;
m:= nops({x, y, z});
if m=3 then V[r]:= V[r]+6
elif m=2 then V[r]:= V[r]+3
else V[r]:= V[r]+1
fi
fi
od od od:
0, op(ListTools:-PartialSums(convert(V, list))); # Robert Israel, Nov 30 2017
MATHEMATICA
a[n_] := Sum[Boole[GCD[i, j, k] == 1], {i, 0, Sqrt[n]}, {j, 0, Sqrt[n - i^2]}, {k, 0, Sqrt[n - i^2 - j^2]}];
Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Jul 07 2018, after Andrew Howroyd *)
PROG
(PARI) a(n) = {sum(i=0, sqrtint(n), sum(j=0, sqrtint(n-i^2), sum(k=0, sqrtint(n-i^2-j^2), gcd([i, j, k]) == 1)))} \\ Andrew Howroyd, Dec 12 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Nov 29 2017
STATUS
approved