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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A146333 Numbers k such that continued fraction of (1+Sqrt[k])/2 has period 8 2
31, 40, 46, 71, 76, 88, 91, 92, 96, 104, 108, 152, 153, 155, 176, 188, 192, 200, 206, 207, 234, 238, 261, 266, 276, 279, 280, 282, 320, 328, 335, 336, 348, 366, 378, 383, 386, 392, 408, 414, 450, 476, 477, 480, 488, 501, 503, 504, 505, 540, 542, 555, 558, 581 (list; graph; refs; listen; history; internal format)
OFFSET

1,1

COMMENTS

For primes in this sequence see A146353.

EXAMPLE

a(1) = 31 because continued fraction of (1+Sqrt[31])/2 = 3, 3, 1, 1, 10, 1, 1, 3, 5, 3, 1, 1, 10, 1, 1, 3, 5, 3, 1, 1, 10, 1, ...

has period (3, 1, 1, 10, 1, 1, 3, 5) length 8

MAPLE

A146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: isA146333 := proc(n) RETURN(A146326(n) = 8) ; end: for n from 2 to 700 do if isA146333(n) then printf("%d, ", n) ; fi; od: [From R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Sep 06 2009]

MATHEMATICA

s = 10; aa = {}; Do[k = ContinuedFraction[(1 + Sqrt[n])/2, 1000]; If[Length[k] < 190, AppendTo[aa, 0], 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]], 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]], 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]], m++ ]; AppendTo[aa, m]], {n, 1, 500}]; bb = {}; Do[If[aa[[n]] == 8, AppendTo[bb, n]], {n, 1, Length[aa]}]; bb (*Artur Jasinski*)

CROSSREFS

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

Sequence in context: A039378 A043201 A043981 * A109645 A156298 A099180

Adjacent sequences:  A146330 A146331 A146332 * A146334 A146335 A146336

KEYWORD

nonn

AUTHOR

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

EXTENSIONS

155 and 279 etc. added, 311 etc. removed by R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Sep 06 2009

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 16 20:01 EST 2012. Contains 205955 sequences.