A033677 -id:A033677 - cp's OEIS frontend

A248835 a(n) = n + A033677(n), where A033677(n) is the smallest divisor of n >= sqrt(n).

2, 4, 6, 6, 10, 9, 14, 12, 12, 15, 22, 16, 26, 21, 20, 20, 34, 24, 38, 25, 28, 33, 46, 30, 30, 39, 36, 35, 58, 36, 62, 40, 44, 51, 42, 42, 74, 57, 52, 48, 82, 49, 86, 55, 54, 69, 94, 56, 56, 60, 68, 65, 106, 63, 66, 64, 76, 87, 118, 70, 122, 93

Offset: 1

Bob Selcoe, Oct 15 2014

nonn

Conjecture: For n > 1 there is at least one prime in [n, a(n)] exclusive.
a(n) = 2*n when n is prime.
When n = A002620(m), then a(n) = A002620(m+1), i.e., quarter-squares. Oppermann's conjecture states that there is at least one prime in [A002620(m), A002620(m+1)] exclusive.
When n is square, repeated values for a(n) occur at n-1 and n. These values are A002378(sqrt(n)), i.e., oblong numbers.
When n = A002378(m), then a(n) = (m+1)^2.

			When n = 40, the smallest divisor of 40 that is greater than or equal to sqrt(40) is 8 so a(40)=48.

Wikipedia, Oppermann's conjecture

Cf. A002378, A002620, A033677.

a248835[n_Integer] := n + Min[Select[Divisors[n], # >= Sqrt[n] &]]; a248835 /@ Range[120] (* Michael De Vlieger, Nov 10 2014 *)

a(n)=fordiv(n,d,if(d^2>=n,return(n+d))) \\ Charles R Greathouse IV, Oct 21 2014

[n+min([x for x in divisors(n) if x>=sqrt(n)]) for n in [1..100]] # Tom Edgar, Oct 15 2014

A095787 Numbers n such that A095163(n) is different from A033677(n).

Original entry on oeis.org

20, 30, 42, 48, 54, 56, 72, 80, 88, 90, 96, 99, 104, 108, 110, 117, 120, 126, 130, 132, 140, 143, 150, 154, 156, 160, 168, 180, 182, 192, 195, 204, 208, 210, 216, 221, 224, 228, 234, 238, 240, 252, 255, 264, 266, 270, 272, 280, 285, 288, 294, 300, 304, 306

Offset: 1

Klaus Brockhaus, Jun 05 2004

nonn

A095163(n) <= A033677(n).

			A095163(54) = 6, A033677(54) = 9, so 54 is in the sequence.
For n = 12 we have divisors 1, 2, 3, 4, 6, 12; 1 occurs earlier once, 2 occurs earlier twice, 3 occurs earlier 3 times, but 4 occurs earlier only once, hence a(12) = 4.

Cf. A095163, A033677.

{m=310;v=vector(m);for(n=1,m,d=divisors(n);j=1;while(v[d[j]]>=d[j],j++);x=d[j];v[d[j]]=v[d[j]]+1;r=sqrt(n);j=1;while(d[j]

{m=80;v=vector(m);for(n=1,m,d=divisors(n);j=1;while(v[d[j]]>=d[j],j++);print1(d[j],",");v[d[j]]=v[d[j]]+1)}

A101257 Remainder when the least divisor of n greater than or equal to the square root of n (A033677(n)) is divided by the greatest divisor of n less than or equal to the square root of n (A033676(n)).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 2, 0, 0, 0, 0, 1, 1, 1, 0, 2, 0, 1, 0, 3, 0, 1, 0, 0, 2, 1, 2, 0, 0, 1, 1, 3, 0, 1, 0, 3, 4, 1, 0, 2, 0, 0, 2, 1, 0, 3, 1, 1, 1, 1, 0, 4, 0, 1, 2, 0, 3, 5, 0, 1, 2, 3, 0, 1, 0, 1, 0, 3, 4, 1, 0, 2, 0, 1, 0, 5, 2, 1, 2, 3, 0, 1, 6, 3, 1, 1, 4, 4, 0, 0, 2, 0, 0, 5, 0, 5, 1

Offset: 1

Joseph Biberstine (jrbibers(AT)indiana.edu), Dec 17 2004

Given n points, sort them into the most-square rectangular point lattice possible. Now sort the points into square point lattices of dimension equal to the lesser dimension of the earlier rectangle. a(n) is the number of points left over. a(n) is trivially 0 for prime numbers n (the most-square and only rectangular point lattice on a prime number of points is a linear point lattice). a(n) != 0 iff n is a member of A080363.

			a(6)=1 because the least divisor of 6 greater than sqrt(6) is 3, the greater divisor of 6 less than sqrt(6) is 2 and 3 mod 2 = 1

Antti Karttunen, Table of n, a(n) for n = 1..65537
Eric Weisstein's World of Mathematics, Point Lattice.
Eric Weisstein's World of Mathematics, Divisor.

Cf. A033676, A033677, A080363.

num[n_] := If[OddQ[DivisorSigma[0, n]], Sqrt[n], Divisors[n][[DivisorSigma[0, n]/2 + 1]]] den[n_] := If[OddQ[DivisorSigma[0, n]], Sqrt[n], Divisors[n][[DivisorSigma[0, n]/2]]] Table[Mod[num[n], den[n]], {n, 1, 128}]

A033676(n) = if(n<2, 1, my(d=divisors(n)); d[(length(d)+1)\2]); \\ From A033676
A033677(n) = (n/A033676(n));
A101257(n) = (A033677(n)%A033676(n)); \\ Antti Karttunen, Sep 23 2018

Definition corrected by Antti Karttunen, Sep 23 2018

A090680 Number of primes p such that n < p < n+A033677(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 3, 2, 2, 2, 4, 2, 4, 1, 1, 3, 5, 1, 1, 3, 2, 2, 6, 1, 7, 1, 3, 4, 2, 2, 9, 4, 3, 3, 9, 2, 9, 2, 2, 5, 9, 1, 1, 2, 4, 3, 11, 2, 2, 2, 5, 7, 13, 2, 12, 6, 2, 2, 3, 1, 13, 4, 5, 3, 14, 2, 13, 8, 3, 3, 2, 3, 15, 2, 2, 8, 15, 1, 3, 8, 7, 2, 16, 1, 3, 6, 6, 10, 6, 4, 19, 4, 4, 4, 20, 4, 19, 3, 3, 10, 19, 2, 18, 1, 5, 1, 18, 2, 3, 4, 1, 10, 2, 2, 2, 12, 8, 6, 5, 4, 23, 3, 8, 3, 23, 2, 4, 14

Offset: 2

Philippe Deléham, Dec 18 2003

Conjecture: no a(n) is zero. Example: n=15: A033677(15) = 5; 2 primes between 15 and 15+5=20 (17 and 19), so a(15) = 2.

Cf. A033677.

More terms from Benoit Cloitre, Dec 23 2003

A355498 a(n) = A000217(A033676(n)) * A000217(A033677(n)).

Original entry on oeis.org

1, 3, 6, 9, 15, 18, 28, 30, 36, 45, 66, 60, 91, 84, 90, 100, 153, 126, 190, 150, 168, 198, 276, 210, 225, 273, 270, 280, 435, 315, 496, 360, 396, 459, 420, 441, 703, 570, 546, 540, 861, 588, 946, 660, 675, 828, 1128, 756, 784, 825, 918, 910, 1431, 945, 990, 1008, 1140, 1305, 1770

Offset: 1

Steven Lu, Jul 04 2022

nonn

When n squares are arranged in a rectangular grid which is as nearly square as possible, a(n) represents the count of rectangles in the grid. The whole grid itself must be a rectangle too.

			For n=10, the grid as nearly square as possible is 2*5. Thus a(10)=3*15=45 is the number of rectangles in this grid.

Cf. A000217, A033676, A033677.

Table[(# (# + 1)/2) &[
     First[Select[Divisors[n], # >= Sqrt[n] &]]] (# (# + 1)/2) &[
  Last[Select[Divisors[n], # <= Sqrt[n] &]]], {n, 80}]

t(n) = n*(n+1)/2; \\ A000217
largdiv(n) = if(n<2, 1, my(d=divisors(n)); d[(length(d)+1)\2]); \\ A033676
a(n) = my(d=largdiv(n)); t(d)*t(n/d); \\ Michel Marcus, Jul 18 2022

A033676 Largest divisor of n <= sqrt(n).

Original entry on oeis.org

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

Offset: 1

N. J. A. Sloane

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
a(n^4+n^2+1) = n^2-n+1: suppose that n^2-n+k divides n^4+n^2+1 = (n^2-n+k)*(n^2+n-k+2) - (k-1)*(2*n+1-k) for 2 <= k <= 2*n, then (k-1)*(2*n+1-k) >= n^2-n+k, or n^2 - (2*k-1)*n + (k^2-k+1) = (n-k+1/2)^2 + 3/4 < 0, which is impossible. Hence the next smallest divisor of n^4+n^2+1 than n^2-n+1 is at least n^2-n+(2*n+1) = n^2+n+1 > sqrt(n^4+n^2+1). - Jianing Song, Oct 23 2022

G. Tenenbaum, pp. 268 ff, in: R. L. Graham et al., eds., Mathematics of Paul Erdős I.

T. D. Noe, Table of n, a(n) for n = 1..10000

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)).
The strictly inferior case is A060775. Cf. also A140271.
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).

a033676 n = last $ takeWhile (<= a000196 n) $ a027750_row n
-- Reinhard Zumkeller, Jun 04 2012

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

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 *)

A033676(n) = {local(d);if(n<2,1,d=divisors(n);d[(length(d)+1)\2])} \\ Michael B. Porter, Jan 30 2010

from sympy import divisors
def A033676(n):
    d = divisors(n)
    return d[(len(d)-1)//2]  # Chai Wah Wu, Apr 05 2021

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

A060775 The greatest divisor d|n such that d < n/d, with a(1) = 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 3, 1, 2, 3, 2, 1, 3, 1, 4, 3, 2, 1, 4, 1, 2, 3, 4, 1, 5, 1, 4, 3, 2, 5, 4, 1, 2, 3, 5, 1, 6, 1, 4, 5, 2, 1, 6, 1, 5, 3, 4, 1, 6, 5, 7, 3, 2, 1, 6, 1, 2, 7, 4, 5, 6, 1, 4, 3, 7, 1, 8, 1, 2, 5, 4, 7, 6, 1, 8, 3, 2, 1, 7, 5, 2, 3

Offset: 1

Labos Elemer, Apr 26 2001

Also: Largest divisor of n which is less than sqrt(n).
If n is not a square, then a(n) = A033676(n), else a(n) is strictly smaller than A033676(n) = sqrt(n) (except for a(1) = 1).  - M. F. Hasler, Sep 20 2011
Record values occur for n = k * (k+1), for which a(n) = k. - Franklin T. Adams-Watters, May 01 2015
If we define a divisor d|n to be strictly inferior if d < n/d, then strictly inferior divisors are counted by A056924 and listed by A341674. This sequence gives the greatest strictly inferior divisor, which may differ from the lower central divisor A033676. Central divisors are listed by A207375. - Gus Wiseman, Feb 28 2021

			n = 252, D = {1, 2, 3, 4, 6, 7, 9, 12, 14, 18, 21, 28, 36, 42, 63, 84, 126, 252}, 18 divisors, the 9th is 14, so a(252) = 14.
From _Gus Wiseman_, Feb 28 2021: (Start)
The strictly inferior divisors of selected n:
n = 1  2  6  12  20  30  42  56  72  90  110  132  156  182  210  240
    -----------------------------------------------------------------
    {} 1  1  1   1   1   1   1   1   1   1    1    1    1    1    1
          2  2   2   2   2   2   2   2   2    2    2    2    2    2
             3   4   3   3   4   3   3   5    3    3    7    3    3
                     5   6   7   4   5   10   4    4    13   5    4
                                 6   6        6    6         6    5
                                 8   9        11   12        7    6
                                                             10   8
                                                             14   10
                                                                  12
                                                                  15
(End)

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (terms n = 2..1000 from Harry J. Smith)

Cf. A033677, A000196, A000005, A000142, A027423, A055226, A060776, A060777, A002378.
The weakly inferior version is A033676.
Positions of first appearances are A180291.
These are the row-maxima of A341674.
A038548 counts superior (or inferior) divisors.
A056924 counts strictly superior (or strictly inferior) divisors.
A070039 adds up strictly inferior divisors.
A207375 lists central divisors.
A333805 counts strictly inferior odd divisors.
A333806 counts strictly inferior prime divisors.
A341596 counts strictly inferior squarefree divisors.
A341677 counts strictly inferior prime-power divisors.
- Inferior: A063962, A066839, A069288, A161906, A217581, A333749, A333750.
- Superior: A051283, A059172, A063538, A063539, A070038, A161908, A341591.
- Strictly Superior: A048098, A064052, A140271, A238535, A341642, A341673.
Cf. A000203, A001248, A006530, A020639, A112798, A161901.

with(numtheory):
a:= n-> max(select(d-> is(d=1 or dAlois P. Heinz, Jan 29 2018

Table[Part[Divisors[w], Floor[DivisorSigma[0, w]/2]], {w, 1, 256}]
Table[If[n==1,1,Max[Select[Divisors[n],#Gus Wiseman, Feb 28 2021 *)

A060775(n)=if(n>1,divisors(n)[numdiv(n)\2],1) \\ M. F. Hasler, Sep 21 2011

a(n) = max { d: d|n and d < sqrt(n) or d = 1 }, where "|" means "divides". [Corrected by M. F. Hasler, Apr 03 2019]

a(1) = 1 added (to preserve the relation a(n) | n) by Franklin T. Adams-Watters, Jan 27 2018
Edited by M. F. Hasler, Apr 03 2019
Name changed by Gus Wiseman, Feb 28 2021 (was: Lower central (median) divisor of n, with a(1) = 1.)

A063539 Numbers n that are sqrt(n-1)-smooth: largest prime factor of n (=A006530(n)) < sqrt(n).

Original entry on oeis.org

1, 8, 12, 16, 18, 24, 27, 30, 32, 36, 40, 45, 48, 50, 54, 56, 60, 63, 64, 70, 72, 75, 80, 81, 84, 90, 96, 98, 100, 105, 108, 112, 120, 125, 126, 128, 132, 135, 140, 144, 147, 150, 154, 160, 162, 165, 168, 175, 176, 180, 182, 189, 192, 195, 196

Offset: 1

N. J. A. Sloane, Aug 14 2001

nonn

Sometimes (Weisstein) called the "usual numbers" as opposed to what Greene and Knuth define as "unusual numbers" (A063538), which turn out to not be so unusual after all (Greene and Knuth 1990, Finch 2001). - Jonathan Vos Post, Sep 11 2010
If we define a divisor d|n to be superior if d >= n/d, then superior divisors are counted by A038548 and listed by A161908. This sequence lists numbers without a superior prime divisor, which is unique (A341676) when it exists. For example, the set of superior prime divisors of each n starts: {},{2},{3},{2},{5},{3},{7},{},{3},{5},{11},{},{13},{7}. The positions of empty sets give the sequence. - Gus Wiseman, Feb 24 2021
As Jonathan Vos Post's comment suggests, the sqrt(n-1)-smooth numbers are asymptotically less dense than their "unusual" complement. This is part of a larger picture of "typical" relative sizes of a number's prime factors: see, for example, the medians of the n-th smallest prime factors of the positive integers in A281889. - Peter Munn, Mar 03 2021

			a(100) = 360; a(1000) = 3744; a(10000) = 37665; a(100000)=375084;
a(10^6) = 3697669; a(10^7) = 36519633; a(10^8) = 360856296;
a(10^9) = 3571942311; a(10^10) = 35410325861; a(10^11) = 351498917129. - _Giovanni Resta_, Apr 12 2020

Greene, D. H. and Knuth, D. E., Mathematics for the Analysis of Algorithms, 3rd ed. Boston, MA: Birkhäuser, pp. 95-98, 1990.

N. J. A. Sloane, Table of a(n) for n = 1..10622 [Extending and correcting earlier b-files from T. D. Noe and Marius A. Burtea]
M. Beeler, R. W. Gosper and R. Schroeppel, HAKMEM, ITEM 29
Steven Finch, "RE: Unusual Numbers.", Aug 27 2001.
Hugo Pfoertner, Illustration of deviation from linear fit 3.7642*n, (2020).
Hugo Pfoertner, Illustration of relative error of asymptotic approximation, (2020).
Project Euler, Problem 668: Square root smooth numbers
V. Ramaswami, On the number of positive integers less than x and free of prime divisors greater than x^c, Bull. Amer. Math. Soc. 55 (1949), 1122-1127.
Eric W. Weisstein, Rough Number. [From _Jonathan Vos Post_, Sep 11 2010]
Wikipedia, Dickman function

Set difference of A048098 and A001248.
Complement of A063538.
Cf. A006530.
The following are all different versions of sqrt(n)-smooth numbers: A048098, A063539, A064775, A295084, A333535, A333536.
Positions of zeros in A341591.
A001221 counts prime divisors, with sum A001414.
A001222 counts prime-power divisors.
A033677 selects the smallest superior divisor.
A038548 counts superior (or inferior) divisors.
A051283 lists numbers without a superior prime-power divisor.
A056924 counts strictly superior (or strictly inferior) divisors.
A059172 lists numbers without a superior squarefree divisor.
A063962 counts inferior prime divisors.
A116882/A116883 list numbers with/without a superior odd divisor.
A161908 lists superior divisors.
A207375 lists central divisors.
A217581 selects the greatest inferior prime divisor.
A341642 counts strictly superior prime divisors.
A341676 gives unique superior prime divisors, with strict case A341643.
- Inferior: A033676, A066839, A069288, A161906, A333749, A333750.
- Superior: A070038, A072500, A341592, A341593, A341675.
- Strictly inferior: A060775, A070039, A333805, A333806, A341596, A341674.
- Strictly superior: A064052, A140271, A238535, A341594, A341595, A341644, A341645, A341646, A341673.
Cf. A000005, A000203, A020639, A112798, A281889.

[1] cat [m:m in [2..200]| Max(PrimeFactors(m)) lt Sqrt(m) ]; // Marius A. Burtea, May 08 2019

N:= 1000: # to get all terms <= N
Primes:= select(isprime, [2, seq(2*i+1, i=1..floor((N-1)/2))]):
S:= {$1..N} minus {seq(seq(m*p, m = 1 .. min(p, N/p)), p=Primes)}:
sort(convert(S, list)); # Robert Israel, Sep 02 2015

Prepend[Select[Range[192], FactorInteger[#][[-1, 1]] < Sqrt[#] &], 1] (* Ivan Neretin, Sep 02 2015 *)

from math import isqrt
from sympy import primepi
def A063539(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return int(n+primepi(x//(y:=isqrt(x)))+sum(primepi(x//i)-primepi(i) for i in range(1,y)))
    return bisection(f,n,n) # Chai Wah Wu, Oct 05 2024

From Hugo Pfoertner, Apr 02 - Apr 12 2020: (Start)
For small n (e.g. n < 10000) a(n) can apparently be approximated by 3.7642*n.
Asymptotically, the number of sqrt(n)-smooth numbers < x is known to be (1-log(2))*x + O(x/log(x)), see Ramaswami (1949).
n = (1-log(2))*a(n) - 0.59436*a(n)/log(a(n)) is a fitted approximation. (End)
However, it is known that this fit only leads to an increase of accuracy in the range up to a(10^11). The improvement in accuracy suggested by the plot of the relative error for even larger n does not occur. For larger n the behavior of the error term O(x/log(x)) is not known. - Hugo Pfoertner, Nov 12 2023

A161906 Triangle read by rows in which row n lists the divisors of n that are <= sqrt(n).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 3, 1, 2, 1, 1, 2, 3, 1, 1, 2, 1, 3, 1, 2, 4, 1, 1, 2, 3, 1, 1, 2, 4, 1, 3, 1, 2, 1, 1, 2, 3, 4, 1, 5, 1, 2, 1, 3, 1, 2, 4, 1, 1, 2, 3, 5, 1, 1, 2, 4, 1, 3, 1, 2, 1, 5, 1, 2, 3, 4, 6, 1, 1, 2, 1, 3, 1, 2, 4, 5, 1, 1, 2, 3, 6, 1, 1, 2, 4, 1, 3, 5, 1, 2, 1, 1, 2, 3

Offset: 1

Omar E. Pol, Jun 27 2009

If we define a divisor d|n to be inferior if d <= n/d, then inferior divisors are counted by A038548 and listed by this sequence. - Gus Wiseman, Mar 08 2021

			Triangle begins:
   1....... 1;
   2....... 1;
   3....... 1;
   4..... 1,2;
   5....... 1;
   6..... 1,2;
   7....... 1;
   8..... 1,2;
   9..... 1,3;
  10..... 1,2;
  11....... 1;
  12... 1,2,3;
  13....... 1;
  14..... 1,2;
  15..... 1,3;
  16... 1,2,4;

Reinhard Zumkeller, Rows n = 1..1000 of triangle, flattened

Initial terms are A000012.
Final terms are A033676.
Row lengths are A038548 (number of inferior divisors).
Row sums are A066839 (sum of inferior divisors).
The prime terms are counted by A063962.
The odd terms are counted by A069288.
Row products are A072499.
Row LCMs are A072504.
The superior version is A161908.
The squarefree terms are counted by A333749.
The prime-power terms are counted by A333750.
The strictly superior version is A341673.
The strictly inferior version is A341674.
A001221 counts prime divisors, with sum A001414.
A000005 counts divisors, listed by A027750 with sum A000203.
A056924 count strictly superior (or strictly inferior divisors).
A207375 lists central divisors.
- Inferior: A217581.
- Superior: A033677, A051283, A059172, A063538, A063539, A070038, A116882, A116883, A341591, A341592, A341593, A341675, A341676.
- Strictly Inferior: A060775, A070039, A333805, A333806, A341596, A341677.
- Strictly Superior: A048098, A064052, A140271, A238535, A341594, A341595, A341642, A341643, A341644, A341645, A341646.
Cf. A000196, A001055, A001248, A006530, A020639, A050320, A068101, A161901.

a161906 n k = a161906_tabf !! (n-1) !! (k-1)
a161906_row n = a161906_tabf !! (n-1)
a161906_tabf = zipWith (\m ds -> takeWhile ((<= m) . (^ 2)) ds)
                       [1..] a027750_tabf'
-- Reinhard Zumkeller, Jun 24 2015, Mar 08 2013

div[n_] := Select[Divisors[n], # <= Sqrt[n] &]; div /@ Range[48] // Flatten (* Amiram Eldar, Nov 13 2020 *)

row(n) = select(x->(x<=sqrt(n)), divisors(n)); \\ Michel Marcus, Nov 13 2020

More terms from Sean A. Irvine, Nov 29 2010

A063655 Smallest semiperimeter of integral rectangle with area n.

Original entry on oeis.org

2, 3, 4, 4, 6, 5, 8, 6, 6, 7, 12, 7, 14, 9, 8, 8, 18, 9, 20, 9, 10, 13, 24, 10, 10, 15, 12, 11, 30, 11, 32, 12, 14, 19, 12, 12, 38, 21, 16, 13, 42, 13, 44, 15, 14, 25, 48, 14, 14, 15, 20, 17, 54, 15, 16, 15, 22, 31, 60, 16, 62, 33, 16, 16, 18, 17, 68, 21, 26

Offset: 1

Floor van Lamoen, Jul 24 2001

nonn

Similar to A027709, which is minimal perimeter of polyomino of n cells, or equivalently, minimal perimeter of rectangle of area at least n and with integer sides. Present sequence is minimal semiperimeter of rectangle with area exactly n and with integer sides. - Winston C. Yang (winston(AT)cs.wisc.edu), Feb 03 2002
Semiperimeter b+d, d >= b, of squarest (smallest d-b) integral rectangle with area bd = n. That is, b = largest divisor of n <= sqrt(n), d = smallest divisor of n >= sqrt(n). a(n) = n+1 iff n is noncomposite (1 or prime). - Daniel Forgues, Nov 22 2009
From Juhani Heino, Feb 05 2019: (Start)
Basis for any thickness "frames" around the minimal area. Perimeter can be thought as the 0-thick frame, it is obviously 2a(n). Thickness 1 is achieved by laying unit tiles around the area, there are 2(a(n)+2) of them. Thickness 2 comes from the second such layer, now there are 4(a(n)+4) and so on. They all depend only on a(n), so they share this structure:
Every n > 1 is included. (For different thicknesses, every integer that can be derived from these with the respective formula. So, the perimeter has every even n > 2.)
For each square n > 1, a(n) = a(n-1).
a(1), a(2) and a(6) are the only unique values - the others appear multiple times.
(End)
Gives a discrete Uncertainty Principle. A complex function on an abelian group of order n and its Discrete Fourier Transform must have at least a(n) nonzero entries between them. This bound is achieved by the indicator function on a subgroup of size closest to sqrt(n). - Oscar Cunningham, Oct 10 2021
Also two times the median divisor of n, where the median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length). The version for mean instead of median is A057020/A057021. Other doubled medians of multisets are: A360005 (prime indices), A360457 (distinct prime indices), A360458 (distinct prime factors), A360459 (prime factors), A360460 (prime multiplicities), A360555 (0-prepended differences). - Gus Wiseman, Mar 18 2023

			Since 15 = 1*15 = 3*5 and the 3*5 rectangle gives smallest semiperimeter 8, we have a(15)=8.

T. D. Noe, Table of n, a(n) for n = 1..10000
Roy Meshulam, An uncertainty inequality for finite abelian groups, arXiv:math/0312407 [math.CO], 2003.
Roy Meshulam, An uncertainty inequality for finite abelian groups, European Journal of Combinatorics, 27 (2006) 63-67.

Cf. A027709, A033676, A033677, A092510, A162348.
Positions of odd terms are A139710.
Positions of even terms are A139711.
A000005 counts divisors, listed by A027750.
A000975 counts subsets with integer median.
Cf. A003601, A026424, A057020/A057021, A325347, A359893, A359901.

A063655 := proc(n)
    local i,j;
    for i from floor(sqrt(n)) to 1 by -1 do
        j := floor(n/i) ;
        if i*j = n then
            return i+j;
        end if;
    end do:
end proc:
seq(A063655(n), n=1..80); # Winston C. Yang, Feb 03 2002

Table[d = Divisors[n]; len = Length[d]; If[OddQ[len], 2*Sqrt[n], d[[len/2]] + d[[1 + len/2]]], {n, 100}] (* T. D. Noe, Mar 06 2012 *)
Table[2*Median[Divisors[n]],{n,100}] (* Gus Wiseman, Mar 18 2023 *)

A063655(n) = { my(c=1); fordiv(n,d,if((d*d)>=n,if((d*d)==n,return(2*d),return(c+d))); c=d); (0); }; \\ Antti Karttunen, Oct 20 2017

from sympy import divisors
def A063655(n):
    d = divisors(n)
    l = len(d)
    return d[(l-1)//2] + d[l//2] # Chai Wah Wu, Jun 14 2019

a(n) = A033676(n) + A033677(n).
a(n) = A162348(2n-1) + A162348(2n). - Daniel Forgues, Sep 29 2014
a(n) = Min_{d|n} (n/d + d). - Ridouane Oudra, Mar 17 2024

Corrected and extended by Larry Reeves (larryr(AT)acm.org) and Dean Hickerson, Jul 26 2001

cp's OEIS Frontend

A248835 a(n) = n + A033677(n), where A033677(n) is the smallest divisor of n >= sqrt(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Sage

A095787 Numbers n such that A095163(n) is different from A033677(n).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

PARI

A101257 Remainder when the least divisor of n greater than or equal to the square root of n (A033677(n)) is divided by the greatest divisor of n less than or equal to the square root of n (A033676(n)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A090680 Number of primes p such that n < p < n+A033677(n).

Views

Author

Keywords

Comments

Crossrefs

Extensions

A355498 a(n) = A000217(A033676(n)) * A000217(A033677(n)).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

A033676 Largest divisor of n <= sqrt(n).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

A060775 The greatest divisor d|n such that d < n/d, with a(1) = 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula