OFFSET
0,17
LINKS
R. J. Mathar, Table of n, a(n) for n = 0..1000
EXAMPLE
a(16)=2 because we have 130=11^2+3^2=9^2+7^2. a(2)=0 because although 18=3^2+3^2, these components are not mutually prime.
MAPLE
A106602 := proc(n)
local a, x, y, fourn;
fourn := 8*n+2 ;
a := 0 ;
for x from 1 do
if x^2 >= fourn then
return a;
else
y := fourn-x^2 ;
if issqr(y) then
y := sqrt(y) ;
if y <= x and igcd(x, y) = 1 then
a := a+1 ;
end if;
end if;
end if:
end do:
end proc: # R. J. Mathar, Sep 21 2013
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Colin Mallows, May 10 2005
STATUS
approved
