login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A146355 Primes p such that continued fraction of (1+Sqrt[p])/2 has period 10 : primes in A146335. 1
43, 67, 563, 827, 1787, 1811, 2099, 2459, 5107, 7643, 8363, 9323, 9371, 9467 (list; graph; refs; listen; history; internal format)
OFFSET

1,1

MATHEMATICA

$MaxExtraPrecision = 4000; s = 10; aa = {}; Do[k = ContinuedFraction[(1 + Sqrt[Prime[n]])/2, 3000]; m = 1; While[k[[s]] != k[[s + m]] || k[[s + m]] != k[[s + 2 m]] || k[[s + 2 m]] != k[[s + 3 m]] || k[[s + 3 m]] != k[[s + 4 m]] || k[[s + 4 m]] != k[[s + 5 m]], m++ ]; s = s + 1; While[k[[s]] != k[[s + m]] || k[[s + m]] != k[[s + 2 m]] || k[[s + 2 m]] != k[[s + 3 m]] || k[[s + 3 m]] != k[[s + 4 m]] || k[[s + 4 m]] != k[[s + 5 m]], m++ ]; s = s + 1; While[k[[s]] != k[[s + m]] || k[[s + m]] != k[[s + 2 m]] || k[[s + 2 m]] != k[[s + 3 m]] || k[[s + 3 m]] != k[[s + 4 m]] || k[[s + 4 m]] != k[[s + 5 m]]; AppendTo[aa, m]], {n, 1, 1495}]; bb = {}; Do[If[aa[[n]] == 10, AppendTo[bb, Prime[n]]], {n, 1, Length[aa]}]; bb (*Artur Jasinski*)

CROSSREFS

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

Sequence in context: A033229 A139875 A174812 * A175730 A087765 A141971

Adjacent sequences:  A146352 A146353 A146354 * A146356 A146357 A146358

KEYWORD

nonn

AUTHOR

Artur Jasinski (grafix(AT)csl.pl), Oct 30 2008

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 17 11:46 EST 2012. Contains 206011 sequences.