|
| |
|
|
A096491
|
|
a(n) = sqrt(n) of n if n is a perfect square, otherwise a(n) = largest term in period of continued fraction expansion of square root of n.
|
|
6
| |
|
|
1, 2, 2, 2, 4, 4, 4, 4, 3, 6, 6, 6, 6, 6, 6, 4, 8, 8, 8, 8, 8, 8, 8, 8, 5, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 6, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 7, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 8, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,2
|
|
|
EXAMPLE
| n=127: the period={3,1,2,2,7,11,7,2,2,1,3,22}, max=a[127]=22;
|
|
|
MAPLE
| (Maple code from R. J. Mathar, Mar 18 2010)
A096491 := proc(n)
if issqr(n) then
sqrt(n) ;
else
numtheory[cfrac](sqrt(n), 'periodic', 'quotients') ;
%[2] ;
max(op(%)) ;
end if;
end proc:
|
|
|
MATHEMATICA
| u=1; Do[s=Max[Last[ContinuedFraction[n^(1/2)]]]; tc[[u]]=s; u=u+1, {n, 1, m}]
|
|
|
CROSSREFS
| Cf. A003285, A013646.
Sequence in context: A130872 A087627 A195051 * A106160 A007614 A113402
Adjacent sequences: A096488 A096489 A096490 * A096492 A096493 A096494
|
|
|
KEYWORD
| cofr,nonn
|
|
|
AUTHOR
| Labos E. (labos(AT)ana.sote.hu), Jun 29 2004
|
|
|
EXTENSIONS
| Definition revised by N. J. A. Sloane, Mar 18 2010
|
| |
|
|
