OFFSET
1,3
MAPLE
a:=proc(n)
local x, y, z, count, num;
count:=0;
for x from 0 to n do
for y from 0 to n do
for z from -n to n do
if x^2+y^2+z^2 > n^2 then next; fi;
if x>0 and x*y > z^2 then count:=count+1; fi;
od:
od:
od:
count;
end:
CROSSREFS
KEYWORD
nonn
AUTHOR
W. Edwin Clark, Nov 26 2012
STATUS
approved