|
|
EXAMPLE
| a(1) = 521 because continued fraction of (1+Sqrt[521])/2 = 11, 1, 10, 2, 5, 4, 2, 1, 1, 1, 1, 2, 4, 5, 2, 10, 1, 21, 1, 10, 2, 5, 4, 2, 1, 1, 1, 1, 2, 4, 5, 2, 10, 1, 21, 1, 10, 2, 5,...
has period (1, 10, 2, 5, 4, 2, 1, 1, 1, 1, 2, 4, 5, 2, 10, 1, 21) length 17
|
|
|
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++ ]; AppendTo[aa, m]], {n, 1, 500}]; bb = {}; Do[If[aa[[n]] == 17, AppendTo[bb, n]], {n, 1, Length[aa]}]; bb (*Artur Jasinski*)
|
