|
|
A065082
|
|
Numbers k such that (period of the continued fraction for sqrt(k))^3 = k.
|
|
0
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Next term, if it exists, exceeds 6600000. - Sam Handler (sam_5_5_5_0(AT)yahoo.com), Nov 14 2004
|
|
LINKS
|
|
|
MATHEMATICA
|
Do[ If[ Length[ Last[ ContinuedFraction[ Sqrt[n]]]]^3 == n, Print[n]], {n, 1, 10^5} ]
|
|
PROG
|
(PARI) f(n) = {if (issquare(n), 0, my(p=100); while (! isokq(n, p), p+=100); localprec(p); my(cf = contfrac(sqrt(n))); for (k=2, #cf, if (cf[k] == 2*cf[1], return (k-1))); ); } \\ A003285
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|