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.
%I A331805 #82 Mar 18 2024 11:49:55
%S A331805 42,1316,131080256
%N A331805 Integers k such that k is equal to the sum of the nonprime proper divisors of k.
%C A331805 The number 37778715690312487141376 is also in the sequence. - _Daniel Suteu_, Jan 27 2020
%C A331805 The first 3 terms have the form (2^p-1)*(2^(p-1))*((2^p-1)^2-2), i.e., a Perfect number times a Carol prime. - _G. L. Honaker, Jr._, Jan 27 2020
%C A331805 In other words, the values of p are given by the intersection of A091515 and A000043. Currently, only four such values of p are known: {2, 3, 7, 19}. - _Daniel Suteu_, Jan 27 2020
%C A331805 From _Bernard Schott_, Jan 28 2020: (Start)
%C A331805 Proposition: If a number N_p is of the form Q_p * C_p where Q_p = (2^(p-1)) * (2^p - 1) is a perfect number and C_p = (2^p-1)^2-2 is a Carol prime then, the sum of the nonprime proper divisors of N_p called S_p(N_p) is equal to N_p.
%C A331805 Proof:
%C A331805 The sum of the nonprime proper divisors of N_p is:
%C A331805 S_p(N_p) = (2* Q_p - 2 - (2^p-1)) + ((Q_p - 1) * C_p).
%C A331805 In the first parenthesis, there is the sum of the nonprime proper divisors of N_p coming only from the perfect number Q_p, then in the second parenthesis, there is the sum of the nonprime proper divisors of N_p coming from C_p.
%C A331805 Then, this sum of the nonprime proper divisors of N_p, S_p(N_p) is indeed equal to N_p = (2^(p-1)) * (2^p-1) * ((2^p-1)^2-2).
%C A331805 Hence, (2^19-1)*(2^(19-1))*((2^19-1)^2-2) = 37778715690312487141376 is a term. (End)
%C A331805 10^13 < a(4) <= 72872313094554244192 = 2^5 * 109 * 151 * 65837 * 2101546957. - _Giovanni Resta_, Jan 28 2020
%H A331805 Chris K. Caldwell and G. L. Honaker, Jr., <a href="https://primes.utm.edu/curios/page.php?curio_id=37395">Prime Curio for 42</a>
%H A331805 Wikipedia, <a href="https://en.wikipedia.org/wiki/42_(number)">42 (number)</a>
%e A331805 42 is a term because 42 = 1 + 6 + 14 + 21.
%e A331805 1316 is a term because 1316 = 1 + 4 + 14 + 28 + 94 + 188 + 329 + 658.
%t A331805 fun[p_, e_] := (p^(e+1) - 1)/(p - 1); npsigma[n_] := Times @@ fun @@@ (f = FactorInteger[n]) - Plus @@ First /@ f;; Select[Range[2, 1500], npsigma[#] == 2# &] (* _Amiram Eldar_, Jan 26 2020 *)
%o A331805 (PARI) isok(n) = sigma(n) - n - vecsum(factor(n)[,1]) == n; \\ _Daniel Suteu_, Jan 27 2020
%Y A331805 Cf. A001065, A018252, A023890, A331858.
%Y A331805 Cf. A000043, A091515, A091516 (Carol primes).
%K A331805 nonn,bref,more
%O A331805 1,1
%A A331805 _G. L. Honaker, Jr._, Jan 26 2020
%E A331805 a(2) from Chuck Gaydos
%E A331805 a(3) from _Amiram Eldar_, Jan 26 2020