OFFSET
0,4
LINKS
MAPLE
N:= 100: A:= Array(0..N):
for a from 1 to N do
for b from 0 to a-1 do
r:= ceil(sqrt(a^2+b^2));
if r > N then break fi;
A[r]:= A[r]+1;
od od:
seq(A[i], i=0..N); # Robert Israel, Dec 18 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved