A146358 Primes p such that continued fraction of (1 + sqrt(p))/2 has period 13: primes in A333640.

421, 757, 1021, 1097, 1117, 1301, 1553, 1973, 2069, 2237, 2273, 2789, 2861, 3373, 3461, 3517, 3877, 3917, 4133, 4397, 4481, 5521, 5573, 5717, 6221, 6317, 6637, 6997, 7253, 7517, 8741, 9049, 9173, 9437, 10181, 10949, 11597, 11789, 12497, 15473, 15797, 16141, 18353

Offset: 1

Author

Artur Jasinski, Oct 30 2008

Keywords

• nonn

Links

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

Crossrefs

Cf. A000290, A050950-A050969, A078370, A146326-A146345, A333640, A146348-A146360.

Programs

Mathematica

Select[Range[2*10^4], PrimeQ[#] && Length[ContinuedFraction[(1+Sqrt[#])/2][[2]]] == 13 &] (* Amiram Eldar, Mar 30 2020 *)
Select[Prime[Range[2500]],Length[ContinuedFraction[(1+Sqrt[#])/2][[2]]]==13&] (* Harvey P. Dale, Mar 05 2023 *)

Extensions

Definition corrected, 3 terms added. - R. J. Mathar, Nov 08 2008

More terms from Amiram Eldar, Mar 30 2020