OFFSET
1,2
COMMENTS
If p is prime p^2 is in the sequence since the continued fraction for Sum_{d|n} 1/d is [1, p-1, p+1] and there are 3 divisors for p^2.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
aQ[n_] := ! PrimeQ[n] && Length@ContinuedFraction[DivisorSigma[1, n]/n] == DivisorSigma[0, n]; Select[Range[488], aQ] (* Amiram Eldar, Aug 30 2019 *)
PROG
(PARI) for(n=1, 1000, if(length(contfrac(sumdiv(n, d, 1/d)))==numdiv(n)*(1-isprime(n)), print1(n, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Jun 09 2002
STATUS
approved