OFFSET
1,1
COMMENTS
For primes in this sequence see A146354.
Superset of A146354. - R. J. Mathar, Nov 05 2008
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 73 because continued fraction of (1+sqrt(73))/2 = 4, 1, 3, 2, 1, 1, 2, 3, 1, 7, 1, 3, 2, 1, 1, 2, 3, 1, 7, 1, 3, 2, 1, 1, 2, 3, 1, 7, 1, 3, ... has period (1, 3, 2, 1, 1, 2, 3, 1, 7) length 9 .
MAPLE
isA143577 := proc(k) local c; try c := numtheory[cfrac](1/2+sqrt(k)/2, 'periodic', 'quotients') ; if nops(c[2]) = 9 then RETURN(true) ; else RETURN(false) ; fi; catch: RETURN(false) ; end try; end: for k from 2 to 80000 do if isA143577(k) then printf("%d, ", k) ; fi; od: # R. J. Mathar, Nov 05 2008
MATHEMATICA
Select[Range[1000], !IntegerQ @ Sqrt[#] && Length[ContinuedFraction[(1 + Sqrt[#])/2][[2]]] == 9 &] (* Amiram Eldar, Mar 19 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
Extended by R. J. Mathar, Nov 05 2008
More terms from Amiram Eldar, Mar 19 2020
STATUS
approved