OFFSET
0,1
LINKS
Piotr Miska, Maciej Ulas, On consecutive 1's in continued fractions expansions of square roots of prime numbers, arXiv:1904.03404 [math.NT], 2019. See Table 1 p. 15.
FORMULA
Limit_{n->infinity} (sqrt(a(n)) - floor(sqrt(a(n)))) = A094214. - Daniel Suteu, Apr 09 2019
EXAMPLE
PROG
(PARI) isok(p, n) = {my(c=contfrac(sqrt(p))); for (k=2, n+1, if (c[k] != 1, return (0)); ); return(c[n+2] != 1); }
a(n) = {my(p=2); while (! isok(p, n), p = nextprime(p+1)); p; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Apr 09 2019
EXTENSIONS
a(21)-a(29) from Daniel Suteu, Apr 09 2019
a(0) added by Chai Wah Wu, Apr 09 2019
STATUS
approved