Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A146350
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A146350 Primes p such that continued fraction of (1+Sqrt[p])/2 has period 5 : primes in A146330. +0
2
41, 149, 157, 181, 269, 397, 761, 941, 1013, 2081, 2153, 2477, 2693, 3181, 3221, 3533, 4253, 4409, 5273, 5297, 5741, 6949, 8069, 8501, 8597, 9293, 10301, 10357, 10957, 11321, 12281 (list; graph; listen)
OFFSET

1,1

COMMENT

A050954 is subset of this sequence.

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]] == 5, AppendTo[bb, Prime[n]]], {n, 1, Length[aa]}]; bb (*Artur Jasinski*)

CROSSREFS

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

Sequence in context: A044373 A044754 A146330 this_sequence A050954 A141957 A108016

Adjacent sequences: A146347 A146348 A146349 this_sequence A146351 A146352 A146353

KEYWORD

nonn

AUTHOR

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

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 20:09 EST 2009. Contains 167514 sequences.


AT&T Labs Research