OFFSET
1,1
COMMENTS
Also values of y 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
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 100: # to get values corresponding to primes <= 4*N+1
P:= select(isprime, [seq(4*i+1, i=0..N)]):
F:= proc(p) local f; f:= GaussInt:-GIfactors(p)[2][1][1]; abs(Re(f))+abs(Im(f)) end proc:
map(F, P); # Robert Israel, Jul 07 2014
MATHEMATICA
pp = Select[ Range[200] // Prime, Mod[#, 4] == 1 &]; f[p_] := x + y /. ToRules[ Reduce[0 <= x <= y && p == x^2 + y^2, {x, y}, Integers]]; A079886 = f /@ pp (* Jean-François Alcover, Jan 15 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Jan 13 2003
STATUS
approved