OFFSET
0,2
COMMENTS
Equivalently, integers k such that (1+i)*(2+i)*...*(k+i) is not in the same quadrant of the complex plane that (1+i)*(2+i)*...*(k-1+i) is in (if one of these numbers lies on the real or imaginary axis, it is taken to be in the quadrant immediately clockwise from it).
The only time that (1+i)*(2+i)*...*(k+i) lies on the real or imaginary axis is when k = 3, which follows from a result of Cilleruelo (see links). - Nathaniel Johnston, Dec 27 2013
The ratio between successive terms quickly approaches exp(Pi/2), which can be proved using the Taylor series of the arctangent function and the (basic) definition of Euler's constant.
LINKS
Javier Cilleruelo, Squares in (1^2 + 1)*...*(n^2 + 1), Journal of Number Theory 128:8 (2008), pp. 2488-2491.
PROG
(PARI)
{
a=1; s=0; S=Pi/2;
while(1, s+=atan(1/a); if(s>S,
S+=Pi/2; print(a)); a++)
}
CROSSREFS
KEYWORD
nonn,more
AUTHOR
James G. Merickel, Oct 10 2013
EXTENSIONS
a(17) added by James G. Merickel, Oct 14 2013
STATUS
approved