A008578 -id:A008578 - cp's OEIS frontend

A006005 The odd prime numbers together with 1.

1, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271

Offset: 1

N. J. A. Sloane

The odd noncomposite numbers. Also odd primes at the beginning of the 20th century. - Omar E. Pol, Mar 19 2008
Indices at which records occur in A002322. - Artur Jasinski, Apr 05 2008
Odd numbers n such that their largest divisor <= sqrt(n) equals 1. (See A161344). - Omar E. Pol, Aug 03 2009
All k for which cos((k-1)!*Pi/k) is negative. - Gerry Martens, Jun 01 2018

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972, p. 870 [alternative scanned copy].
S. M. Losanitsch, Die Isomerie-Arten bei den Homologen der Paraffin-Reihe, Chem. Ber. 30 (1897), 1917-1926. (Annotated scanned copy)

Cf. A000040, A008578, A065091, A002322.

a = {}; max = 0; Do[w = CarmichaelLambda[k]; If[w > max, AppendTo[a, k]; max = k], {k, 1, 200}]; a (* Artur Jasinski, Apr 05 2008 *)
Join[{1},Prime[Range[2,60]]] (* Harvey P. Dale, Apr 15 2019 *)

prime(n)-(n==1) \\ Charles R Greathouse IV, Aug 26 2011

a(n) = A000040(n) - A000007(n) = A000040(n) - ((-1)^A000040(n)+1)/2. - Juri-Stepan Gerasimov, Oct 25 2009

a(n) = A175524(n) for n <= 30. - Reinhard Zumkeller, Jul 12 2011

A034891 Number of different products of partitions of n; number of partitions of n into prime parts (1 included); number of distinct orders of Abelian subgroups of symmetric group S_n.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 8, 11, 14, 18, 23, 29, 36, 45, 55, 67, 81, 98, 117, 140, 166, 196, 231, 271, 317, 369, 429, 496, 573, 660, 758, 869, 993, 1133, 1290, 1465, 1662, 1881, 2125, 2397, 2699, 3035, 3407, 3820, 4276, 4780, 5337, 5951, 6628, 7372, 8191, 9090

Offset: 0

Wouter Meeussen

a(n) = length of n-th row in A212721. - Reinhard Zumkeller, Jun 14 2012

Number of partitions of n into noncomposite parts. - Omar E. Pol, Jun 23 2022

Alois P. Heinz, Table of n, a(n) for n = 0..10000 (terms n = 1..1000 from T. D. Noe)
Index entries for sequences related to partitions

Cf. A000041, A000792, A000793, A009490, A008578, A140436, A353188.

a034891 = length . a212721_row  -- Reinhard Zumkeller, Jun 14 2012

b:= proc(n, i) option remember; `if`(n=0, 1, (p->
      `if`(i<0, 0, b(n, i-1)+ `if`(p>n, 0,
         b(n-p, i))))(`if`(i<1, 1, ithprime(i))))
    end:
a:= n-> b(n, numtheory[pi](n)):
seq(a(n), n=0..100);  # Alois P. Heinz, Feb 15 2013

Table[ Length[ Union[ Apply[ Times, Partitions[ n], 1]]], {n, 30}]
CoefficientList[ Series[ (1/(1 - x)) Product[1/(1 - x^Prime[i]), {i, 100}], {x, 0, 50}], x] (* Robert G. Wilson v, Aug 17 2013 *)
b[n_, i_] := b[n, i] = Module[{p}, p = If[i<1, 1, Prime[i]]; If[n == 0, 1, If[i<0, 0, b[n, i-1] + If[p>n, 0, b[n-p, i]]]]]; a[n_] := b[n, PrimePi[n] ]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Nov 05 2015, after Alois P. Heinz *)

[Partitions(n, parts_in=(prime_range(n+1)+[1])).cardinality() for n in xsrange(1000)] # Giuseppe Coppoletta, Jul 11 2016

G.f.: (1/(1-x))*(1/Product_{k>0} (1-x^prime(k))). a(n) = (1/n)*Sum_{k=1..n} A074372(k)*a(n-k). Partial sums of A000607. - Vladeta Jovovic, Sep 19 2002

a(n) = A000041(n) - A353188(n). - Omar E. Pol, Jun 23 2022

More terms from Vladeta Jovovic

a(0)=1 from Michael Somos, Feb 05 2011

A161344 Numbers k with A033676(k)=2, where A033676 is the largest divisor <= sqrt(k).

Original entry on oeis.org

4, 6, 8, 10, 14, 22, 26, 34, 38, 46, 58, 62, 74, 82, 86, 94, 106, 118, 122, 134, 142, 146, 158, 166, 178, 194, 202, 206, 214, 218, 226, 254, 262, 274, 278, 298, 302, 314, 326, 334, 346, 358, 362, 382, 386, 394, 398, 422, 446, 454, 458, 466, 478, 482, 502, 514

Offset: 1

Omar E. Pol, Jun 20 2009

Define a sieve operation with parameter s that eliminates integers of the form s^2 + s*i (i >= 0) from the set A000027 of natural numbers. The sequence lists those natural numbers that are eliminated by the sieve s=2 and cannot be eliminated by any sieve s >= 3. - R. J. Mathar, Jun 24 2009
After a(3)=8 all terms are 2*prime; for n > 3, a(n) = 2*prime(n-1) = 2*A000040(n-1). - Zak Seidov, Jul 18 2009
From Omar E. Pol, Jul 18 2009: (Start)
A classification of the natural numbers A000027.
=============================================================
Numbers k whose largest divisor <= sqrt(k) equals j
=============================================================
j       Sequence     Comment
=============================================================
1 ..... A008578      1 together with the prime numbers
2 ..... A161344      This sequence
3 ..... A161345
4 ..... A161424
5 ..... A161835
6 ..... A162526
7 ..... A162527
8 ..... A162528
9 ..... A162529
10 .... A162530
11 .... A162531
12 .... A162532
... And so on. (End)
The numbers k whose largest divisor <= sqrt(k) is j are exactly those numbers j*m where m is either a prime >= k or one of the numbers in row j of A163925. - Franklin T. Adams-Watters, Aug 06 2009
See also A163280, the main entry for this sequence. - Omar E. Pol, Oct 24 2009
Also A100484 UNION 8. - Omar E. Pol, Nov 29 2012 (after Seidov and Hasler)
Is this the union of {4} and A073582? - R. J. Mathar, May 30 2025

Omar E. Pol, Determinacion geometrica de los numeros primos y perfectos
Omar E. Pol, Illustration: Divisors and pi(x)
Omar E. Pol, Illustration of initial terms
Omar E. Pol. Illustration of initial terms of A008578, A161344, A161345, A161424

Cf. A000005, A018253, A160811, A160812, A161205, A161346, A033676, A008578, A161345, A161424, A161835, A162526, A162527, A162528, A162529, A162530, A162531, A162532, A163925.

Second column of array in A163280. Also, second row of array in A163990.

isA := proc(n,s) if n mod s <> 0 then RETURN(false); fi; if n/s-s >= 0 then RETURN(true); else RETURN(false); fi; end: isA161344 := proc(n) for s from 3 to n do if isA(n,s) then RETURN(false); fi; od: isA(n,2) ; end: for n from 1 to 3000 do if isA161344(n) then printf("%d,",n) ; fi; od; # R. J. Mathar, Jun 24 2009

a[n_] := If[n <= 3, 2n+2, 2*Prime[n-1]]; Table[a[n], {n, 1, 56}] (* Jean-François Alcover, Nov 26 2012, after Zak Seidov *)

a(n)=if(n>3,prime(n-1),n+1)*2 \\ M. F. Hasler, Nov 27 2012

Equals 2*A000040 union {8}. - M. F. Hasler, Nov 27 2012

a(n) = 2*A046022(n+1) = 2*A175787(n). - Omar E. Pol, Nov 27 2012

More terms from R. J. Mathar, Jun 24 2009

Definition added by R. J. Mathar, Jun 28 2009

A033286 a(n) = n * prime(n).

Original entry on oeis.org

2, 6, 15, 28, 55, 78, 119, 152, 207, 290, 341, 444, 533, 602, 705, 848, 1003, 1098, 1273, 1420, 1533, 1738, 1909, 2136, 2425, 2626, 2781, 2996, 3161, 3390, 3937, 4192, 4521, 4726, 5215, 5436, 5809, 6194, 6513, 6920, 7339, 7602, 8213, 8492, 8865, 9154, 9917

Offset: 1

N. J. A. Sloane

Does an n exist such that n*prime(n)/(n+prime(n)) is an integer? - Ctibor O. Zizka, Mar 04 2008. The answer to Zizka's question is easily seen to be No: such an integer k would be positive and less than prime(n), but then k*(n + prime(n)) = prime(n)*n would be impossible. - Robert Israel, Apr 20 2015
Sums of rows of the triangle in A005145. - Reinhard Zumkeller, Aug 05 2009
Complement of A171520(n). - Jaroslav Krizek, Dec 13 2009
Partial sums of A090942. - Omar E. Pol, Apr 20 2015

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Albert Frank, International Contest Of Logical Sequences, 2002 - 2003. Item 1.
Albert Frank, Solutions of International Contest Of Logical Sequences, 2002 - 2003.

Cf. A000040, A007504, A014689, A090942, A124012, A141042, A152535.

Cf. A005145 (primes repeated), A171520 (complement), A076146 (iterated).

a033286 n = a000040 n * n  -- Reinhard Zumkeller, Jul 24 2013

[ n*NthPrime(n): n in [1..47] ]; // Klaus Brockhaus, Sep 09 2009

A033286 := proc(n) n*ithprime(n) ; end proc:
seq(A033286(n),n=1..20) ; # R. J. Mathar, Mar 21 2011

Table[Prime[n]*n, {n, 38}] (* Alonso del Arte *)

ithprime(i)*i $ i = 1..47 // Zerinvary Lajos, Feb 26 2007

a(n)=n*prime(n) \\ Charles R Greathouse IV, Jul 01 2013

a(n) = n * A000040(n) = n * A008578(n+1) = n * A158611(n+2). - Jaroslav Krizek, Aug 31 2009
a(n) = A007504(n) + A152535(n). - Omar E. Pol, Aug 09 2012
Sum_{n>=1} 1/a(n) = A124012. - Amiram Eldar, Oct 15 2020

Correction for change of offset in A158611 and A008578 in Aug 2009 from Jaroslav Krizek, Jan 27 2010

A065387 a(n) = sigma(n) + phi(n).

Original entry on oeis.org

2, 4, 6, 9, 10, 14, 14, 19, 19, 22, 22, 32, 26, 30, 32, 39, 34, 45, 38, 50, 44, 46, 46, 68, 51, 54, 58, 68, 58, 80, 62, 79, 68, 70, 72, 103, 74, 78, 80, 106, 82, 108, 86, 104, 102, 94, 94, 140, 99, 113, 104, 122, 106, 138, 112, 144, 116, 118, 118, 184, 122, 126, 140

Offset: 1

Labos Elemer, Nov 05 2001

a(n) = 2n for n listed in A008578, the prime numbers at the beginning of the 20th century. When a(n) = a(n + 1), n is probably listed in A066198, numbers n where phi changes as fast as sigma (the only exceptions below 10000 are 2 and 854). - Alonso del Arte, Nov 16 2005
A. Makowski proved that n is prime if and only if a(n) = n * d(n), where d is A000005. - Charles R Greathouse IV, Mar 19 2012
If n is semiprime, a(n) = 2n+1+ceiling(sqrt(n))-floor(sqrt(n)). - Wesley Ivan Hurt, May 05 2015
Atanassov proves that a(n) >= n + A001414(n). - Charles R Greathouse IV, Dec 06 2016
a(n) = 2*n+1 iff n is square of prime (A001248), a(n) = 2*(n+1) iff n is squarefree semiprime (A006881). - Bernard Schott, Feb 09 2020

			a(10) = 22 because there are 4 coprimes to 10 below 10, the divisors of 10 add up to 18, and 4 + 18 = 22.

K. Atanassov, New integer functions, related to ψ and σ functions. IV., Bull. Number Theory Related Topics 12 (1988), pp. 31-35.
Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004. See Section B41, p. 149.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 162.

N. J. A. Sloane, Table of n, a(n) for n = 1..10000 (First 1000 terms from T. D. Noe.)
A. Makowski, Aufgaben 339, Elemente der Mathematik 15 (1960), pp. 39-40.

Cf. A000010, A000203, A065388, A015702, A051709, A011774.

See A292768 for partial sums, A051612 for sigma - phi.

[DivisorSigma(1,k)+EulerPhi(k):k in [1..65]]; // Marius A. Burtea, Feb 09 2020

with(numtheory); A065387:=n->phi(n) + sigma(n); seq(A065387(n), n=1..100); # Wesley Ivan Hurt, Apr 08 2014

Table[EulerPhi[n] + DivisorSigma[1,n], {n, 65}] (* Alonso del Arte *)
a[n_] := SeriesCoefficient[Sum[(1+MoebiusMu[k])*x^k/(1-x^k)^2, {k, 1, n}], {x, 0, n}]; Array[a, 63] (* Jean-François Alcover, Sep 29 2017, after Ilya Gutkovskiy *)

a(n) = sigma(n) + eulerphi(n) \\ Harry J. Smith, Oct 17 2009

[sigma(n,1)+euler_phi(n) for n in range(1, 64)] # Stefano Spezia, Jul 20 2025

a(n) = A000203(n) + A000010(n).
a(n) = A051709(n) + 2n. - N. J. A. Sloane, Jun 12 2004
G.f.: Sum_{k>=1} (mu(k) + 1)*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Sep 29 2017

A207375 Irregular array read by rows in which row n lists the (one or two) central divisors of n in increasing order.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 1, 5, 2, 3, 1, 7, 2, 4, 3, 2, 5, 1, 11, 3, 4, 1, 13, 2, 7, 3, 5, 4, 1, 17, 3, 6, 1, 19, 4, 5, 3, 7, 2, 11, 1, 23, 4, 6, 5, 2, 13, 3, 9, 4, 7, 1, 29, 5, 6, 1, 31, 4, 8, 3, 11, 2, 17, 5, 7, 6, 1, 37, 2, 19, 3, 13, 5, 8, 1, 41, 6, 7, 1, 43

Offset: 1

Omar E. Pol, Feb 23 2012

If n is a square then row n lists only the square root of n because the squares (A000290) have only one central divisor.
If n is not a square then row n lists the pair (j, k) of divisors of n, nearest to the square root of n, such that j*k = n.
Conjecture 1: It appears that the n-th record in this sequence is the last member of row A008578(n).
Column 1 gives A033676. Right border gives A033677. - Omar E. Pol, Feb 26 2019
The conjecture 1 follows from Bertrand's Postulate. - Charles R Greathouse IV, Feb 11 2022
Row products give A097448. - Omar E. Pol, Feb 17 2022

			Array begins:
  1;
  1,  2;
  1,  3;
  2;
  1,  5;
  2,  3;
  1,  7;
  2,  4;
  3;
  2,  5;
  1, 11;
  3,  4;
  1, 13;
...

Alois P. Heinz, Rows n = 1..5000
Omar E. Pol, Illustration of the divisors of the first 12 positive integers

Row n has length A169695(n).
Row sums give A207376.
Cf. A000005, A000290, A008578, A027750, A033676, A033677, A097448, A161901, A161904.

A207375row[n_] := ArrayPad[#, -Floor[(Length[#] - 1)/2]] & [Divisors[n]];
Array[A207375row, 50] (* Paolo Xausa, Apr 07 2025 *)

A001747 2 together with primes multiplied by 2.

Original entry on oeis.org

2, 4, 6, 10, 14, 22, 26, 34, 38, 46, 58, 62, 74, 82, 86, 94, 106, 118, 122, 134, 142, 146, 158, 166, 178, 194, 202, 206, 214, 218, 226, 254, 262, 274, 278, 298, 302, 314, 326, 334, 346, 358, 362, 382, 386, 394, 398, 422, 446, 454, 458, 466, 478, 482, 502

Offset: 1

N. J. A. Sloane

When supplemented with 8, may be considered the "even primes", since these are the even numbers n = 2k which are divisible just by 1, 2, k and 2k. - Louis Zuckerman (louis(AT)trapezoid.com), Sep 12 2000
Sequence gives solutions of sigma(n) - phi(n) = n + tau(n) where tau(n) is the number of divisors of n.
Numbers n such that sigma(n) = 3*(n - phi(n)).
Except for 2, orders of non-cyclic groups k (in A060679(n)) such that x^k==1 (mod k) has only 1 solution 2<=x<=k. - Benoit Cloitre, May 10 2002
Numbers n such that A092673(n) = 2. - Jon Perry, Mar 02 2004
Except for initial terms, this sequence = A073582 = A074845 = A077017. Starting with the term 10, they are identical. - Robert G. Wilson v, Jun 15 2004
Together with 8 and 16, even numbers n such that n^2 does not divide (n/2)!. - Arkadiusz Wesolowski, Jul 16 2011
Twice noncomposite numbers. - Omar E. Pol, Jan 30 2012

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

Cf. A060679, A009530, A098764.

Equals {2} UNION {A100484}.

Concatenation([2], List([1..60], n-> 2*Primes[n])); # G. C. Greubel, May 18 2019

[2] cat [2*NthPrime(n): n in [1..60]]; // G. C. Greubel, May 18 2019

Join[{2},2*Prime[Range[60]]] (* Harvey P. Dale, Jul 23 2013 *)

print1(2);forprime(p=2,97,print1(", "2*p)) \\ Charles R Greathouse IV, Jan 31 2012

[2]+[2*nth_prime(n) for n in (1..60)] # G. C. Greubel, May 18 2019

a(n) = A001043(n) - A001223(n+1), except for initial term.
a(n) = A116366(n-2,n-2) for n>2. - Reinhard Zumkeller, Feb 06 2006
A006093(n) = A143201(a(n+1)) for n>1. - Reinhard Zumkeller, Aug 12 2008
a(n) = 2*A008578(n). - Omar E. Pol, Jan 30 2012, and Reinhard Zumkeller, Feb 16 2012

A036497 Number of partitions of n into distinct primes (counting 1 as a prime).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 2, 3, 2, 3, 3, 3, 4, 4, 4, 5, 5, 6, 7, 7, 8, 8, 9, 10, 10, 11, 11, 11, 13, 13, 15, 16, 16, 18, 18, 20, 22, 22, 24, 25, 26, 29, 30, 32, 33, 34, 37, 39, 41, 44, 45, 47, 51, 53, 57, 59, 61, 64, 67, 72, 76, 79, 82, 86, 89, 95, 100, 103, 108, 112, 118

Offset: 0

Wouter Meeussen, Dec 17 1998

nonn

Honsberger shows that the primes-including-1 are a complete sequence and therefore all numbers in this sequence exceed zero. - Ron Knott, Aug 27 2016

Number of partitions of n into distinct noncomposite numbers. - Omar E. Pol, Dec 14 2024

			a(11) = 3 since 11 = 1+2+3+5=1+3+7 has 3 partitions of distinct primes-including-1. - _Ron Knott_, Aug 27 2016

Ross Honsberger, Mathematical Gems III, The Mathematical Association of America, 1985, pages 127-128.

Alois P. Heinz, Table of n, a(n) for n = 0..10000

Cf. A000586, A008578.

s:= proc(n) option remember;
      `if`(n<1, n+1, ithprime(n)+s(n-1))
    end:
b:= proc(n, i) option remember; (p-> `if`(n=0, 1,
      `if`(n>s(i), 0, b(n, i-1)+ `if`(p>n, 0,
       b(n-p, i-1)))))(`if`(i<1, 1, ithprime(i)))
    end:
a:= n-> b(n, numtheory[pi](n)):
seq(a(n), n=0..100);  # Alois P. Heinz, Aug 27 2016

myprime[ n_ ] := If[ n===0, 1, Prime[ n ] ]; ta1=Table[ Product[ 1+z^myprime[ k ], {k, 0, n} ]~CoefficientList~z, {n, 31, 32} ]; leveled=Count[ Take[ Last@ta1, Length@ta1[ [ -2 ] ] ]-ta1[ [ -2 ] ], 0 ]; Take[ Last@ta1, leveled ]
Table[Length@ DeleteCases[DeleteCases[IntegerPartitions@ n, {_, a_, _} /; CompositeQ@ a], w_ /; MemberQ[Differences@ w, 0]], {n, 0, 60}] (* Michael De Vlieger, Aug 27 2016 *)

G.f.: (1 + x)*Product_{k>=1} (1 + x^prime(k)). - Ilya Gutkovskiy, Dec 31 2016

A161345 Numbers k whose largest divisor <= sqrt(k) is 3.

Original entry on oeis.org

9, 12, 15, 18, 21, 27, 33, 39, 51, 57, 69, 87, 93, 111, 123, 129, 141, 159, 177, 183, 201, 213, 219, 237, 249, 267, 291, 303, 309, 321, 327, 339, 381, 393, 411, 417, 447, 453, 471, 489, 501, 519, 537, 543, 573, 579, 591, 597, 633, 669, 681, 687, 699, 717, 723

Offset: 1

Omar E. Pol, Jun 20 2009

Define a sieve operation with parameter s that eliminates integers of the form s^2+s*i (i >= 0) from the set A000027 of natural numbers. The sequence lists those natural numbers that are eliminated by the sieve s=3 and cannot be eliminated by any sieve s >= 4. - R. J. Mathar, Jun 24 2009
See A161344 for more information. - Omar E. Pol, Jul 05 2009
See also the array in A163280, the main entry for this sequence. - Omar E. Pol, Oct 24 2009
Union of {12, 18, 27} and all the numbers of the form 3*p, where p is an odd prime. - Amiram Eldar, Apr 17 2024

Omar E. Pol, Determinacion geometrica de los numeros primos y perfectos.
Omar E. Pol, Illustration: Divisors and pi(x).
Omar E. Pol, Illustration for A008578, A161344, A161345 and A161424. [From _Omar E. Pol_, Oct 24 2009]

Cf. A000005, A018253, A160811, A160812, A161205, A161344, A161346, A033676, A008578, A161424, A161835, A162526, A162527, A162528, A162529, A162530, A162531, A162532.

Third column of the array in A163280. Also, third row of array in A163990. - Omar E. Pol, Oct 24 2009

isA := proc(n,s) if n mod s <> 0 then RETURN(false); fi; if n/s-s >= 0 then RETURN(true); else RETURN(false); fi; end: isA161345 := proc(n) for s from 4 to n do if isA(n,s) then RETURN(false); fi; od: isA(n,3) ; end: for n from 1 to 3000 do if isA161345(n) then printf("%d,",n) ; fi; od; # R. J. Mathar, Jun 24 2009

md3Q[n_]:=Max[Select[Divisors[n],#<=Sqrt[n]&]]==3; Select[Range[800],md3Q] (* Harvey P. Dale, Aug 12 2013 *)

Numbers k such that A033676(k)=3. - Omar E. Pol, Jul 05 2009

Terms beyond a(10) from R. J. Mathar, Jun 24 2009

Definition added by R. J. Mathar, Jun 28 2009

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

cp's OEIS Frontend

A006005 The odd prime numbers together with 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A034891 Number of different products of partitions of n; number of partitions of n into prime parts (1 included); number of distinct orders of Abelian subgroups of symmetric group S_n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Sage

Formula

Extensions

A161344 Numbers k with A033676(k)=2, where A033676 is the largest divisor <= sqrt(k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A033286 a(n) = n * prime(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

MuPAD

PARI

Formula

Extensions

A065387 a(n) = sigma(n) + phi(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A207375 Irregular array read by rows in which row n lists the (one or two) central divisors of n in increasing order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs