|
| |
|
|
A033676
|
|
Largest divisor of n <= sqrt(n).
|
|
122
|
|
|
|
1, 1, 1, 2, 1, 2, 1, 2, 3, 2, 1, 3, 1, 2, 3, 4, 1, 3, 1, 4, 3, 2, 1, 4, 5, 2, 3, 4, 1, 5, 1, 4, 3, 2, 5, 6, 1, 2, 3, 5, 1, 6, 1, 4, 5, 2, 1, 6, 7, 5, 3, 4, 1, 6, 5, 7, 3, 2, 1, 6, 1, 2, 7, 8, 5, 6, 1, 4, 3, 7, 1, 8, 1, 2, 5, 4, 7, 6, 1, 8, 9, 2, 1, 7, 5, 2, 3
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,4
|
|
|
COMMENTS
|
a(n) = sqrt(n) is a new record if and only if n is a square. - Zak Seidov, Jul 17 2009
a(n) = A060775(n) unless n is a square, when a(n) = A033677(n) = sqrt(n) is strictly larger than A060775(n). It would be nice to have an efficient algorithm to calculate these terms when n has a large number of divisors, as for example in A060776, A060777 and related problems such as A182987. - M. F. Hasler, Sep 20 2011
a(n) = 1 when n = 1 or n is prime. - Alonso del Arte, Nov 25 2012
a(n) is the smallest central divisor of n. Column 1 of A207375. - Omar E. Pol, Feb 26 2019
|
|
|
REFERENCES
|
G. Tenenbaum, pp. 268 ff, in: R. L. Graham et al., eds., Mathematics of Paul Erdős I.
|
|
|
LINKS
|
T. D. Noe, Table of n, a(n) for n = 1..10000
|
|
|
FORMULA
|
a(n) = n / A033677(n).
a(n) = A161906(n,A038548(n)). - Reinhard Zumkeller, Mar 08 2013
a(n) = A162348(2n-1). - Daniel Forgues, Sep 29 2014
|
|
|
MAPLE
|
A033676 := proc(n) local a, d; a := 0 ; for d in numtheory[divisors](n) do if d^2 <= n then a := max(a, d) ; end if; end do: a; end proc: # R. J. Mathar, Aug 09 2009
|
|
|
MATHEMATICA
|
largestDivisorLEQR[n_Integer] := Module[{dvs = Divisors[n]}, dvs[[Ceiling[Length@dvs/2]]]]; largestDivisorLEQR /@ Range[100] (* Borislav Stanimirov, Mar 28 2010 *)
Table[Last[Select[Divisors[n], #<=Sqrt[n]&]], {n, 100}] (* Harvey P. Dale, Mar 17 2017 *)
|
|
|
PROG
|
(PARI) A033676(n) = {local(d); if(n<2, 1, d=divisors(n); d[(length(d)+1)\2])} \\ Michael B. Porter, Jan 30 2010
(Haskell)
a033676 n = last $ takeWhile (<= a000196 n) $ a027750_row n
-- Reinhard Zumkeller, Jun 04 2012
(Python)
from sympy import divisors
def A033676(n):
d = divisors(n)
return d[(len(d)-1)//2] # Chai Wah Wu, Apr 05 2021
|
|
|
CROSSREFS
|
Cf. A033677 (n/a(n)), A000196 (sqrt), A027750 (list of divisors), A056737 (n/a(n) - a(n)), A219695 (half of this for odd numbers), A207375 (list the central divisor(s)).
Indices of given values: A008578 (1 and the prime numbers: a(n) = 1), A161344 (a(n) = 2), A161345 (a(n) = 3), A161424 (4), A161835 (5), A162526 (6), A162527 (7), A162528 (8), A162529 (9), A162530 (10), A162531 (11), A162532 (12), A282668 (indices of primes).
Sequence in context: A348581 A124044 A059981 * A095165 A046805 A034880
Adjacent sequences: A033673 A033674 A033675 * A033677 A033678 A033679
|
|
|
KEYWORD
|
nonn,easy,nice
|
|
|
AUTHOR
|
N. J. A. Sloane
|
|
|
STATUS
|
approved
|
| |
|
|
