OFFSET
1,2
MATHEMATICA
A[n_] := Sum[If[GCD[a^2 + b^2, n] == a^2 + b^2, 1, 0], {a, 0, n - 1}, {b, 0, n - 1}]; Array[A, 100]
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Mar 23 2014
STATUS
approved