OFFSET
0,3
LINKS
MAPLE
N:= 100: V:= Vector(N):
for a from 1 to N do
for b from 0 to a-1 do
nmin:= ceil(sqrt(a^2+b^2)-1/2);
if nmin > N then break fi;
V[nmin..N]:= V[nmin..N]+~1
od od:
0, seq(V[i], i=1..N); # Robert Israel, Jul 14 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved