OFFSET
1,1
COMMENTS
For primes in this sequence see A146356.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
a(4) = 661 because continued fraction of (1+sqrt(661))/2 = 13, 2, 1, 4, 2, 8, 8, 2, 4, 1, 2, 25, 2, 1, 4, 2, 8, 8, 2, 4, 1, 2, 25, 2, 1, 4, 2, 8, 8, 2, 4, 1, 2, 25, 2, 1, 4, 2, 8, 8 ... has period (2, 1, 4, 2, 8, 8, 2, 4, 1, 2, 25) length 11.
MAPLE
A146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: isA146335 := proc(n) RETURN(A146326(n) = 11) ; end: for n from 2 to 2000 do if isA146335(n) then printf("%d, ", n) ; fi; od: # R. J. Mathar, Sep 06 2009
MATHEMATICA
Select[Range[10^4], !IntegerQ @ Sqrt[#] && Length[ContinuedFraction[(1 + Sqrt[#])/2][[2]]] == 11 &] (* Amiram Eldar, Mar 31 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
916 removed by R. J. Mathar, Sep 06 2009
More terms from Amiram Eldar, Mar 31 2020
STATUS
approved