A070824 -id:A070824 - cp's OEIS frontend

A354697 a(n) is the least integer that can be written in two or more ways as the product of the integers in two subsets of its A070824(a(n)) nontrivial divisors, each of size n and with empty intersection.

12, 120, 720, 10080, 110880, 1814400, 26611200, 518918400, 10378368000, 261534873600, 5928123801600, 168951528345600, 4505374089216000, 152056375511040000, 4663062182338560000, 167870238564188160000, 6463004184721244160000, 249902828475888107520000, 10495918795987300515840000

Offset: 2

Hugo Pfoertner, Jun 03 2022

a(6) <= 110880 = 2*3*6*10*14*22 = 4*5*7*8*9*11.

			a(2) = 12 = 2*6 = 3*4,
a(3) = 120 = 2*3*20 = 4*5*6,
a(4) = 720 = 2*4*9*10 = 3*5*6*8,
a(5) = 10080 = 2*3*6*10*28 = 4*5*7*8*9.
a(6) = 110880 = 2*3*6*10*14*22 = 4*5*7*8*9*11.
a(7) = 1814400 = 2*3*4*14*15*18*20 = 5*6*7*8*9*10*12.

Zhao Hui Du, Table of n, a(n) for n = 2..25

Cf. A002182, A025487, A354457.

A070824(a(n)) >= 2*n.

a(6) confirmed by and a(7)-a(13) from David A. Corneth, Jun 04 2022

a(14) onwards from Zhao Hui Du, May 12 2024

A307118 a(1) = 0; for n>1, a(n) = dr(n-1) + dr(n) + dr(n+1), where dr(n) is the number of nontrivial divisors of n (A070824).

Original entry on oeis.org

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

Offset: 1

Todor Szimeonov, Mar 25 2019

nonn

Real divisibility of n's one-area (or 1-area). This is the first step to examine the divisibility of n's k-area. n's k-area is the set of m for which |n-m| is less than or equal to k (where n, k, m are natural numbers). 1's 1-area is {1,2}, 5's 1-area {4,5,6}, 3's 2-area {1,2,3,4,5}. We could call this natural area, and still talk about nonnegative or integer areas, etc.

Cf. A070824, A307119, A307120.

{0}~Join~MapAt[# + 1 &, Total /@ Partition[DivisorSigma[0, Range@ 82] - 2, 3, 1], 1] (* Michael De Vlieger, Jun 06 2019 *)

dr(n) = if (n<2, 0, numdiv(n)-2);
a(n) = if (n==1, 0, dr(n-1) + dr(n) + dr(n+1)); \\ Michel Marcus, Apr 11 2019

A147588 Duplicate of A070824.

Original entry on oeis.org

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

Offset: 1

dead

A010051 Characteristic function of primes: 1 if n is prime, else 0.

Original entry on oeis.org

0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0

Offset: 1

N. J. A. Sloane

The following sequences all have the same parity (with an extra zero term at the start of a(n)): a(n), A061007, A035026, A069754, A071574. - Jeremy Gardiner, Aug 09 2002
Hardy and Wright prove that the real number 0.011010100010... is irrational. See Nasehpour link. - Michel Marcus, Jun 21 2018
The spectral components (excluding the zero frequency) of the Fourier transform of the partial sequences {a(j)} with j=1..n and n an even number, exhibit a remarkable symmetry with respect to the central frequency component at position 1 + n/4. See the Fourier spectrum of the first 2^20 terms in Links, Comments in A289777, and Conjectures in A001223 of Sep 01 2019. It also appears that the symmetry grows with n. - Andres Cicuttin, Aug 23 2020

J.-P. Allouche and J. Shallit, Automatic Sequences, Cambridge Univ. Press, 2003, p. 3.
V. Brun, Über das Goldbachsche Gesetz und die Anzahl der Primzahlpaare, Arch. Mat. Natur. B, 34, no. 8, 1915.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Oxford University Press, London, 1975.
Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 65.
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 132.

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 10000 terms from N. J. A. Sloane)
Andres Cicuttin, Fourier spectrum of the first 2^20 terms of the characteristic function of primes.
William Craig, Jan-Willem van Ittersum and Ken Ono, Integer partitions detect the primes, PNAS, Vol. 121, No. 39 (2024), e2409417121.
Yoichi Motohashi, An overview of the Sieve Method and its History, arXiv:math/0505521 [math.NT], 2005-2006.
Peyman Nasehpour, A Computational Criterion for the Irrationality of Some Real Numbers, Journal of Algorithms and Computation, Vol. 52, No. 1 (2020), pp. 97-104, preprint, arXiv:1806.07560 [math.AC], 2018.
José L. Ramírez and Gustavo N. Rubiano, Properties and Generalizations of the Fibonacci Word Fractal, The Mathematica Journal, Vol. 16 (2014).
Eric Weisstein's World of Mathematics, Prime Number.
Eric Weisstein's World of Mathematics, Prime Constant.
Eric Weisstein's World of Mathematics, Prime zeta function primezeta(s).
Index entries for characteristic functions

Cf. A051006 (constant 0.4146825... (base 10) = 0.01101010001010001010... (base 2)), A001221 (inverse Moebius transform), A143519, A156660, A156659, A156657, A059500, A053176, A059456, A072762.
First differences of A000720, so A000720 gives partial sums.
Column k=1 of A117278.
Characteristic function of A000040.
Cf. A008683.

import Data.List (unfoldr)
a010051 :: Integer -> Int
a010051 n = a010051_list !! (fromInteger n-1)
a010051_list = unfoldr ch (1, a000040_list) where
   ch (i, ps'@(p:ps)) = Just (fromEnum (i == p),
                              (i + 1, if i == p then ps else ps'))
-- Reinhard Zumkeller, Apr 17 2012, Sep 15 2011

s:=[]; for n in [1..100] do if IsPrime(n) then s:=Append(s,1); else s:=Append(s,0); end if; end for; s;

[IsPrime(n) select 1 else 0: n in [1..100]];  // Bruno Berselli, Mar 02 2011

A010051:= n -> if isprime(n) then 1 else 0 fi;

Table[ If[ PrimeQ[n], 1, 0], {n, 105}] (* Robert G. Wilson v, Jan 15 2005 *)
Table[Boole[PrimeQ[n]], {n, 105}] (* Alonso del Arte, Aug 09 2011 *)
Table[PrimePi[n] - PrimePi[n-1], {n,50}] (* G. C. Greubel, Jan 05 2017 *)

a(n)=isprime(n) \\ Charles R Greathouse IV, Apr 16 2011

from sympy import isprime
def A010051(n): return int(isprime(n)) # Chai Wah Wu, Jan 20 2022

a(n) = floor(cos(Pi*((n-1)! + 1)/n)^2) for n >= 2. - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Nov 07 2002
Let M(n) be the n X n matrix m(i, j) = 0 if n divides ij + 1, m(i, j) = 1 otherwise; then for n > 0 a(n) = -det(M(n)). - Benoit Cloitre, Jan 17 2003
n >= 2, a(n) = floor(phi(n)/(n - 1)) = floor(A000010(n)/(n - 1)). - Benoit Cloitre, Apr 11 2003
a(n) = Sum_{d|gcd(n, A034386(n))} mu(d). [Brun]
a(m*n) = a(m)*0^(n - 1) + a(n)*0^(m - 1). - Reinhard Zumkeller, Nov 25 2004
a(n) = 1 if n has no divisors other than 1 and n, and 0 otherwise. - Jon Perry, Jul 02 2005
Dirichlet generating function: Sum_{n >= 1} a(n)/n^s = primezeta(s), where primezeta is the prime zeta function. - Franklin T. Adams-Watters, Sep 11 2005
a(n) = (n-1)!^2 mod n. - Franz Vrabec, Jun 24 2006
a(n) = A047886(n, 1). - Reinhard Zumkeller, Apr 15 2008
Equals A051731 (the inverse Möbius transform) * A143519. - Gary W. Adamson, Aug 22 2008
a(n) = A051731((n + 1)! + 1, n) from Wilson's theorem: n is prime if and only if (n + 1)! is congruent to -1 mod n. - N-E. Fahssi, Jan 20 2009, Jan 29 2009
a(n) = A166260/A001477. - Mats Granvik, Oct 10 2009
a(n) = 0^A070824, where 0^0=1. - Mats Granvik, Gary W. Adamson, Feb 21 2010
It appears that a(n) = (H(n)*H(n + 1)) mod n, where H(n) = n!*Sum_{k=1..n} 1/k = A000254(n). - Gary Detlefs, Sep 12 2010
Dirichlet generating function: log( Sum_{n >= 1} 1/(A112624(n)*n^s) ). - Mats Granvik, Apr 13 2011
a(n) = A100995(n) - sqrt(A100995(n)*A193056(n)). - Mats Granvik, Jul 15 2011
a(n) * (2 - n mod 4) = A151763(n). - Reinhard Zumkeller, Oct 06 2011
(n - 1)*a(n) = ( (2*n + 1)!! * Sum_{k=1..n}(1/(2*k + 1))) mod n, n > 2. - Gary Detlefs, Oct 07 2011
For n > 1, a(n) = floor(1/A001222(n)). - Enrique Pérez Herrero, Feb 23 2012
a(n) = mu(n) * Sum_{d|n} mu(d)*omega(d), where mu is A008683 and omega A001222 or A001221 indistinctly. - Enrique Pérez Herrero, Jun 06 2012
a(n) = A003418(n+1)/A003418(n) - A217863(n+1)/A217863(n) = A014963(n) - A072211(n). - Eric Desbiaux, Nov 25 2012
For n > 1, a(n) = floor(A014963(n)/n). - Eric Desbiaux, Jan 08 2013
a(n) = ((abs(n-2))! mod n) mod 2. - Timothy Hopper, May 25 2015
a(n) = abs(F(n)) - abs(F(n)-1/2) - abs(F(n)-1) + abs(f(n)-3/2), where F(n) = Sum_{m=2..n+1} (abs(1 - (n mod m)) - abs(1/2 - (n mod m)) + 1/2), n > 0. F(n) = 1 if n is prime, > 1 otherwise, except F(1) = 0. a(n) = 1 if F(n) = 1, 0 otherwise. - Timothy Hopper, Jun 16 2015
For n > 4, a(n) = (n-2)! mod n. - Thomas Ordowski, Jul 24 2016
From Ilya Gutkovskiy, Jul 24 2016: (Start)
G.f.: A(x) = Sum_{n>=1} x^A000040(n) = B(x)*(1 - x), where B(x) is the g.f. for A000720.
a(n) = floor(2/A000005(n)), for n>1. (End)
a(n) = pi(n) - pi(n-1) = A000720(n) - A000720(n-1), for n>=1. - G. C. Greubel, Jan 05 2017
Decimal expansion of Sum_{k>=1} (1/10)^prime(k) = 9 * Sum_{k>=1} pi(k)/10^(k+1), where pi(k) = A000720(k). - Amiram Eldar, Aug 11 2020
a(n) = 1 - ceiling((2/n) * Sum_{k=2..floor(sqrt(n))} floor(n/k)-floor((n-1)/k)), n>1. - Gary Detlefs, Sep 08 2023
a(n) = Sum_{d|n} mu(d)*omega(n/d), where mu = A008683 and omega = A001221. - Ridouane Oudra, Apr 12 2025
a(n) = 0 if (n^2 - 3*n + 2) * A000203(n) - 8 * A002127(n) > 0 else 1 (n>2, see Craig link). - Bill McEachen, Jul 04 2025

A032741 a(0) = 0; for n > 0, a(n) = number of proper divisors of n (divisors of n which are less than n).

Original entry on oeis.org

0, 0, 1, 1, 2, 1, 3, 1, 3, 2, 3, 1, 5, 1, 3, 3, 4, 1, 5, 1, 5, 3, 3, 1, 7, 2, 3, 3, 5, 1, 7, 1, 5, 3, 3, 3, 8, 1, 3, 3, 7, 1, 7, 1, 5, 5, 3, 1, 9, 2, 5, 3, 5, 1, 7, 3, 7, 3, 3, 1, 11, 1, 3, 5, 6, 3, 7, 1, 5, 3, 7, 1, 11, 1, 3, 5, 5, 3, 7, 1, 9, 4, 3, 1, 11, 3, 3, 3, 7, 1, 11, 3, 5, 3, 3, 3, 11, 1, 5, 5

Offset: 0

Patrick De Geest, May 15 1998

Number of d < n which divide n.
Call an integer k between 1 and n a "semi-divisor" of n if n leaves a remainder of 1 when divided by k, i.e., n == 1 (mod k). a(n) gives the number of semi-divisors of n+1. - Joseph L. Pe, Sep 11 2002
a(n+1) is also the number of k, 0 <= k <= n-1, such that C(n,k) divides C(n,k+1). - Benoit Cloitre, Oct 17 2002
a(n+1) is also the number of factors of the n-th degree polynomial x^n + x^(n-1) + x^(n-2) + ... + x^2 + x + 1. Example: 1 + x + x^2 + x^3 = (1+x)(1+x^2) implies a(4)=2.
a(n) is also the number of factors of the n-th Fibonacci polynomial. - T. D. Noe, Mar 09 2006
Number of partitions of n into 2 parts with the second dividing the first. - Franklin T. Adams-Watters, Sep 20 2006
Number of partitions of n+1 into exactly one q and at least one q+1. Example: a(12)=5; indeed, we have 13 = 7 + 6 = 5 + 4 + 4 = 4 + 3 + 3 + 3 = 3 + 2 + 2 + 2 + 2 + 2 = 2 + 11*1.
Differences of A002541. - George Beck, Feb 12 2012
For n > 1: number of ones in row n+1 of triangle A051778. - Reinhard Zumkeller, Dec 03 2014
For n > 0, a(n) is the number of strong divisors of n. - Omar E. Pol, May 03 2015
a(n) is also the number of factors of the (n-1)-th degree polynomial ((x+1)^n-1)/x. Example: for n=6, ((x+1)^6-1)/x = x^5 + 6*x^4 + 15*x^3 + 20*x^2 + 15*x + 6 = (2+x)(1+x+x^2)(3+3x+x^2) implies a(6)=3. - Federico Provvedi, Oct 09 2018
Consider the polynomial P(n,z) = Sum_{i=1..q} d(i)*z^(i-1) where d(1), d(2), ..., d(q) are are the q ordered divisors of n. The sequence lists the numbers of zeros of P(n,z) strictly inside the unit circle. - Michel Lagneau, Apr 06 2025

			a(6) = 3 since the proper divisors of 6 are 1, 2, 3.

André Weil, Number Theory, An approach through history, From Hammurapi to Legendre, Birkhäuser, 1984, page 5.

T. D. Noe, Table of n, a(n) for n = 0..10000
Eric Weisstein's World of Mathematics, Proper divisors.

Column 2 of A122934.
Cf. A003238, A001065, A027749, A027751 (list of proper divisors).
Cf. A001620, A039653, A051778, A070824, A091818.

Concatenation([0],List([1..100],n->Tau(n)-1)); # Muniru A Asiru, Oct 09 2018

a032741 n = if n == 0 then 0 else a000005 n - 1
-- Reinhard Zumkeller, Jul 31 2014

A032741 := proc(n)
    if n = 0 then
        0 ;
    else
        numtheory[tau](n)-1 ;
    end if;
end proc: # R. J. Mathar, Feb 03 2013

Prepend[DivisorSigma[0, Range[99]]-1, 0] (* Jayanta Basu, May 25 2013 *)

PARI
```
a(n) = if(n<1,0,numdiv(n)-1)
```

{a(n)=polcoeff(2*sum(m=1,n\2+1,sumdiv(m,d,log(1-x^(m/d) +x*O(x^n) )^(2*d)/(2*d)!)), n)} \\ Paul D. Hanna, Aug 21 2014

from sympy import divisor_count
def A032741(n): return divisor_count(n)-1 if n else 0 # Chai Wah Wu, Mar 14 2023

a(n) = tau(n)-1 = A000005(n)-1. Cf. A039653.
G.f.: Sum_{n>=1} x^(2*n)/(1-x^n). - Michael Somos, Apr 29 2003
G.f.: Sum_{i>=1} (1-x^i+x^(2*i))/(1-x^i). - Jon Perry, Jul 03 2004
a(n) = Sum_{k=1..floor(n/2)} A051731(n-k,k). - Reinhard Zumkeller, Nov 01 2009
G.f.: 2*Sum_{n>=1} Sum_{d|n} log(1 - x^(n/d))^(2*d) / (2*d)!. - Paul D. Hanna, Aug 21 2014
Dirichlet g.f.: zeta(s)*(zeta(s)-1). - Geoffrey Critzer, Dec 06 2014
a(n) = Sum_{k=1..n-1} binomial((n-1) mod k, k-1). - Wesley Ivan Hurt, Sep 26 2016
a(n) = Sum_{i=1..n-1} floor(n/i)-floor((n-1)/i). - Wesley Ivan Hurt, Nov 15 2017
a(n) = Sum_{i=1..n-1} 1-sign(i mod (n-i)). - Wesley Ivan Hurt, Sep 27 2018
Sum_{k=1..n} a(k) ~ n*log(n) + 2*(gamma - 1)*n, where gamma is Euler's constant (A001620). - Amiram Eldar, Nov 27 2022

Typos in definition corrected by Omar E. Pol, Dec 13 2008

A049820 a(n) = n - d(n), where d(n) is the number of divisors of n (A000005).

Original entry on oeis.org

0, 0, 1, 1, 3, 2, 5, 4, 6, 6, 9, 6, 11, 10, 11, 11, 15, 12, 17, 14, 17, 18, 21, 16, 22, 22, 23, 22, 27, 22, 29, 26, 29, 30, 31, 27, 35, 34, 35, 32, 39, 34, 41, 38, 39, 42, 45, 38, 46, 44, 47, 46, 51, 46, 51, 48, 53, 54, 57, 48, 59, 58, 57, 57, 61, 58, 65

Offset: 1

Clark Kimberling

a(n) is the number of non-divisors of n in 1..n. - Jaroslav Krizek, Nov 14 2009
Also equal to the number of partitions p of n such that max(p)-min(p) = 1. The number of partitions of n with max(p)-min(p) <= 1 is n; there is one with k parts for each 1 <= k <= n. max(p)-min(p) = 0 iff k divides n, leaving n-d(n) with a difference of 1. It is easiest to see this by looking at fixed k with increasing n: for k=3, starting with n=3 the partitions are [1,1,1], [2,1,1], [2,2,1], [2,2,2], [3,2,2], etc. - Giovanni Resta, Feb 06 2006 and Franklin T. Adams-Watters, Jan 30 2011
Number of positive numbers in n-th row of array T given by A049816.
Number of proper non-divisors of n. - Omar E. Pol, May 25 2010
a(n+2) is the sum of the n-th antidiagonal of A225145. - Richard R. Forberg, May 02 2013
For n > 2, number of nonzero terms in n-th row of triangle A051778. - Reinhard Zumkeller, Dec 03 2014
Number of partitions of n of the form [j,j,...,j,i] (j > i). Example: a(7)=5 because we have [6,1], [5,2], [4,3], [3,3,1], and [2,2,2,1]. - Emeric Deutsch, Sep 22 2016

			a(7) = 5; the 5 non-divisors of 7 in 1..7 are 2, 3, 4, 5, and 6.
The 5 partitions of 7 with max(p) - min(p) = 1 are [4,3], [3,2,2], [2,2,2,1], [2,2,1,1,1] and [2,1,1,1,1,1]. - _Emeric Deutsch_, Mar 01 2006

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
G. E. Andrews, M. Beck, N. Robbins, Partitions with fixed differences between largest and smallest parts, arXiv preprint arXiv:1406.3374 [math.NT], 2014-2015.

Cf. A000005.
One less than A062968, two less than A059292.
Cf. A161664 (partial sums).
Cf. A060990 (number of solutions to a(x) = n).
Cf. A045765 (numbers not occurring in this sequence).
Cf. A236561 (same sequence sorted into ascending order), A236562 (with also duplicates removed), A236565, A262901 and A262903.
Cf. A262511 (numbers that occur only once).
Cf. A055927 (positions of repeated terms).
Cf. A245388 (positions of squares).
Cf. A155043 (number of steps needed to reach zero when iterating a(n)), A262680 (number of nonzero squares encountered).
Cf. A259934 (an infinite trunk of the tree defined by edge-relation a(child) = parent, conjectured to be unique).
Cf. tables and arrays A047916, A051731, A051778, A173540, A173541.
Cf. also arrays A225145, A262898, A263255 and tables A263265, A263267.
Other related sequences: A006218, A006590, A051953, A070824, A094181, A062249, A067391, A076627, A128508, A131187, A134156, A140826, A161886, A177235, A177236, A227874, A228453, A230653, A230654, A231167, A245197, A253473.

List([1..80],n->n-Tau(n)); # Muniru A Asiru, Sep 28 2018

a049820 n = n - a000005 n  -- Reinhard Zumkeller, Feb 06 2012

A049820 := n->n-numtheory[tau](n):
seq(A049820(n), n=1..100);

Table[n - DivisorSigma[0, n], {n, 100}] (* Wesley Ivan Hurt, Nov 19 2014 *)
Array[(# - DivisorSigma[0, #])&, 70] (* Vincenzo Librandi, Dec 29 2015 *)

PARI
```
a(n)=n-numdiv(n)
```

(define (A049820 n) (- n (A000005 n))) ;; Antti Karttunen, Nov 27 2015

a(n) = Sum_{k=1..n} ceiling(n/k)-floor(n/k). - Benoit Cloitre, May 11 2003
G.f.: Sum_{k>0} x^(2*k+1)/(1-x^k)/(1-x^(k+1)). - Emeric Deutsch, Mar 01 2006
a(n) = A006590(n) - A006218(n) = A161886(n) - A000005(n) - A006218(n) + 1 for n >= 1. - Jaroslav Krizek, Nov 14 2009
a(n) = Sum_{k=1..n} A000007(A051731(n,k)). - Reinhard Zumkeller, Mar 09 2010
a(n) = A076627(n) / A000005(n). - Reinhard Zumkeller, Feb 06 2012
For n >= 2, a(n) = A094181(n) / A051953(n). - Antti Karttunen, Nov 27 2015
a(n) = Sum_{k=1..n} ((n mod k) + (-n mod k))/k. - Wesley Ivan Hurt, Dec 28 2015
G.f.: Sum_{j>=2} (x^(j+1)*(1-x^(j-1))/(1-x^j))/(1-x). - Emeric Deutsch, Sep 22 2016
Dirichlet g.f.: zeta(s-1)- zeta(s)^2. - Ilya Gutkovskiy, Apr 12 2017
a(n) = Sum_{i=1..n-1} sign(i mod n-i). - Wesley Ivan Hurt, Sep 27 2018

Edited by Franklin T. Adams-Watters, Jan 30 2012

A356224 Number of divisors of n whose prime indices cover an initial interval of positive integers.

Original entry on oeis.org

1, 2, 1, 3, 1, 3, 1, 4, 1, 2, 1, 5, 1, 2, 1, 5, 1, 4, 1, 3, 1, 2, 1, 7, 1, 2, 1, 3, 1, 4, 1, 6, 1, 2, 1, 7, 1, 2, 1, 4, 1, 3, 1, 3, 1, 2, 1, 9, 1, 2, 1, 3, 1, 5, 1, 4, 1, 2, 1, 7, 1, 2, 1, 7, 1, 3, 1, 3, 1, 2, 1, 10, 1, 2, 1, 3, 1, 3, 1, 5, 1, 2, 1, 5, 1, 2, 1

Offset: 1

Gus Wiseman, Aug 04 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The a(n) gapless divisors of n = 1..24:
  1  2  1  4  1  6  1  8  1  2  1  12  1  2  1  16  1  18  1  4  1  2  1  24
     1     2     2     4     1     6      1     8      6      2     1     12
           1     1     2           4            4      2      1           8
                       1           2            2      1                  6
                                   1            1                         4
                                                                          2
                                                                          1
For example, the divisors of 12 are {1,2,3,4,6,12}, of which {1,2,4,6,12} belong to A055932, so a(12) = 5.

These divisors belong to A055932, a subset of A073491 (complement A073492).
The complement is A356225.
A001223 lists the prime gaps.
A328338 has third-largest divisor prime.
A356226 gives the lengths of maximal gapless intervals of prime indices.
Cf. A000005, A001222, A028334, A029709, A055874, A056239, A070824, A112798, A119313, A137921, A287170, A289509, A356223.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
normQ[m_]:=m=={}||Union[m]==Range[Max[m]];
Table[Length[Select[Divisors[n],normQ[primeMS[#]]&]],{n,100}]

A251683 Irregular triangular array: T(n,k) is the number of ordered factorizations of n with exactly k factors, n >= 1, 1 <= k <= A086436(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 2, 1, 2, 1, 3, 3, 1, 1, 1, 4, 3, 1, 1, 4, 3, 1, 2, 1, 2, 1, 1, 6, 9, 4, 1, 1, 1, 2, 1, 2, 1, 1, 4, 3, 1, 1, 6, 6, 1, 1, 4, 6, 4, 1, 1, 2, 1, 2, 1, 2, 1, 7, 12, 6, 1, 1, 2, 1, 2, 1, 6, 9, 4

Offset: 1

Geoffrey Critzer, Dec 06 2014

Row sums = A074206.
Row lengths give A086436.
T(n,2) = A070824(n).
T(n,3) = A200221(n).
Sum_{k>=1} k*T(n,k) = A254577.
For all n > 1,  Sum_{k=1..A086436(n)} (-1)^k*T(n,k) = A008683(n). - Geoffrey Critzer, May 25 2018
From Gus Wiseman, Aug 21 2020: (Start)
Also the number of strict length k + 1 chains of divisors from n to 1. For example, row n = 24 counts the following chains:
  24/1  24/2/1   24/4/2/1   24/8/4/2/1
        24/3/1   24/6/2/1   24/12/4/2/1
        24/4/1   24/6/3/1   24/12/6/2/1
        24/6/1   24/8/2/1   24/12/6/3/1
        24/8/1   24/8/4/1
        24/12/1  24/12/2/1
                 24/12/3/1
                 24/12/4/1
                 24/12/6/1
(End)

			Triangle T(n,k) begins:
  1;
  1;
  1;
  1, 1;
  1;
  1, 2;
  1;
  1, 2, 1;
  1, 1;
  1, 2;
  1;
  1, 4, 3;
  1;
  1, 2;
  1, 2;
  ...
There are 8 ordered factorizations of the integer 12: 12, 6*2, 4*3, 3*4, 2*6, 3*2*2, 2*3*2, 2*2*3.  So T(12,1)=1, T(12,2)=4, and T(12,3)=3.

Alois P. Heinz, Rows n = 1..4000, flattened
Jeffery Kline, On the eigenstructure of sparse matrices related to the prime number theorem, Linear Algebra and its Applications (2020) Vol. 584, 409-430.
Arnold Knopfmacher and Michael Mays, Ordered and Unordered Factorizations of Integers, The Mathematica Journal, Vol 10 (1).
Eric Weisstein's World of Mathematics, Ordered Factorization

Cf. A008683, A070824, A200221, A254577.
A008480 gives rows ends.
A086436 gives row lengths.
A124433 is the same except for signs and zeros.
A334996 is the same except for zeros.
A337107 is the restriction to factorial numbers (but with zeros).
A000005 counts divisors.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A074206 counts strict chains of divisors from n to 1.
A067824 counts strict chains of divisors starting with n.
A122651 counts strict chains of divisors summing to n.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict nonempty chains of divisors of n.
A337071 counts strict chains of divisors starting with n!.
A337256 counts strict chains of divisors of n.
Cf. A001221, A002033, A124010, A167865, A337070, A337105.

with(numtheory):
b:= proc(n) option remember; expand(x*(1+
      add(b(n/d), d=divisors(n) minus {1, n})))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=1..degree(p)))(b(n)):
seq(T(n), n=1..100);  # Alois P. Heinz, Dec 07 2014

f[1] = {{}};
f[n_] := f[n] =
  Level[Table[
    Map[Prepend[#, d] &, f[n/d]], {d, Rest[Divisors[n]]}], {2}];
Prepend[Map[Select[#, # > 0 &] &,
  Drop[Transpose[
    Table[Map[Count[#, k] &,
      Map[Length, Table[f[n], {n, 1, 40}], {2}]], {k, 1, 10}]],
   1]],{1}] // Grid
(* Second program: *)
b[n_] := b[n] = x(1+Sum[b[n/d], {d, Divisors[n]~Complement~{1, n}}]);
T[n_] := CoefficientList[b[n]/x, x];
Array[T, 100] // Flatten (* Jean-François Alcover, Nov 17 2020, after Alois P. Heinz *)

Dirichlet g.f.: 1/(1 - y*(zeta(x)-1)).

A130130 a(0)=0, a(1)=1, a(n)=2 for n >= 2.

Original entry on oeis.org

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

Offset: 0

Paul Curtz, Aug 01 2007

a(n) is also total number of positive integers below 10^(n+1) requiring 9 positive cubes in their representation as sum of cubes (cf. Dickson, 1939).
A061439(n) + A181375(n) + A181377(n) + A181379(n) + A181381(n) + A181400(n) + A181402(n) + A181404(n) + a(n) = A002283(n).
a(n) = number of obvious divisors of n. The obvious divisors of n are the numbers 1 and n. - Jaroslav Krizek, Mar 02 2009
Number of colors needed to paint n adjacent segments on a line. - Jaume Oliver Lafont, Mar 20 2009
a(n) = ceiling(n-th nonprimes/n) = ceiling(A018252(n)/A000027(n)) for n >= 1. Numerators of (A018252(n)/A000027(n)) in A171529(n), denominators of (A018252(n)/A000027(n)) in A171530(n). a(n) = A171624(n) + 1 for n >= 5. - Jaroslav Krizek, Dec 13 2009
a(n) is also the continued fraction for sqrt(1/2). - Enrique Pérez Herrero, Jul 12 2010
For n >= 1, a(n) = minimal number of divisors of any n-digit number. See A066150 for maximal number of divisors of any n-digit number. - Jaroslav Krizek, Jul 18 2010
Central terms in the triangle A051010. - Reinhard Zumkeller, Jun 27 2013
Decimal expansion of 11/900. - Elmo R. Oliveira, May 05 2024

Leonard Eugene Dickson, All integers except 23 and 239 are sums of eight cubes, Bulletin of the American Mathematical Society 45 (1939), p. 588-591.
Eric Weisstein's World of Mathematics, Waring's Problem.
Index entries for linear recurrences with constant coefficients, signature (1).

Cf. A000005, A000027, A002283, A018252, A051010, A061439, A066150, A070824, A158411, A171529, A171530, A171624, A181375, A181377, A181379, A181381, A181400, A181402, A181404.

a130130 = min 2
a130130_list = 0 : 1 : repeat 2  -- Reinhard Zumkeller, Jun 27 2013

A130130[0]:=0; A130130[1]:=1; A130130[n_]:=2; (* Enrique Pérez Herrero, Jul 12 2010 *)
A130130[n_]:=ContinuedFraction[Sqrt[1/2],n+1][[n+1]] (* Enrique Pérez Herrero, Jul 12 2010 *)
Join[{0, 1},LinearRecurrence[{1},{2},96]] (* Ray Chandler, Sep 23 2015 *)
PadRight[{0,1},120,{2}] (* Harvey P. Dale, Sep 15 2022 *)

a(n)=min(n,2) \\ Charles R Greathouse IV, Jun 01 2011

G.f.: x*(1+x)/(1-x) = x*(1-x^2)/(1-x)^2. - Jaume Oliver Lafont, Mar 20 2009
a(n) = A000005(n) - A070824(n) for n >= 1.
E.g.f.: 2*exp(x) - x - 2. - Stefano Spezia, May 19 2024

A356225 Number of divisors of n whose prime indices do not cover an initial interval of positive integers.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 13 2022

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The a(70) = 6 divisors: 5, 7, 10, 14, 35, 70.

These divisors belong to the complement of A055932, a subset of A073491.
These divisors belong to A080259, a superset of A073492.
The complement is counted by A356224.
A001223 lists the prime gaps.
A328338 has third-largest divisor prime, smallest A119313.
A356226 gives the lengths of maximal gapless intervals of prime indices.
Cf. A000005, A001222, A028334, A055874, A056239, A070824, A112798, A137921, A287170, A356233-A356237.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
normQ[m_]:=m=={}||Union[m]==Range[Max[m]];
Table[Length[Select[Divisors[n],!normQ[primeMS[#]]&]],{n,100}]

a(n) = A000005(n) - A356224(n).

cp's OEIS Frontend

A354697 a(n) is the least integer that can be written in two or more ways as the product of the integers in two subsets of its A070824(a(n)) nontrivial divisors, each of size n and with empty intersection.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A307118 a(1) = 0; for n>1, a(n) = dr(n-1) + dr(n) + dr(n+1), where dr(n) is the number of nontrivial divisors of n (A070824).

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

PARI

A147588 Duplicate of A070824.

Views

Author

Keywords

A010051 Characteristic function of primes: 1 if n is prime, else 0.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

Python

Formula

A032741 a(0) = 0; for n > 0, a(n) = number of proper divisors of n (divisors of n which are less than n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A049820 a(n) = n - d(n), where d(n) is the number of divisors of n (A000005).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Haskell

Maple

Mathematica

PARI

Scheme

Formula

Extensions

A356224 Number of divisors of n whose prime indices cover an initial interval of positive integers.

Views

Author