%I #15 Mar 30 2020 05:44:48
%S 541,593,661,701,857,1061,1109,1217,1237,1709,1733,1949,2333,2557,
%T 2957,3229,3677,3701,4373,5081,5237,5309,6133,7013,8693,9533,10333,
%U 10853,12437,14197,19213,20693,21101,23173,29753,30949,33797,36677,37781,37993,41813
%N Primes p such that continued fraction of (1 + sqrt(p))/2 has period 11: primes in A146335.
%H Amiram Eldar, <a href="/A146356/b146356.txt">Table of n, a(n) for n = 1..2000</a> (terms 1..250 from Harvey P. Dale)
%p A146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic','quotients') ; nops(%[2]) ; else 0 ; fi; end: isA146356 := proc(n) RETURN(isprime(n) and A146326(n) = 11) ; end: for n from 2 to 30000 do if isA146356(n) then printf("%d,\n",n) ; fi; od: # _R. J. Mathar_, Sep 06 2009
%t Select[Prime[Range[5000]],Length[ContinuedFraction[(1+Sqrt[#])/2][[2]]] == 11&] (* _Harvey P. Dale_, Apr 27 2016 *)
%Y Cf. A000290, A050950-A050969, A078370, A146326-A146345, A146348-A146360.
%K nonn
%O 1,1
%A _Artur Jasinski_, Oct 30 2008
%E 1721 and 6491 removed by _R. J. Mathar_, Sep 06 2009
%E More terms from _Harvey P. Dale_, Apr 27 2016