OFFSET
1,1
COMMENTS
For distinct primes p,q, includes p^a*q^b if (a+1)*(b+1)=p+q. - Robert Israel, Jan 05 2018
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
MAPLE
filter:= n -> numtheory:-tau(n)=convert(numtheory:-factorset(n), `+`):
select(filter, [$1..10^5]); # Robert Israel, Jan 05 2018
MATHEMATICA
Select[Range[20000], Plus @@ ((f = FactorInteger[#])[[;; , 1]]) == Times @@ (f[[;; , 2]] + 1) &] (* Amiram Eldar, Jun 06 2022 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 13 2002
STATUS
approved