OFFSET
2,1
MATHEMATICA
Do[Print[Length[Union[Last[ContinuedFraction[((-1+10^n)/9)^(1/2)]]]]], {n, 2, 18}]
PROG
(Python)
from sympy.ntheory.continued_fraction import continued_fraction
from sympy import sqrt
def A096488(n): return len(set(continued_fraction(sqrt((10**n-1)//9))[-1])) # Chai Wah Wu, Mar 30 2021
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Labos Elemer, Jun 25 2004
EXTENSIONS
Name edited by Michel Marcus, Aug 22 2019
a(19)-a(22) from Daniel Suteu, Aug 22 2019
a(23)-a(28) from Chai Wah Wu, Apr 05 2021
STATUS
approved