Fact-checked by Grok 2 weeks ago

Aliquot sequence

An aliquot sequence is a sequence of positive integers defined by starting with an initial number n > 0 and iteratively applying the aliquot sum function s(k), which computes the sum of the proper divisors of k (all positive divisors excluding k itself), given by s(k) = \sigma(k) - k where \sigma(k) is the sum of all positive divisors of k. Thus, the sequence is a_1 = n, a_{k+1} = s(a_k) for k \geq 1. Aliquot sequences trace their origins to ancient , particularly the Pythagorean study of perfect numbers, where s(n) = n, forming a cycle of length 1, with the smallest example being 6 since s(6) = 1 + 2 + 3 = 6. They also encompass amicable pairs, such as and 284, where s(220) = 284 and s(284) = 220, creating a cycle of length 2, and sociable numbers that form longer cycles, like the 28-term cycle starting at 14316. Many sequences terminate by reaching a prime p (where s(p) = 1), followed by s(1) = 0, effectively ending at 0. The behavior of aliquot sequences remains a central topic in , with possibilities including termination, periodic cycling, or unbounded growth. The –Dickson conjecture, proposed by in 1888 and extended by Leonard Eugene Dickson in 1913, asserts that every aliquot sequence either terminates or enters a , but this is widely doubted following computational evidence of diverging sequences, such as those starting at 276, which grow without bound. In 1975, H. W. Lenstra Jr. proved the existence of strictly increasing aliquot sequences of any finite length k, while showed in 1976 that the set of numbers n for which the aliquot sequence strictly increases for k steps has asymptotic density zero. Statistical analyses indicate that, on average, terms in sequences starting from even numbers tend to decrease slightly ( ratio \approx e^{-0.033}), but for multiples of 4, they grow moderately (\approx e^{0.175}), supporting the likelihood of unbounded growth in certain cases. Open questions abound, including whether there are infinitely many aliquot cycles beyond known lengths, the existence of cycles of length 3, and the asymptotic density of sociable numbers, which is conjectured to be zero. As of 2025, the longest known sociable cycle remains the 28-term one, and no cycles of length 3 have been found. Computational efforts, such as the Aliquot Sequences Project, have explored sequences up to enormous terms (e.g., over 100 digits for starting value 3630, which terminates), revealing no new cycles but confirming divergences. These sequences connect to broader themes in divisor theory, including untouchable numbers (those not equal to s(m) for any m), whose density is heuristically around 0.17.

Fundamentals

Definition

An aliquot sequence begins with a positive n > 0 and is constructed iteratively by applying the aliquot sum function, which computes the sum of the proper divisors of the previous term. The proper divisors of a number are all positive divisors excluding the number itself. The foundation of this construction relies on the \sigma(n), which gives the sum of all positive divisors of n. For example, the divisors of 10 are 1, 2, 5, and 10, so \sigma(10) = 18. The aliquot sum s(n) is then defined as s(n) = \sigma(n) - n, yielding s(10) = 8 in this case. Formally, the aliquot sequence starting from n is the sequence of terms s^k(n) for k = 0, 1, 2, \dots, where s^0(n) = n, s^1(n) = s(n), and s^k(n) = s(s^{k-1}(n)) for k \geq 2. By convention, the sequence terminates at 0 upon reaching 1, since the proper divisors of 1 are none, so s(1) = 0.

Basic properties

The terms in an aliquot sequence are classified based on their abundance relative to the aliquot sum function s(n): a term n is deficient if s(n) < n, perfect if s(n) = n, or abundant if s(n) > n. This classification influences the sequence's progression, with sequences often exhibiting alternations between these categories depending on the primality or prime factorization of the terms; for instance, prime numbers and prime powers tend to produce deficient terms leading to decreases, while highly composite numbers can yield abundant terms that cause increases. Aliquot sequences can be conceptualized as directed paths in an infinite whose nodes are the positive s and whose edges connect each n to s(n). This graphical representation highlights the structure of the sequences as trajectories traversing the , potentially merging with other paths or entering . A fundamental property of aliquot sequences starting from even s n > 0 is that they eventually reach an or enter a (with known being ). This occurs because s(n) \equiv n \pmod{2} unless n is a square or twice a square, in which case the parity changes; thus, even sequences remain even until such a form is encountered, transitioning to . For n, s(n) remains unless n is an odd square, in which case s(n) is even. For powers of 2, the aliquot sum admits a closed form: s(2^k) = 2^k - 1 for k \geq 1. To derive this, note that the proper divisors of $2^k are $1, 2, 4, \dots, 2^{k-1}, forming a whose sum is $1 + 2 + \cdots + 2^{k-1} = 2^k - 1. Thus, the sequence starting at $2^k proceeds as $2^k, 2^k - 1, s(2^k - 1), \dots, where $2^k - 1 (a Mersenne number) is odd and deficient, leading to a rapid descent to 1 and then 0, ensuring termination.

Examples and Classifications

Terminating sequences

Terminating aliquot sequences are those that eventually reach 1 and then 0, effectively ending the process. This occurs when the sequence encounters a prime number p, for which the sum of proper divisors s(p) = 1, followed by s(1) = 0. All prime numbers and 1 terminate immediately in this manner, as their proper divisor sums are trivially 1 and 0, respectively. A concrete example is the aliquot sequence starting at 10: s(10) = 1 + 2 + 5 = 8, s(8) = 1 + 2 + 4 = 7, s(7) = 1, and s(1) = 0. Here, the sequence descends through composite numbers until hitting the prime 7, after which it terminates rapidly. This illustrates the typical descent in terminating sequences, where each step reduces the value until a prime is reached. Certain numbers, known as or nonaliquot numbers, cannot appear as terms in any aliquot sequence except possibly as starting points, because no has them as its of proper divisors. Examples include 2, 5, and ; despite this property, most starting values, including many , lead to terminating sequences. Early computations of aliquot sequences, beginning with tabulations by Eugene Dickson in 1913 and extensive work by Paul Poulet in 1918, revealed that the vast majority of small starting values terminate. For instance, powers of 2 terminate via their Mersenne number counterparts: s(2^k) = 2^k - 1, which either is prime (leading directly to ) or factors into smaller terms that continue descending to a prime. Sequences originating from deficient numbers—those with abundance measure d(n) = s(n)/n < 1—often terminate quickly, as the terms strictly decrease toward smaller values until a prime is encountered. This descent is a key characteristic driving termination in such cases.

Cyclic sequences

In an aliquot sequence, a cycle occurs when the sequence returns to a previously encountered term, thereby forming a periodic loop of length k \geq 1. Such cycles represent bounded behaviors distinct from sequences that terminate at 0. Cycles are classified by their length. For k=1, the sequence is fixed at a perfect number n, where the sum of proper divisors s(n) = n. The smallest such number is 6, whose proper divisors are 1, 2, and 3, summing to 6; thus, the sequence is 6 → 6 → ⋯. All known perfect numbers are even and of the form $2^{p-1}(2^p - 1) for prime p where $2^p - 1 is a ; 52 such numbers are known as of 2025. For k=2, the cycle consists of an amicable pair (m, n) with m \neq n, s(m) = n, and s(n) = m. The smallest pair is (220, 284), discovered in antiquity and attributed to around 300 CE. The proper divisors of 220 are 1, 2, 4, 5, 10, 11, 20, 22, 44, 55, 110, summing to 284; those of 284 are 1, 2, 4, 71, and 142, summing to 220. Thus, the sequence cycles as 220 → 284 → 220 → ⋯. Over 1.2 billion amicable pairs are known as of 2023, with ongoing searches yielding more. For k > 2, the numbers are sociable, forming cycles longer than two. The smallest known sociable cycle, of length 4, was discovered by H. J. J. te Riele in (initially attributed to ) and starts at 1264460:
1264460 → 1547860 → 1727636 → 1305184 → 1264460 → ⋯,
where each term is the sum of proper divisors of the previous. As of July 2025, 5433 sociable cycles of length greater than 2 are known, with 5421 of length 4, one of length 5, five of length 6, four of length 8, one of length 9, and one of length 28 (discovered in 1995). These cycles are predominantly even and consist of a mix of abundant and deficient numbers, though no cycles of length 3 are known.

Conjectures and Open Problems

Catalan-Dickson conjecture

The Catalan-Dickson conjecture posits that every aliquot sequence, starting from any positive integer n, either terminates by reaching 0 or eventually enters a cycle of the form of a , an amicable pair, or a sociable cycle of longer length. This statement was first proposed by around 1888 based on computations for small values of n, where he observed that sequences appeared to end at 1 (equivalent to terminating at 0 in modern notation) or a . Leonard Eugene Dickson extended the conjecture in the early to account for amicable pairs and longer cycles, formalizing it in his comprehensive treatment of divisor functions and related sequences. The implications of the are profound for the study of aliquot sequences, as it asserts the absence of unbounded sequences that grow indefinitely without repeating. Under this hypothesis, all aliquot sequences are eventually periodic or finite, aligning with the behavior observed in terminating cases (where the sum of proper divisors reaches a prime, leading to 1 and then 0) or cycling cases (where the sequence returns to a previous term after a fixed ). This boundedness would imply a complete of aliquot in terms of known or discoverable cycles, contrasting with earlier, more limited views that overlooked sociable numbers beyond pairs. While no definitive counterexamples have been identified, the growth observed in some open sequences has led to doubts about the conjecture's validity, despite extensive computational verification as of the latest available data in . As of , all aliquot sequences initiated by starting values n \leq 10^6 have either terminated at 0, entered a known , or progressed to sufficiently large terms without proven , leaving approximately 18,361 sequences open. These results, derived from efforts, highlight ongoing challenges for deeper exploration.

Lehmer five and unbounded sequences

The Lehmer five refers to the five smallest starting values for aliquot sequences whose ultimate behavior remains unresolved: 276, 552, 564, 660, and 966. These were identified in by mathematician Lehmer during his manual computations of aliquot sequences for all starting numbers up to , where he determined that these particular even, abundant numbers did not terminate or enter a known cycle within the limits of his calculations. As of , extensive computational efforts continue on these sequences without resolution. For instance, the sequence starting at 276 has been extended to 2157, with the most recent term exceeding 200 digits, yet showing no sign of cycling or terminating. The other four sequences exhibit similar open statuses, having reached indices in the thousands and terms with over 180 digits each, with suspicions of eventual periodic merging into larger known trajectories based on patterns observed in related computations. No resolutions have occurred for any of the Lehmer five between and , and they form part of a broader set of unresolved sequences starting below 10,000, such as 1074. These sequences are characterized by a "" pattern of overall growth punctuated by occasional sharp drops, typically when a is fully factored and reveals smaller sums. Numerical evidence from their trajectories suggests phases of increase in size, raising the that, should the Catalan-Dickson prove false—predicting all sequences must eventually or terminate—these could diverge unbounded to .

Computational Exploration

Search methods

Early explorations of aliquot sequences relied on manual computations, as pioneered by in the late , who examined small starting values by hand to observe patterns in their behavior. These efforts were limited to sequences with modest terms, often terminating quickly or cycling in amicable pairs, due to the labor-intensive process of listing divisors without computational aids. By the 1970s, the advent of electronic computers enabled more systematic investigations, such as the aliquot initiated by and John Selfridge, which used early programming to extend sequences beyond manual reach and test conjectures on their boundedness. Contemporary search methods model aliquot sequences as paths in an infinite directed graph, known as the aliquot graph, where each node represents a positive integer and a directed edge connects n to s(n), the sum of its proper divisors. Exploration typically employs breadth-first search (BFS) or depth-first search (DFS) traversals starting from seed values, prioritizing even starting numbers since odd sequences often terminate rapidly at 1. This graph-theoretic framework allows detection of merges, where distinct sequences converge to a common term, and cycles, by tracking visited nodes to avoid redundant computation. The core of any search is efficient computation of s(n) = \sigma(n) - n, where \sigma(n) is the sum-of-divisors function, which is multiplicative and thus computed from the prime of n as \sigma(n) = \prod_{p^k \| n} \frac{p^{k+1} - 1}{p - 1}. For small n (up to around 10^12), trial division suffices to find factors by checking divisibility up to \sqrt{n}. Larger terms, common in extended sequences, require advanced factorization algorithms: the elliptic curve method () excels at discovering factors of 20-60 digits, while efficiently handles composites up to about 100 digits by simulating a pseudorandom walk in the modulo n. For numbers exceeding 100 digits, the general number field sieve (GNFS) is employed, though its high computational cost limits its use to critical terms. To scale searches across vast parameter spaces, parallelization leverages among volunteers, coordinated through forums like mersenneforum.org, where participants reserve specific sequences (e.g., starting below 10^7) to extend using custom software such as the alique program. Projects like those hosted on rechenkraft.net track progress for over 3 million starting values, integrating with databases for collaborative effort, while tools like YAFU facilitate on-demand in BOINC-like environments. This distributed model has enabled exploration of sequences up to terms with 10^7 starting indices. Optimizations are crucial for feasibility, including caching precomputed s(n) values in databases like factordb.com, which stores factorizations and parts for billions of integers, allowing instant retrieval and preventing recomputation when sequences merge into known paths. Searches are bounded by term size, typically halting at 10^6 digits if growth persists, to manage resources, with heuristics like "" analysis (tracking the 2-adic valuation) used to predict and prune unstable branches.

Known results and records

By the year 2000, all aliquot sequences starting from positive integers less than 1000 had been fully resolved except for the five known as the Lehmer five (276, 552, 564, 660, and 966), which remain open and are conjectured to be unbounded. A notable computational achievement in this domain is the termination of the sequence starting at 3630, which reached a maximum of 100 digits at index 1263 before converging to 1 at index 2624; this result, initially reported in 2002, was confirmed and highlighted as a record in subsequent analyses around 2007. Regarding cycle discoveries, 5,433 sociable cycles (of length greater than 2) are known as of July 2025, primarily of length 4, with one each of lengths 5, 9, and 28, and a few of lengths 6 and 8. In recent years, extensive computational searches have continued to uncover new cycles, particularly of short lengths. Record lengths for computed sequences include the paired sequences starting at 314718 and 4788, which merged at index 6466 (with 314718 reaching index 6466 equivalent to 4788 at index 6), establishing a pre-2021 benchmark exceeding 6000 terms before stalling at large composites. More recently, the sequence for 276 has been extended to index 2157 as of November 2025, reaching more than 215 digits, underscoring the computational challenges for the Lehmer five. Ongoing projects such as rechenkraft.net systematically track and compute sequences starting below 4 million, with status updates provided quarterly through integration with the Factor Database and community reservations. Similarly, aliquot.de maintains detailed records for the Lehmer five, with the last major extension reported in 2014, advancing these sequences to over 200 digits in some cases. In 2025, sequences like 3326040 have been confirmed to terminate, highlighting continued advancements in the field. In recent years (2021–2025), numerical analyses have focused on growth statistics, such as average term expansion rates and distribution of sequence behaviors across large samples, as detailed in a 2021 study examining iterations of the aliquot sum function. Explorations of sequences starting from powers of small primes have also seen updates in 2025, revealing patterns in early growth but no resolutions for open cases.

References

  1. [1]
    [2110.14136] Numerical and Statistical Analysis of Aliquot Sequences
    Oct 27, 2021 · Abstract:We present a variety of numerical data related to the growth of terms in aliquot sequences, iterations of the function s(n) ...
  2. [2]
    [PDF] Aliquot Sequences - Dartmouth Mathematics
    Oct 2, 2020 · Lenstra (1975) showed that for every k there is an increasing aliquot sequence of length k, and the next year, Erd˝os showed that this commonly ...
  3. [3]
    Aliquot sequence 3630 ends after reaching 100 digits - Project Euclid
    Abstract. In this paper we present a new computational record: the aliquot sequence starting at 3630 converges to 1 after reaching a hundred decimal digits.
  4. [4]
    Aliquot Sequence -- from Wolfram MathWorld
    The sequence of numbers s^0(n)=n,s^1(n)=s(n),s^2(n)=s(s(n)),... is called an aliquot sequence. If the sequence for a given
  5. [5]
    Aliquot sequences - OeisWiki
    An aliquot sequence (or aliquot trajectory) is a recurrence in which each term is the sum of proper divisors of the previous term, where the initial term ...
  6. [6]
    [PDF] Illustration of Aliquot Sequence Mergers - OEIS
    An aliquot sequence is an iterative mapping of positive integers n → σ(n) − n, from n to the sum of its divisors that are less than n [5, 1, 8, 7][6, Sec. B5]. ...Missing: definition | Show results with:definition
  7. [7]
    What Drives an Aliquot Sequence ? - American Mathematical Society
    WHAT DRIVES AN ALIQUOT SEQUENCE? 105. If the prime 7 is missing, this product ... GUY & J. L. SELFRIDGE, "Interim report on aliquot series," Proc.
  8. [8]
    [PDF] CWI Scanprofile/PDF/300
    A tabulation of the s function was given by Dickson (1913), and Poulet (1929) computed several long aliquot series until a term increased beyond his practical ...
  9. [9]
    Aliquot Cycle -- from Wolfram MathWorld
    If an aliquot sequence {s^0(n),s(n),s^2(n),...} for a given n is bounded, it either ends at s(1)=0 or becomes periodic.
  10. [10]
    [PDF] On the distribution of amicable numbers - Dartmouth Mathematics
    Amicable numbers have a very old history dating back at least to Pythagoras who was aware that 220 and 284 form an amicable pair. We now know of more than ...
  11. [11]
    [PDF] Sociable numbers - Dartmouth Mathematics
    Apr 22, 2009 · (The smallest sociable number of order. 4 was found by Cohen in 1970; it is 1,264,460.) We know 46 perfect numbers and about 12 million amicable.
  12. [12]
    A list of currently known aliquot cycles of length greater than 2 - djm.cc
    There are 5433 cycles in all: 5421 of length 4, 1 of length 5, 5 of length 6, 4 of length 8, 1 of length 9, and 1 of length 28. P. Moews and I have ...
  13. [13]
    Catalan's Aliquot Sequence Conjecture -- from Wolfram MathWorld
    The conjecture proposed by Catalan in 1888 and extended by E. Dickson that each aliquot sequence ends in a prime, a perfect number, or a set of sociable ...
  14. [14]
    History of the Theory of Numbers: Divisibility and Primality, Volume 1
    Jan 27, 2012 · This 1st volume in the series History of the Theory of Numbers presents the material related to the subjects of divisibility and primality.
  15. [15]
    Aliquot sequence - Prime-Wiki
    Feb 21, 2023 · Definiton. An aliquot sequence is a sequence of numbers generated from an initial number using the sigma σ ( n ) function.Missing: theory | Show results with:theory
  16. [16]
    [PDF] aliquot sequence 3630 ends after reaching 100 digits
    Abstract. In this paper we present a new computational record: the aliquot sequence starting at 3630 converges to 1 after reaching a hundred decimal digits.
  17. [17]
    factordb.com
    - **Status**: Unchecked
  18. [18]
    The aliquot project - ACM Digital Library
    INTRODUCTION. This paper is divided into two parts. Part i first presents an old and charming number theoretic conjecture which has been gen-.
  19. [19]
    [PDF] Aliquot sequences
    This infinite graph has several topological properties and as regards the. "lengths" of its branches. These properties are listed on this page and can not be ...
  20. [20]
    [PDF] Aliquot sequences with small starting values - Radboud Universiteit
    As a rather naive indication for the truth of the Catalan-Dickson conjecture, we have also calculated a first order, linear, approximation (−0.019·x+ 18.93) to ...
  21. [21]
    Intro to Aliquot Sequences Analysis - Rechenkraft.net
    This is a nearly-elementary introduction to the theory and analysis of aliquot sequences, assuming only a very little bit of basic number theory.
  22. [22]
    Current status of aliquot sequences with start term below 4 million
    The current settings of a 90 day max update period, and 200 sequences per hour, is suspected to be more than sufficient to support a table size of at least five ...Missing: Lehmer | Show results with:Lehmer
  23. [23]
  24. [24]
    factordb.com
    **Summary of factordb.com Usage in Computing Aliquot Sequences:**
  25. [25]
    Aliquot Sequence 3630 Ends After - Reaching 100 Digits
    Mar 4, 2021 · Aliquot sequences whose end is in doubt. go up to 200 or 1000 digit numbers for n = 276, which would make this sequence out of reach.
  26. [26]
  27. [27]
    A292217 - OEIS
    Up to the known 1593 sociable number cycles, 96.1% of the sociable number cycles satisfy this condition (up to the first 10 sociable number cycles: 40%; up ...
  28. [28]
    Aliquot Sequences from the Trenches - Rechenkraft.net
    Serge Batalov got his second terminating sequence in 134856 after it had reached a maximum height of 102 digits. His computations began where Wolfgang ...
  29. [29]
    Lehmer Five - Aliquot-Seiten und Westsahara
    Die ersten fünf Offenendketten werden "Lehmer Five" genannt. Der Begriff wurde ursprünglich für die offenen Sequenzen mit den Startzahlen 276, 552, 564, 660, ...Top · data