OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..2000
MAPLE
A146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: isA146352 := proc(n) RETURN(isprime(n) and A146326(n) = 7) ; end: for n from 2 to 13000 do if isA146352(n) then printf("%d, \n", n) ; fi; od: # R. J. Mathar, Sep 06 2009
MATHEMATICA
Select[Range[2*10^4], PrimeQ[#] && Length[ContinuedFraction[(1+Sqrt[#])/2][[2]]] == 7 &] (* Amiram Eldar, Mar 30 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
607 removed, 797 inserted by R. J. Mathar, Sep 06 2009
More terms from Amiram Eldar, Mar 30 2020
STATUS
approved