A146351 Primes p such that continued fraction of (1 + sqrt(p))/2 has period 6: primes in A146331.

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

Author

Artur Jasinski, Oct 30 2008

Keywords

• nonn

Links

• Amiram Eldar, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A000290, A078370, A146326-A146345, A146348-A146360.

Programs

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 *)

Extensions

797 removed by R. J. Mathar, Sep 06 2009

More terms from Harvey P. Dale, Dec 22 2013