|
Search: id:A146345
|
|
| |
|
| 2, 6, 18, 31, 43, 94, 106, 151, 211, 331, 394, 526, 694, 751, 886, 919, 1114, 1324
(list; graph; listen)
|
|
|
OFFSET
|
1,1
|
|
|
MAPLE
|
A146326 := 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..1400)] ; RECORDS(a26)[2] ; [From R. J. Mathar (mathar(AT)strw.leidenuniv.nl), 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, k]; m = aa[[k]]], {k, 1, Length[aa]}]; bb (*Artur Jasinski*)
|
|
CROSSREFS
|
A000290, A078370, A146326-A146345, A146348-A146360, A146363.
Sequence in context: A032649 A066286 A034881 this_sequence A064842 A101695 A014741
Adjacent sequences: A146342 A146343 A146344 this_sequence A146346 A146347 A146348
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Artur Jasinski (grafix(AT)csl.pl), Oct 30 2008
|
|
EXTENSIONS
|
19 replaced by 18. 394 inserted. 4 more terms added - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Sep 06 2009
|
|
|
Search completed in 0.002 seconds
|
