OFFSET
1,3
COMMENTS
Also values of x where p runs through the primes of form 4k+1 and 2*p=x^2+y^2, 0<=x<y. - Colin Barker, Jul 07 2014
FORMULA
MATHEMATICA
pp = Select[ Range[200] // Prime, Mod[#, 4] == 1 &]; f[p_] := y - x /. ToRules[ Reduce[0 <= x <= y && p == x^2 + y^2, {x, y}, Integers]]; A079887 = f /@ pp (* Jean-François Alcover, Jan 15 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Jan 13 2003
STATUS
approved