Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A146331
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A146331 Numbers k such that continued fraction of (1+Sqrt[k])/2 has period 6 +0
3
18, 19, 22, 38, 39, 44, 54, 57, 58, 59, 66, 68, 70, 74, 86, 102, 105, 107, 111, 112, 114, 115, 130, 131, 146, 147, 148, 150, 159, 164, 175, 178, 183, 186, 187, 198, 203, 253, 258, 260, 264, 267, 273, 275, 278, 294, 303, 308, 309, 326, 327, 330, 333, 341, 346 (list; graph; listen)
OFFSET

1,1

COMMENT

For primes in this sequence see A146351.

EXAMPLE

a(2) = 19 because continued fraction of (1+Sqrt[19])/2 = 2, 1, 2, 8, 2, 1, 3, 1, 2, 8, 2, 1, 3, 1, 2, 8, 2, 1, 3, 1, 2, 8, 2, 1, 3, 1, 2, 8, 2, 1 ...

has period (1, 2, 8, 2, 1, 3) length 6

MAPLE

A146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: isA146331 := proc(n) RETURN(A146326(n) = 6) ; end: for n from 2 to 380 do if isA146331(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]] == 6, AppendTo[bb, n]], {n, 1, Length[aa]}]; bb (*Artur Jasinski*)

CROSSREFS

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

Sequence in context: A004507 A018823 A025144 this_sequence A151900 A031956 A095393

Adjacent sequences: A146328 A146329 A146330 this_sequence A146332 A146333 A146334

KEYWORD

nonn

AUTHOR

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

EXTENSIONS

Added 39, 68, 150, 203 etc. - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Sep 06 2009

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 December 16 17:18 EST 2009. Contains 170825 sequences.


AT&T Labs Research