login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A065082
Numbers k such that (period of the continued fraction for sqrt(k))^3 = k.
0
8, 125, 216, 1728
OFFSET
1,1
COMMENTS
Next term, if it exists, exceeds 6600000. - Sam Handler (sam_5_5_5_0(AT)yahoo.com), Nov 14 2004
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
isok(m) = f(m)^3 == m; \\ Michel Marcus, Jul 12 2021
CROSSREFS
Cf. A003285.
Sequence in context: A024282 A340888 A231638 * A053058 A050803 A061103
KEYWORD
nonn,more
AUTHOR
Naohiro Nomoto, Nov 09 2001
STATUS
approved