

A079886


Values of x+y where p runs through the primes of form 4k+1 and p=x^2+y^2, 0<=x<=y.


4



3, 5, 5, 7, 7, 9, 9, 11, 11, 13, 13, 11, 13, 15, 15, 17, 17, 15, 19, 19, 15, 17, 21, 19, 17, 23, 23, 21, 19, 25, 25, 25, 23, 25, 25, 27, 25, 21, 23, 29, 29, 27, 25, 29, 27, 31, 31, 33, 33, 25, 31, 29, 35, 35, 29, 35, 31, 35, 27, 31, 37, 29, 35, 39, 37, 39, 37, 33, 39, 37, 41, 33
(list;
graph;
refs;
listen;
history;
text;
internal format)



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



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:


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 (* JeanFrançois Alcover, Jan 15 2015 *)


CROSSREFS



KEYWORD

nonn


AUTHOR



STATUS

approved



