%I #3 Oct 01 2013 17:57:43
%S 1,1,3,2,1,3,2,1,4,3,7,2,5,1,3,9,2,5,4,1,3,7,6,11,2,4,1,7,3,6,5,2,9,4,
%T 7,1,3,5,9,2,8,4,7,1,6,3,10,9,5,8,2,4,12,11,1,3,9,5,8,19,2,12,7,4,10,
%U 6,1,9,3,14,5,12,21,2,7,4,10,6,9,1,3,8,5,11,23,2,10,4,15,14,6,13,1,3,8,5,11
%N Consider numbers x = 4n + 1 which are of the form x = y^2 + z^2 with y <= z; sequence gives smallest value of y.
%e x=45, y=3, z=6. 45-9 = 36 = 6^2
%o (PARI) fourkp1(n) = { forstep(x=1,n,4, for(y=1,floor(sqrt(x)), z=x-y^2; if(issquare(z),print1(y" "); break; ) ) ) }
%Y Cf. A087819.
%K nonn
%O 1,3
%A _Cino Hilliard_, Oct 06 2003