OFFSET
1,1
COMMENTS
For primes in this sequence see A146350.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 41 because continued fraction of (1+sqrt(41))/2 = 3, 1, 2, 2, 1, 5, 1, 2, 2, 1, 5, 1, 2, 2, 1, 5, 1, 2, ... has period (1,2,2,1,5) length 5.
MAPLE
isA146330 := proc(n) RETURN(A146326(n) = 5) ; end:
for n from 2 to 2000 do if isA146330(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]]] == 5 &] (* Amiram Eldar, Mar 31 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
259 and 1026 removed by R. J. Mathar, Sep 06 2009
More terms from Amiram Eldar, Mar 31 2020
STATUS
approved