login
A146356
Primes p such that continued fraction of (1 + sqrt(p))/2 has period 11: primes in A146335.
3
541, 593, 661, 701, 857, 1061, 1109, 1217, 1237, 1709, 1733, 1949, 2333, 2557, 2957, 3229, 3677, 3701, 4373, 5081, 5237, 5309, 6133, 7013, 8693, 9533, 10333, 10853, 12437, 14197, 19213, 20693, 21101, 23173, 29753, 30949, 33797, 36677, 37781, 37993, 41813
OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..2000 (terms 1..250 from Harvey P. Dale)
MAPLE
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
MATHEMATICA
Select[Prime[Range[5000]], Length[ContinuedFraction[(1+Sqrt[#])/2][[2]]] == 11&] (* Harvey P. Dale, Apr 27 2016 *)
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 30 2008
EXTENSIONS
1721 and 6491 removed by R. J. Mathar, Sep 06 2009
More terms from Harvey P. Dale, Apr 27 2016
STATUS
approved