|
| |
|
|
A146348
|
|
Primes p such that continued fraction of (1+sqrt(p))/2 has period 3.
|
|
37
| |
|
|
17, 37, 61, 101, 197, 257, 317, 401, 461, 557, 577, 677, 773, 1129, 1297, 1429, 1601, 1877, 1901, 2917, 3137, 4357, 4597, 5417, 5477, 6053, 7057, 8101, 8761, 8837, 10733, 11621, 12101, 13457, 13877, 14401, 15277, 15377, 15877, 16333, 16901, 17737, 17957, 18329, 21317, 22501, 23593, 24337, 25601, 28901, 30137, 30977, 32401, 33857, 41453, 41617, 42437, 44101
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,1
|
|
|
COMMENTS
| Primes in A146328. Finite A050952 is subset of this sequence.
|
|
|
LINKS
| Zak Seidov, Table of n, a(n) for n = 1..200
|
|
|
MAPLE
| A146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: isA146348 := proc(n) RETURN(isprime(n) and A146326(n) = 3) ; end: for n from 2 to 4000 do if isA146348(n) then printf("%d, \n", n) ; fi; od: [From R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Sep 06 2009]
|
|
|
MATHEMATICA
| okQ[n_] := Length[ContinuedFraction[(1 + Sqrt[n])/2][[2]]] == 3; Select[Prime[Range[100]], okQ]
|
|
|
CROSSREFS
| Cf. A000290, A050952, A078370, A146326-A146345, A146348-A146360, A146363.
Sequence in context: A059425 A146328 A161549 * A050952 A200865 A093930
Adjacent sequences: A146345 A146346 A146347 * A146349 A146350 A146351
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Artur Jasinski (grafix(AT)csl.pl), Oct 30 2008
|
|
|
EXTENSIONS
| 1019 removed; more terms added - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Sep 06 2009
More terms from Zak Seidov (zakseidov(AT)yahoo.com), Mar 09 2011
|
| |
|
|
