A054973 -id:A054973 - cp's OEIS frontend

A159953 Values in A054973 larger than 1.

2, 2, 3, 2, 2, 3, 3, 2, 2, 3, 5, 2, 3, 3, 4, 2, 2, 2, 2, 4, 2, 2, 2, 2, 2, 5, 2, 2, 6, 4, 2, 2, 5, 2, 5, 3, 3, 3, 7, 3, 6, 2, 3, 2, 2, 6, 3, 2, 4, 2, 3, 8, 2, 9, 4, 2, 6, 2, 2, 2, 2, 2, 2, 4, 8, 4, 2, 2, 2, 3, 4, 3, 9, 2, 10, 2, 3, 2, 4, 4, 3, 4, 2, 2, 11, 5, 2, 5, 2, 3, 4, 2, 2, 3, 5, 3, 8, 7, 4, 15, 2, 4, 7, 8

Offset: 1

Jaroslav Krizek, Apr 27 2009

nonn

This is a survey of how many solutions the equation sigma(x)=k has for k in A159886, or about the lengths of the plateaus in A007609.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1095 from Jean-François Alcover)
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).

Cf. A000203, A007609, A054973, A159886.

read("transforms3") ; a054973 := BFILETOLIST("b054973.txt") ;
for i from 1 to 1000 do if op(i,a054973) > 1 then printf("%d,", op(i,a054973)) ; fi; od: # R. J. Mathar, May 22 2009

b[n_] := Sum[Boole[DivisorSigma[1, k] == n], {k, 1, n}];
Select[Array[b, 1000], # > 1&] (* Jean-François Alcover, Apr 06 2020 *)

list(lim) = {my(s); for(k = 1, lim, s = invsigmaNum(k); if(s > 1, print1(s, ", ")));} \\ Amiram Eldar, Dec 25 2024, using Max Alekseyev's invphi.gp

Edited and extended by R. J. Mathar, May 22 2009

A258968 a(n) is the least positive integer x with A054973(x) = n and A258913(x) < x.

Original entry on oeis.org

2, 3, 124, 10714158

Offset: 0

Jeppe Stig Nielsen, Jun 15 2015

Is a(n) well-defined for all n?
For n > 1, a(n) belongs to sequence A258931.
The terms a(k), for 4 <= k <= 100, if they exist, are larger than 3*10^10. - Giovanni Resta, Jun 15 2015

			For n=2, sigma(x)=124 for x=48 and 75, and 48+75 = 123 < 124.
For n=3, sigma(x)=10714158 for x=3031200, 3417300, and 3987450; and their sum is 10435950 (<10714158).

Cf. A054973, A258913, A258931, A007368, A258914.

a(n)=x=0;while(x++,u=List();for(i=1,x,if(sigma(i)==x,listput(u,i)));if(#u==n&vecsum(Vec(u))

A000203 a(n) = sigma(n), the sum of the divisors of n. Also called sigma_1(n).

Original entry on oeis.org

1, 3, 4, 7, 6, 12, 8, 15, 13, 18, 12, 28, 14, 24, 24, 31, 18, 39, 20, 42, 32, 36, 24, 60, 31, 42, 40, 56, 30, 72, 32, 63, 48, 54, 48, 91, 38, 60, 56, 90, 42, 96, 44, 84, 78, 72, 48, 124, 57, 93, 72, 98, 54, 120, 72, 120, 80, 90, 60, 168, 62, 96, 104, 127, 84, 144, 68, 126, 96, 144

Offset: 1

N. J. A. Sloane

Multiplicative: If the canonical factorization of n into prime powers is the product of p^e(p) then sigma_k(n) = Product_p ((p^((e(p)+1)*k))-1)/(p^k-1).
Sum_{d|n} 1/d^k is equal to sigma_k(n)/n^k. So sequences A017665-A017712 also give the numerators and denominators of sigma_k(n)/n^k for k = 1..24. The power sums sigma_k(n) are in sequences A000203 (this sequence) (k=1), A001157-A001160 (k=2,3,4,5), A013954-A013972 for k = 6,7,...,24. - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 05 2001
A number n is abundant if sigma(n) > 2n (cf. A005101), perfect if sigma(n) = 2n (cf. A000396), deficient if sigma(n) < 2n (cf. A005100).
a(n) is the number of sublattices of index n in a generic 2-dimensional lattice. - Avi Peretz (njk(AT)netvision.net.il), Jan 29 2001 [In the language of group theory, a(n) is the number of index-n subgroups of Z x Z. - Jianing Song, Nov 05 2022]
The sublattices of index n are in one-to-one correspondence with matrices [a b; 0 d] with a>0, ad=n, b in [0..d-1]. The number of these is Sum_{d|n} d = sigma(n), which is a(n). A sublattice is primitive if gcd(a,b,d) = 1; the number of these is n * Product_{p|n} (1+1/p), which is A001615. [Cf. Grady reference.]
Sum of number of common divisors of n and m, where m runs from 1 to n. - Naohiro Nomoto, Jan 10 2004
a(n) is the cardinality of all extensions over Q_p with degree n in the algebraic closure of Q_p, where p>n. - Volker Schmitt (clamsi(AT)gmx.net), Nov 24 2004. Cf. A100976, A100977, A100978 (p-adic extensions).
Let s(n) = a(n-1) + a(n-2) - a(n-5) - a(n-7) + a(n-12) + a(n-15) - a(n-22) - a(n-26) + ..., then a(n) = s(n) if n is not pentagonal, i.e., n != (3 j^2 +- j)/2 (cf. A001318), and a(n) is instead s(n) - ((-1)^j)*n if n is pentagonal. - Gary W. Adamson, Oct 05 2008 [corrected Apr 27 2012 by William J. Keith based on Ewell and by Andrey Zabolotskiy, Apr 08 2022]
Write n as 2^k * d, where d is odd. Then a(n) is odd if and only if d is a square. - Jon Perry, Nov 08 2012
Also total number of parts in the partitions of n into equal parts. - Omar E. Pol, Jan 16 2013
Note that sigma(3^4) = 11^2. On the other hand, Kanold (1947) shows that the equation sigma(q^(p-1)) = b^p has no solutions b > 2, q prime, p odd prime. - N. J. A. Sloane, Dec 21 2013, based on postings to the Number Theory Mailing List by Vladimir Letsko and Luis H. Gallardo
Limit_{m->infinity} (Sum_{n=1..prime(m)} a(n)) / prime(m)^2 = zeta(2)/2 = Pi^2/12 (A072691). See more at A244583. - Richard R. Forberg, Jan 04 2015
a(n) + A000005(n) is an odd number iff n = 2m^2, m>=1. - Richard R. Forberg, Jan 15 2015
a(n) = a(n+1) for n = 14, 206, 957, 1334, 1364 (A002961). - Zak Seidov, May 03 2016
Equivalent to the Riemann hypothesis: a(n) < H(n) + exp(H(n))*log(H(n)), for all n>1, where H(n) is the n-th harmonic number (Jeffrey Lagarias). See A057641 for more details. - Ilya Gutkovskiy, Jul 05 2016
a(n) is the total number of even parts in the partitions of 2*n into equal parts. More generally, a(n) is the total number of parts congruent to 0 mod k in the partitions of k*n into equal parts (the comment dated Jan 16 2013 is the case for k = 1). - Omar E. Pol, Nov 18 2019
From Jianing Song, Nov 05 2022: (Start)
a(n) is also the number of order-n subgroups of C_n X C_n, where C_n is the cyclic group of order n. Proof: by the correspondence theorem in the group theory, there is a one-to-one correspondence between the order-n subgroups of C_n X C_n = (Z x Z)/(nZ x nZ) and the index-n subgroups of Z x Z containing nZ x nZ. But an index-n normal subgroup of a (multiplicative) group G contains {g^n : n in G} automatically. The desired result follows from the comment from Naohiro Nomoto above.
The number of subgroups of C_n X C_n that are isomorphic to C_n is A001615(n). (End)

			For example, 6 is divisible by 1, 2, 3 and 6, so sigma(6) = 1 + 2 + 3 + 6 = 12.
Let L = <V,W> be a 2-dimensional lattice. The 7 sublattices of index 4 are generated by <4V,W>, <V,4W>, <4V,W+-V>, <2V,2W>, <2V+W,2W>, <2V,2W+V>. Compare A001615.

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 840.
T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 38.
A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 116ff.
Florian Cajori, A History of Mathematical Notations, Dover edition (2012), par. 407.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 162, #16, (6), 2nd formula.
G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, AMS Chelsea Publishing, Providence, Rhode Island, 2002, pp. 141, 166.
H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fifth Edition, Clarendon Press, Oxford, 2003.
Ross Honsberger, "Mathematical Gems, Number One," The Dolciani Mathematical Expositions, Published and Distributed by The Mathematical Association of America, page 116.
Kanold, Hans Joachim, Kreisteilungspolynome und ungerade vollkommene Zahlen. (German), Ber. Math.-Tagung Tübingen 1946, (1947). pp. 84-87.
M. Krasner, Le nombre des surcorps primitifs d'un degré donné et le nombre des surcorps métagaloisiens d'un degré donné d'un corps de nombres p-adiques. Comptes Rendus Hebdomadaires, Académie des Sciences, Paris 254, 255, 1962.
A. Lubotzky, Counting subgroups of finite index, Proceedings of the St. Andrews/Galway 93 group theory meeting, Th. 2.1. LMS Lecture Notes Series no. 212 Cambridge University Press 1995.
D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, Section III.1, page 77.
G. Pólya, Induction and Analogy in Mathematics, vol. 1 of Mathematics and Plausible Reasoning, Princeton Univ Press 1954, page 92.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 91, 395.
Robert M. Young, Excursions in Calculus, The Mathematical Association of America, 1992 p. 361.

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 20000 terms from N. J. A. Sloane)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
B. Apostol and L. Petrescu, Extremal Orders of Certain Functions Associated with Regular Integers (mod n), Journal of Integer Sequences, 2013, # 13.7.5.
Joerg Arndt, On computing the generalized Lambert series, arXiv:1202.6525v3 [math.CA], (2012).
M. Baake and U. Grimm, Quasicrystalline combinatorics
Henry Bottomley, Illustration of initial terms
C. K. Caldwell, The Prime Glossary, sigma function
Imanuel Chen and Michael Z. Spivey, Integral Generalized Binomial Coefficients of Multiplicative Functions, Preprint 2015; Summer Research Paper 238, Univ. Puget Sound.
D. Christopher and T. Nadu, Partitions with Fixed Number of Sizes, Journal of Integer Sequences, 15 (2015), #15.11.5.
J. N. Cooper and A. W. N. Riasanovsky, On the Reciprocal of the Binary Generating Function for the Sum of Divisors, 2012. - _N. J. A. Sloane_, Dec 25 2012
R. L. Duncan, Note on the Divisors of a Number, The American Mathematical Monthly, 68(4), 356-359, (1961).
Jason Earls, The Smarandache sum of composites between factors function, in Smarandache Notions Journal (2004), Vol. 14.1, page 243.
L. Euler, Discovery of a most extraordinary law of numbers, relating to the sum of their divisors
L. Euler, Observatio de summis divisorum
L. Euler, An observation on the sums of divisors, arXiv:math/0411587 [math.HO], 2004-2009.
J. A. Ewell, Recurrences for the sum of divisors, Proc. Amer. Math. Soc. 64 (2) 1977.
Farideh Firoozbakht and M. F. Hasler, Variations on Euclid's formula for Perfect Numbers, JIS 13 (2010) #10.3.1.
Daniele A. Gewurz and Francesca Merola, Sequences realized as Parker vectors ..., J. Integer Seqs., Vol. 6, 2003.
Johan Gielis and Ilia Tavkhelidze, The general case of cutting of GML surfaces and bodies, arXiv:1904.01414 [math.GM], 2019.
J. W. L. Glaisher, On the function chi(n), Quarterly Journal of Pure and Applied Mathematics, 20 (1884), 97-167.
J. W. L. Glaisher, On the function chi(n), Quarterly Journal of Pure and Applied Mathematics, 20 (1884), 97-167. [Annotated scanned copy]
M. J. Grady, A group theoretic approach to a famous partition formula, Amer. Math. Monthly, 112 (No. 7, 2005), 645-651.
Masazumi Honda and Takuya Yoda, String theory, N = 4 SYM and Riemann hypothesis, arXiv:2203.17091 [hep-th], 2022.
Douglas E. Iannucci, On sums of the small divisors of a natural number, arXiv:1910.11835 [math.NT], 2019.
Antti Karttunen, Scheme program for computing this sequence.
P. A. MacMahon, Divisors of numbers and their continuations in the theory of partitions, Proc. London Math. Soc., 19 (1921), 75-113.
M. Maia and M. Mendez, On the arithmetic product of combinatorial species, arXiv:math.CO/0503436, 2005.
K. Matthews, Factorizing n and calculating phi(n),omega(n),d(n),sigma(n) and mu(n)
Walter Nissen, Abundancy : Some Resources
P. Pollack and C. Pomerance, Some problems of Erdős on the sum-of-divisors function, For Richard Guy on his 99th birthday: May his sequence be unbounded, Trans. Amer. Math. Soc. Ser. B 3 (2016), 1-26; errata.
Carl Pomerance and Hee-Sung Yang, Variant of a theorem of Erdős on the sum-of-proper-divisors function, Math. Comp. 83 (2014), 1903-1913.
John S. Rutherford, The enumeration and symmetry-significant properties of derivative lattices, Act. Cryst. (1992) A48, 500-508. - _N. J. A. Sloane_, Mar 14 2009
John S. Rutherford, The enumeration and symmetry-significant properties of derivative lattices II, Acta Cryst. A49 (1993), 293-300. - _N. J. A. Sloane_, Mar 14 2009
John S. Rutherford, Sublattice enumeration. IV. Equivalence classes of plane sublattices by parent Patterson symmetry and colour lattice group type, Acta Cryst. (2009). A65, 156-163. [See Table 1]. - _N. J. A. Sloane_, Feb 23 2009
N. J. A. Sloane, "A Handbook of Integer Sequences" Fifty Years Later, arXiv:2301.03149 [math.NT], 2023, p. 3.
Eric Weisstein's World of Mathematics, Divisor Function
Wikipedia, Divisor function
Index entries for sequences related to sublattices
Index entries for sequences related to sigma(n)
Index entries for "core" sequences

See A034885, A002093 for records. Bisections give A008438, A062731. Values taken are listed in A007609. A054973 is an inverse function.
For partial sums see A024916.
Row sums of A127093.
sigma_i (i=0..25): A000005, A000203, A001157, A001158, A001159, A001160, A013954, A013955, A013956, A013957, A013958, A013959, A013960, A013961, A013962, A013963, A013964, A013965, A013966, A013967, A013968, A013969, A013970, A013971, A013972, A281959.
Cf. A144736, A158951, A158902, A174740, A147843, A001158, A001160, A001065, A002192, A001001, A001615 (primitive sublattices), A039653, A088580, A074400, A083728, A006352, A002659, A083238, A000593, A050449, A050452, A051731, A027748, A124010, A069192, A057641, A001318.
Cf. A009194, A082062 (gcd(a(n),n) and its largest prime factor), A179931, A192795 (gcd(a(n),A001157(n)) and largest prime factor).
Cf. also A034448 (sum of unitary divisors).
Cf. A007955 (products of divisors).
Cf. A144613, A144614, A144615, A146076.
A001227, A000593 and this sequence have the same parity: A053866. - Omar E. Pol, May 14 2016
Cf. A001221, A054533.

A000203:=List([1..10^2],n->Sigma(n)); # Muniru A Asiru, Oct 01 2017

a000203 n = product $ zipWith (\p e -> (p^(e+1)-1) `div` (p-1)) (a027748_row n) (a124010_row n)
-- Reinhard Zumkeller, May 07 2012

Magma
```
[SumOfDivisors(n): n in [1..70]];
```

[DivisorSigma(1,n): n in [1..70]]; // Bruno Berselli, Sep 09 2015

with(numtheory): A000203 := n->sigma(n); seq(A000203(n), n=1..100);

Table[ DivisorSigma[1, n], {n, 100}]
a[ n_] := SeriesCoefficient[ QPolyGamma[ 1, 1, q] / Log[q]^2, {q, 0, n}]; (* Michael Somos, Apr 25 2013 *)

makelist(divsum(n),n,1,1000); /* Emanuele Munarini, Mar 26 2011 */

numlib::sigma(n)$ n=1..81 // Zerinvary Lajos, May 13 2008

PARI
```
{a(n) = if( n<1, 0, sigma(n))};
```

{a(n) = if( n<1, 0, direuler( p=2, n, 1 / (1 - X) /(1 - p*X))[n])};

{a(n) = if( n<1, 0, polcoeff( sum( k=1, n, x^k / (1 - x^k)^2, x * O(x^n)), n))}; /* Michael Somos, Jan 29 2005 */

max_n = 30; ser = - sum(k=1,max_n,log(1-x^k)); a(n) = polcoeff(ser,n)*n \\ Gottfried Helms, Aug 10 2009

from sympy import divisor_sigma
def a(n): return divisor_sigma(n, 1)
print([a(n) for n in range(1, 71)]) # Michael S. Branicky, Jan 03 2021

from math import prod
from sympy import factorint
def a(n): return prod((p**(e+1)-1)//(p-1) for p, e in factorint(n).items())
print([a(n) for n in range(1, 51)]) # Michael S. Branicky, Feb 25 2024
(APL, Dyalog dialect) A000203 ← +/{ð←⍵{(0=⍵|⍺)/⍵}⍳⌊⍵*÷2 ⋄ 1=⍵:ð ⋄ ð,(⍵∘÷)¨(⍵=(⌊⍵*÷2)*2)↓⌽ð} ⍝ Antti Karttunen, Feb 20 2024

[sigma(n, 1) for n in range(1, 71)]  # Zerinvary Lajos, Jun 04 2009

(definec (A000203 n) (if (= 1 n) n (let ((p (A020639 n)) (e (A067029 n))) (* (/ (- (expt p (+ 1 e)) 1) (- p 1)) (A000203 (A028234 n)))))) ;; Uses macro definec from http://oeis.org/wiki/Memoization#Scheme - Antti Karttunen, Nov 25 2017

(define (A000203 n) (let ((r (sqrt n))) (let loop ((i (inexact->exact (floor r))) (s (if (integer? r) (- r) 0))) (cond ((zero? i) s) ((zero? (modulo n i)) (loop (- i 1) (+ s i (/ n i)))) (else (loop (- i 1) s)))))) ;; (Stand-alone program) - Antti Karttunen, Feb 20 2024

Multiplicative with a(p^e) = (p^(e+1)-1)/(p-1). - David W. Wilson, Aug 01 2001
For the following bounds and many others, see Mitrinovic et al. - N. J. A. Sloane, Oct 02 2017
If n is composite, a(n) > n + sqrt(n).
a(n) < n*sqrt(n) for all n.
a(n) < (6/Pi^2)*n^(3/2) for n > 12.
G.f.: -x*deriv(eta(x))/eta(x) where eta(x) = Product_{n>=1} (1-x^n). - Joerg Arndt, Mar 14 2010
L.g.f.: -log(Product_{j>=1} (1-x^j)) = Sum_{n>=1} a(n)/n*x^n. - Joerg Arndt, Feb 04 2011
Dirichlet convolution of phi(n) and tau(n), i.e., a(n) = sum_{d|n} phi(n/d)*tau(d), cf. A000010, A000005.
a(n) is odd iff n is a square or twice a square. - Robert G. Wilson v, Oct 03 2001
a(n) = a(n*prime(n)) - prime(n)*a(n). - Labos Elemer, Aug 14 2003 (Clarified by Omar E. Pol, Apr 27 2016)
a(n) = n*A000041(n) - Sum_{i=1..n-1} a(i)*A000041(n-i). - Jon Perry, Sep 11 2003
a(n) = -A010815(n)*n - Sum_{k=1..n-1} A010815(k)*a(n-k). - Reinhard Zumkeller, Nov 30 2003
a(n) = f(n, 1, 1, 1), where f(n, i, x, s) = if n = 1 then s*x else if p(i)|n then f(n/p(i), i, 1+p(i)*x, s) else f(n, i+1, 1, s*x) with p(i) = i-th prime (A000040). - Reinhard Zumkeller, Nov 17 2004
Recurrence: n^2*(n-1)*a(n) = 12*Sum_{k=1..n-1} (5*k*(n-k) - n^2)*a(k)*a(n-k), if n>1. - Dominique Giard (dominique.giard(AT)gmail.com), Jan 11 2005
G.f.: Sum_{k>0} k * x^k / (1 - x^k) = Sum_{k>0} x^k / (1 - x^k)^2. Dirichlet g.f.: zeta(s)*zeta(s-1). - Michael Somos, Apr 05 2003. See the Hardy-Wright reference, p. 312. first equation, and p. 250, Theorem 290. - Wolfdieter Lang, Dec 09 2016
For odd n, a(n) = A000593(n). For even n, a(n) = A000593(n) + A074400(n/2). - Jonathan Vos Post, Mar 26 2006
Equals the inverse Moebius transform of the natural numbers. Equals row sums of A127093. - Gary W. Adamson, May 20 2007
A127093 * [1/1, 1/2, 1/3, ...] = [1/1, 3/2, 4/3, 7/4, 6/5, 12/6, 8/7, ...]. Row sums of triangle A135539. - Gary W. Adamson, Oct 31 2007
a(n) = A054785(2*n) - A000593(2*n). - Reinhard Zumkeller, Apr 23 2008
a(n) = n*Sum_{k=1..n} A060642(n,k)/k*(-1)^(k+1). - Vladimir Kruchinin, Aug 10 2010
Dirichlet convolution of A037213 and A034448. - R. J. Mathar, Apr 13 2011
G.f.: A(x) = x/(1-x)*(1 - 2*x*(1-x)/(G(0) - 2*x^2 + 2*x)); G(k) = -2*x - 1 - (1+x)*k + (2*k+3)*(x^(k+2)) - x*(k+1)*(k+3)*((-1 + (x^(k+2)))^2)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Dec 06 2011
a(n) = A001065(n) + n. - Mats Granvik, May 20 2012
a(n) = A006128(n) - A220477(n). - Omar E. Pol, Jan 17 2013
a(n) = Sum_{k=1..A003056(n)} (-1)^(k-1)*A196020(n,k). - conjectured by Omar E. Pol, Feb 02 2013, and proved by Max Alekseyev, Nov 17 2013
a(n) = Sum_{k=1..A003056(n)} (-1)^(k-1)*A000330(k)*A000716(n-A000217(k)). - Mircea Merca, Mar 05 2014
a(n) = A240698(n, A000005(n)). - Reinhard Zumkeller, Apr 10 2014
a(n) = Sum_{d^2|n} A001615(n/d^2) = Sum_{d^3|n} A254981(n/d^3). - Álvar Ibeas, Mar 06 2015
a(3*n) = A144613(n). a(3*n + 1) = A144614(n). a(3*n + 2) = A144615(n). - Michael Somos, Jul 19 2015
a(n) = Sum{i=1..n} Sum{j=1..i} cos((2*Pi*n*j)/i). - Michel Lagneau, Oct 14 2015
a(n) = A000593(n) + A146076(n). - Omar E. Pol, Apr 05 2016
a(n) = A065475(n) + A048050(n). - Omar E. Pol, Nov 28 2016
a(n) = (Pi^2*n/6)*Sum_{q>=1} c_q(n)/q^2, with the Ramanujan sums c_q(n) given in A054533 as a c_n(k) table. See the Hardy reference, p. 141, or Hardy-Wright, Theorem 293, p. 251. - Wolfdieter Lang, Jan 06 2017
G.f. also (1 - E_2(q))/24, with the g.f. E_2 of A006352. See e.g., Hardy, p. 166, eq. (10.5.5). - Wolfdieter Lang, Jan 31 2017
From Antti Karttunen, Nov 25 2017: (Start)
a(n) = A048250(n) + A162296(n).
a(n) = A092261(n) * A295294(n). [This can be further expanded, see comment in A291750.] (End)
a(n) = A000593(n) * A038712(n). - Ivan N. Ianakiev and Omar E. Pol, Nov 26 2017
a(n) = Sum_{q=1..n} c_q(n) * floor(n/q), where c_q(n) is the Ramanujan's sum function given in A054533. - Daniel Suteu, Jun 14 2018
a(n) = Sum_{k=1..n} gcd(n, k) / phi(n / gcd(n, k)), where phi(k) is the Euler totient function. - Daniel Suteu, Jun 21 2018
a(n) = (2^(1 + (A000005(n) - A001227(n))/(A000005(n) - A183063(n))) - 1)*A000593(n) = (2^(1 + (A183063(n)/A001227(n))) - 1)*A000593(n). - Omar E. Pol, Nov 03 2018
a(n) = Sum_{i=1..n} tau(gcd(n, i)). - Ridouane Oudra, Oct 15 2019
From Peter Bala, Jan 19 2021: (Start)
G.f.: A(x) = Sum_{n >= 1} x^(n^2)*(x^n + n*(1 - x^(2*n)))/(1 - x^n)^2 - differentiate equation 5 in Arndt w.r.t. x, and set x = 1.
A(x) = F(x) + G(x), where F(x) is the g.f. of A079667 and G(x) is the g.f. of A117004. (End)
a(n) = Sum_{k=1..n} tau(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). - Richard L. Ollerton, May 07 2021
With the convention that a(n) = 0 for n <= 0 we have the recurrence a(n) = t(n) + Sum_{k >= 1} (-1)^(k+1)*(2*k + 1)*a(n - k*(k + 1)/2), where t(n) = (-1)^(m+1)*(2*m+1)*n/3 if n = m*(m + 1)/2, with m positive, is a triangular number else t(n) = 0. For example, n = 10 = (4*5)/2 is a triangular number, t(10) = -30, and so a(10) = -30 + 3*a(9) - 5*a(7) + 7*a(4) = -30 + 39 - 40 + 49 = 18. - Peter Bala, Apr 06 2022
Recurrence: a(p^x) = p*a(p^(x-1)) + 1, if p is prime and for any integer x. E.g., a(5^3) = 5*a(5^2) + 1 = 5*31 + 1 = 156. - Jules Beauchamp, Nov 11 2022
Sum_{n>=1} a(n)/exp(2*Pi*n) = 1/24 - 1/(8*Pi) = A319462. - Vaclav Kotesovec, May 07 2023
a(n) < (7n*A001221(n) + 10*n)/6 [Duncan, 1961] (see Duncan and Tattersall). - Stefano Spezia, Jul 13 2025

A002191 Possible values for sum of divisors of n.

Original entry on oeis.org

1, 3, 4, 6, 7, 8, 12, 13, 14, 15, 18, 20, 24, 28, 30, 31, 32, 36, 38, 39, 40, 42, 44, 48, 54, 56, 57, 60, 62, 63, 68, 72, 74, 78, 80, 84, 90, 91, 93, 96, 98, 102, 104, 108, 110, 112, 114, 120, 121, 124, 126, 127, 128, 132, 133, 138, 140, 144, 150, 152, 156

Offset: 1

N. J. A. Sloane

nonn

Distinct values attained by the sigma(n) function, in ascending order.

The asymptotic density of this sequence is 0 (Niven, 1951, Rao and Murty, 1979). - Amiram Eldar, Jul 23 2020

			a(100) = 272, a(10^3) = 3696, a(10^4) = 44496, a(10^5) = 510356, a(10^6) = 5691216. - _M. F. Hasler_, Nov 22 2019

J. W. L. Glaisher, Number-Divisor Tables. British Assoc. Math. Tables, Vol. 8, Camb. Univ. Press, 1940, p. 85.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972, p. 840.
Ivan Niven, The asymptotic density of sequences, Bull. Amer. Math. Soc., Vol. 57 (1951), pp. 420-434.
R. Sita Rama Chandra Rao and G. Sri Rama Chandra Murty, On a theorem of Niven, Canadian Mathematical Bulletin, Vol 22, No. 1 (1979), pp. 113-115.

Cf. A000203, A007609.
Complement of A007369. A175192(a(n)) = 1, A054973(a(n)) >= 1. - Jaroslav Krizek, Mar 01 2010
See A083531 for the gaps, i.e., first differences. - M. F. Hasler, Mar 12 2018
Subsequence of A211347.

N:= 1000: # to get all entries <= N
select(`<=`,{seq(numtheory[sigma](i),i=1..N)},N); # Robert Israel, Jun 16 2014

lim=1000; Select[Union[DivisorSigma[1,Range[lim]]], #<=lim &] (* T. D. Noe, May 06 2010 *)

list(lim)=select(n->n<=lim,Set(vector(lim\=1,n,sigma(n)))) \\ Charles R Greathouse IV, Nov 12 2013

A002191_upto(N,M=N\1+1)=Set(apply(t->min(sigma(t),M), [1..N\1-1]))[^-1] \\ Needs big stack for N >= 10^6; slower alternative: {A002191_upto(N)= my(L=List(1),s); for(n=2,N\=1,N<(s=sigma(n))||listput(L,s));Set(L)}
A2191=A002191_upto(1e4); A002191(n)={#A2191A002191_upto(n*logint(n,10)+n); A2191[n]} \\ - M. F. Hasler, Nov 22 2019

a(n)/n < log_10(n) + O(1) with O(1) <= 1 for all n. - M. F. Hasler, Nov 22 2019

A007369 Numbers n such that sigma(x) = n has no solution.

Original entry on oeis.org

2, 5, 9, 10, 11, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 33, 34, 35, 37, 41, 43, 45, 46, 47, 49, 50, 51, 52, 53, 55, 58, 59, 61, 64, 65, 66, 67, 69, 70, 71, 73, 75, 76, 77, 79, 81, 82, 83, 85, 86, 87, 88, 89, 92, 94, 95, 97, 99, 100, 101, 103, 105, 106, 107, 109, 111, 113

Offset: 1

N. J. A. Sloane, Mira Bernstein, Robert G. Wilson v

nonn

With an initial 1, may be constructed inductively in stages from the list L = {1,2,3,....} by the following sieve procedure. Stage 1. Add 1 as the first term of the sequence a(n) and strike off 1 from L. Stage n+1. Add the first (i.e. leftmost) term k of L as a new term of the sequence a(n) and strike off k, sigma(k), sigma(sigma(k)),.... from L. - Joseph L. Pe, May 08 2002
This sieve is a special case of a more general sieve. Let D be a subset of N and let f be an injection on D satisfying f(n) > n. Define the sieve process as follows: 1. Start with the empty sequence S and let E = D. 2. Append the smallest element s of E to S. 3. Remove s, f(s), f(f(s)), f(f(f(s))), ... from E. 4. Go to step 2. After this sieving process, S = D - f(D). To get the current sequence, take f = sigma and D = {n | n >= 2}. - Max Alekseyev, Aug 08 2005
By analogy with the untouchable numbers (A005114), these numbers could be named "sigma-untouchable". - Daniel Lignon, Mar 28 2014
The asymptotic density of this sequence is 1 (Niven, 1951, Rao and Murty, 1979). - Amiram Eldar, Jul 23 2020

			a(4) = 10 because there is no x < 10 whose sigma(x) = 10.

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 840.
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 (first 1000 terms from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Ivan Niven, The asymptotic density of sequences, Bull. Amer. Math. Soc., Vol. 57 (1951), pp. 420-434.
R. Sita Rama Chandra Rao and G. Sri Rama Chandra Murty, On a theorem of Niven, Canadian Mathematical Bulletin, Vol 22, No. 1 (1979), pp. 113-115.
R. G. Wilson, V, Letter to N. J. A. Sloane, Jul. 1992

Complement of A002191.
See A083532 for the gaps, i.e., first differences.
See A048995 for the missed sums of nontrivial divisors.

a = {}; Do[s = DivisorSigma[1, n]; a = Append[a, s], {n, 1, 115} ]; Complement[ Table[ n, {n, 1, 115} ], Union[a] ]

list(lim)=my(v=List(),u=vectorsmall(lim\1),t); for(n=1,lim, t=sigma(n); if(t<=lim, u[t]=1)); for(n=2,lim, if(u[n]==0, listput(v,n))); Vec(v) \\ Charles R Greathouse IV, Mar 09 2017

A007369_list(LIM,m=0,L=List(),s)={for(n=2,LIM,(s=sigma(n-1))>LIM || bittest(m,s) || m+=1<M. F. Hasler, Mar 12 2018

A175192(a(n)) = 0, A054973(a(n)) = 0. - Jaroslav Krizek, Mar 01 2010

a(n) < 2n + sqrt(8n). - Charles R Greathouse IV, Oct 23 2015

More terms from David W. Wilson

A007368 Smallest k such that sigma(x) = k has exactly n solutions.

Original entry on oeis.org

2, 1, 12, 24, 96, 72, 168, 240, 336, 360, 504, 576, 1512, 1080, 1008, 720, 2304, 3600, 5376, 2520, 2160, 1440, 10416, 13392, 3360, 4032, 3024, 7056, 6720, 2880, 6480, 10800, 13104, 5040, 6048, 4320, 13440, 5760, 18720, 20736, 19152, 22680, 43680

Offset: 0

N. J. A. Sloane, Mira Bernstein, Robert G. Wilson v

nonn

It's not obvious that a(n) exists for all n; I'd like to see a proof. - David Wasserman, Jun 07 2002
Note that k-1 is frequently prime. See A115374 for the least prime. For each n, it appears that there are an infinite number of k such that sigma(x)=k has exactly n solutions. - T. D. Noe, Jan 21 2006
According to Sierpiński, H. J. Kanold proved that there is a k such that sigma(x)=k has n or more solutions. Sierpiński states that Erdős proved that if, for some k, sigma(x)=k has exactly n solutions, then there are an infinite number of such k. - T. D. Noe, Oct 18 2006
Index of the first occurrence of n in A054973. - Jaroslav Krizek, Apr 25 2009

			a(10) = 504; {204, 220, 224, 246, 284, 286, 334, 415, 451, 503} is the set of x such that sigma(x) = 504.

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 840.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe and Donovan Johnson, Table of n, a(n) for n = 0..5000 (terms up to a(429) from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Kevin Ford, Sergei Konyagin, On two conjectures of Sierpiński concerning the arithmetic functions σ and ϕ, arXiv:1910.08452 [math.NT], 2019.
W. Sierpiński, Elementary Theory of Numbers, Warszawa 1964, page 166.
R. G. Wilson, V, Letter to N. J. A. Sloane, Jul. 1992

Cf. A000203, A054973, A002191, A007609.
Cf. A115374 (least prime p such that sigma(x)=sigma(p) has exactly n solutions).
Cf. A007369, A007370, A007371, A007372 (n such that sigma(x)=k has 0, 1, 2 and 3 solutions).
Cf. A184393, A184394, A201915 (smallest solution, largest solution, triangle of solutions for sigma(x)=a(n)).

Needs["Statistics`DataManipulation`"]; s=DivisorSigma[1, Range[10^5]]; f=Frequencies[s]; fs=Sort[f]; tfs=Transpose[fs][[1]]; utfs=Union[tfs]; firstMissing=First[Complement[Range[Last[utfs]], utfs]]; pos=1; Table[While[tfs[[pos]]T. D. Noe *)
terms = 100; cnt = DivisorSigma[1, Range[terms^3]] // Tally // Sort; a[0] = 2; a[n_] := SelectFirst[cnt, #[[2]] == n&][[1]]; Table[a[n], {n, 0, terms - 1}] (* Jean-François Alcover, Jul 18 2017 *)

More terms from David W. Wilson

A007609 Values taken by the sigma function A000203, listed with multiplicity and in ascending order.

Original entry on oeis.org

1, 3, 4, 6, 7, 8, 12, 12, 13, 14, 15, 18, 18, 20, 24, 24, 24, 28, 30, 31, 31, 32, 32, 36, 38, 39, 40, 42, 42, 42, 44, 48, 48, 48, 54, 54, 56, 56, 57, 60, 60, 60, 62, 63, 68, 72, 72, 72, 72, 72, 74, 78, 80, 80, 84, 84, 84, 90, 90, 90, 91, 93, 96, 96, 96, 96, 98, 98

Offset: 1

Walter Nissen

A175192(a(n)) = 1, A054973(a(n)) >= 1. - Jaroslav Krizek, Mar 01 2010

a(n) is the median of the values of A000203(m) from m=1 to m=2n-1. (This needs confirmation as it relies on the assumption that A000203(n) is always bigger than the median of the values A000203(x) from x=1 to x=n.) - Chayim Lowen, May 27 2015

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

Cf. A000203, A002191 (duplicates removed), A007368, A085790.

sort(select(`<=`,map(numtheory:-sigma,[$1..1000]),1001)); # Robert Israel, Jun 04 2015

terms = 68; ClearAll[t]; t[k_] := t[k] = Sort[ Table[ DivisorSigma[1, n], {n, 1, k*terms}]][[1 ;; terms]]; t[k = 2]; While[t[k] != t[k-1], k++]; t[k] (* Jean-François Alcover, Nov 21 2012 *)
With[{nn=80},Take[Sort[DivisorSigma[1,Range[nn*100]]],nn]] (* Harvey P. Dale, Mar 09 2016 *)

list(lim)=select(k->k<=lim,Set(apply(sigma,[1..lim\1]))) \\ Charles R Greathouse IV, Mar 09 2014

a(n) = sigma(A085790(n)). - Jinyuan Wang, Apr 15 2020

A085790 Integers sorted by the sum of their divisors.

Original entry on oeis.org

1, 2, 3, 5, 4, 7, 6, 11, 9, 13, 8, 10, 17, 19, 14, 15, 23, 12, 29, 16, 25, 21, 31, 22, 37, 18, 27, 20, 26, 41, 43, 33, 35, 47, 34, 53, 28, 39, 49, 24, 38, 59, 61, 32, 67, 30, 46, 51, 55, 71, 73, 45, 57, 79, 44, 65, 83, 40, 58, 89, 36, 50, 42, 62, 69, 77, 52, 97, 101, 63, 103, 85

Offset: 1

Hugo Pfoertner, Jul 23 2003

Integers having the same sum of divisors are sorted in ascending order, e.g., sigma(14)=sigma(15)=sigma(23)=24 -> a(15)=14, a(16)=15, a(17)=23.
Also an irregular triangle where the k-th row consists of all numbers with divisor sum k. See A054973(k) for the k-th row length. - Jeppe Stig Nielsen, Jan 29 2015
By definition this is a permutation of the positive integers. Also positive integers of A299762. - Omar E. Pol, Mar 14 2018

			a(9) = 9, a(10) = 13, a(11) = 8 because sigma(9) = 9 + 3 + 1 = 13, sigma(13) = 13 + 1 = 14, sigma(8) = 8 + 4 + 2 + 1 = 15 and there are no other numbers with those sigma values.
Irregular triangle starts: (row numbers to the left are not part of the sequence)
   n : row(n)
   1 : 1,
   2 :
   3 : 2,
   4 : 3,
   5 :
   6 : 5,
   7 : 4,
   8 : 7,
   9 :
  10 :
  11 :
  12 : 6, 11,
  13 : 9,
  14 : 13,
  15 : 8,
  16 :
  17 :
  18 : 10, 17,
  19 :
  20 : 19,
  21 :
  22 :
  23 :
  24 : 14, 15, 23,
  25 :
- _Jeppe Stig Nielsen_, Feb 02 2015, edited by _M. F. Hasler_, Nov 21 2019

Hugo Pfoertner, Table of n, a(n) for n=1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems: invphi.gp, Oct. 2005
Jeppe Stig Nielsen, First 10000 rows of the triangle for a(n)
Index entries for sequences that are permutations of the natural numbers

Cf. A000203 (sigma), A007609 (values taken by sigma, with multiplicity), A002191 (possible values for sigma), A002192 (first column).

Cf. A152454 (similar sequence for proper divisors only (aliquot parts)).

SortBy[Table[{n,DivisorSigma[1,n]},{n,120}],Last][[;;,1]] (* Harvey P. Dale, Sep 10 2024 *)

A085790_row(n)=invsigma(n) \\ Cf. Alekseyev link for invsigma(). - M. F. Hasler, Nov 21 2019

A066938 Primes of the form p*q+p+q, where p and q are primes.

Original entry on oeis.org

11, 17, 23, 31, 41, 47, 53, 59, 71, 79, 83, 89, 107, 113, 127, 131, 151, 167, 179, 191, 227, 239, 251, 263, 269, 271, 293, 311, 359, 383, 419, 431, 439, 443, 449, 479, 491, 503, 521, 587, 593, 599, 607, 631, 647, 659, 683, 701, 719, 727, 743, 773, 809, 827

Offset: 1

Reinhard Zumkeller, Jan 24 2002

nonn

For p not equal to q, either p*q or p+q is odd, so their sum is odd.
The representation is ambiguous, e.g. 2*7+2+7 = 23 = 3*5+3+5.
Complement of A198273 with respect to A000040. - Reinhard Zumkeller, Oct 23 2011
None of these primes are in A158913 since if p*q+p+q is a prime, then sigma(p*q+p+q) = sigma(p*q). - Amiram Eldar, Nov 15 2021

			59 is in the sequence because 59 = 2 * 19 + 2 + 19.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A054973, A067432, A072668, A072673, A088709, A158913, A198273, A198277.

a066938 n = a066938_list !! (n-1)
a066938_list = map a000040 $ filter ((> 0) . a067432) [1..]
-- Reinhard Zumkeller, Oct 23 2011

nn = 1000; n2 = PrimePi[nn/3]; Select[Union[Flatten[Table[(Prime[i] + 1) (Prime[j] + 1) - 1, {i, n2}, {j, n2}]]], # <= nn && PrimeQ[#] &]

is(n)=fordiv(n+1,d,my(p=d-1,q=(n+1)/d-1); if(isprime(p) && isprime(q), return(isprime(n)))); 0 \\ Charles R Greathouse IV, Jul 23 2013

A067432(A049084(a(n))) > 0. - Reinhard Zumkeller, Oct 23 2011

A054973(a(n)+1) >= 2. - Amiram Eldar, Nov 15 2021

Edited by Robert G. Wilson v, Feb 01 2002

A066075 Number of solutions x to prime(n) = sigma(x) - 1, where prime(n) is the n-th prime.

Original entry on oeis.org

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

Offset: 1

Labos Elemer, Dec 03 2001

nonn

prime(n) itself is always the largest solution, but often composite solutions also occur.

If a(n) = 1, then the single solution is prime(n).

			For n = 96, prime(96) = 503, 503 = sigma(x)-1 has 10 solutions together with 503: {204, 220, 224, 246, 284, 286, 334, 415, 451, 503}, so a(96) = 10.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).

Number of solutions to A000040(n) = A000203(x) - 1.

Cf. A000040, A000203, A054973, A058340, A066071, A066072, A066073, A066074, A066075, A066076, A066077, A066080.

{ for (n=1, 1000, a=1; for (x=1, prime(n) - 1, if (prime(n) == (sigma(x) - 1), a++)); write("b066075.txt", n, " ", a) ) } \\ Harry J. Smith, Nov 10 2009

a(n) = invsigmaNum(prime(n)+1); \\ Amiram Eldar, Dec 16 2024, using Max Alekseyev's invphi.gp

a(n) = A054973(prime(n)+1). - Amiram Eldar, Dec 16 2024

cp's OEIS Frontend

A159953 Values in A054973 larger than 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A258968 a(n) is the least positive integer x with A054973(x) = n and A258913(x) < x.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

A000203 a(n) = sigma(n), the sum of the divisors of n. Also called sigma_1(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Magma

Maple

Mathematica

Maxima

MuPAD

PARI

PARI

PARI

PARI

Python

Python

SageMath

Scheme

Scheme

Formula

A002191 Possible values for sum of divisors of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

A007369 Numbers n such that sigma(x) = n has no solution.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions