A161908 -id:A161908 - cp's OEIS frontend

A038548 Number of divisors of n that are at most sqrt(n).

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

Offset: 1

Tom Verhoeff, N. J. A. Sloane

Number of ways to arrange n identical objects in a rectangle, modulo rotation.
Number of unordered solutions of x*y = n. - Colin Mallows, Jan 26 2002
Number of ways to write n-1 as n-1 = x*y + x + y, 0 <= x <= y <= n. - Benoit Cloitre, Jun 23 2002
Also number of values for x where x+2n and x-2n are both squares (e.g., if n=9, then 18+18 and 18-18 are both squares, as are 82+18 and 82-18 so a(9)=2); this is because a(n) is the number of solutions to n=k(k+r) in which case if x=r^2+2n then x+2n=(r+2k)^2 and x-2n=r^2 (cf. A061408). - Henry Bottomley, May 03 2001
Also number of sums of sequences of consecutive odd numbers or consecutive even numbers including sequences of length 1 (e.g., 12 = 5+7 or 2+4+6 or 12 so a(12)=3). - Naohiro Nomoto, Feb 26 2002
Number of partitions whose consecutive parts differ by exactly two.
a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24=2^3*3 and 375=3*5^3 both have prime signature (3,1). - Christian G. Bower, Jun 06 2005
Also number of partitions of n such that if k is the largest part, then each of the parts 1,2,...,k-1 occurs exactly twice. Example: a(12)=3 because we have [3,3,2,2,1,1],[2,2,2,2,2,1,1] and [1,1,1,1,1,1,1,1,1,1,1,1]. - Emeric Deutsch, Mar 07 2006
a(n) is also the number of nonnegative integer solutions of the Diophantine equation 4*x^2 - y^2 = 16*n. For example, a(24)=4 because there are 4 solutions: (x,y) = (10,4), (11,10), (14,20), (25,46). - N-E. Fahssi, Feb 27 2008
a(n) is the number of even divisors of 2*n that are <= sqrt(2*n). - Joerg Arndt, Mar 04 2010
First differences of A094820. - John W. Layman, Feb 21 2012
a(n) = #{k: A027750(n,k) <= A000196(n)}; a(A008578(n)) = 1; a(A002808(n)) > 1. - Reinhard Zumkeller, Dec 26 2012
Row lengths of the tables in A161906 and A161908. - Reinhard Zumkeller, Mar 08 2013
Number of positive integers in the sequence defined by x_0 = n, x_(k+1) = (k+1)*(x_k-2)/(k+2) or equivalently by x_k = n/(k+1) - k. - Luc Rousseau, Mar 03 2018
Expanding the first comment: Number of rectangles with area n and integer side lengths, modulo rotation. Also number of 2D grids of n congruent squares, in a rectangle, modulo rotation (cf. A000005 for rectangles instead of squares; cf. A034836 for the 3D case). - Manfred Boergens, Jun 08 2021
Number of divisors of n that have an even number of prime divisors (counted with multiplicity), or in other words, number of terms of A028260 that divide n. - Antti Karttunen, Apr 17 2022

			a(4) = 2 since 4 = 2 * 2 = 4 * 1. Also A034178(4*4) = 2 since 16 = 4^2 - 0^2 = 5^2 - 3^2. - _Michael Somos_, May 11 2011
x + x^2 + x^3 + 2*x^4 + x^5 + 2*x^6 + x^7 + 2*x^8 + 2*x^9 + 2*x^10 + x^11 + ...

George E. Andrews and Kimmo Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, page 18, exer. 21, 22.

T. D. Noe, Table of n, a(n) for n = 1..10000
Cristina Ballantine and Mircea Merca, New convolutions for the number of divisors, Journal of Number Theory, 2016, vol. 170, pp. 17-34.
Christopher Briggs, Y. Hirano, and H. Tsutsui, Positive Solutions to Some Systems of Diophantine Equations, Journal of Integer Sequences, Vol. 19 (2016), Article 16.8.4.
S.-H. Cha, E. G. DuCasse, and L. V. Quintas, Graph invariants based on the divides relation and ordered by prime signatures, arXiv:1405.5283 [math.NT], 2014, (2.27).
Madeline Locus Dawsey, Matthew Just and Robert Schneider, A "supernormal" partition statistic, arXiv:2107.14284 [math.NT], 2021. See Table 2 p. 21.
T. Verhoeff, Rectangular and Trapezoidal Arrangements, J. Integer Sequences, Vol. 2 (1999), Article 99.1.6.
Index entries for sequences computed from exponents in factorization of n.

Different from A068108. Records give A038549, A004778, A086921.
Cf. A000005, A001620, A028260, A056924, A072670, A094820, A161841, A108504.
Cf. A066839, A033676, row sums of A303300.
Inverse Möbius transform of A065043.
Cf. A244664 (Dgf at s=2), A244665 (Dgf at s=3).

a038548 n = length $ takeWhile (<= a000196 n) $ a027750_row n
-- Reinhard Zumkeller, Dec 26 2012
(C#)
int A038548(int n)  {
    System.Numerics.BigInteger erg = 0, i;
    for (i = 1; i * i <= n; i++)
        if (n % i == 0) erg++;
    return (int)erg;
} // Frank Hollstein, Jan 08 2023

with(numtheory): A038548 := n->ceil(sigma[0](n)/2);

Table[ Floor[ (DivisorSigma[0, n] + 1)/2], {n, 105}] (* Robert G. Wilson v, Mar 02 2009 *)
Table[Count[Divisors[n],?(#<=Sqrt[n]&)],{n,110}] (* _Harvey P. Dale, Jul 10 2021 *)
Table[Sum[If[n > k*(k-1), 1, 0], {k, Divisors[n]}], {n, 1, 100}] (* Vaclav Kotesovec, Oct 22 2024 *)

{a(n) = if( n<1, 0, sumdiv(n, d, d*d <= n))} /* Michael Somos, Jan 25 2005 */

a(n)=ceil(numdiv(n)/2) \\ Charles R Greathouse IV, Sep 28 2012

from sympy import divisor_count
def A038548(n): return divisor_count(n)+1>>1 # Chai Wah Wu, Dec 19 2023

a(n) = ceiling(d(n)/2), where d(n) = number of divisors of n (A000005).
a(2k) = A034178(2k) + A001227(k). a(2k+1) = A034178(2k+1). - Naohiro Nomoto, Feb 26 2002
G.f.: Sum_{k>=1} x^(k^2)/(1-x^k). - Jon Perry, Sep 10 2004
Dirichlet g.f.: (zeta(s)^2 + zeta(2*s))/2. - Christian G. Bower, Jun 06 2005 [corrected by Vaclav Kotesovec, Aug 19 2019]
a(n) = (A000005(n) + A010052(n))/2. - Omar E. Pol, Jun 23 2009
a(n) = A034178(4*n). - Michael Somos, May 11 2011
2*a(n) = A161841(n). - R. J. Mathar, Mar 07 2021
a(n) = A000005(n) - A056924(n) = A056924(n) + A010052(n) = Sum_{d|n} A065043(d). - Antti Karttunen, Apr 17 2022
Sum_{k=1..n} a(k) ~ n*log(n)/2 + (gamma - 1/2)*n, where gamma is Euler's constant (A001620). - Amiram Eldar, Nov 27 2022

A033677 Smallest divisor of n >= sqrt(n).

Original entry on oeis.org

1, 2, 3, 2, 5, 3, 7, 4, 3, 5, 11, 4, 13, 7, 5, 4, 17, 6, 19, 5, 7, 11, 23, 6, 5, 13, 9, 7, 29, 6, 31, 8, 11, 17, 7, 6, 37, 19, 13, 8, 41, 7, 43, 11, 9, 23, 47, 8, 7, 10, 17, 13, 53, 9, 11, 8, 19, 29, 59, 10, 61, 31, 9, 8, 13, 11, 67, 17, 23, 10, 71, 9, 73, 37, 15, 19, 11, 13, 79, 10

Offset: 1

N. J. A. Sloane

a(n) is the smallest k such that n appears in the k X k multiplication table and A027424(k) is the number of n with a(n) <= k.
a(n) is the largest central divisor of n. Right border of A207375. - Omar E. Pol, Feb 26 2019
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 selects the smallest superior divisor of n. - Gus Wiseman, Feb 19 2021
a(p) = p for p a prime or 1, these are also the record high points in this sequence. - Charles Kusniec, Aug 26 2022
a(n^4+n^2+1) = n^2+n+1 (see A033676). - Jianing Song, Oct 23 2022

			From _Gus Wiseman_, Feb 19 2021: (Start)
The divisors of 36 are {1,2,3,4,6,9,12,18,36}. Of these {1,2,3,4,6} are inferior and {6,9,12,18,36} are superior, so a(36) = 6.
The divisors of 40 are {1,2,4,5,8,10,20,40}. Of these {1,2,4,5} are inferior and {8,10,20,40} are superior, so a(40) = 8.
(End)

G. Tenenbaum, pp. 268ff of 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. A027424, A056737, A060775, A219695.
The lower central divisor is A033676.
The strictly superior case is A140271.
Leftmost column of A161908 (superior divisors).
Rightmost column of A207375 (central divisors).
A038548 counts superior (or inferior) divisors.
A056924 counts strictly superior (or strictly inferior) divisors.
A063538/A063539 list numbers with/without a superior prime divisor.
A070038 adds up superior divisors.
A341676 selects the unique superior prime divisor.
- Inferior: A063962, A066839, A069288, A161906, A217581, A333749, A333750.
- Superior: A059172, A116882, A116883, A341591, A341592, A341593, A341675.
- Strictly Inferior: A070039, A333805, A333806, A341596, A341674, A341677.
- Strictly Superior: A048098, A064052, A238535, A341594, A341595, A341642, A341643, A341644, A341645, A341646, A341673.
Cf. A000005, A000203, A001248, A001221, A001222, A020639, A051283.

a033677 n = head $
   dropWhile ((< n) . (^ 2)) [d | d <- [1..n], mod n d == 0]
-- Reinhard Zumkeller, Oct 20 2011

A033677 := proc(n)
    n/A033676(n) ;
end proc:

Table[Select[Divisors[n], # >= Sqrt[n] &, 1] // First, {n, 80}]  (* Jean-François Alcover, Apr 01 2011 *)

A033677(n) = {local(d); d=divisors(n); d[length(d)\2+1]} \\ Michael B. Porter, Feb 26 2010

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

a(n) = n/A033676(n).

a(n) = A162348(2n). - 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.)

A005245 The (Mahler-Popken) complexity of n: minimal number of 1's required to build n using + and *.

Original entry on oeis.org

1, 2, 3, 4, 5, 5, 6, 6, 6, 7, 8, 7, 8, 8, 8, 8, 9, 8, 9, 9, 9, 10, 11, 9, 10, 10, 9, 10, 11, 10, 11, 10, 11, 11, 11, 10, 11, 11, 11, 11, 12, 11, 12, 12, 11, 12, 13, 11, 12, 12, 12, 12, 13, 11, 12, 12, 12, 13, 14, 12, 13, 13, 12, 12, 13, 13, 14, 13, 14, 13, 14, 12, 13, 13, 13, 13, 14, 13, 14

Offset: 1

N. J. A. Sloane

The complexity of an integer n is the least number of 1's needed to represent it using only additions, multiplications and parentheses. This does not allow juxtaposition of 1's to form larger integers, so for example, 2 = 1+1 has complexity 2, but 11 does not ("pasting together" two 1's is not an allowed operation).
The complexity of a number has been defined in several different ways by different authors. See the Index to the OEIS for other definitions.
Guy asks if a(p) = a(p-1) + 1 for prime p. Martin Fuller found the least counterexample p = 353942783 in 2008, see Fuller link. - Charles R Greathouse IV, Oct 04 2012
It appears that this sequence is lower than A348262 {1,+,^} only a finite number of times. - Gordon Hamilton and Brad Ballinger, May 23 2022
The second Altman links proves that {a(n) - 3*log_3(n)} is a well-ordered subset of the reals whose intersection with [0,k) has order type omega^k for each positive integer k, so this set itself has order type omega^omega. - Jianing Song, Apr 13 2024

			From _Lekraj Beedassy_, Jul 04 2009: (Start)
   n.........minimal expression........ a(n) = number of 1's
   1..................1...................1
   2.................1+1..................2
   3................1+1+1.................3
   4.............(1+1)*(1+1)..............4
   5............(1+1)*(1+1)+1.............5
   6............(1+1)*(1+1+1).............5
   7...........(1+1)*(1+1+1)+1............6
   8..........(1+1)*(1+1)*(1+1)...........6
   9...........(1+1+1)*(1+1+1)............6
  10..........(1+1+1)*(1+1+1)+1...........7
  11.........(1+1+1)*(1+1+1)+1+1..........8
  12.........(1+1)*(1+1)*(1+1+1)..........7
  13........(1+1)*(1+1)*(1+1+1)+1.........8
  14.......{(1+1)*(1+1+1)+1}*(1+1)........8
  15.......{(1+1)*(1+1)+1}*(1+1+1)........8
  16.......(1+1)*(1+1)*(1+1)*(1+1)........8
  17......(1+1)*(1+1)*(1+1)*(1+1)+1.......9
  18........(1+1)*(1+1+1)*(1+1+1).........8
  19.......(1+1)*(1+1+1)*(1+1+1)+1........9
  20......{(1+1+1)*(1+1+1)+1}*(1+1).......9
  21......{(1+1)*(1+1+1)+1}*(1+1+1).......9
  22.....{(1+1)*(1+1+1)+1}*(1+1+1)+1.....10
  23....{(1+1)*(1+1+1)+1}*(1+1+1)+1+1....11
  24......(1+1)*(1+1)*(1+1)*(1+1+1).......9
  25.....(1+1)*(1+1)*(1+1)*(1+1+1)+1.....10
  26....{(1+1)*(1+1)*(1+1+1)+1}*(1+1)....10
  27.......(1+1+1)*(1+1+1)*(1+1+1)........9
  28......(1+1+1)*(1+1+1)*(1+1+1)+1......10
  29.....(1+1+1)*(1+1+1)*(1+1+1)+1+1.....11
  30.....{(1+1+1)*(1+1+1)+1}*(1+1+1).....10
  31....{(1+1+1)*(1+1+1)+1}*(1+1+1)+1....11
  32....(1+1)*(1+1)*(1+1)*(1+1)*(1+1)....10
  33...(1+1)*(1+1)*(1+1)*(1+1)*(1+1)+1...11
  34..{(1+1)*(1+1)*(1+1)*(1+1)+1}*(1+1)..11
  .........................................
(End)

M. Criton, "Les uns de Germain", Problem No. 4, pp. 13 ; 68 in '7 x 7 Enigmes Et Défis Mathématiques pour tous', vol. 25, Editions POLE, Paris 2001.
R. K. Guy, Unsolved Problems in Number Theory, Sect. F26.
K. Mahler and J. Popken, Over een Maximumprobleem uit de Rekenkunde (Dutch: On a maximum problem in arithmetic), Nieuw Arch. Wiskunde, (3) 1 (1953), pp. 1-15.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

M. F. Hasler, Table of n, a(n) for n = 1..10000
Harry Altman, Highest few sums and products of ones, Preprint, undated but before 2013.
Harry Altman, Integer Complexity and Well-Ordering, arXiv:1310.2894v1 [math.NT] , Oct 10, 2013.
Harry Altman and Joshua Zelinsky, Numbers with integer complexity close to the lower bound, arXiv:1207.4841 [math.NT], 2012.
J. Arias de Reyna, Complejidad de los números naturales, Gaceta R. Soc. Mat. Esp., 3 (2000), 230-250.
J. Arias de Reyna, Complejidad de los números naturales, Gaceta R. Soc. Mat. Esp., 3 (2000), 230-250. [Cached copy, with permission]
J. Arias de Reyna, Complexity of natural numbers, arXiv:2111.03345 [math.NT], 2021.
J. Arias de Reyna and J. van de Lune, The question "How many 1's are needed?" revisited, arXiv preprint arXiv:1404.1850 [math.NT], 2014.
J. Arias de Reyna and J. van de Lune, Algorithms for determining integer complexity, arXiv preprint arXiv:1404.2183 [math.NT], 2014.
Juris Cernenoks, Janis Iraids, Martins Opmanis, Rihards Opmanis and Karlis Podnieks, Integer Complexity: Experimental and Analytical Results II, arxiv:1409.0446 [math.NT], 2014.
Katherine Cordwell, Alyssa Epstein, Anand Hemmady, Steven J. Miller, Eyvindur A. Palsson, Aaditya Sharma, Stefan Steinerberger and Yen Nhi Truong Vu, On algorithms to calculate integer complexity, arXiv:1706.08424 [math.NT], 2017.
Martin N. Fuller, C program
R. K. Guy, Some suspiciously simple sequences, Amer. Math. Monthly 93 (1986), 186-190; 94 (1987), 965; 96 (1989), 905.
Qizheng He, Improved Algorithms for Integer Complexity, arXiv:2308.10301 [cs.DS], 2023.
Janis Iraids, Online calculator of a(n) for n < 10^12
Jānis Iraids, Kaspars Balodis, Juris Čerņenoks, Mārtiņš Opmanis, Rihards Opmanis, and Kārlis Podnieks, Integer complexity: experimental and analytical results, arXiv:1203.6462 [math.NT], 2012.
Daniel A. Rawsthorne, How many 1's are needed?, Fibonacci Quarterly 27 (1989), pp. 14-17.
Srinivas, Vivek V. and Shankar, B. R., Integer Complexity: Breaking the Theta(n^2) barrier, World Academy of Science 41 (2008), pp. 690-691.
Stefan Steinerberger, A Short Note on Integer Complexity, Contributions to Discrete Mathematics, Vol. 9.1 (2014).
Stefan Steinerberger, A short note on integer complexity, Contributions to Discrete Mathematics, Vol. 9.1 (2014). [Cached copy, with permission.]
Joshua Zelinsky, Upper Bounds on Integer Complexity, arXiv preprint (2022). arXiv:2211.02995 [math.NT]
Venecia Wang, A counterexample to the prime conjecture of expressing numbers using just ones, YouTube video, 2012.
Venecia Wang, A counterexample to the prime conjecture of expressing numbers using just ones, Journal of Number Theory, Volume 133, Issue 2, February 2013, Pages 391-397.
Eric Weisstein's World of Mathematics, Integer Complexity
Dimitri Zucker, The Most Underrated Concept in Number Theory, Combo Class Youtube video that mentions this sequence.
Index to sequences related to the complexity of n

Cf. A025280 (variant using +, *, and ^), A091333 (using +, -, and *), A091334 (using +, -, *, and ^), A348089 (using +, -, *, / and ^), A348262 (using + and ^).

Cf. A000792 (largest integer of given complexity), A003313, A076142, A076091, A061373, A005421, A064097, A005520, A003037, A161906, A161908, A244743.

import Data.List (genericIndex)
a005245 n = a005245_list `genericIndex` (n-1)
a005245_list = 1 : f 2 [1] where
   f x ys = y : f (x + 1) (y : ys) where
     y = minimum $
         (zipWith (+) (take (x `div` 2) ys) (reverse ys)) ++
         (zipWith (+) (map a005245 $ tail $ a161906_row x)
                      (map a005245 $ reverse $ init $ a161908_row x))
-- Reinhard Zumkeller, Mar 08 2013

with(numtheory):
a:= proc(n) option remember;
      `if`(n=1, 1, min(seq(a(i)+a(n-i), i=1..n/2),
       seq(a(d)+a(n/d), d=divisors(n) minus {1, n})))
    end:
seq(a(n), n=1..100); # Alois P. Heinz, Apr 18 2012

a[n_] := a[n] = If[n == 1, 1,
   Min[Table[a[i] + a[n-i], {i, 1, n/2}] ~Join~
   Table[a[d] + a[n/d], {d, Divisors[n] ~Complement~ {1, n}}]]];
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Dec 12 2013, after Alois P. Heinz *)

A005245(n /* start by calling this with the largest needed n */, lim/* see below */) = { local(d); n<6 && return(n);
if(n<=#A005245, A005245[n]&return(A005245[n]) /* return memoized result if available */,
A005245=vector(n) /* create vector if needed - should better reuse existing data if available */);
for(i=1, n-1, A005245[i] || A005245[i]=A005245(i,lim)); /* compute all previous elements */
A005245[n]=min( vecmin(vector(min(n\2,if(lim>0,lim,n)), k, A005245[k]+A005245[n-k])) /* additive possibilities - if lim>0 is given, consider a(k)+a(n-k) only for k<=lim - we know it is save to use lim=1 up to n=2e7 */, if( isprime(n), n, vecmin(vector((-1+#d=divisors(n))\2, i, A005245[d[i+1]]+A005245[d[ #d-i]]))/* multiplicative possibilities */))}
\\ See also the Python program by Tim Peters at A005421.
\\ M. F. Hasler, Jan 30 2008

from functools import lru_cache
from itertools import takewhile
from sympy import divisors
@lru_cache(maxsize=None)
def A005245(n): return min(min(A005245(a)+A005245(n-a) for a in range(1,(n>>1)+1)),min((A005245(d)+A005245(n//d) for d in takewhile(lambda d:d*d<=n,divisors(n)) if d>1),default=n)) if n>1 else 1 # Chai Wah Wu, Apr 29 2023

It is known that a(n) <= A061373(n) but I think 0 <= A061373(n) - a(n) <= 1 also holds. - Benoit Cloitre, Nov 23 2003 [That's false: the numbers {46, 235, 649, 1081, 7849, 31669, 61993} require {1, 2, 3, 4, 5, 6, 7} fewer 1's in A005245 than in A061373. - Ed Pegg Jr, Apr 13 2004]
It is known from the work of Selfridge and Coppersmith that 3 log_3 n <= a(n) <= 3 log_2 n = 4.754... log_3 n for all n > 1. [Guy, Unsolved Problems in Number Theory, Sect. F26.] - Charles R Greathouse IV, Apr 19 2012 [Comment revised by N. J. A. Sloane, Jul 17 2016]
Zelinsky (2022) improves the upper bound to a(n) <= A*log n where A = 41/log(55296) = 3.754422.... This compares to the constant 2.7307176... of the lower bound. - Charles R Greathouse IV, Dec 11 2022
a(n) >= A007600(n) is a very accurate lower bound. - Mehmet Sarraç, Dec 18 2022

More terms from David W. Wilson, May 15 1997

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

A116882 A number k is included if (highest odd divisor of k)^2 <= k.

Original entry on oeis.org

1, 2, 4, 8, 12, 16, 24, 32, 40, 48, 56, 64, 80, 96, 112, 128, 144, 160, 176, 192, 208, 224, 240, 256, 288, 320, 352, 384, 416, 448, 480, 512, 544, 576, 608, 640, 672, 704, 736, 768, 800, 832, 864, 896, 928, 960, 992, 1024, 1088, 1152, 1216, 1280, 1344, 1408

Offset: 1

Leroy Quet, Feb 24 2006

Also k is included if (and only if) the greatest power of 2 dividing k is >= the highest odd divisor of k. All terms of the sequence are even besides the 1.
Equivalently, positive integers of the form k*2^m, where odd k <= 2^m. - Thomas Ordowski, Oct 19 2014
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 consists of 1 and all numbers without a superior odd divisor. - Gus Wiseman, Feb 18 2021
Numbers k such that A006519(k) >= A000265(k), with equality only when k = 1. - Amiram Eldar, Jan 24 2023

			40 = 8 * 5, where 8 is highest power of 2 dividing 40 and 5 is the highest odd dividing 40. 8 is >= 5 (so 5^2 <= 40), so 40 is in the sequence.

G. C. Greubel, Table of n, a(n) for n = 1..1000
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Identities and periodic oscillations of divide-and-conquer recurrences splitting at half, arXiv:2210.10968 [cs.DS], 2022, p. 38.
Milan Janjic and Boris Petkovic, A Counting Function, arXiv 1301.4550 [math.CO], 2013.
Amelia Carolina Sparavigna, Discussion of the groupoid of Proth numbers (OEIS A080075), Politecnico di Torino, Italy (2019).

Cf. A000265, A006519, A080075, A112714.
The complement is A116883.
Positions of zeros (and 1) in A341675.
A051283 = numbers without a superior prime-power divisor (zeros of A341593).
A059172 = numbers without a superior squarefree divisor (zeros of A341592).
A063539 = numbers without a superior prime divisor (zeros of A341591).
A333805 counts strictly inferior odd divisors.
A341594 counts strictly superior odd divisors.
- Inferior: A033676, A038548, A063962, A066839, A069288, A161906, A217581.
- Superior: A033677, A063538, A070038, A072500, A161908, A341676.
- Strictly Inferior: A056924, A060775, A070039, A333806, A341596, A341674.
- Strictly Superior: A056924, A064052/A048098, A140271, A238535, A341642, A341643, A341673.
Cf. A000005, A000203, A001248, A006530, A020639, A026804, A027193, A340101, A340854 (zeros of A340832), A001008, A002805.
Subsequence of A082662, {1} U A363122.

f[n_] := Select[Divisors[n], OddQ[ # ] &][[ -1]]; Insert[Select[Range[2, 1500], 2^FactorInteger[ # ][[1]][[2]] > f[ # ] &], 1, 1] (* Stefan Steinerberger, Apr 10 2006 *)
q[n_] := 2^(2*IntegerExponent[n, 2]) >= n; Select[Range[1500], q] (* Amiram Eldar, Jan 24 2023 *)

isok(n) = vecmax(select(x->((x % 2)==1), divisors(n)))^2 <= n; \\ Michel Marcus, Sep 06 2016

isok(n) = 2^(valuation(n,2)*2) >= n \\ Jeppe Stig Nielsen, Feb 19 2019

from itertools import count, islice
def A116882_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:(n&-n)**2>=n,count(max(startvalue,1)))
A116882_list = list(islice(A116882_gen(),20)) # Chai Wah Wu, May 17 2023

a(n) = A080075(n-1)-1. - Klaus Brockhaus, Georgi Guninski and M. F. Hasler, Aug 16 2010
a(n) ~ n^2/2. - Thomas Ordowski, Oct 19 2014
Sum_{n>=1} 1/a(n) = 1 + (3/4) * Sum_{k>=1} H(2^k-1)/2^k = 2.3388865091..., where H(k) = A001008(k)/A002805(k) is the k-th harmonic number. - Amiram Eldar, Jan 24 2023

More terms from Stefan Steinerberger, Apr 10 2006

A342094 Number of integer partitions of n with no adjacent parts having quotient > 2.

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 9, 13, 16, 21, 27, 37, 44, 59, 75, 94, 117, 153, 186, 238, 296, 369, 458, 573, 701, 870, 1068, 1312, 1601, 1964, 2384, 2907, 3523, 4270, 5159, 6235, 7491, 9021, 10819, 12964, 15494, 18517, 22049, 26260, 31195, 37020, 43851, 51906, 61290

Offset: 1

Gus Wiseman, Mar 02 2021

nonn

The decapitation of such a partition (delete the greatest part) is term-wise greater than or equal to its negated first-differences.

			The a(1) = 1 through a(8) = 13 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (21)   (22)    (32)     (33)      (43)       (44)
             (111)  (211)   (221)    (42)      (322)      (53)
                    (1111)  (2111)   (222)     (421)      (332)
                            (11111)  (321)     (2221)     (422)
                                     (2211)    (3211)     (2222)
                                     (21111)   (22111)    (3221)
                                     (111111)  (211111)   (4211)
                                               (1111111)  (22211)
                                                          (32111)
                                                          (221111)
                                                          (2111111)
                                                          (11111111)

Fausto A. C. Cariboni, Table of n, a(n) for n = 1..250

The version with no adjacent parts having quotient < 2 is A000929.
The case of equality (all adjacent parts having quotient 2) is A154402.
A strong multiplicative version is A342083 or A342084.
The multiplicative version is A342085, with reciprocal version A337135.
The strict case is A342095.
The version with all adjacent parts having quotient < 2 is A342096, with strict case A342097.
The version with all adjacent parts having quotient > 2 is A342098.
The Heinz numbers of these partitions are listed by A342191.
A000009 counts strict partitions.
A003114 counts partitions with adjacent parts differing by more than 1.
A034296 counts partitions with adjacent parts differing by at most 1.
A038548 counts inferior (or superior) divisors, listed by A161906.
A161908 lists superior prime divisors.
Cf. A001055, A001227, A003242, A027193, A167606, A178470, A253784, A341591.

Table[Length[Select[IntegerPartitions[n],And@@Thread[Differences[-#]<=Rest[#]]&]],{n,30}]

A063962 Number of distinct prime divisors of n that are <= sqrt(n).

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 2, 0, 1, 1, 1, 0, 2, 0, 1, 1, 1, 0, 2, 1, 1, 1, 1, 0, 3, 0, 1, 1, 1, 1, 2, 0, 1, 1, 2, 0, 2, 0, 1, 2, 1, 0, 2, 1, 2, 1, 1, 0, 2, 1, 2, 1, 1, 0, 3, 0, 1, 2, 1, 1, 2, 0, 1, 1, 3, 0, 2, 0, 1, 2, 1, 1, 2, 0, 2, 1, 1, 0, 3, 1, 1, 1, 1, 0, 3, 1, 1, 1, 1, 1, 2, 0, 2, 1, 2, 0, 2, 0, 1, 3

Offset: 1

Reinhard Zumkeller, Sep 04 2001

nonn

For all primes p: a(p) = 0 (not marked) and for k > 1 a(p^k) = 1.
a(1) = 0 and for n > 0 a(n) is the number of marks when applying the sieve of Eratosthenes where a stage for prime p starts at p^2.
If we define a divisor d|n to be inferior if d <= n/d, then inferior divisors are counted by A038548 and listed by A161906. This sequence counts inferior prime divisors. - Gus Wiseman, Feb 25 2021

			a(33) = a(3*11) = 1, as 3^2 = 9 < 33 and 11^2 = 121 > 33.
From _Gus Wiseman_, Feb 25 2021: (Start)
The a(n) inferior prime divisors (columns) for selected n:
n =  3  8  24  3660  390  3570 87780
   ---------------------------------
    {}  2   2     2    2     2     2
            3     3    3     3     3
                  5    5     5     5
                      13     7     7
                            17    11
                                  19
(End)

Harry J. Smith, Table of n, a(n) for n = 1..1000

Cf. A055399, A001221.
Cf. A027748, A063962.
Zeros are at indices A008578.
The divisors are listed by A161906 and add up to A097974.
Dominates A333806 (the strictly inferior version).
The superior version is A341591.
The strictly superior version is A341642.
A001221 counts prime divisors, with sum A001414.
A033677 selects the smallest superior divisor.
A038548 counts inferior divisors.
A063538/A063539 have/lack a superior prime divisor.
A161908 lists superior divisors.
A207375 lists central divisors.
A217581 selects the greatest inferior prime divisor.
A341676 lists the unique superior prime divisors.
- Inferior: A033676, A066839, A069288, A072499, A333749, A333750.
- Superior: A051283, A059172, A070038, A072500, A116883, A341592, A341675.
- Strictly Inferior: A056924, A060775, A070039, A333805, A341596, A341674.
- Strictly Superior: A056924, A140271, A238535, A341594, A341595, A341673.
Cf. A000005, A001248, A005117, A006530, A020639, A048098, A064052, A341643.

a063962 n = length [p | p <- a027748_row n, p ^ 2 <= n]
-- Reinhard Zumkeller, Apr 05 2012

with(numtheory): a:=proc(n) local c,F,f,i: c:=0: F:=factorset(n): f:=nops(F): for i from 1 to f do if F[i]^2<=n then c:=c+1 else c:=c: fi od: c; end: seq(a(n),n=1..105); # Emeric Deutsch

Join[{0},Table[Count[Transpose[FactorInteger[n]][[1]],?(#<=Sqrt[n]&)],{n,2,110}]] (* _Harvey P. Dale, Mar 26 2015 *)

{ for (n=1, 1000, f=factor(n)~; a=0; for (i=1, length(f), if (f[1, i]^2<=n, a++, break)); write("b063962.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 04 2009

G.f.: Sum_{k>=1} x^(prime(k)^2) / (1 - x^prime(k)). - Ilya Gutkovskiy, Apr 04 2020

a(A002110(n)) = n for n > 2. - Gus Wiseman, Feb 25 2021

Revised definition from Emeric Deutsch, Jan 31 2006

A063538 Numbers n that are not sqrt(n-1)-smooth: largest prime factor of n (=A006530(n)) >= sqrt(n).

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 19, 20, 21, 22, 23, 25, 26, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 46, 47, 49, 51, 52, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 68, 69, 71, 73, 74, 76, 77, 78, 79, 82, 83, 85, 86, 87, 88, 89, 91

Offset: 1

N. J. A. Sloane, Aug 14 2001

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 all numbers with a superior prime divisor, which is unique (A341676) when it exists. For example, 42 is in the sequence because it has a prime divisor 7 which is greater than the quotient 42/7 = 6. - Gus Wiseman, Feb 19 2021

D. H. Greene and D. E. Knuth, Mathematics for the Analysis of Algorithms; see pp. 95-98.

Robert Israel, Table of n, a(n) for n = 1..10000
Beeler, M., Gosper, R. W. and Schroeppel, R., HAKMEM, ITEM 29

Cf. A006530, A063762.
Complement of A063539. Supersequence of A001358 (semiprimes).
The strictly superior version is A064052 (complement: A048098), with associated unique prime divisor A341643.
The case of odd instead of prime divisors is A116883 (complement: A116882).
Also nonzeros of A341591 (number of superior prime divisors).
The unique superior prime divisors of the terms are A341676.
A001221 counts prime divisors, with sum A001414.
A033677 selects the smallest superior divisor.
A038548 counts superior (also inferior) divisors.
A161908 lists superior divisors.
- Inferior: A033676, A063962, A066839, A161906, A217581, A333749, A333750.
- Superior: A051283, A059172, A070038, A072500, A341592, A341593, A341675.
- Strictly Inferior: A056924, A060775, A070039, A333805, A333806, A341596, A341674, A341677.
- Strictly Superior: A140271, A238535, A341594, A341595, A341642, A341644, A341645/A341646, A341673.
Cf. A000005, A001055, A001222, A001248, A020639, A056239, A112798.

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

Select[Range[2, 91], FactorInteger[#][[-1, 1]] >= Sqrt[#] &] (* Ivan Neretin, Aug 30 2015 *)

from math import isqrt
from sympy import primepi
def A063538(n):
    def f(x): return int(n+x-primepi(x//(y:=isqrt(x)))-sum(primepi(x//i)-primepi(i) for i in range(1,y)))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Oct 05 2024

Union of A001248 and A064052. - Gus Wiseman, Feb 24 2021

cp's OEIS Frontend

A038548 Number of divisors of n that are at most sqrt(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Formula

A033677 Smallest divisor of n >= sqrt(n).

Views

Author

Keywords

Comments

Examples

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

Extensions

A005245 The (Mahler-Popken) complexity of n: minimal number of 1's required to build n using + and *.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Python

Formula