OFFSET
1,1
COMMENTS
The period of the continued fraction expansion of sqrt(a(n)) = A003285(a(n)) is 10, so the a(n) are a subset of A020349. The periodic part of the continued fraction of sqrt(a(n)) is (D-d)/(2d),1,1,(D-1)/2,2D/d,(D-1)/2,1,1,(D-d)/(2d),2D. See the Bernstein paper.
a(n) seems to be always congruent 5 (mod 8).
LINKS
Leon Bernstein, Fundamental units and cycles in the period of real quadratic number fields, I. Pacific J. Math 63 (1976): 37-61.
PROG
(PARI) list(n)=for(i=1, n, D=2*i+1; fordiv(D, d, if(d>1&&d<D, print1(D^2+4*d, ", "))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Ralf Stephan, Feb 10 2014
STATUS
approved