login
A146351
Primes p such that continued fraction of (1 + sqrt(p))/2 has period 6: primes in A146331.
3
19, 59, 107, 131, 499, 659, 1627, 1907, 2251, 2467, 3803, 4139, 4283, 5827, 6779, 9539, 10067, 11491, 12619, 13763, 16987, 18587, 18803, 19507, 22003, 23003, 23819, 24859, 28643, 30859, 37507, 40939, 42083, 42299, 43403, 43867, 44563, 52747, 53507, 55339
OFFSET
1,1
LINKS
MAPLE
A146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: isA146351 := proc(n) RETURN(isprime(n) and A146326(n) = 6) ; end: for n from 2 to 13000 do if isA146351(n) then printf("%d, \n", n) ; fi; od: # R. J. Mathar, Sep 06 2009
MATHEMATICA
Select[Prime[Range[10000]], Length[ContinuedFraction[(1+Sqrt[#])/2][[2]]] == 6&] (* Harvey P. Dale, Dec 22 2013 *)
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
797 removed by R. J. Mathar, Sep 06 2009
More terms from Harvey P. Dale, Dec 22 2013
STATUS
approved