OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..2000
MAPLE
N:= 100: # for a(1)..a(N)
V:= Vector(N):
for x from 1 to N/2 do
for y from 0 to x do
v:= ceil(2*sqrt(x^2+y^2));
if v > N then break fi;
if y=x or y=0 then V[v]:= V[v]+4 else V[v]:= V[v]+8 fi
od od:
ListTools:-PartialSums(convert(V, list)); # Robert Israel, Mar 29 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Mar 27 2020
STATUS
approved