OFFSET
1,1
COMMENTS
For primes in this sequence see A028871 - {2}.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
a(2) = 7 because continued fraction of (1 + sqrt(7))/2 = 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 4, ... has period (1,1,4,1) length 4.
MAPLE
isA146329 := proc(n) RETURN(A146326(n) = 4) ; end:
for n from 2 to 400 do if isA146329(n) then printf("%d, ", n) ; fi; od: # R. J. Mathar, Sep 06 2009
MATHEMATICA
cf4Q[n_]:=Module[{s=(1+Sqrt[n])/2}, If[IntegerQ[s], 1, Length[ ContinuedFraction[ s][[2]]]]==4]; Select[Range[300], cf4Q] (* Harvey P. Dale, Dec 14 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
39, 68, 150, 155, etc. removed by R. J. Mathar, Sep 06 2009
STATUS
approved