OFFSET
1,1
COMMENTS
Subsequence of A146326 (length of period continued fraction of (1 + sqrt(n))/2).
LINKS
Zak Seidov, Table of n, a(n) for n = 1..1000
MAPLE
A146478 := proc(p) local c; c := numtheory[cfrac](1/2+sqrt(p)/2, 'periodic', 'quotients') ; nops(c[2]) ; end: for n from 1 to 100 do printf("%d, ", A146478(ithprime(n))) ; od: # R. J. Mathar, Nov 05 2008
MATHEMATICA
Table[Length[ContinuedFraction[(1+Sqrt[Prime[n]])/2][[2]]], {n, 100}] (* Zak Seidov, Mar 22 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
a(59) changed from 7 to 9 by R. J. Mathar, Nov 05 2008
STATUS
approved