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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A146344 Records in A146326. 1
2, 4, 6, 8, 10, 12, 18, 20, 26, 34, 42, 48, 50, 52, 54, 60 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

Table of n, a(n) for n=1..16.

MAPLE

146326 := proc(n) if not issqr(n) then numtheory[cfrac]( (1+sqrt(n))/2, 'periodic', 'quotients') ; nops(%[2]) ; else 0 ; fi; end: read("transforms") ; a26 := [seq(A146326(n), n=1..1000)] ; RECORDS(a26)[1] ; # R. J. Mathar, Sep 06 2009

MATHEMATICA

$MaxExtraPrecision = 300; 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++ ]; AppendTo[aa, m]], {n, 1, 1024}]; Print[aa]; bb = {}; m = 0; Do[If[aa[[k]] > m, AppendTo[bb, aa[[k]]]; m = aa[[k]]], {k, 1, Length[aa]}]; bb (* Artur Jasinski *)

CROSSREFS

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

Sequence in context: A258137 A101814 A034090 * A162763 A113242 A309716

Adjacent sequences:  A146341 A146342 A146343 * A146345 A146346 A146347

KEYWORD

nonn

AUTHOR

Artur Jasinski, Oct 30 2008

EXTENSIONS

42 inserted by R. J. Mathar, Sep 06 2009

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 18 22:26 EDT 2019. Contains 328211 sequences. (Running on oeis4.)