OFFSET
1,1
COMMENTS
Banks and Luca (2007) showed that the number of terms <= x, N(x) <= x * exp(-((1/3)*(log(8))^(1/3) + o(1))*(log(x))^(1/3)*(log(log(x)))^(1/3)) as x -> infinity, and that under Dickson's conjecture this sequence is infinite, since for each positive integer m, if p = 5m + 1 and q = 20m + 13 are primes, then p*q is a term. - Amiram Eldar, Apr 13 2020
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..10000
William D. Banks and Florian Luca, When the sum of aliquots divides the totient, Proceedings of the Edinburgh Mathematical Society, Vol. 50, No. 3 (2007), pp. 563-569.
Douglas E. Iannucci, On the Equation sigma(n) = n + phi(n), Journal of Integer Sequences, Vol. 20 (2017), Article 17.6.2.
EXAMPLE
k = 133 = 7*19: phi(133)=108, sigma(133)-133 = 1+7+19 = 27, q = 4.
MATHEMATICA
Do[s=EulerPhi[n]/(DivisorSigma[1, n]-n); If[ !PrimeQ[n]&&IntegerQ[s], Print[n]], {n, 2, 1000000}]
Select[Range[54000], CompositeQ[#]&&IntegerQ[EulerPhi[#]/(DivisorSigma[ 1, #]-#)]&] (* Harvey P. Dale, Nov 16 2021 *)
PROG
(PARI) is(n)=!isprime(n) && n>1 && eulerphi(n)%(sigma(n)-n)==0 \\ Charles R Greathouse IV, Jan 02 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jul 22 2000
STATUS
approved