login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A120423 a(n) = maximum value among all k where 1<=k<=n of gcd(k,floor(n/k)). 1
1, 1, 1, 2, 2, 1, 1, 2, 3, 3, 3, 2, 2, 2, 2, 4, 4, 4, 4, 3, 3, 3, 3, 2, 5, 5, 5, 5, 5, 3, 3, 4, 4, 4, 4, 6, 6, 6, 6, 6, 6, 3, 3, 3, 3, 3, 3, 4, 7, 7, 7, 7, 7, 7, 7, 5, 5, 5, 5, 3, 3, 3, 3, 8, 8, 8, 8, 8, 8, 8, 8, 6, 6, 6, 6, 6, 6, 6, 6, 6, 9, 9, 9, 9, 9, 9, 9, 9, 9, 4, 4, 4, 4, 4, 4, 4, 4, 7, 7, 10, 10, 10, 10 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,4
COMMENTS
a(n) is the greatest m such that floor(n/m^2) > n/(m^2+m). - Robert Israel, Jun 04 2019
LINKS
EXAMPLE
For n = 10, we have the pairs {k,floor(n/k)} of {1,10},{2,5},{3,3},{4,2},{5,2},{6,1},{7,1},{8,1},{9,1},{10,1}. The GCD's of these 10 pairs are 1,1,3,2,1,1,1,1,1,1. Of these, 3 is the largest. So a(10) = 3.
MAPLE
a:=n->max(seq(gcd(k, floor(n/k)), k=1..n)): seq(a(n), n=1..112); # Emeric Deutsch, Jul 24 2006
# Alternative:
f:= proc(n) local m, a;
for m from floor(sqrt(n)) by -1 do
a:= floor(n/m^2);
if n < a*(m^2+m) then return m fi
od
end proc:
map(f, [$1..200]); # Robert Israel, Jun 04 2019
MATHEMATICA
Table[Max[Table[GCD[k, Floor[n/k]], {k, 1, n}]], {n, 1, 100}] (* Stefan Steinerberger, Jul 22 2006 *)
CROSSREFS
Sequence in context: A203776 A343559 A242357 * A113137 A220603 A238404
KEYWORD
nonn,look
AUTHOR
Leroy Quet, Jul 11 2006
EXTENSIONS
More terms from Stefan Steinerberger and Emeric Deutsch, Jul 22 2006
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 19 04:58 EDT 2024. Contains 370952 sequences. (Running on oeis4.)