%I #12 Nov 10 2021 07:06:42
%S 0,0,0,0,0,0,2,0,0,0,0,0,3,1,0,0,0,0,2,0,3,2,1,0,0,0,0,1,2,0,4,2,1,1,
%T 0,0,0,0,0,0,1,0,6,5,3,7,1,0,0,0,0,2,2,2,2,0,3,5,2,1,6,1,0,0,0,0,5,0,
%U 4,3,4,0,4,3,2,7,2,1,1,0,0,0,0,0,2,6,0,3,2,0,5,4,6,10,1,1,8,1,0,0,0,0
%N Index of redundancy in period of continued fraction for sqrt(n).
%H Amiram Eldar, <a href="/A028833/b028833.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A003285(n) - A028832(n).
%t a[n_] := If[IntegerQ[Sqrt[n]], 0, Length[p = ContinuedFraction[Sqrt[n]] // Last] - Length[p // Union]]; Table[a[n], {n, 1, 102}] (* _Jean-François Alcover_, Dec 26 2012 *)
%Y Cf. A028832, A003285.
%K nonn,nice
%O 1,7
%A _Olivier Gérard_
|