cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A005114 Untouchable numbers, also called nonaliquot numbers: impossible values for the sum of aliquot parts function (A001065).

Original entry on oeis.org

2, 5, 52, 88, 96, 120, 124, 146, 162, 188, 206, 210, 216, 238, 246, 248, 262, 268, 276, 288, 290, 292, 304, 306, 322, 324, 326, 336, 342, 372, 406, 408, 426, 430, 448, 472, 474, 498, 516, 518, 520, 530, 540, 552, 556, 562, 576, 584, 612, 624, 626, 628, 658
Offset: 1

Views

Author

Keywords

Comments

Complement of A078923. - Lekraj Beedassy, Jul 19 2005
Chen & Zhao show that the lower density of this sequence is at least 0.06, improving on te Riele. - Charles R Greathouse IV, Dec 28 2013
Numbers k such that A048138(k) = 0. A048138(k) measures how "touchable" k is. - Jeppe Stig Nielsen, Jan 12 2020
From Amiram Eldar, Feb 13 2021: (Start)
The term "untouchable number" was coined by Alanen (1972). He found the 570 terms below 5000.
Erdős (1973) proved that the lower asymptotic density of untouchable numbers is positive, te Riele (1976) proved that it is > 0.0324, and Banks and Luca (2004, 2005) proved that it is > 1/48.
Pollack and Pomerance (2016) conjectured that the asymptotic density is ~ 0.17. (End)
The upper asymptotic density is less than 1/2 by the 'almost all' binary Goldbach conjecture, independently proved by Nikolai Chudakov, Johannes van der Corput, and Theodor Estermann. (In this context, this shows that the density of the odd numbers of this form is 0 (consider A001065(p*q) for prime p, q); full Goldbach would prove that 5 is the only odd number in this sequence.) - Charles R Greathouse IV, Dec 05 2022

References

  • Richard K. Guy, Unsolved Problems in Number Theory, 3rd ed., Springer, 2004, section B10, pp. 100-101.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 65.
  • József Sándor, Dragoslav S. Mitrinovic, Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, page 93.
  • 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, page 147.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 125.

Crossrefs

Programs

  • Mathematica
    untouchableQ[n_] := Catch[ Do[ If[n == DivisorSigma[1, k]-k, Throw[True]], {k, 0, (n-1)^2}]] === Null; Reap[ Table[ If[ untouchableQ[n], Print[n]; Sow[n]], {n, 2, 700}]][[2, 1]] (* Jean-François Alcover, Jun 29 2012, after Benoit Cloitre *)
  • PARI
    isA078923(n)=if(n==0 || n==1, return(1)); for(m=1,(n-1)^2, if( sigma(m)-m == n, return(1))); 0
    isA005114(n)=!isA078923(n)
    for(n=1,700, if (isA005114(n), print(n))) \\ R. J. Mathar, Aug 10 2006
    
  • PARI
    is(n)=if(n%2 && n<4e18, return(n==5)); forfactored(m=1,(n-1)^2, if(sigma(m)-m[1]==n, return(0))); 1 \\ Charles R Greathouse IV, Dec 05 2022
    
  • Python
    from sympy import divisor_sigma as sigma
    from functools import cache
    @cache
    def f(m): return sigma(m)-m
    def okA005114(n):
        if n < 2: return 0
        return not any(f(m) == n for m in range(1, (n-1)**2+1))
    print([k for k in range(289) if okA005114(k)]) # Michael S. Branicky, Nov 16 2024
    
  • Python
    # faster for intial segment of sequence
    from itertools import count, islice
    from sympy import divisor_sigma as sigma
    def agen(): # generator of terms
        n, touchable, t = 2, {0, 1}, 1
        for m in count(2):
            touchable.add(sigma(m)-m)
            while m > t:
                if n not in touchable:
                    yield n
                else:
                    touchable.discard(n)
                n += 1
                t = (n-1)**2
    print(list(islice(agen(), 20))) # Michael S. Branicky, Nov 16 2024

Extensions

More terms from David W. Wilson